HLA-G Gene Variability Is Associated with Papillary Thyroid Carcinoma Morbidity and the HLA-G Protein Profile.

Human leukocyte antigen (HLA)-G is an immune checkpoint molecule that is highly expressed in papillary thyroid carcinoma (PTC). The HLA-G gene presents several functional polymorphisms distributed across the coding and regulatory regions (5'URR: 5' upstream regulatory region and 3'UTR: 3' untranslated region) and some of them may impact HLA-G expression and human malignancy. To understand the contribution of the HLA-G genetic background in PTC, we studied the HLA-G gene variability in PTC patients in association with tumor morbidity, HLA-G tissue expression, and plasma soluble (sHLA-G) levels. We evaluated 185 PTC patients and 154 healthy controls. Polymorphic sites defining coding, regulatory and extended haplotypes were characterized by sequencing analyses. HLA-G tissue expression and plasma soluble HLA-G levels were evaluated by immunohistochemistry and ELISA, respectively. Compared to the controls, the G0104a(5'URR)G*01:04:04(coding)UTR-03(3'UTR) extended haplotype was underrepresented in the PTC patients, while G0104a(5'URR)G*01:04:01(coding)UTR-03(3'UTR) was less frequent in patients with metastatic and multifocal tumors. Decreased HLA-G tissue expression and undetectable plasma sHLA-G were associated with the G010102a(5'URR)G*01:01:02:01(coding)UTR-02(3'UTR) extended haplotype. We concluded that the HLA-G variability was associated with PTC development and morbidity, as well as the magnitude of the encoded protein expression at local and systemic levels.

BRAF, TERT and HLA-G Status in the Papillary Thyroid Carcinoma: A Clinicopathological Association Study.

As BRAF, TERT, HLA-G, and microRNAs have been individually associated with papillary thyroid carcinoma (PTC), we aimed to evaluate the individual and collaborative role of these markers in PTC in the same patient cohort. HLA-G and BRAF tumor expression was evaluated by immunohistochemistry. Using molecular methods, BRAFV600E and TERT promoter mutations were evaluated in thyroid fine needle aspirates. MicroRNA tumor profiling was investigated using massively parallel sequencing. We observed strong HLA-G (67.96%) while BRAF (62.43%) staining was observed in PTC specimens. BRAF overexpression was associated with poor response to therapy. The BRAFV600E (52.9%) and TERTC228T (13%) mutations were associated with extrathyroidal extension, advanced-age, and advanced-stage cancer. The TERT rs2853669 CC+TC genotypes (38%) were overrepresented in metastatic tumors. Nine modulated microRNAs targeting the BRAF, TERT, and/or HLA-G genes were observed in PTC and involved with cancer-related signaling pathways. The markers were individually associated with PTC features, emphasizing the synergistic effect of BRAFV600E and TERTC228T; however, their collaborative role on PTC outcome was not fully demonstrated. The differentially expressed miRNAs targeting the BRAF and/or HLA-G genes may explain their increased expression in the tumor milieu.

HLA-G liver expression and HLA-G extended haplotypes are associated with chronic hepatitis C in HIV-negative and HIV-coinfected patients.

Chronic hepatitis C virus (HCV) infection induces liver damage and the HCV/Human Immunodeficiency Virus (HIV)-coinfection may further contribute to its progression. The HLA-G molecule inhibits innate and adaptive immunity and may be deleterious for chronically virus-infected cells. Thus we studied 204 HCV-mono-infected patients, 142 HCV/HIV-coinfected patients, 104 HIV-mono-infected patients and 163 healthy subjects. HLA-G liver expression was similarly induced in HCV and HCV/HIV specimens, increasing with advanced fibrosis and necroinflammatory activity, and with increased levels of liver function-related enzymes. Plasma soluble HLA-G (sHLA-G) levels were higher in HCV/HIV patients compared to HCV, HIV and to healthy individuals. sHLA-G continued to be higher in coinfected patients even after stratification of samples according to degree of liver fibrosis and necroinflammatory activity when compared to mono-infected patients. Some HLA-G gene haplotypes differentiated patient groups and presented few associations with liver and plasma HLA-G expression. HLA-G thus may help to distinguish patient groups.

Origin and dynamics of admixture in Brazilians and its effect on the pattern of deleterious mutations.

While South Americans are underrepresented in human genomic diversity studies, Brazil has been a classical model for population genetics studies on admixture. We present the results of the EPIGEN Brazil Initiative, the most comprehensive up-to-date genomic analysis of any Latin-American population. A population-based genome-wide analysis of 6,487 individuals was performed in the context of worldwide genomic diversity to elucidate how ancestry, kinship, and inbreeding interact in three populations with different histories from the Northeast (African ancestry: 50%), Southeast, and South (both with European ancestry >70%) of Brazil. We showed that ancestry-positive assortative mating permeated Brazilian history. We traced European ancestry in the Southeast/South to a wider European/Middle Eastern region with respect to the Northeast, where ancestry seems restricted to Iberia. By developing an approximate Bayesian computation framework, we infer more recent European immigration to the Southeast/South than to the Northeast. Also, the observed low Native-American ancestry (6-8%) was mostly introduced in different regions of Brazil soon after the European Conquest. We broadened our understanding of the African diaspora, the major destination of which was Brazil, by revealing that Brazilians display two within-Africa ancestry components: one associated with non-Bantu/western Africans (more evident in the Northeast and African Americans) and one associated with Bantu/eastern Africans (more present in the Southeast/South). Furthermore, the whole-genome analysis of 30 individuals (42-fold deep coverage) shows that continental admixture rather than local post-Columbian history is the main and complex determinant of the individual amount of deleterious genotypes.

Global and local ancestry estimation in a captive baboon colony.

The last couple of decades have highlighted the importance of studying hybridization, particularly among primate species, as it allows us to better understand our own evolutionary trajectory. Here, we report on genetic ancestry estimates using dense, full genome data from 881 olive (Papio anubus), yellow (Papio cynocephalus), or olive-yellow crossed captive baboons from the Southwest National Primate Research Center. We calculated global and local ancestry information, imputed low coverage genomes (n = 830) to improve marker quality, and updated the genetic resources of baboons available to assist future studies. We found evidence of historical admixture in some putatively purebred animals and identified errors within the Southwest National Primate Research Center pedigree. We also compared the outputs between two different phasing and imputation pipelines along with two different global ancestry estimation software. There was good agreement between the global ancestry estimation software, with R2 > 0.88, while evidence of phase switch errors increased depending on what phasing and imputation pipeline was used. We also generated updated genetic maps and created a concise set of ancestry informative markers (n = 1,747) to accurately obtain global ancestry estimates.

Characterization of pharmacogenomic variants in a Brazilian admixed cohort of elderly individuals based on whole-genome sequencing data.

Introduction: Research in the field of pharmacogenomics (PGx) aims to identify genetic variants that modulate response to drugs, through alterations in their pharmacokinetics (PK) or pharmacodynamics (PD). The distribution of PGx variants differs considerably among populations, and whole-genome sequencing (WGS) plays a major role as a comprehensive approach to detect both common and rare variants. This study evaluated the frequency of PGx markers in the context of the Brazilian population, using data from a population-based admixed cohort from Sao Paulo, Brazil, which includes variants from WGS of 1,171 unrelated, elderly individuals. Methods: The Stargazer tool was used to call star alleles and structural variants (SVs) from 38 pharmacogenes. Clinically relevant variants were investigated, and the predicted drug response phenotype was analyzed in combination with the medication record to assess individuals potentially at high-risk of gene-drug interaction. Results: In total, 352 unique star alleles or haplotypes were observed, of which 255 and 199 had a frequency < 0.05 and < 0.01, respectively. For star alleles with frequency > 5% (n = 97), decreased, loss-of-function and unknown function accounted for 13.4%, 8.2% and 27.8% of alleles or haplotypes, respectively. Structural variants (SVs) were identified in 35 genes for at least one individual, and occurred with frequencies >5% for CYP2D6, CYP2A6, GSTM1, and UGT2B17. Overall 98.0% of the individuals carried at least one high risk genotype-predicted phenotype in pharmacogenes with PharmGKB level of evidence 1A for drug interaction. The Electronic Health Record (EHR) Priority Result Notation and the cohort medication registry were combined to assess high-risk gene-drug interactions. In general, 42.0% of the cohort used at least one PharmGKB evidence level 1A drug, and 18.9% of individuals who used PharmGKB evidence level 1A drugs had a genotype-predicted phenotype of high-risk gene-drug interaction. Conclusion: This study described the applicability of next-generation sequencing (NGS) techniques for translating PGx variants into clinically relevant phenotypes on a large scale in the Brazilian population and explores the feasibility of systematic adoption of PGx testing in Brazil.

Identification of ancestry proportions in admixed groups across the Americas using clinical pharmacogenomic SNP panels.

We evaluated the performance of three PGx panels to estimate biogeographical ancestry: the DMET panel, and the VIP and Preemptive PGx panels described in the literature. Our analysis indicate that the three panels capture quite well the individual variation in admixture proportions observed in recently admixed populations throughout the Americas, with the Preemptive PGx and DMET panels performing better than the VIP panel. We show that these panels provide reliable information about biogeographic ancestry and can be used to guide the implementation of PGx clinical decision-support (CDS) tools. We also report that using these panels it is possible to control for the effects of population stratification in association studies in recently admixed populations, as exemplified with a warfarin dosing GWA study in a sample from Brazil.

HLA-E gene polymorphisms in chronic hepatitis C: Impact on HLA-E liver expression and disease severity.

Hepatitis C virus usually produces chronic infection and liver damage. Considering that: i) the human leukocyte antigen-E (HLA-E) molecule may modulate the immune response, and ii) little is known about the role of HLA-E gene variability on chronic hepatitis C, we studied the impact of HLA-E polymorphisms on the magnitude of HLA-E liver expression and severity of hepatitis C. HLA-E variability was evaluated in terms of: i) single nucleotide polymorphism (SNP) alleles and genotypes along the gene (beginning of the promoter region, coding region and 3'UTR), and ii) ensemble of SNPs that defines the coding region alleles, considered individually or as genotypes. The comparisons of the HLA-E variation sites between patients and controls revealed no significant results. The HLA-E + 424 T > C (rs1059510), +756 G > A (rs1264457) and + 3777 G > A (rs1059655) variation sites and the HLA-E*01:01:01:01 and HLA-E*01:03:02:01 alleles, considered at single or double doses, were associated with the magnitude of HLA-E liver expression in Kupfer cell, steatosis, inflammatory activity and liver fibrosis. Although these associations were lost after corrections for multiple comparisons, these variable sites may propitiate biological clues for the understanding of the mechanisms associated with hepatitis C severity.

Novel insights on demographic history of tribal and caste groups from West Maharashtra (India) using genome-wide data.

The South Asian subcontinent is characterized by a complex history of human migrations and population interactions. In this study, we used genome-wide data to provide novel insights on the demographic history and population relationships of six Indo-European populations from the Indian State of West Maharashtra. The samples correspond to two castes (Deshastha Brahmins and Kunbi Marathas) and four tribal groups (Kokana, Warli, Bhil and Pawara). We show that tribal groups have had much smaller effective population sizes than castes, and that genetic drift has had a higher impact in tribal populations. We also show clear affinities between the Bhil and Pawara tribes, and to a lesser extent, between the Warli and Kokana tribes. Our comparisons with available modern and ancient DNA datasets from South Asia indicate that the Brahmin caste has higher Ancient Iranian and Steppe pastoralist contributions than the Kunbi Marathas caste. Additionally, in contrast to the two castes, tribal groups have very high Ancient Ancestral South Indian (AASI) contributions. Indo-European tribal groups tend to have higher Steppe contributions than Dravidian tribal groups, providing further support for the hypothesis that Steppe pastoralists were the source of Indo-European languages in South Asia, as well as Europe.

Insights on hair, skin and eye color of ancient and contemporary Native Americans.

Over the past few years, tools capable of predicting pigmentation phenotypes have been developed aiming to contribute for criminal and anthropological investigations. In this study, we used eight genetic systems to infer eye, hair, and skin color of ancient and contemporary Native Americans. To achieve this goal, we retrieved 61 SNPs from 42 samples available in free online repositories of DNA sequences. We performed pigmentation predictions using two freely available tools, HIrisPlex-S and Snipper, in addition to two other published models. This workflow made possible to predict all three phenotypes with at least one tool for 29 out of the 42 samples. Considering these 29 individuals, predictions for eye, hair, and skin color were obtained with HIrisPlex-S for 27, 28 and 27 individuals, respectively, while 24, 25 and 25 individuals had such predictions with Snipper. In general, ancient and contemporary Native Americans were predicted to have intermediate/brown eyes, black hair, and intermediate/darker skin pigmentation.

HLA-G Polymorphisms Are Associated with Non-segmental Vitiligo among Brazilians.

(1) Background: Vitiligo is characterized by white patches on the skin caused by loss of melanocyte activity or the absence of these cells. The available treatments minimize the symptoms by retarding the process of skin depigmentation or re-pigmenting the affected regions. New studies are required for a better comprehension of the mechanisms that trigger the disease and for the development of more efficient treatments. Studies have suggested an autoimmune feature for vitiligo, based on the occurrence of other autoimmune diseases in vitiligo patients and their relatives, and on the involvement of genes related to the immune response. (2) Methods: We evaluated, by massive parallel sequencing, polymorphisms of the HLA-G gene in vitiligo patients and control samples, to verify if variants of this gene could influence the susceptibility to vitiligo. (3) Results: We detected an association with non-segmental vitiligo regarding the haplotype Distal-010101a/G*01:01:01:01/UTR-1, adjusting for population stratification by using ancestry-informative markers (AIMs). (4) Conclusions: It remains unclear whether the HLA-G variants associated with vitiligo were detected because of the high linkage disequilibrium (LD) with HLA-A*02, or if the HLA-A variants previously reported as associated with vitiligo were detected because of the high LD with HLA-G*01:01:01:01/UTR-1, or if both genes jointly contribute to vitiligo susceptibility.

Applicability of the SNPforID 52-plex panel for human identification and ancestry evaluation in a Brazilian population sample by next-generation sequencing.

SNP analysis is of paramount importance in forensic genetics. The development of new technologies in next-generation sequencing allowed processing a large number of markers in various samples simultaneously. Although SNPs are less informative than STRs, they present lower mutation rates and perform better when using degraded samples. Some SNP systems were developed for forensic usage, such as the SNPforID 52-plex, from the SNPforID Consortium, containing 52 bi-allelic SNPs for human identification. In this paper we evaluated the informativeness of this system in a Brazilian population sample (n = 340). DNA libraries were prepared using a customized HaloPlex Target Enrichment System kit (Agilent Technologies, Inc.) and sequenced in the MiSeq Personal Sequencer platform (Illumina Inc.). The methodology presented here allowed the analysis of 51 out of 52 SNPforID markers. Allele frequencies and forensic parameters were estimated, revealing high informativeness: the combined match probability and power of exclusion were 6.48 × 10-21 and 0.9997, respectively. Population admixture analysis indicates high European contribution (more than 70%) and low Amerindian contribution (less than 10%) in our population, while individual admixture analyses were consistent with the majority of individuals presenting high European contribution. This study demonstrates that the 52-plex kit is suitable for forensic cases in a Brazilian population, presenting results comparable with those obtained using a 16 STR panel.

Extended HLA-G genetic diversity and ancestry composition in a Brazilian admixed population sample: Implications for HLA-G transcriptional control and for case-control association studies.

Human leukocyte antigen-G (HLA-G) is a nonclassical Major Histocompatibility Complex (MHC) molecule with immunomodulatory function and restricted tissue expression. The genetic diversity of HLA-G has been extensively studied in several populations, however, the segment located upstream -1406 has not yet been evaluated. We characterized the nucleotide variation and haplotype structure of an extended distal region (-2635), all exons and the 3'UTR segment of HLA-G by next-generation sequencing (NGS) in a sample of 335 Brazilian individuals. We detected 29 variants at the HLA-G distal promoter region, arranged into 19 haplotypes, among which we identified sites that may influence transcription factor targeting. Although the variation pattern in the distal region resembled the one observed in the conventional promoter segment, molecular signature for balancing selection was observed in the promoter segment from -1406 to -1 (Tajima's D = 2.315, P = 0.017), but not in this distal segment (D = 1.049, P = 0.118). Furthermore, the ancestry composition of this Brazilian population sample was determined by the analysis of SNPforID 34-plex ancestry informative marker (AIM) SNP panel. The distribution of HLA-G haplotypes was ancestry-dependent, corroborating previous findings and emphasizing the importance of considering the ancestry information in association studies.

Cardiometabolic risk markers of normal weight and excess body weight in Brazilian adolescents.

Excess body weight leads to a variety of metabolic changes and increases the risk for cardiovascular diseases (CVD) in adulthood. The objective of this study was to investigate the presence of risk markers for CVD among Brazilian adolescents of normal weight and with excess body weight. The markers included blood pressure, C-reactive protein, homocysteine, tumor necrosis factor alpha, fibrinogen, fasting insulin and glucose, homeostasis model assessment of insulin resistance (HOMA-IR), leptin, total cholesterol, low-density lipoprotein cholesterol (LDL-c), high-density lipoprotein cholesterol (HDL-c), and triglycerides. We calculated odds ratios (OR) using logistic regression and adjusted for potential confounders such as age, sex, physical activity, and socioeconomic background. Compared with normal weight subjects, overweight/obese adolescents were more likely to have higher systolic blood pressure (OR = 3.49, p < 0.001), fasting insulin (OR = 8.03, p < 0.001), HOMA-IR (OR = 8.03, p < 0.001), leptin (OR = 5.55, p < 0.001), and LDL-c (OR = 5.50, p < 0.001) and lower serum HDL-c concentrations (OR = 2.76, p = 0.004). After adjustment for confounders, the estimates did not change substantially, except for leptin for which the risk associated with overweight increased to 11.09 (95% CI: 4.05-30.35). In conclusion, excess body weight in adolescents exhibits strong associations with several markers that are established as causes of CVD in adults. This observation stresses the importance of primary prevention and of maintaining a healthy body weight throughout adolescence to reduce the global burden of CVD.

IL28B gene polymorphisms in mono- and HIV-coinfected chronic hepatitis C patients.

INTRODUCTION: Single-nucleotide polymorphisms (SNPs) associated with hepatitis C virus (HCV) clearance were identified near the IL28B gene. Coinfection by the human immunodeficiency virus (HIV) influences the course of HCV contributing to liver damage. Nevertheless, little is known about the relationship between these SNPs and HCV/HIV coinfection. Our aim was to estimate the frequencies of the allelic and genotypic variants of the IL28B polymorphisms rs12979860 (C/T) and rs8099917 (T/G) and their possible association with the establishment of HCV infection. METHODOLOGY: A total of 199 non-infected controls and 230 patients with chronic hepatitis C, including 53 coinfected with HIV, participated in the study. Genotyping consisted of polymerase chain reaction and subsequent analysis of the restriction patterns resulting from exposure to endonucleases. RESULTS: Among the controls with established results, 47.4% (90/190) exhibited the rs12979860 CC genotype, 43.7 CT, and 8.9% TT, whereas 29.1% (66/227), 51.5%, and 19.4% of the patients exhibited the CC, CT, and TT genotypes, respectively. With respect to rs8099917, 66.8% (133/199) of the controls exhibited the TT genotype, 31.2% TG, and 2.0% GG, whereas 56.1% (129/230), 40.9%, and 3.0% of the patients exhibited the TT, TG, and GG genotypes, respectively. CONCLUSION: The frequencies of the rs12979860 C allele and CC genotype and of the rs8099917 T allele and TT genotype were significantly higher among controls compared with patients, thus confirming the suggested protective effect against HCV infection. No significant difference was observed in the genotype and allelic distributions between the mono- and coinfected patients.