Paths and timings of the peopling of Polynesia inferred from genomic networks.

Polynesia was settled in a series of extraordinary voyages across an ocean spanning one third of the Earth1, but the sequences of islands settled remain unknown and their timings disputed. Currently, several centuries separate the dates suggested by different archaeological surveys2-4. Here, using genome-wide data from merely 430 modern individuals from 21 key Pacific island populations and novel ancestry-specific computational analyses, we unravel the detailed genetic history of this vast, dispersed island network. Our reconstruction of the branching Polynesian migration sequence reveals a serial founder expansion, characterized by directional loss of variants, that originated in Samoa and spread first through the Cook Islands (Rarotonga), then to the Society (Totaiete ma) Islands (11th century), the western Austral (Tuha'a Pae) Islands and Tuamotu Archipelago (12th century), and finally to the widely separated, but genetically connected, megalithic statue-building cultures of the Marquesas (Te Henua 'Enana) Islands in the north, Raivavae in the south, and Easter Island (Rapa Nui), the easternmost of the Polynesian islands, settled in approximately AD 1200 via Mangareva.

Native American gene flow into Polynesia predating Easter Island settlement.

The possibility of voyaging contact between prehistoric Polynesian and Native American populations has long intrigued researchers. Proponents have pointed to the existence of New World crops, such as the sweet potato and bottle gourd, in the Polynesian archaeological record, but nowhere else outside the pre-Columbian Americas1-6, while critics have argued that these botanical dispersals need not have been human mediated7. The Norwegian explorer Thor Heyerdahl controversially suggested that prehistoric South American populations had an important role in the settlement of east Polynesia and particularly of Easter Island (Rapa Nui)2. Several limited molecular genetic studies have reached opposing conclusions, and the possibility continues to be as hotly contested today as it was when first suggested8-12. Here we analyse genome-wide variation in individuals from islands across Polynesia for signs of Native American admixture, analysing 807 individuals from 17 island populations and 15 Pacific coast Native American groups. We find conclusive evidence for prehistoric contact of Polynesian individuals with Native American individuals (around AD 1200) contemporaneous with the settlement of remote Oceania13-15. Our analyses suggest strongly that a single contact event occurred in eastern Polynesia, before the settlement of Rapa Nui, between Polynesian individuals and a Native American group most closely related to the indigenous inhabitants of present-day Colombia.

Genomic insights into the origin and diversification of late maritime hunter-gatherers from the Chilean Patagonia.

Patagonia was the last region of the Americas reached by humans who entered the continent from Siberia ∼15,000-20,000 y ago. Despite recent genomic approaches to reconstruct the continental evolutionary history, regional characterization of ancient and modern genomes remains understudied. Exploring the genomic diversity within Patagonia is not just a valuable strategy to gain a better understanding of the history and diversification of human populations in the southernmost tip of the Americas, but it would also improve the representation of Native American diversity in global databases of human variation. Here, we present genome data from four modern populations from Central Southern Chile and Patagonia (n = 61) and four ancient maritime individuals from Patagonia (∼1,000 y old). Both the modern and ancient individuals studied in this work have a greater genetic affinity with other modern Native Americans than to any non-American population, showing within South America a clear structure between major geographical regions. Native Patagonian Kawéskar and Yámana showed the highest genetic affinity with the ancient individuals, indicating genetic continuity in the region during the past 1,000 y before present, together with an important agreement between the ethnic affiliation and historical distribution of both groups. Lastly, the ancient maritime individuals were genetically equidistant to a ∼200-y-old terrestrial hunter-gatherer from Tierra del Fuego, which supports a model with an initial separation of a common ancestral group to both maritime populations from a terrestrial population, with a later diversification of the maritime groups.

Human-SARS-CoV-2 interactome and human genetic diversity: TMPRSS2-rs2070788, associated with severe influenza, and its population genetics caveats in Native Americans.

For human/SARS-CoV-2 interactome genes ACE2, TMPRSS2 and BSG, there is a convincing evidence of association in Asians with influenza-induced SARS for TMPRSS2-rs2070788, tag-SNP of the eQTL rs383510. This case illustrates the importance of population genetics and of sequencing data in the design of genetic association studies in different human populations: the high linkage disequilibrium (LD) between rs2070788 and rs383510 is Asian-specific. Leveraging on a combination of genotyping and sequencing data for Native Americans (neglected in genetic studies), we show that while their frequencies of the Asian tag-SNP rs2070788 is, surprisingly, the highest worldwide, it is not in LD with the eQTL rs383510, that therefore, should be directly genotyped in genetic association studies of SARS in populations with Native American ancestry.

Polymorphisms PSCA rs2294008, IL-4 rs2243250 and MUC1 rs4072037 are associated with gastric cancer in a high risk population.

Genetic variants are considered risk factors for gastric cancer. To date, 61 polymorphisms have been identified as associated with this disease. The aim of the present study was to analyze the association of some of those polymorphisms with GC in Chile. We performed a case-control study including 310 gastric cancer cases and 311 controls to assess the association of 36 single-nucleotide polymorphisms genotyped by Global Screening Array (GSA). Three polymorphisms was significantly associated: PSCA rs2294008 (allele model, OR = 1.49, 95%CI 1.17-1.88, P = 1.08 × 10-3), IL-4 rs2243250 (allele model, OR = 1.28, 95%CI 1.01-1.62, P = 0.04), and MUC1 rs4072037 (allele model, OR = 0.78, 95%CI 0.61-0.99, P = 0.04).PSCA rs2294008, IL-4 rs2243250 and MUC1 rs4072037 are associated with gastric cancer in Chile. It suggests that those polymorphisms could be used as biomarkers to assess the genetic risk for this cancer outside of the previously studied populations, not only for East Asians and Caucasians populations.

Development of a small panel of SNPs to infer ancestry in Chileans that distinguishes Aymara and Mapuche components.

BACKGROUND: Current South American populations trace their origins mainly to three continental ancestries, i.e. European, Amerindian and African. Individual variation in relative proportions of each of these ancestries may be confounded with socio-economic factors due to population stratification. Therefore, ancestry is a potential confounder variable that should be considered in epidemiologic studies and in public health plans. However, there are few studies that have assessed the ancestry of the current admixed Chilean population. This is partly due to the high cost of genome-scale technologies commonly used to estimate ancestry. In this study we have designed a small panel of SNPs to accurately assess ancestry in the largest sampling to date of the Chilean mestizo population (n = 3349) from eight cities. Our panel is also able to distinguish between the two main Amerindian components of Chileans: Aymara from the north and Mapuche from the south. RESULTS: A panel of 150 ancestry-informative markers (AIMs) of SNP type was selected to maximize ancestry informativeness and genome coverage. Of these, 147 were successfully genotyped by KASPar assays in 2843 samples, with an average missing rate of 0.012, and a 0.95 concordance with microarray data. The ancestries estimated with the panel of AIMs had relative high correlations (0.88 for European, 0.91 for Amerindian, 0.70 for Aymara, and 0.68 for Mapuche components) with those obtained with AXIOM LAT1 array. The country's average ancestry was 0.53 ± 0.14 European, 0.04 ± 0.04 African, and 0.42 ± 0.14 Amerindian, disaggregated into 0.18 ± 0.15 Aymara and 0.25 ± 0.13 Mapuche. However, Mapuche ancestry was highest in the south (40.03%) and Aymara in the north (35.61%) as expected from the historical location of these ethnic groups. We make our results available through an online app and demonstrate how it can be used to adjust for ancestry when testing association between incidence of a disease and nongenetic risk factors. CONCLUSIONS: We have conducted the most extensive sampling, across many different cities, of current Chilean population. Ancestry varied significantly by latitude and human development. The panel of AIMs is available to the community for estimating ancestry at low cost in Chileans and other populations with similar ancestry.

HIST1H1E heterozygous protein-truncating variants cause a recognizable syndrome with intellectual disability and distinctive facial gestalt: A study to clarify the HIST1H1E syndrome phenotype in 30 individuals.

Histone Gene Cluster 1 Member E, HIST1H1E, encodes Histone H1.4, is one of a family of epigenetic regulator genes, acts as a linker histone protein, and is responsible for higher order chromatin structure. HIST1H1E syndrome (also known as Rahman syndrome, OMIM #617537) is a recently described intellectual disability (ID) syndrome. Since the initial description of five unrelated individuals with three different heterozygous protein-truncating variants (PTVs) in the HIST1H1E gene in 2017, we have recruited 30 patients, all with HIST1H1E PTVs that result in the same shift in frame and that cluster to a 94-base pair region in the HIST1H1E carboxy terminal domain. The identification of 30 patients with HIST1H1E variants has allowed the clarification of the HIST1H1E syndrome phenotype. Major findings include an ID and a recognizable facial appearance. ID was reported in all patients and is most frequently of moderate severity. The facial gestalt consists of a high frontal hairline and full lower cheeks in early childhood and, in later childhood and adulthood, affected individuals have a strikingly high frontal hairline, frontal bossing, and deep-set eyes. Other associated clinical features include hypothyroidism, abnormal dentition, behavioral issues, cryptorchidism, skeletal anomalies, and cardiac anomalies. Brain magnetic resonance imaging (MRI) is frequently abnormal with a slender corpus callosum a frequent finding.

ADAR1-mediated RNA-editing of 3'UTRs in breast cancer.

BACKGROUND: Whole transcriptome RNA variant analyses have shown that adenosine deaminases acting on RNA (ADAR) enzymes modify a large proportion of cellular RNAs, contributing to transcriptome diversity and cancer evolution. Despite the advances in the understanding of ADAR function in breast cancer, ADAR RNA editing functional consequences are not fully addressed. RESULTS: We characterized A to G(I) mRNA editing in 81 breast cell lines, showing increased editing at 3'UTR and exonic regions in breast cancer cells compared to immortalized non-malignant cell lines. In addition, tumors from the BRCA TCGA cohort show a 24% increase in editing over normal breast samples when looking at 571 well-characterized UTRs targeted by ADAR1. Basal-like subtype breast cancer patients with high level of ADAR1 mRNA expression shows a worse clinical outcome and increased editing in their 3'UTRs. Interestingly, editing was particularly increased in the 3'UTRs of ATM, GINS4 and POLH transcripts in tumors, which correlated with their mRNA expression. We confirmed the role of ADAR1 in this regulation using a shRNA in a breast cancer cell line (ZR-75-1). CONCLUSIONS: Altogether, these results revealed a significant association between the mRNA editing in genes related to cancer-relevant pathways and clinical outcomes, suggesting an important role of ADAR1 expression and function in breast cancer.

chromPlot: visualization of genomic data in chromosomal context.

UNLABELLED: : Visualizing genomic data in chromosomal context can help detecting errors in data processing and may suggest new hypotheses to be tested. Here, we report a new tool for displaying large and diverse genomic data along chromosomes. The software is implemented in R so that visualization can be easily integrated with its numerous packages for processing genomic data. It supports simultaneous visualization of multiple tracks of data. Large genomic regions such as QTLs or synteny tracts may be shown along histograms of number of genes, genetic variants, or any other type of genomic element. Tracks can also contain values for continuous or categorical variables and the user can choose among points, connected lines, colored segments, or histograms for representing data. chromPlot takes data from tables in data.frame in GRanges formats. The information necessary to draw chromosomes for mouse and human is included with the package. For other organisms, chromPlot can read Gap and cytoBandIdeo tables from the UCSC Genome Browser. We present common use cases here, and a full tutorial is included as the package's vignette. AVAILABILITY AND IMPLEMENTATION: chromPlot is distributed under a GLP2 licence at http://www.bioconductor.org CONTACT: raverdugo@u.uchile.cl SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Overlapping mouse subcongenic strains successfully separate two linked body fat QTL on distal MMU 2.

BACKGROUND: Mouse chromosome 2 is linked to growth and body fat phenotypes in many mouse crosses. With the goal to identify the underlying genes regulating growth and body fat on mouse chromosome 2, we developed five overlapping subcongenic strains that contained CAST/EiJ donor regions in a C57BL/6J (hg/hg) background (hg is a spontaneous deletion of 500 Kb on mouse chromosome 10). To fine map QTL on distal mouse chromosome 2 a total of 1,712 F2 mice from the five subcongenic strains, plus 278 F2 mice from the HG2D founder congenic strain were phenotyped and analyzed. Interval mapping (IM) and composite IM (CIM) were performed on body weight and body fat traits on a combination of SNP and microsatellite markers, which generated a high-density genotyping panel. RESULTS: Phenotypic analysis and interval mapping of total fat mass identified two QTL on distal mouse chromosome 2. One QTL between 150 and 161 Mb, Fatq2a, and the second between 173.3 and 175.6 Mb, Fatq2b. The two QTL reside in different congenic strains with significant total fat differences between homozygous cast/cast and b6/b6 littermates. Both of these QTL were previously identified only as a single QTL affecting body fat, Fatq2. Furthermore, through a novel approach referred here as replicated CIM, Fatq2b was mapped to the Gnas imprinted locus. CONCLUSIONS: The integration of subcongenic strains, high-density genotyping, and CIM succesfully partitioned two previously linked QTL 20 Mb apart, and the strongest QTL, Fatq2b, was fine mapped to a ~2.3 Mb region interval encompassing the Gnas imprinted locus.

Genetic analysis of complex traits in the emerging Collaborative Cross.

The Collaborative Cross (CC) is a mouse recombinant inbred strain panel that is being developed as a resource for mammalian systems genetics. Here we describe an experiment that uses partially inbred CC lines to evaluate the genetic properties and utility of this emerging resource. Genome-wide analysis of the incipient strains reveals high genetic diversity, balanced allele frequencies, and dense, evenly distributed recombination sites-all ideal qualities for a systems genetics resource. We map discrete, complex, and biomolecular traits and contrast two quantitative trait locus (QTL) mapping approaches. Analysis based on inferred haplotypes improves power, reduces false discovery, and provides information to identify and prioritize candidate genes that is unique to multifounder crosses like the CC. The number of expression QTLs discovered here exceeds all previous efforts at eQTL mapping in mice, and we map local eQTL at 1-Mb resolution. We demonstrate that the genetic diversity of the CC, which derives from random mixing of eight founder strains, results in high phenotypic diversity and enhances our ability to map causative loci underlying complex disease-related traits.

Early mutational signatures and transmissibility of SARS-CoV-2 Gamma and Lambda variants in Chile.

Genomic surveillance (GS) programmes were crucial in identifying and quantifying the mutating patterns of SARS-CoV-2 during the COVID-19 pandemic. In this work, we develop a Bayesian framework to quantify the relative transmissibility of different variants tailored for regions with limited GS. We use it to study the relative transmissibility of SARS-CoV-2 variants in Chile. Among the 3443 SARS-CoV-2 genomes collected between January and June 2021, where sampling was designed to be representative, the Gamma (P.1), Lambda (C.37), Alpha (B.1.1.7), B.1.1.348, and B.1.1 lineages were predominant. We found that Lambda and Gamma variants' reproduction numbers were 5% (95% CI: [1%, 14%]) and 16% (95% CI: [11%, 21%]) larger than Alpha's, respectively. Besides, we observed a systematic mutation enrichment in the Spike gene for all circulating variants, which strongly correlated with variants' transmissibility during the studied period (r = 0.93, p-value = 0.025). We also characterised the mutational signatures of local samples and their evolution over time and with the progress of vaccination, comparing them with those of samples collected in other regions worldwide. Altogether, our work provides a reliable method for quantifying variant transmissibility under subsampling and emphasises the importance of continuous genomic surveillance.

Early-life nutrition interacts with developmental genes to shape the brain and sleep behavior in Drosophila melanogaster.

The mechanisms by which the genotype interacts with nutrition during development to contribute to the variation of complex behaviors and brain morphology of adults are not well understood. Here we use the Drosophila Genetic Reference Panel to identify genes and pathways underlying these interactions in sleep behavior and mushroom body morphology. We show that early-life nutritional restriction effects on sleep behavior and brain morphology depends on the genotype. We mapped genes associated with sleep sensitivity to early-life nutrition, which were enriched for protein-protein interactions responsible for translation, endocytosis regulation, ubiquitination, lipid metabolism, and neural development. By manipulating the expression of candidate genes in the mushroom bodies (MBs) and all neurons, we confirm that genes regulating neural development, translation and insulin signaling contribute to the variable response of sleep and brain morphology to early-life nutrition. We show that the interaction between differential expression of candidate genes with nutritional restriction in early life resides in the MBs or other neurons and that these effects are sex-specific. Natural variations in genes that control the systemic response to nutrition and brain development and function interact with early-life nutrition in different types of neurons to contribute to the variation of brain morphology and adult sleep behavior.

Clinical and pulmonary function analysis in long-COVID revealed that long-term pulmonary dysfunction is associated with vascular inflammation pathways and metabolic syndrome.

Introduction: Long-term pulmonary dysfunction (L-TPD) is one of the most critical manifestations of long-COVID. This lung affection has been associated with disease severity during the acute phase and the presence of previous comorbidities, however, the clinical manifestations, the concomitant consequences and the molecular pathways supporting this clinical condition remain unknown. The aim of this study was to identify and characterize L-TPD in patients with long-COVID and elucidate the main pathways and long-term consequences attributed to this condition by analyzing clinical parameters and functional tests supported by machine learning and serum proteome profiling. Methods: Patients with L-TPD were classified according to the results of their computer-tomography (CT) scan and diffusing capacity of the lungs for carbon monoxide adjusted for hemoglobin (DLCOc) tests at 4 and 12-months post-infection. Results: Regarding the acute phase, our data showed that L-TPD was favored in elderly patients with hypertension or insulin resistance, supported by pathways associated with vascular inflammation and chemotaxis of phagocytes, according to computer proteomics. Then, at 4-months post-infection, clinical and functional tests revealed that L-TPD patients exhibited a restrictive lung condition, impaired aerobic capacity and reduced muscular strength. At this time point, high circulating levels of platelets and CXCL9, and an inhibited FCgamma-receptor-mediated-phagocytosis due to reduced FcγRIII (CD16) expression in CD14+ monocytes was observed in patients with L-TPD. Finally, 1-year post infection, patients with L-TPD worsened metabolic syndrome and augmented body mass index in comparison with other patient groups. Discussion: Overall, our data demonstrated that CT scan and DLCOc identified patients with L-TPD after COVID-19. This condition was associated with vascular inflammation and impair phagocytosis of virus-antibody immune complexes by reduced FcγRIII expression. In addition, we conclude that COVID-19 survivors required a personalized follow-up and adequate intervention to reduce long-term sequelae and the appearance of further metabolic diseases.

Using bioinformatics and systems genetics to dissect HDL-cholesterol genetics in an MRL/MpJ x SM/J intercross.

A higher incidence of coronary artery disease is associated with a lower level of HDL-cholesterol. We searched for genetic loci influencing HDL-cholesterol in F2 mice from a cross between MRL/MpJ and SM/J mice. Quantitative trait loci (QTL) mapping revealed one significant HDL QTL (Apoa2 locus), four suggestive QTL on chromosomes 10, 11, 13, and 18 and four additional QTL on chromosomes 1 proximal, 3, 4, and 7 after adjusting HDL for the strong Apoa2 locus. A novel nonsynonymous polymorphism supports Lipg as the QTL gene for the chromosome 18 QTL, and a difference in Abca1 expression in liver tissue supports it as the QTL gene for the chromosome 4 QTL. Using weighted gene co-expression network analysis, we identified a module that after adjustment for Apoa2, correlated with HDL, was genetically determined by a QTL on chromosome 11, and overlapped with the HDL QTL. A combination of bioinformatics tools and systems genetics helped identify several candidate genes for both the chromosome 11 HDL and module QTL based on differential expression between the parental strains, cis regulation of expression, and causality modeling. We conclude that integrating systems genetics to a more-traditional genetics approach improves the power of complex trait gene identification.

Integration of QTL and bioinformatic tools to identify candidate genes for triglycerides in mice.

To identify genetic loci influencing lipid levels, we performed quantitative trait loci (QTL) analysis between inbred mouse strains MRL/MpJ and SM/J, measuring triglyceride levels at 8 weeks of age in F2 mice fed a chow diet. We identified one significant QTL on chromosome (Chr) 15 and three suggestive QTL on Chrs 2, 7, and 17. We also carried out microarray analysis on the livers of parental strains of 282 F2 mice and used these data to find cis-regulated expression QTL. We then narrowed the list of candidate genes under significant QTL using a "toolbox" of bioinformatic resources, including haplotype analysis; parental strain comparison for gene expression differences and nonsynonymous coding single nucleotide polymorphisms (SNP); cis-regulated eQTL in livers of F2 mice; correlation between gene expression and phenotype; and conditioning of expression on the phenotype. We suggest Slc25a7 as a candidate gene for the Chr 7 QTL and, based on expression differences, five genes (Polr3 h, Cyp2d22, Cyp2d26, Tspo, and Ttll12) as candidate genes for Chr 15 QTL. This study shows how bioinformatics can be used effectively to reduce candidate gene lists for QTL related to complex traits.

Importance of randomization in microarray experimental designs with Illumina platforms.

Measurements of gene expression from microarray experiments are highly dependent on experimental design. Systematic noise can be introduced into the data at numerous steps. On Illumina BeadChips, multiple samples are assayed in an ordered series of arrays. Two experiments were performed using the same samples but different hybridization designs. An experiment confounding genotype with BeadChip and treatment with array position was compared to another experiment in which these factors were randomized to BeadChip and array position. An ordinal effect of array position on intensity values was observed in both experiments. We demonstrate that there is increased rate of false-positive results in the confounded design and that attempts to correct for confounded effects by statistical modeling reduce power of detection for true differential expression. Simple analysis models without post hoc corrections provide the best results possible for a given experimental design. Normalization improved differential expression testing in both experiments but randomization was the most important factor for establishing accurate results. We conclude that lack of randomization cannot be corrected by normalization or by analytical methods. Proper randomization is essential for successful microarray experiments.

Serious limitations of the QTL/microarray approach for QTL gene discovery.

BACKGROUND: It has been proposed that the use of gene expression microarrays in nonrecombinant parental or congenic strains can accelerate the process of isolating individual genes underlying quantitative trait loci (QTL). However, the effectiveness of this approach has not been assessed. RESULTS: Thirty-seven studies that have implemented the QTL/microarray approach in rodents were reviewed. About 30% of studies showed enrichment for QTL candidates, mostly in comparisons between congenic and background strains. Three studies led to the identification of an underlying QTL gene. To complement the literature results, a microarray experiment was performed using three mouse congenic strains isolating the effects of at least 25 biometric QTL. Results show that genes in the congenic donor regions were preferentially selected. However, within donor regions, the distribution of differentially expressed genes was homogeneous once gene density was accounted for. Genes within identical-by-descent (IBD) regions were less likely to be differentially expressed in chromosome 2, but not in chromosomes 11 and 17. Furthermore, expression of QTL regulated in cis (cis eQTL) showed higher expression in the background genotype, which was partially explained by the presence of single nucleotide polymorphisms (SNP). CONCLUSIONS: The literature shows limited successes from the QTL/microarray approach to identify QTL genes. Our own results from microarray profiling of three congenic strains revealed a strong tendency to select cis-eQTL over trans-eQTL. IBD regions had little effect on rate of differential expression, and we provide several reasons why IBD should not be used to discard eQTL candidates. In addition, mismatch probes produced false cis-eQTL that could not be completely removed with the current strains genotypes and low probe density microarrays. The reviewed studies did not account for lack of coverage from the platforms used and therefore removed genes that were not tested. Together, our results explain the tendency to report QTL candidates as differentially expressed and indicate that the utility of the QTL/microarray as currently implemented is limited. Alternatives are proposed that make use of microarray data from multiple experiments to overcome the outlined limitations.

The Association between Fasting Glucose and Sugar Sweetened Beverages Intake Is Greater in Latin Americans with a High Polygenic Risk Score for Type 2 Diabetes Mellitus.

Chile is one of the largest consumers of sugar-sweetened beverages (SSB) world-wide. However, it is unknown whether the effects from this highly industrialized food will mimic those reported in industrialized countries or whether they will be modified by local lifestyle or population genetics. Our goal is to evaluate the interaction effect between SSB intake and T2D susceptibility on fasting glucose. We calculated a weighted genetic risk score (GRSw) based on 16 T2D risk SNPs in 2828 non-diabetic participants of the MAUCO cohort. SSB intake was categorized in four levels using a food frequency questionnaire. Log-fasting glucose was regressed on SSB and GRSw tertiles while accounting for socio-demography, lifestyle, obesity, and Amerindian ancestry. Fasting glucose increased systematically per unit of GRSw (ß = 0.02 ± 0.006, p = 0.00002) and by SSB intake (ß[cat4] = 0.04 ± 0.01, p = 0.0001), showing a significant interaction, where the strongest effect was observed in the highest GRSw-tertile and in the highest SSB consumption category (ß = 0.05 ± 0.02, p = 0.02). SNP-wise, SSB interacted with additive effects of rs7903146 (TCF7L2) (ß = 0.05 ± 0.01, p = 0.002) and with the G/G genotype of rs10830963 (MTNRB1B) (ß = 0.19 ± 0.05, p = 0.001). Conclusions: The association between SSB intake and fasting glucose in the Chilean population without diabetes is modified by T2D genetic susceptibility.

Inflammatory profiles in Chilean Mapuche and non-Mapuche women with gallstones at risk of developing gallbladder cancer.

Chile has high incidence rates of gallbladder cancer globally, particularly among Amerindian women, who also have a high prevalence of gallstones. We examined differences in inflammatory biomarkers between Mapuche and non-Mapuche women from the Chile Biliary Longitudinal Study, a cohort of women with ultrasound-detected gallstones. We randomly selected 200 Mapuche women frequency matched to non-Mapuche women on age and statin use Inflammatory biomarkers were analyzed using a multiplex assay and linear regression to assess associations of a priori markers (CCL20, CXCL10, IL-6, and IL-8) with ethnicity. Novel biomarkers were analyzed using exploratory factor analysis (EFA) and sufficient dimension reduction (SDR) to identify correlated marker groups, followed by linear regression to examine their association with ethnicity. The mean values of IL-8 were higher in Mapuche than non-Mapuche women (P = 0.04), while CCL20, CXCL10, and IL-6 did not differ significantly by ethnicity. EFA revealed two marker groups associated with ethnicity (P = 0.03 and P < 0.001). SDR analysis confirmed correlation between the biomarkers and ethnicity. We found higher IL-8 levels among Mapuche than non-Mapuche women. Novel inflammatory biomarkers were correlated with ethnicity and should be studied further for their role in gallbladder disease. These findings may elucidate underlying ethnic disparities in gallstones and carcinogenesis among Amerindians.