Genomic context sensitizes regulatory elements to genetic disruption.

Genomic context critically modulates regulatory function but is difficult to manipulate systematically. The murine insulin-like growth factor 2 (Igf2)/H19 locus is a paradigmatic model of enhancer selectivity, whereby CTCF occupancy at an imprinting control region directs downstream enhancers to activate either H19 or Igf2. We used synthetic regulatory genomics to repeatedly replace the native locus with 157-kb payloads, and we systematically dissected its architecture. Enhancer deletion and ectopic delivery revealed previously uncharacterized long-range regulatory dependencies at the native locus. Exchanging the H19 enhancer cluster with the Sox2 locus control region (LCR) showed that the H19 enhancers relied on their native surroundings while the Sox2 LCR functioned autonomously. Analysis of regulatory DNA actuation across cell types revealed that these enhancer clusters typify broader classes of context sensitivity genome wide. These results show that unexpected dependencies influence even well-studied loci, and our approach permits large-scale manipulation of complete loci to investigate the relationship between regulatory architecture and function.

Packaging and containerization of computational methods.

Methods for analyzing the full complement of a biomolecule type, e.g., proteomics or metabolomics, generate large amounts of complex data. The software tools used to analyze omics data have reshaped the landscape of modern biology and become an essential component of biomedical research. These tools are themselves quite complex and often require the installation of other supporting software, libraries and/or databases. A researcher may also be using multiple different tools that require different versions of the same supporting materials. The increasing dependence of biomedical scientists on these powerful tools creates a need for easier installation and greater usability. Packaging and containerization are different approaches to satisfy this need by delivering omics tools already wrapped in additional software that makes the tools easier to install and use. In this systematic review, we describe and compare the features of prominent packaging and containerization platforms. We outline the challenges, advantages and limitations of each approach and some of the most widely used platforms from the perspectives of users, software developers and system administrators. We also propose principles to make the distribution of omics software more sustainable and robust to increase the reproducibility of biomedical and life science research.

Genomic context sensitizes regulatory elements to genetic disruption.

Enhancer function is frequently investigated piecemeal using truncated reporter assays or single deletion analysis. Thus it remains unclear to what extent enhancer function at native loci relies on surrounding genomic context. Using the Big-IN technology for targeted integration of large DNAs, we analyzed the regulatory architecture of the murine Igf2/H19 locus, a paradigmatic model of enhancer selectivity. We assembled payloads containing a 157-kb functional Igf2/H19 locus and engineered mutations to genetically direct CTCF occupancy at the imprinting control region (ICR) that switches the target gene of the H19 enhancer cluster. Contrasting activity of payloads delivered at the endogenous Igf2/H19 locus or ectopically at Hprt revealed that the Igf2/H19 locus includes additional, previously unknown long-range regulatory elements. Exchanging components of the Igf2/H19 locus with the well-studied Sox2 locus showed that the H19 enhancer cluster functioned poorly out of context, and required its native surroundings to activate Sox2 expression. Conversely, the Sox2 locus control region (LCR) could activate both Igf2 and H19 outside its native context, but its activity was only partially modulated by CTCF occupancy at the ICR. Analysis of regulatory DNA actuation across different cell types revealed that, while the H19 enhancers are tightly coordinated within their native locus, the Sox2 LCR acts more independently. We show that these enhancer clusters typify broader classes of loci genome-wide. Our results show that unexpected dependencies may influence even the most studied functional elements, and our synthetic regulatory genomics approach permits large-scale manipulation of complete loci to investigate the relationship between locus architecture and function.

Mouse genome rewriting and tailoring of three important disease loci.

Genetically engineered mouse models (GEMMs) help us to understand human pathologies and develop new therapies, yet faithfully recapitulating human diseases in mice is challenging. Advances in genomics have highlighted the importance of non-coding regulatory genome sequences, which control spatiotemporal gene expression patterns and splicing in many human diseases1,2. Including regulatory extensive genomic regions, which requires large-scale genome engineering, should enhance the quality of disease modelling. Existing methods set limits on the size and efficiency of DNA delivery, hampering the routine creation of highly informative models that we call genomically rewritten and tailored GEMMs (GREAT-GEMMs). Here we describe 'mammalian switching antibiotic resistance markers progressively for integration' (mSwAP-In), a method for efficient genome rewriting in mouse embryonic stem cells. We demonstrate the use of mSwAP-In for iterative genome rewriting of up to 115 kb of a tailored Trp53 locus, as well as for humanization of mice using 116 kb and 180 kb human ACE2 loci. The ACE2 model recapitulated human ACE2 expression patterns and splicing, and notably, presented milder symptoms when challenged with SARS-CoV-2 compared with the existing K18-hACE2 model, thus representing a more human-like model of infection. Finally, we demonstrated serial genome writing by humanizing mouse Tmprss2 biallelically in the ACE2 GREAT-GEMM, highlighting the versatility of mSwAP-In in genome writing.

Synthetic regulatory genomics uncovers enhancer context dependence at the Sox2 locus.

Sox2 expression in mouse embryonic stem cells (mESCs) depends on a distal cluster of DNase I hypersensitive sites (DHSs), but their individual contributions and degree of interdependence remain a mystery. We analyzed the endogenous Sox2 locus using Big-IN to scarlessly integrate large DNA payloads incorporating deletions, rearrangements, and inversions affecting single or multiple DHSs, as well as surgical alterations to transcription factor (TF) recognition sequences. Multiple mESC clones were derived for each payload, sequence-verified, and analyzed for Sox2 expression. We found that two DHSs comprising a handful of key TF recognition sequences were each sufficient for long-range activation of Sox2 expression. By contrast, three nearby DHSs were entirely context dependent, showing no activity alone but dramatically augmenting the activity of the autonomous DHSs. Our results highlight the role of context in modulating genomic regulatory element function, and our synthetic regulatory genomics approach provides a roadmap for the dissection of other genomic loci.

Triple-Negative Breast Cancer circRNAome Reveals Hsa_circ_0072309 as a Potential Risk Biomarker.

Circular RNAs (circRNAs) are a class of long non-coding RNAs that have the ability to sponge RNA-Binding Proteins (RBPs). Triple-negative breast cancer (TNBC) has very aggressive behavior and poor prognosis for the patient. Here, we aimed to characterize the global expression profile of circRNAs in TNBC, in order to identify potential risk biomarkers. For that, we obtained RNA-Seq data from TNBC and control samples and performed validation experiments using FFPE and frozen tissues of TNBC patients and controls, followed by in silico analyses to explore circRNA-RBP interactions. We found 16 differentially expressed circRNAs between TNBC patients and controls. Next, we mapped the RBPs that interact with the top five downregulated circRNAs (hsa_circ_0072309, circ_0004365, circ_0006677, circ_0008599, and circ_0009043) and hsa_circ_0000479, resulting in a total of 16 RBPs, most of them being enriched to pathways related to cancer and gene regulation (e.g., AGO1/2, EIF4A3, ELAVL1, and PTBP1). Among the six circRNAs, hsa_circ_0072309 was the one that presented the most confidence results, being able to distinguish TNBC patients from controls with an AUC of 0.78 and 0.81, respectively. This circRNA may be interacting with some RBPs involved in important cancer-related pathways and is a novel potential risk biomarker of TNBC.

Identification of Genomic Variants Associated with the Risk of Acute Lymphoblastic Leukemia in Native Americans from Brazilian Amazonia.

A number of genomic variants related to native American ancestry may be associated with an increased risk of developing Acute Lymphoblastic Leukemia (ALL), which means that Latin American and hispanic populations from the New World may be relatively susceptible to this disease. However, there has not yet been any comprehensive investigation of the variants associated with susceptibility to ALL in traditional Amerindian populations from Brazilian Amazonia. We investigated the exomes of the 18 principal genes associated with susceptibility to ALL in samples of 64 Amerindians from this region, including cancer-free individuals and patients with ALL. We compared the findings with the data on populations representing five continents available in the 1000 Genomes database. The variation in the allele frequencies found between the different groups was evaluated using Fisher's exact test. The analyses of the exomes of the Brazilian Amerindians identified 125 variants, seven of which were new. The comparison of the allele frequencies between the two Amerindian groups analyzed in the present study (ALL patients vs. cancer-free individuals) identified six variants (rs11515, rs2765997, rs1053454, rs8068981, rs3764342, and rs2304465) that may be associated with susceptibility to ALL. These findings contribute to the identification of genetic variants that represent a potential risk for ALL in Amazonian Amerindian populations and might favor precision oncology measures.

Genomic context sensitivity of insulator function.

The specificity of interactions between genomic regulatory elements and potential target genes is influenced by the binding of insulator proteins such as CTCF, which can act as potent enhancer blockers when interposed between an enhancer and a promoter in a reporter assay. But not all CTCF sites genome-wide function as insulator elements, depending on cellular and genomic context. To dissect the influence of genomic context on enhancer blocker activity, we integrated reporter constructs with promoter-only, promoter and enhancer, and enhancer blocker configurations at hundreds of thousands of genomic sites using the Sleeping Beauty transposase. Deconvolution of reporter activity by genomic position reveals distinct expression patterns subject to genomic context, including a compartment of enhancer blocker reporter integrations with robust expression. The high density of integration sites permits quantitative delineation of characteristic genomic context sensitivity profiles and their decomposition into sensitivity to both local and distant DNase I hypersensitive sites. Furthermore, using a single-cell expression approach to test the effect of integrated reporters for differential expression of nearby endogenous genes reveals that CTCF insulator elements do not completely abrogate reporter effects on endogenous gene expression. Collectively, our results lend new insight into genomic regulatory compartmentalization and its influence on the determinants of promoter-enhancer specificity.

Comprehensive analysis of germline mutations in northern Brazil: a panel of 16 genes for hereditary cancer-predisposing syndrome investigation.

BACKGROUND: Next generation sequencing (NGS) has been a handy tool in clinical practice, mainly due to its efficiency and cost-effectiveness. It has been widely used in genetic diagnosis of several inherited diseases, and, in clinical oncology, it may enhance the discovery of new susceptibility genes and enable individualized care of cancer patients. In this context, we explored a pan-cancer panel in the investigation of germline variants in Brazilian patients presenting clinical criteria for hereditary cancer syndromes or familial history. METHODS: Seventy-one individuals diagnosed or with familial history of hereditary cancer syndromes were submitted to custom pan-cancer panel including 16 high and moderate penetrance genes previously associated with hereditary cancer syndromes (APC, BRCA1, BRCA2, CDH1, CDKN2A, CHEK2, MSH2, MSH6, MUTYH, PTEN, RB1, RET, TP53, VHL, XPA and XPC). All pathogenic variants were validated by Sanger sequencing. RESULTS: We identified a total of eight pathogenic variants among 12 of 71 individuals (16.9%). Among the mutation-positive subjects, 50% were diagnosed with breast cancer and had mutations in BRCA1, CDH1 and MUTYH. Notably, 33.3% were individuals diagnosed with polyposis or who had family cases and harbored pathogenic mutations in APC and MUTYH. The remaining individuals (16.7%) were gastric cancer patients with pathogenic variants in CDH1 and MSH2. Overall, 54 (76.05%) individuals presented at least one variant uncertain significance (VUS), totalizing 81 VUS. Of these, seven were predicted to have disease-causing potential. CONCLUSION: Overall, analysis of all these genes in NGS-panel allowed the identification not only of pathogenic variants related to hereditary cancer syndromes but also of some VUS that need further clinical and molecular investigations. The results obtained in this study had a significant impact on patients and their relatives since it allowed genetic counselling and personalized management decisions.

Systematic analysis of binding of transcription factors to noncoding variants.

Many sequence variants have been linked to complex human traits and diseases1, but deciphering their biological functions remains challenging, as most of them reside in noncoding DNA. Here we have systematically assessed the binding of 270 human transcription factors to 95,886 noncoding variants in the human genome using an ultra-high-throughput multiplex protein-DNA binding assay, termed single-nucleotide polymorphism evaluation by systematic evolution of ligands by exponential enrichment (SNP-SELEX). The resulting 828 million measurements of transcription factor-DNA interactions enable estimation of the relative affinity of these transcription factors to each variant in vitro and evaluation of the current methods to predict the effects of noncoding variants on transcription factor binding. We show that the position weight matrices of most transcription factors lack sufficient predictive power, whereas the support vector machine combined with the gapped k-mer representation show much improved performance, when assessed on results from independent SNP-SELEX experiments involving a new set of 61,020 sequence variants. We report highly predictive models for 94 human transcription factors and demonstrate their utility in genome-wide association studies and understanding of the molecular pathways involved in diverse human traits and diseases.

Exome Sequencing of Native Populations From the Amazon Reveals Patterns on the Peopling of South America.

Studies on the peopling of South America have been limited by the paucity of sequence data from Native Americans, especially from the east part of the Amazon region. Here, we investigate the whole exome variation from 58 Native American individuals (eight different populations) from the Amazon region and draw insights into the peopling of South America. By using the sequence data generated here together with data from the public domain, we confirmed a strong genetic distinction between Andean and Amazonian populations. By testing distinct demographic models, our analysis supports a scenario of South America occupation that involves migrations along the Pacific and Atlantic coasts. Occupation of the southeast part of South America would involve migrations from the north, rather than from the west of the continent.

Skin pigmentation polymorphisms associated with increased risk of melanoma in a case-control sample from southern Brazil.

BACKGROUND: Melanoma is the most aggressive type of skin cancer and is associated with environmental and genetic risk factors. It originates in melanocytes, the pigment-producing cells. Single nucleotide polymorphisms (SNPs) in pigmentation genes have been described in melanoma risk modulation, but knowledge in the field is still limited. METHODS: In a case-control approach (107 cases and 119 controls), we investigated the effect of four pigmentation gene SNPs (TYR rs1126809, HERC2 rs1129038, SLC24A5 rs1426654, and SLC45A2 rs16891982) on melanoma risk in individuals from southern Brazil using a multivariate logistic regression model and multifactor dimensionality reduction (MDR) analysis. RESULTS: Two SNPs were associated with an increased risk of melanoma in a dominant model: rs1129038AA and rs1426654AA [OR = 2.094 (95% CI: 1.106-3.966), P = 2.3 10- 2 and OR = 7.126 (95% CI: 1.873-27.110), P = 4.0 10- 3, respectively]. SNP rs16891982CC was associated with a lower risk to melanoma development in a log-additive model when the allele C was inherited [OR = 0.081 (95% CI: 0.008-0.782), P = 3 10- 2]. In addition, MDR analysis showed that the combination of the rs1426654AA and rs16891982GG genotypes was associated with a higher risk for melanoma (P = 3 10- 3), with a redundant effect. CONCLUSIONS: These results contribute to the current knowledge and indicate that epistatic interaction of these SNPs, with an additive or correlational effect, may be involved in modulating the risk of melanoma in individuals from a geographic region with a high incidence of the disease.

Sequencing identifies multiple early introductions of SARS-CoV-2 to the New York City region.

Effective public response to a pandemic relies upon accurate measurement of the extent and dynamics of an outbreak. Viral genome sequencing has emerged as a powerful approach to link seemingly unrelated cases, and large-scale sequencing surveillance can inform on critical epidemiological parameters. Here, we report the analysis of 864 SARS-CoV-2 sequences from cases in the New York City metropolitan area during the COVID-19 outbreak in spring 2020. The majority of cases had no recent travel history or known exposure, and genetically linked cases were spread throughout the region. Comparison to global viral sequences showed that early transmission was most linked to cases from Europe. Our data are consistent with numerous seeds from multiple sources and a prolonged period of unrecognized community spreading. This work highlights the complementary role of genomic surveillance in addition to traditional epidemiological indicators.

Global Analyses of Expressed Piwi-Interacting RNAs in Gastric Cancer.

Gastric cancer (GC) represents a notable amount of morbidity and mortality worldwide. Understanding the molecular basis of CG will offer insight into its pathogenesis in an attempt to identify new molecular biomarkers to early diagnose this disease. Therefore, studies involving small non-coding RNAs have been widely explored. Among these, PIWI-interacting RNAs (piRNAs) are an emergent class that can play important roles in carcinogenesis. In this study, small-RNA sequencing was used to identify the global piRNAs expression profile (piRNome) of gastric cancer patients. We found 698 piRNAs in gastric tissues, 14 of which were differentially expressed (DE) between gastric cancer (GC), adjacent to gastric cancer (ADJ), and non-cancer tissues (NC). Moreover, three of these DE piRNAs (piR-48966*, piR-49145, piR-31335*) were differently expressed in both GC and ADJ samples in comparison to NC samples, indicating that the tumor-adjacent tissue was molecularly altered and should not be considered as a normal control. These three piRNAs are potential risk biomarkers for GC, especially piR-48966* and piR-31335*. Furthermore, an in-silico search for mRNAs targeted by the differentially expressed piRNAs revealed that these piRNAs may regulate genes that participate in cancer-related pathways, suggesting that these small non-coding RNAs may be directly and indirectly involved in gastric carcinogenesis.

Sequencing identifies multiple early introductions of SARS-CoV-2 to the New York City Region.

Effective public response to a pandemic relies upon accurate measurement of the extent and dynamics of an outbreak. Viral genome sequencing has emerged as a powerful approach to link seemingly unrelated cases, and large-scale sequencing surveillance can inform on critical epidemiological parameters. Here, we report the analysis of 864 SARS-CoV-2 sequences from cases in the New York City metropolitan area during the COVID-19 outbreak in Spring 2020. The majority of cases had no recent travel history or known exposure, and genetically linked cases were spread throughout the region. Comparison to global viral sequences showed that early transmission was most linked to cases from Europe. Our data are consistent with numerous seeds from multiple sources and a prolonged period of unrecognized community spreading. This work highlights the complementary role of genomic surveillance in addition to traditional epidemiological indicators.

miRNome Reveals New Insights Into the Molecular Biology of Field Cancerization in Gastric Cancer.

BACKGROUND: MicroRNAs (miRNAs) play an important role in gastric carcinogenesis and have been associated with gastric field cancerization; however, their role is not fully understood in this process. We performed the miRNome sequencing of non-cancerous, adjacent to tumor and gastric cancer samples to understand the involvement of these small RNAs in gastric field cancerization. METHODS: We analyzed samples of patients without cancer as control (non-cancerous gastric samples) and adjacent to cancer and gastric cancer paired samples, and considered miRNAs with |log2(fold change)| > 2 and Padj < 0.05 to be statistically significant. The identification of target genes, functional analysis and enrichment in KEGG pathways were realized in the TargetCompare, miRTargetLink, and DAVID tools. We also performed receiver operating characteristic (ROC) curves and miRNAs that had an AUC > 0.85 were considered to be potential biomarkers. RESULTS: We found 14 miRNAs exclusively deregulated in gastric cancer, of which six have potential diagnostic value for advanced disease. Nine miRNAs with known tumor suppressor activities (TS-miRs) were deregulated exclusively in adjacent tissue. Of these, five have potential diagnostic value for the early stages of gastric cancer. Functional analysis of these TS-miRs revealed that they regulate important cellular signaling pathways (PI3K-Akt, HIF-1, Ras, Rap1, ErbB, and MAPK signaling pathways), that are involved in gastric carcinogenesis. Seven miRNAs were differentially expressed in both gastric cancer and adjacent regarding to non-cancerous tissues; among them, hsa-miR-200a-3p and hsa-miR-873-5p have potential diagnostic value for early and advanced stages of the disease. Only hsa-miR-196a-5p was differentially expressed between adjacent to cancer and gastric cancer tissues. In addition, the other miRNAs identified in this study were not differentially expressed between adjacent to cancer and gastric cancer, suggesting that these tissues are very similar and that share these molecular changes. CONCLUSION: Our results show that gastric cancer and adjacent tissues have a similar miRNA expression profile, indicating that studied miRNAs are intimately associated with field cancerization in gastric cancer. The overexpression of TS-miRs in adjacent tissues may be a barrier against tumorigenesis within these pre-cancerous conditions prior to the eventual formation or relapse of a tumor. Additionally, these miRNAs have a great accuracy in discriminating non-cancerous from adjacent to tumor and cancer tissues and can be potentially useful as biomarkers for gastric cancer.

Whole Genome Sequencing of the Pirarucu (Arapaima gigas) Supports Independent Emergence of Major Teleost Clades.

The Pirarucu (Arapaima gigas) is one of the world's largest freshwater fishes and member of the superorder Osteoglossomorpha (bonytongues), one of the oldest lineages of ray-finned fishes. This species is an obligate air-breather found in the basin of the Amazon River with an attractive potential for aquaculture. Its phylogenetic position among bony fishes makes the Pirarucu a relevant subject for evolutionary studies of early teleost diversification. Here, we present, for the first time, a draft genome version of the A. gigas genome, providing useful information for further functional and evolutionary studies. The A. gigas genome was assembled with 103-Gb raw reads sequenced in an Illumina platform. The final draft genome assembly was ∼661 Mb, with a contig N50 equal to 51.23 kb and scaffold N50 of 668 kb. Repeat sequences accounted for 21.69% of the whole genome, and a total of 24,655 protein-coding genes were predicted from the genome assembly, with an average of nine exons per gene. Phylogenomic analysis based on 24 fish species supported the postulation that Osteoglossomorpha and Elopomorpha (eels, tarpons, and bonefishes) are sister groups, both forming a sister lineage with respect to Clupeocephala (remaining teleosts). Divergence time estimations suggested that Osteoglossomorpha and Elopomorpha lineages emerged independently in a period of ∼30 Myr in the Jurassic. The draft genome of A. gigas provides a valuable genetic resource for further investigations of evolutionary studies and may also offer a valuable data for economic applications.

The comprehensive expression analysis of circular RNAs in gastric cancer and its association with field cancerization.

Circular RNAs comprise a new class of long noncoding RNAs characterized by their 5' and 3' ends covalently joined. Previous studies have demonstrated that some circular RNAs act as microRNA sponges, and are associated with cellular proliferation in cancer. We were the first to analyze the global expression of circular RNAs in samples of patients without gastric cancer, gastric cancer, and matched tumor-adjacent gastric tissue. Among the samples, we identified 736 previously annotated circular RNAs by RNA-Seq. The tumor-adjacent tissue presented the higher abundance of circular RNAs and could not be considered as a normal tissue, reinforcing the notion of field effect in gastric cancer. We identified five differentially expressed circular RNAs that may be potential biomarkers of this type of cancer. We also predicted candidate microRNAs targets of the highest expressed circular RNAs in gastric tissues and found five miRNAs. Overall, our results support the hypothesis of circular RNAs representing a novel factor in the dynamic epigenetic network of gene regulation, which involves the microRNAs, its mRNAs targets, and the circular RNAs-derived genes. Further studies are needed to elucidate the roles and the functional relevance of the circular RNAs in human diseases.

Whole genome sequencing of an African American family highlights toll like receptor 6 variants in Kawasaki disease susceptibility.

Kawasaki disease (KD) is the most common acquired pediatric heart disease. We analyzed Whole Genome Sequences (WGS) from a 6-member African American family in which KD affected two of four children. We sought rare, potentially causative genotypes by sequentially applying the following WGS filters: sequence quality scores, inheritance model (recessive homozygous and compound heterozygous), predicted deleteriousness, allele frequency, genes in KD-associated pathways or with significant associations in published KD genome-wide association studies (GWAS), and with differential expression in KD blood transcriptomes. Biologically plausible genotypes were identified in twelve variants in six genes in the two affected children. The affected siblings were compound heterozygous for the rare variants p.Leu194Pro and p.Arg247Lys in Toll-like receptor 6 (TLR6), which affect TLR6 signaling. The affected children were also homozygous for three common, linked (r2 = 1) intronic single nucleotide variants (SNVs) in TLR6 (rs56245262, rs56083757 and rs7669329), that have previously shown association with KD in cohorts of European descent. Using transcriptome data from pre-treatment whole blood of KD subjects (n = 146), expression quantitative trait loci (eQTL) analyses were performed. Subjects homozygous for the intronic risk allele (A allele of TLR6 rs56245262) had differential expression of Interleukin-6 (IL-6) as a function of genotype (p = 0.0007) and a higher erythrocyte sedimentation rate at diagnosis. TLR6 plays an important role in pathogen-associated molecular pattern recognition, and sequence variations may affect binding affinities that in turn influence KD susceptibility. This integrative genomic approach illustrates how the analysis of WGS in multiplex families with a complex genetic disease allows examination of both the common disease-common variant and common disease-rare variant hypotheses.

Analysis of 12 variants in the development of gastric and colorectal cancers.

AIM: To evaluate the relation between 12 polymorphisms and the development of gastric cancer (GC) and colorectal cancer (CRC). METHODS: In this study, we included 125 individuals with GC diagnosis, 66 individuals with CRC diagnosis and 475 cancer-free individuals. All participants resided in the North region of Brazil and authorized the use of their samples. The 12 polymorphisms (in CASP8, CYP2E1, CYP19A1, IL1A, IL4, MDM2, NFKB1, PAR1, TP53, TYMS, UGT1A1 and XRCC1 genes) were genotyped in a single PCR for each individual, followed by fragment analysis. To avoid misinterpretation due to population substructure, we applied a previously developed set of 61 ancestry-informative markers that can also be genotyped by multiplex PCR. The statistical analyses were performed in Structure v.2.3.4, R environment and SPSS v.20. RESULTS: After statistical analyses with the control of confounding factors, such as genetic ancestry, three markers (rs79071878 in IL4, rs3730485 in MDM2 and rs28362491 in NFKB1) were positively associated with the development of GC. One of these markers (rs28362491) and the marker in the UGT1A1 gene (rs8175347) were positively associated with the development of CRC. Therefore, we investigated whether the joint presence of the deleterious alleles of each marker could affect the development of cancer and we obtained positive results in all analyses. Carriers of the combination of alleles RP1 + DEL (rs79071878 and rs28361491, respectively) are at 10-times greater risk of developing GC than carriers of other combinations. Similarly, carriers of the combination of DEL + RARE (rs283628 and rs8175347) are at about 12-times greater risk of developing CRC than carriers of other combinations. CONCLUSION: These findings are important for the comprehension of gastric and CRC development, particularly in highly admixed populations, such as the Brazilian population.