Differential transcriptomics in sarcoidosis lung and lymph node granulomas with comparisons to pathogen-specific granulomas.

RATIONALE: Despite the availability of multi-"omics" strategies, insights into the etiology and pathogenesis of sarcoidosis have been elusive. This is partly due to the lack of reliable preclinical models and a paucity of validated biomarkers. As granulomas are a key feature of sarcoidosis, we speculate that direct genomic interrogation of sarcoid tissues, may lead to identification of dysregulated gene pathways or biomarker signatures. OBJECTIVE: To facilitate the development sarcoidosis genomic biomarkers by gene expression profiling of sarcoidosis granulomas in lung and lymph node tissues (most commonly affected organs) and comparison to infectious granulomas (coccidiodomycosis and tuberculosis). METHODS: Transcriptomic profiles of immune-related gene from micro-dissected sarcoidosis granulomas within lung and mediastinal lymph node tissues and compared to infectious granulomas from paraffin-embedded blocks. Differentially-expressed genes (DEGs) were profiled, compared among the three granulomatous diseases and analyzed for functional enrichment pathways. RESULTS: Despite histologic similarities, DEGs and pathway enrichment markedly differed in sarcoidosis granulomas from lymph nodes and lung. Lymph nodes showed a clear immunological response, whereas a structural regenerative response was observed in lung. Sarcoidosis granuloma gene expression data corroborated previously reported genomic biomarkers (STAB1, HBEGF, and NOTCH4), excluded others and identified new genomic markers present in lung and lymph nodes, ADAMTS1, NPR1 and CXCL2. Comparisons between sarcoidosis and pathogen granulomas identified pathway divergences and commonalities at gene expression level. CONCLUSION: These findings suggest the importance of tissue and disease-specificity evaluation when exploring sarcoidosis genomic markers. This relevant translational information in sarcoidosis and other two histopathological similar infections provides meaningful specific genomic-derived biomarkers for sarcoidosis diagnosis and prognosis.

Plasma metabolomic profiles enhance precision medicine for volunteers of normal health.

Precision medicine, taking account of human individuality in genes, environment, and lifestyle for early disease diagnosis and individualized therapy, has shown great promise to transform medical care. Nontargeted metabolomics, with the ability to detect broad classes of biochemicals, can provide a comprehensive functional phenotype integrating clinical phenotypes with genetic and nongenetic factors. To test the application of metabolomics in individual diagnosis, we conducted a metabolomics analysis on plasma samples collected from 80 volunteers of normal health with complete medical records and three-generation pedigrees. Using a broad-spectrum metabolomics platform consisting of liquid chromatography and GC coupled with MS, we profiled nearly 600 metabolites covering 72 biochemical pathways in all major branches of biosynthesis, catabolism, gut microbiome activities, and xenobiotics. Statistical analysis revealed a considerable range of variation and potential metabolic abnormalities across the individuals in this cohort. Examination of the convergence of metabolomics profiles with whole-exon sequences (WESs) provided an effective approach to assess and interpret clinical significance of genetic mutations, as shown in a number of cases, including fructose intolerance, xanthinuria, and carnitine deficiency. Metabolic abnormalities consistent with early indications of diabetes, liver dysfunction, and disruption of gut microbiome homeostasis were identified in several volunteers. Additionally, diverse metabolic responses to medications among the volunteers may assist to identify therapeutic effects and sensitivity to toxicity. The results of this study demonstrate that metabolomics could be an effective approach to complement next generation sequencing (NGS) for disease risk analysis, disease monitoring, and drug management in our goal toward precision care.

Defective Store-Operated Calcium Entry Causes Partial Nephrogenic Diabetes Insipidus.

Store-operated calcium entry (SOCE) is the mechanism by which extracellular signals elicit prolonged intracellular calcium elevation to drive changes in fundamental cellular processes. Here, we investigated the role of SOCE in the regulation of renal water reabsorption, using the inbred rat strain SHR-A3 as an animal model with disrupted SOCE. We found that SHR-A3, but not SHR-B2, have a novel truncating mutation in the gene encoding stromal interaction molecule 1 (STIM1), the endoplasmic reticulum calcium (Ca(2+)) sensor that triggers SOCE. Balance studies revealed increased urine volume, hypertonic plasma, polydipsia, and impaired urinary concentrating ability accompanied by elevated circulating arginine vasopressin (AVP) levels in SHR-A3 compared with SHR-B2. Isolated, split-open collecting ducts (CD) from SHR-A3 displayed decreased basal intracellular Ca(2+) levels and a major defect in SOCE. Consequently, AVP failed to induce the sustained intracellular Ca(2+) mobilization that requires SOCE in CD cells from SHR-A3. This effect decreased the abundance of aquaporin 2 and enhanced its intracellular retention, suggesting impaired sensitivity of the CD to AVP in SHR-A3. Stim1 knockdown in cultured mpkCCDc14 cells reduced SOCE and basal intracellular Ca(2+) levels and prevented AVP-induced translocation of aquaporin 2, further suggesting the effects in SHR-A3 result from the expression of truncated STIM1. Overall, these results identify a novel mechanism of nephrogenic diabetes insipidus and uncover a role of SOCE in renal water handling.

Adult genetic risk screening.

Recent advances in genetic analysis especially DNA sequencing technology open a new strategy for adult disease prevention by genetic screening. Physicians presently treat disease pathology with less emphasis on disease risk prevention/reduction. Genetic screening has reduced the incidence of untreatable childhood genetic diseases and improved the care of newborns. The opportunity exists to expand screening programs and reduce the incidence of adult onset diseases via genetic risk identification and disease intervention. This article outlines the approach, challenges, and benefits of such screening for adult genetic disease risks.

Personalized genomic disease risk of volunteers.

Next-generation sequencing (NGS) is commonly used for researching the causes of genetic disorders. However, its usefulness in clinical practice for medical diagnosis is in early development. In this report, we demonstrate the value of NGS for genetic risk assessment and evaluate the limitations and barriers for the adoption of this technology into medical practice. We performed whole exome sequencing (WES) on 81 volunteers, and for each volunteer, we requested personal medical histories, constructed a three-generation pedigree, and required their participation in a comprehensive educational program. We limited our clinical reporting to disease risks based on only rare damaging mutations and known pathogenic variations in genes previously reported to be associated with human disorders. We identified 271 recessive risk alleles (214 genes), 126 dominant risk alleles (101 genes), and 3 X-recessive risk alleles (3 genes). We linked personal disease histories with causative disease genes in 18 volunteers. Furthermore, by incorporating family histories into our genetic analyses, we identified an additional five heritable diseases. Traditional genetic counseling and disease education were provided in verbal and written reports to all volunteers. Our report demonstrates that when genome results are carefully interpreted and integrated with an individual's medical records and pedigree data, NGS is a valuable diagnostic tool for genetic disease risk.

Lymphatic abnormalities are associated with RASA1 gene mutations in mouse and man.

Mutations in gene RASA1 have been historically associated with capillary malformation-arteriovenous malformation, but sporadic reports of lymphatic involvement have yet to be investigated in detail. To investigate the impact of RASA1 mutations in the lymphatic system, we performed investigational near-infrared fluorescence lymphatic imaging and confirmatory radiographic lymphangiography in a Parkes-Weber syndrome (PKWS) patient with suspected RASA1 mutations and correlated the lymphatic abnormalities against that imaged in an inducible Rasa1 knockout mouse. Whole-exome sequencing (WES) analysis and validation by Sanger sequencing of DNA from the patient and unaffected biological parents enabled us to identify an early-frameshift deletion in RASA1 that was shared with the father, who possessed a capillary stain but otherwise no overt disease phenotype. Abnormal lymphatic vasculature was imaged in both affected and unaffected legs of the PKWS subject that transported injected indocyanine green dye to the inguinal lymph node and drained atypically into the abdomen and into dermal lymphocele-like vesicles on the groin. Dermal lymphatic hyperplasia and dilated vessels were observed in Rasa1-deficient mice, with subsequent development of chylous ascites. WES analyses did not identify potential gene modifiers that could explain the variability of penetrance between father and son. Nonetheless, we conclude that the RASA1 mutation is responsible for the aberrant lymphatic architecture and functional abnormalities, as visualized in the PKWS subject and in the animal model. Our unique method to combine investigatory near-infrared fluorescence lymphatic imaging and WES for accurate phenoptyping and unbiased genotyping allows the study of molecular mechanisms of lymphatic involvement of hemovascular disorders.

Somatic mutations affect key pathways in lung adenocarcinoma.

Determining the genetic basis of cancer requires comprehensive analyses of large collections of histopathologically well-classified primary tumours. Here we report the results of a collaborative study to discover somatic mutations in 188 human lung adenocarcinomas. DNA sequencing of 623 genes with known or potential relationships to cancer revealed more than 1,000 somatic mutations across the samples. Our analysis identified 26 genes that are mutated at significantly high frequencies and thus are probably involved in carcinogenesis. The frequently mutated genes include tyrosine kinases, among them the EGFR homologue ERBB4; multiple ephrin receptor genes, notably EPHA3; vascular endothelial growth factor receptor KDR; and NTRK genes. These data provide evidence of somatic mutations in primary lung adenocarcinoma for several tumour suppressor genes involved in other cancers--including NF1, APC, RB1 and ATM--and for sequence changes in PTPRD as well as the frequently deleted gene LRP1B. The observed mutational profiles correlate with clinical features, smoking status and DNA repair defects. These results are reinforced by data integration including single nucleotide polymorphism array and gene expression array. Our findings shed further light on several important signalling pathways involved in lung adenocarcinoma, and suggest new molecular targets for treatment.

Genome-wide association study of a lipedema phenotype among women in the UK Biobank identifies multiple genetic risk factors.

Lipedema is a common disorder characterized by excessive deposition of subcutaneous adipose tissue (SAT) in the legs, hips, and buttocks, mainly occurring in adult women. Although it appears to be heritable, no specific genes have yet been identified. To identify potential genetic risk factors for lipedema, we used bioelectrical impedance analysis and anthropometric data from the UK Biobank to identify women with and without a lipedema phenotype. Specifically, we identified women with both a high percentage of fat in the lower limbs and a relatively small waist, adjusting for hip circumference. We performed a genome-wide association study (GWAS) for this phenotype, and performed multiple sensitivity GWAS. In an independent case/control study of lipedema based on strict clinical criteria, we attempted to replicate our top hits. We identified 18 significant loci (p < 5 × 10-9), several of which have previously been identified in GWAS of waist-to-hip ratio with larger effects in women. Two loci (VEGFA and GRB14-COBLL1) were significantly associated with lipedema in the independent replication study. Follow-up analyses suggest an enrichment of genes expressed in blood vessels and adipose tissue, among other tissues. Our findings provide a starting point towards better understanding the genetic and physiological basis of lipedema.

Examination of eQTL Polymorphisms Associated with Increased Risk of Progressive Complicated Sarcoidosis in European and African Descent Subjects.

Background: A limited pool of SNPs are linked to the development and severity of sarcoidosis, a systemic granulomatous inflammatory disease. By integrating genome-wide association studies (GWAS) data and expression quantitative trait loci (eQTL) single nuclear polymorphisms (SNPs), we aimed to identify novel sarcoidosis SNPs potentially influencing the development of complicated sarcoidosis. Methods: A GWAS (Affymetrix 6.0) involving 209 African-American (AA) and 193 European-American (EA, 75 and 51 complicated cases respectively) and publicly-available GWAS controls (GAIN) was utilized. Annotation of multi-tissue eQTL SNPs present on the GWAS created a pool of ~46,000 eQTL SNPs examined for association with sarcoidosis risk and severity (Logistic Model, Plink). The most significant EA/AA eQTL SNPs were genotyped in a sarcoidosis validation cohort (n=1034) and cross-validated in two independent GWAS cohorts. Results: No single GWAS SNP achieved significance (p<1x10-8), however, analysis of the eQTL/GWAS SNP pool yielded 621 eQTL SNPs (p<10-4) associated with 730 genes that highlighted innate immunity, MHC Class II, and allograft rejection pathways with multiple SNPs validated in an independent sarcoidosis cohort (105 SNPs analyzed) (NOTCH4, IL27RA, BTNL2, ANXA11, HLA-DRB1). These studies confirm significant association of eQTL/GWAS SNPs in EAs and AAs with sarcoidosis risk and severity (complicated sarcoidosis) involving HLA region and innate immunity. Conclusion: Despite the challenge of deciphering the genetic basis for sarcoidosis risk/severity, these results suggest that integrated eQTL/GWAS approaches may identify novel variants/genes and support the contribution of dysregulated innate immune responses to sarcoidosis severity.

The finished DNA sequence of human chromosome 12.

Human chromosome 12 contains more than 1,400 coding genes and 487 loci that have been directly implicated in human disease. The q arm of chromosome 12 contains one of the largest blocks of linkage disequilibrium found in the human genome. Here we present the finished sequence of human chromosome 12, which has been finished to high quality and spans approximately 132 megabases, representing approximately 4.5% of the human genome. Alignment of the human chromosome 12 sequence across vertebrates reveals the origin of individual segments in chicken, and a unique history of rearrangement through rodent and primate lineages. The rate of base substitutions in recent evolutionary history shows an overall slowing in hominids compared with primates and rodents.

MicroRNA and protein-coding gene expression analysis in idiopathic pulmonary fibrosis yields novel biomarker signatures associated to survival.

Idiopathic pulmonary fibrosis (IPF) is a chronic, progressive interstitial lung disease of unknown etiology that poses significant challenges in early diagnosis and prediction of progression. Analyses of microRNA and gene expression in IPF have yielded potentially predictive information. However, the relationship between microRNA/gene expression and quantitative phenotypic value in IPF remains controversial, as is the added value of this approach to current molecular signatures in IPF. To identify biomarkers predictive of survival in IPF via a microRNA-driven strategy. We profiled microRNA and protein-coding gene expression in peripheral blood mononuclear cells from 70 IPF subjects in a discovery cohort. We linked the microRNA/gene expression level with the quantitative phenotypic variation in IPF, including diffusing capacity of the lung for carbon monoxide and the forced vital capacity percent predicted. In silico analyses of expression profiles and quantitative phenotypic data allowed the generation of 2 sets of IPF molecular signatures (unique for microRNAs and protein-coding genes) that predict IPF survival. Each signature performed well in a validation cohort comprised of IPF patients aggregated from distinct patient populations recruited from different sites. Resampling test suggests that the protein-coding gene based signature is comparable and potentially superior to published IPF prognostic gene signatures. In conclusion, these results highlight the utility of microRNA-driven peripheral blood molecular signatures as valuable and novel biomarkers associated to individuals at high survival risk and for potentially facilitating individualized therapies in this enigmatic disorder.

Low Dose Carbon Monoxide Exposure in Idiopathic Pulmonary Fibrosis Produces a CO Signature Comprised of Oxidative Phosphorylation Genes.

Compelling preclinical studies indicate that low-dose carbon monoxide (CO) abrogates experimental lung fibrosis. We recently reported the results of a multicenter, double-blinded, clinical trial of inhaled CO in patients with idiopathic pulmonary fibrosis (IPF). Identifying no significantly changes in metalloproteinase-7 (MMP7) serum concentration, or secondary endpoints of physiologic measurements, hospitalization, death, or patient-reported outcomes. In the present study, we evaluated the effect of low dose CO exposure (100-200 ppm) for 12 weeks on genome-wide gene expression in peripheral blood mononuclear cells (PBMC) derived from these IPF study subjects. We conducted transcriptome profiling on 38 IPF subjects with time points available at 0, 12, and 24 weeks. Total RNA isolated from PBMCs was hybridized onto the Affymetrix Human Gene 2.0 ST Array. We identified 621 genes significantly upregulated in the 24-week CO exposed group compared with the 12-week. Pathway analysis demonstrated association with Oxidative Phosphorylation (adjusted P < 0.05). We identified a clear CO signature dominated with 23 oxidative phosphorylation-related genes (FDR <10%). We confirmed the expression of nine selected gene products using Nanostring's nCounter analysis system. These findings suggest this signature may serve as a potential genomic biomarker for CO exposure and for potential titration of dosage to allow precision testing of therapies in future low dose CO therapeutic studies in IPF.

Whole Exome Sequencing Reveals a Mutation in CRYBB2 in a Large Mexican Family with Autosomal Dominant Pulverulent Cataract.

Congenital cataract, an important cause of reversible blindness, is due to several causes including Mendelian inheritance. Thirty percent of cataracts are hereditary with participation of the gamma crystallin genes. Clinical and genetic heterogeneity is observed in patients with gene mutations and congenital cataract; about 40 genetic loci have been associated with hereditary cataract. In this study, we identified the underlying genetic cause of an autosomal dominant pulverulent cataract (ADPC) in a large Mexican family. Twenty-one affected patients and 20 healthy members of a family with ADPC were included. Genomic DNA was analyzed by whole exome sequencing in the proband, a normal daughter, and in an affected son, whereas DNA Sanger sequencing was performed in all members of the family. After the bioinformatics analysis, all samples were genotyped using Sanger sequencing to eliminate variants that do not cosegregate with the cataract. We observed a perfect cosegregation of a nonsense mutation c.475C>T (p.Q155*) in exon 6 of the CRYBB2 gene with ADPC. We calculated a logarithm of the odds score of 5.5. This mutation was not detected in healthy members of the family and in 100 normal controls. This is the first Mexican family with ADPC associated with a p.Q155* mutation. Interestingly, this specific mutation in the CRYBB2 gene seems to be exclusively associated with pulverulent/cerulean cataract (with some clinical variability) independent of the population's genetic background.

Targeted correction and restored function of the CFTR gene in cystic fibrosis induced pluripotent stem cells.

Recently developed reprogramming and genome editing technologies make possible the derivation of corrected patient-specific pluripotent stem cell sources-potentially useful for the development of new therapeutic approaches. Starting with skin fibroblasts from patients diagnosed with cystic fibrosis, we derived and characterized induced pluripotent stem cell (iPSC) lines. We then utilized zinc-finger nucleases (ZFNs), designed to target the endogenous CFTR gene, to mediate correction of the inherited genetic mutation in these patient-derived lines via homology-directed repair (HDR). We observed an exquisitely sensitive, homology-dependent preference for targeting one CFTR allele versus the other. The corrected cystic fibrosis iPSCs, when induced to differentiate in vitro, expressed the corrected CFTR gene; importantly, CFTR correction resulted in restored expression of the mature CFTR glycoprotein and restoration of CFTR chloride channel function in iPSC-derived epithelial cells.

The road from next-generation sequencing to personalized medicine.

Moving from a traditional medical model of treating pathologies to an individualized predictive and preventive model of personalized medicine promises to reduce the healthcare cost on an overburdened and overwhelmed system. Next-generation sequencing (NGS) has the potential to accelerate the early detection of disorders and the identification of pharmacogenetics markers to customize treatments. This review explains the historical facts that led to the development of NGS along with the strengths and weakness of NGS, with a special emphasis on the analytical aspects used to process NGS data. There are solutions to all the steps necessary for performing NGS in the clinical context where the majority of them are very efficient, but there are some crucial steps in the process that need immediate attention.

Hypertensive renal injury is associated with gene variation affecting immune signaling.

BACKGROUND: The spontaneously hypertensive rat (SHR) strain exists in lines that contrast strongly in susceptibility to renal injury in hypertension. These inbred lines share common ancestry, and only 13% of their genomes arise from different ancestors. METHODS AND RESULTS: We used next gen sequencing to detect natural allelic variation in 5 genes of the immunoreceptor signaling pathway (IgH, Dok3, Src, Syk, and JunD) that arise from different ancestors in the injury-prone SHR-A3 and the resistant SHR-B2 lines. We created an intercross between these lines, and in the F2 progeny, we observed that the inheritance of haplotype blocks containing the SHR-A3 alleles of these 5 genes correlated with increased albuminuria and histological measures of renal injury. To test whether accumulated genetic variation in this pathway may create a therapeutic target in hypertensive renal injury, rats of both lines were treated with the immunosuppressant mycophenolate mofetil (MMF). MMF reduced proteinuria (albumin to creatinine ratio) from 6.6 to 1.2 mg/mg (P<0.001) in SHR-A3. Glomerular injury scores were reduced in MMF-treated SHR-A3 from 1.6 to 1.4 (P<0.002). Tubulo-interstitial injury was reduced in MMF-treated SHR-A3 from 2.62 to 2.0 (P=0.001). MMF treatment also reduced renal fibrosis in SHR-A3 (3.9 versus 2.0; P<0.001). CONCLUSIONS: Polygenic susceptibility to renal injury in hypertension arises in association with genetic variation in genes that participate in immune responses and is dramatically improved by reduction of immune system activity.

Evidence for SH2 domain-containing 5'-inositol phosphatase-2 (SHIP2) contributing to a lymphatic dysfunction.

The lymphatic vasculature plays a critical role in a number of disease conditions of increasing prevalence, such as autoimmune disorders, obesity, blood vascular diseases, and cancer metastases. Yet, unlike the blood vasculature, the tools available to interrogate the molecular basis of lymphatic dysfunction/disease have been lacking. More recently, investigators have reported that dysregulation of the PI3K pathway is involved in syndromic human diseases that involve abnormal lymphatic vasculatures, but there have been few compelling results that show the direct association of this molecular pathway with lymphatic dysfunction in humans. Using near-infrared fluorescence lymphatic imaging (NIRFLI) to phenotype and next generation sequencing (NGS) for unbiased genetic discovery in a family with non-syndromic lymphatic disease, we discovered a rare, novel mutation in INPPL1 that encodes the protein SHIP2, which is a negative regulator of the PI3K pathway, to be associated with lymphatic dysfunction in the family. In vitro interrogation shows that SHIP2 is directly associated with impairment of normal lymphatic endothelial cell (LEC) behavior and that SHIP2 associates with receptors that are associated in lymphedema, implicating its direct involvement in the lymphatic vasculature.

Truncating mutations of MAGEL2 cause Prader-Willi phenotypes and autism.

Prader-Willi syndrome (PWS) is caused by the absence of paternally expressed, maternally silenced genes at 15q11-q13. We report four individuals with truncating mutations on the paternal allele of MAGEL2, a gene within the PWS domain. The first subject was ascertained by whole-genome sequencing analysis for PWS features. Three additional subjects were identified by reviewing the results of exome sequencing of 1,248 cases in a clinical laboratory. All four subjects had autism spectrum disorder (ASD), intellectual disability and a varying degree of clinical and behavioral features of PWS. These findings suggest that MAGEL2 is a new gene causing complex ASD and that MAGEL2 loss of function can contribute to several aspects of the PWS phenotype.

Paclitaxel-dependent mutants have severely reduced microtubule assembly and reduced tubulin synthesis.

A subset of mutant cell lines selected for resistance to the antitumor drug paclitaxel are unable to progress normally through mitosis unless the drug is present in the growth medium. Without paclitaxel the cells form defective spindles, undergo aberrant mitoses, fail to complete cell division and eventually die. Analysis of these drug-dependent cells revealed a low amount of microtubule polymer and less tubulin production than wild-type cells. Ribonuclease protection experiments indicated that the decreased tubulin protein was due to decreased tubulin mRNA. Enhancing microtubule assembly by treating the cells with paclitaxel, restored tubulin to levels comparable with those of paclitaxel-treated wild-type cells, which demonstrated that the drug-dependent cells do not have a permanent impairment in their capacity to synthesize tubulin. Paclitaxel-resistant (but not dependent) cells have a smaller reduction in microtubule polymer with little or no decrease in tubulin production, whereas colcemidresistant cells have increased microtubule assembly but also exhibit little or no change in tubulin production. Finally, a mutant cell line producing an unstable beta-tubulin protein has normal growth as well as normal synthesis and polymerization of tubulin, despite an approximately 30% decrease in steady state tubulin content. These studies establish a lower limit of tubulin assembly needed for cell survival and indicate that tubulin assembly must fall below this point to trigger a significant decrease in tubulin synthesis.

Prevalence in two Mexican cities of human papilomavirus DNA sequences in cervical cancer

El cáncer cérvico uterino (CaCU) ocupa el 30 por ciento de todas las neoplasias de la mujer en México, resultando una de las principales causas de muerte por cáncer. En un estudio anterior en la ciudad de México, nuestro grupo reportó que sólo cinco de 16 muestras analizadas por transferencia Southern contenían secuencias virales integradas del papilomavirus humano tipo 16 (HPV16). En el presente trabajo, hemos extendido esta observación en un estudio comparativo de las ciudades de México, D!F!, y Monterrey, Nuevo León, incluyendo además otro tipo viral, el papilomavirus humano tipo 18 (HPV18), para determinar la frecuencia con la que se presentan en población mexicana los dos tipos de papilomavirus humanos más comunes en CaCU. Se determinó en primera instancia que la prevalencia de secuencias del HPV16 era similar en ambas poblaciones (4 de 14 para la ciudad de México y 6 de 23 para Monterrey); sin embargo, la presencia de secuencias de HPV18 fue aún más baja (una de 14 y una de 10, para las ciudades de México y Monterrey, respectivamente). En todos los tumores, las secuencias virales se encontraron integradas al genoma celular. Nuestros resultados muestran que existen una relativamente baja proporción de HPV16 Y 18 en tumores cervicales de población mexicana, sugiriendo que otros tipos de papilomavirus humanos o la presencia de un nuevo factor de riesgo (p.ej., activación de oncogenes) están involucrados en el desarrollo del CaCU en México.