Prevalence of single-nucleotide variants in twenty-five pharmacogenes from a Cuban sample cohort.

Introduction: The Cuban population is genetically diverse, and information on the prevalence of genetic variants is still limited. As complex admixture processes have occurred, we hypothesized that the frequency of pharmacogenetic variants and drug responses may vary within the country. The aims of the study were to describe the frequency distribution of 43 single-nucleotide variants (SNVs) from 25 genes of pharmacogenetic interest within the Cuba population and in relation to other populations, while taking into consideration some descriptive variables such as place of birth and skin color. Materials and Methods: SNVs were analyzed in 357 unrelated healthy Cuban volunteers. Genotype, allele frequencies, and ancestry proportions were determined, and the pairwise fixation index (FST ) was evaluated. Results: Hardy-Weinberg equilibrium (HWE) deviations in six loci (rs11572103, rs2740574, rs776746, rs3025039, rs861539, and rs1762429) were identified. Minor allele frequencies (MAFs) ranged from 0.00 to 0.15 for variants in genes encoding xenobiotic metabolizing enzymes. They also ranged from 0.01 to 0.21 for variants in DNA repair, growth factors, methyltransferase, and methyl-binding proteins, while they ranged from 0.04 to 0.27 for variants in the O-6-methylguanine-DNA methyltransferase enzyme. Moderate genetic divergence was observed upon comparison to Africans (FST = 0.071 and SD 0.079), with 19 markers exhibiting moderate-to-large genetic differentiation. The average European, African, and Amerindian ancestry proportions were 67.8%, 27.2%, and 5.3%, respectively. Ancestry proportions differed by skin color and birthplace for both African and European components, with the exception of the European component, which showed no significant difference between individuals from Western and Eastern regions. Meanwhile, the statistical significance varied in comparisons by skin color and birthplace within the Amerindian component. Low genetic divergence was observed across geographical regions. We identified 12 variants showing moderate-to-large differentiation between White/Black individuals. Conclusion: Altogether, our results may support national strategies for the introduction of pharmacogenetic tools in clinical practice, contributing to the development of precision medicine in Cuba.

Association between HLA alleles and haplotypes with age at diagnosis of type 1 diabetes in an admixed Brazilian population: A nationwide study.

To investigate the potential relationship between HLA alleles and haplotypes and the age at diagnosis of type 1 diabetes (T1DAgeD) in an admixed Brazilian population. This nationwide study was conducted in public clinics across 12 Brazilian cities. We collected demographic and genetic data from 1,600 patients with T1D. DNA samples were utilised to determine genomic ancestry (GA) and perform HLA typings for DRB1, DQA1 and DQB1. We explored allele and haplotype frequencies and GA in patients grouped by T1DAgeD categories (<6 years, ≥6-<11 years, ≥11-<19 years and ≥19 years) through univariate and multivariate analyses and primary component analyses. Additionally, we considered self-reported colour-race and identified a familiar history of T1D in first-degree relatives. The homozygosity index for DRB1~DQA1~DQB1 haplotypes exhibited the highest variation among T1DAgeD groups, and the percentages of Sub-Saharan African and European ancestries showed opposite trends in principal component analysis (PCA) analyses. Regarding the association of alleles and haplotypes with T1DAgeD, risk alleles such as HLA-DQB1*03:02g, -DQA1*03:01g, -02:01g, DRB1*04:05g and -04:02g were more frequently observed in heterozygosity or homozygosity in T1D patients with an early disease onset. Conversely, alleles such as DRB1*07:01g, -13:03g, DQB1*06:02g and DQA1*02:01 were more prevalent in older T1D patients. The combination DR3/DR4.5 was significantly associated with early disease onset. However, gender, GA, familiar history of T1D and self-reported colour-race identity did not exhibit significant associations with the onset of T1D. It is worth noting that the very common risk haplotype DRB1*03:01g~DQA1*05:01g~DQB1*02:01g did not differentiate between T1DAgeD groups. In the admixed Brazilian population, the high-risk haplotype DRB1*04:05~DQA1*03:01~DQB1*03:02 was more prevalent in individuals diagnosed before 6 years of age. In contrast, the protective alleles DQA1*01:02g, DQB1*06:02g, DRB1*07:01g and DRB1*13:03g and haplotypes DRB1*13:03g~DQA1*05:01g~DQB1*03:01g and DRB1*16:02g~DQA1*01:02g~DQB1*05:02g were more frequently observed in patients diagnosed in adulthood. Notably, these associations were independent of factors such as sex, economic status, GA, familiar history of T1D and region of birth in Brazil. These alleles and haplotypes contribute to our understanding of the disease onset heterogeneity and may have implications for early interventions when detected in association with well-known genomic risk or protection factors for T1D.

A multi-ethnic reference panel to impute HLA classical and non-classical class I alleles in admixed samples: Testing imputation accuracy in an admixed sample from Brazil.

The MHC class I region contains crucial genes for the innate and adaptive immune response, playing a key role in susceptibility to many autoimmune and infectious diseases. Genome-wide association studies have identified numerous disease-associated SNPs within this region. However, these associations do not fully capture the immune-biological relevance of specific HLA alleles. HLA imputation techniques may leverage available SNP arrays by predicting allele genotypes based on the linkage disequilibrium between SNPs and specific HLA alleles. Successful imputation requires diverse and large reference panels, especially for admixed populations. This study employed a bioinformatics approach to call SNPs and HLA alleles in multi-ethnic samples from the 1000 genomes (1KG) dataset and admixed individuals from Brazil (SABE), utilising 30X whole-genome sequencing data. Using HIBAG, we created three reference panels: 1KG (n = 2504), SABE (n = 1171), and the full model (n = 3675) encompassing all samples. In extensive cross-validation of these reference panels, the multi-ethnic 1KG reference exhibited overall superior performance than the reference with only Brazilian samples. However, the best results were achieved with the full model. Additionally, we expanded the scope of imputation by developing reference panels for non-classical, MICA, MICB and HLA-H genes, previously unavailable for multi-ethnic populations. Validation in an independent Brazilian dataset showcased the superiority of our reference panels over the Michigan Imputation Server, particularly in predicting HLA-B alleles among Brazilians. Our investigations underscored the need to enhance or adapt reference panels to encompass the target population's genetic diversity, emphasising the significance of multiethnic references for accurate imputation across different populations.

A pilot study of mitochondrial genomic ancestry in admixed Brazilian patients with type 1 diabetes.

Interactions between multiple genes and environmental factors could be related to the pathogenesis of type 1 diabetes (T1D). The Brazilian population results from different historical miscegenation events, resulting in a highly diverse genetic pool. This study aimed to analyze the mtDNA of patients with T1D and to investigate whether there is a relationship between maternal ancestry, self-reported color and the presence of T1D. The mtDNA control region of 204 patients with T1D residing in three geographic regions of Brazil was sequenced following the International Society for Forensic Genetics (ISFG) recommendations. We obtained a frequency of Native American matrilineal origin (43.6%), African origin (38.2%), and European origin (18.1%). For self-declared color, 42.6% of the patients with diabetes reported that they were White, 50.9% were Brown, and 5.4% were Black. Finally, when we compared the self-declaration data with maternal ancestral origin, we found that for the self-declared White group, there was a greater percentage of haplogroups of Native American origin (50.6%); for the self-declared Black group, there was a greater percentage of African haplogroups (90.9%); and for the Brown group, there was a similar percentage of Native American and African haplogroups (42.3% and 45.2%, respectively). The Brazilian population with diabetic has a maternal heritage of more than 80% Native American and African origin, corroborating the country's colonization history.

Genome-Wide Admixture and Association Study of Serum Selenium Deficiency to Identify Genetic Variants Indirectly Linked to Selenium Regulation in Brazilian Adults.

Blood selenium (Se) concentrations differ substantially by population and could be influenced by genetic variants, increasing Se deficiency-related diseases. We conducted a genome-wide association study (GWAS) to identify single nucleotide polymorphisms (SNPs) associated with serum Se deficiency in 382 adults with admixed ancestry. Genotyping arrays were combined to yield 90,937 SNPs. R packages were applied to quality control and imputation. We also performed the ancestral proportion analysis. The Search Tool for the Retrieval of Interacting Genes was used to interrogate known protein-protein interaction networks (PPIs). Our ancestral proportion analysis estimated 71% of the genome was from Caucasians, 22% was from Africans, and 8% was from East Asians. We identified the SNP rs1561573 in the TraB domain containing 2B (TRABD2B), rs425664 in MAF bZIP transcription factor (MAF), rs10444656 in spermatogenesis-associated 13 (SPATA13), and rs6592284 in heat shock protein nuclear import factor (HIKESHI) genes. The PPI analysis showed functional associations of Se deficiency, thyroid hormone metabolism, NRF2-ARE and the Wnt pathway, and heat stress. Our findings show evidence of a genetic association between Se deficiency and metabolic pathways indirectly linked to Se regulation, reinforcing the complex relationship between Se intake and the endogenous factors affecting the Se requirements for optimal health.

Immune response stability to the SARS-CoV-2 mRNA vaccine booster is influenced by differential splicing of HLA genes.

Many molecular mechanisms that lead to the host antibody response to COVID-19 vaccines remain largely unknown. In this study, we used serum antibody detection combined with whole blood RNA-based transcriptome analysis to investigate variability in vaccine response in healthy recipients of a booster (third) dose schedule of the mRNA BNT162b2 vaccine against COVID-19. The cohort was divided into two groups: (1) low-stable individuals, with antibody concentration anti-SARS-CoV IgG S1 below 0.4 percentile at 180 days after boosting vaccination; and (2) high-stable individuals, with antibody values greater than 0.6 percentile of the range in the same period (median 9525 [185-80,000] AU/mL). Differential gene expression, expressed single nucleotide variants and insertions/deletions, differential splicing events, and allelic imbalance were explored to broaden our understanding of the immune response sustenance. Our analysis revealed a differential expression of genes with immunological functions in individuals with low antibody titers, compared to those with higher antibody titers, underscoring the fundamental importance of the innate immune response for boosting immunity. Our findings also provide new insights into the determinants of the immune response variability to the SARS-CoV-2 mRNA vaccine booster, highlighting the significance of differential splicing regulatory mechanisms, mainly concerning HLA alleles, in delineating vaccine immunogenicity.

Psychiatric Level 1A evidence pharmacogenomics in a Brazilian admixed cohort and global populations.

Purpose: To compare minor allele frequencies (MAFs) of psychiatric drug response variants in a Brazilian admixed cohort with global populations and other Brazilian groups. Methods: PharmGKB MAFs were gathered from publicly available genetic datasets for Brazil and worldwide. Results: Among 146 variants in CYP2D6 and CYP2C19, 41 were present in Brazil, mostly rare (MAF <1%). 11 variants showed significant MAF differences with large effect sizes compared with global populations. CYP2C19*3 (rs4986893), CYP2C19*17 (rs12248560), CYP2D6*17 (rs28371706-A) and CYP2D6*29 (rs61736512) exhibited higher frequencies in Brazil, with the latter three also differing from other Brazilian groups. Conclusion: This study highlights significant pharmacogenomic diversity in Brazil and globally, underscoring the need for more research in personalized psychiatric drug therapy.

LAENALS: epidemiological and clinical features of amyotrophic lateral sclerosis in Latin America.

OBJECTIVE: The Latin American Epidemiologic study of ALS (LAENALS) aims to gather data on ALS epidemiology, phenotype, and risk factors in Cuba, Chile, and Uruguay, to understand the impact of genetic and environmental factors on ALS. METHODS: A harmonized data collection protocol was generated, and a Latin-American Spanish language Register was constructed. Patient data were collected in Uruguay in 2018, in Chile from 2017 to 2019, and in Cuba between 2017 and 2018. Statistical analysis was performed using SPSS 25.0.0 software. Crude cumulative incidence, standardized incidence, and prevalence were calculated in the population aged 15 years and older. RESULTS: During 2017-2019, 90 people with ALS from Uruguay (55.6% men), 219 from Chile (54.6% men), and 49 from Cuba (55.1% men) were included. The cumulative crude incidence in 2018 was 1.73/100,000 persons in Uruguay, 1.08 in Chile and 0.195 in Cuba. Crude prevalence in 2018 was 2.19 per 100,000 persons in Uruguay, 1.39 in Chile and 0.55 in Cuba. Mean age at onset was 61.8 ± 11.96 SD years in Uruguay, 61.9 ± 10.4 SD years in Chile, and 60.21 ± 12.45 SD years in Cuba (p = 0.75). Median survival from onset was 32.43 months (21.93 - 42.36) in Uruguay, 24 months (13.5 - 33.5) in Chile, and 29 months (15 - 42.5) in Cuba (p = 0.006). CONCLUSIONS: These preliminary data from LAENALS confirm the lower incidence and prevalence of ALS in counties with admixed populations. The LAENALS database is now open to other Latin American countries for harmonized prospective data collection.

Idiopathic collapsing glomerulopathy is associated with APOL1 high-risk genotypes or Mendelian variants in most affected individuals in a highly admixed population.

Collapsing glomerulopathy (CG) is most often associated with fast progression to kidney failure with an incidence apparently higher in Brazil than in other countries. However, the reason for this occurrence is unknown. To better understand this, we performed an integrated analysis of clinical, histological, therapeutic, causative genetic and genetic ancestry data in a highly genetically admixed cohort of 70 children and adult patients with idiopathic CG (ICG). The disease onset occurred at 23 (interquartile range: 17-31) years and approximately half of patients progressed to chronic kidney disease requiring kidney replacement therapy (CKD-KRT) 36 months after diagnosis. Causative genetic bases, assessed by targeted-gene panel or whole-exome sequencing, were identified in 58.6% of patients. Among these cases, 80.5% harbored APOL1 high-risk genotypes (HRG) and 19.5% causative Mendelian variants (MV). Self-reported non-White patients more frequently had HRG. MV was an independent risk factor for progression to CKD-KRT by 36 months and the end of follow-up, while remission was an independent protective factor. All patients with HRG manifested CG at 9-44 years of age, whereas in those with APOL1 low-risk genotype, the disease arose throughout life. HRGs were associated with higher proportion of African genetic ancestry. Novel causative MVs were identified in COL4A5, COQ2 and PLCE1 and previously described causative MVs were identified in MYH9, TRPC6, COQ2, COL4A3 and TTC21B. Three patients displayed HRG combined with a variant of uncertain significance (ITGB4, LAMA5 or PTPRO). MVs were associated with worse kidney prognosis. Thus, our data reveal that the genetic status plays a major role in ICG pathogenesis, accounting for more than half of cases in a highly admixed Brazilian population.

Human N-acetyltransferase 2 (NAT2) gene variability in Brazilian populations from different geographical areas.

Introduction: Several polymorphisms altering the NAT2 activity have already been identified. The geographical distribution of NAT2 variants has been extensively studied and has been demonstrated to vary significantly among different ethnic population. Here, we describe the genetic variability of human N-acetyltransferase 2 (NAT2) gene and the predominant genotype-deduced acetylation profiles of Brazilians. Methods: A total of 964 individuals, from five geographical different regions, were genotyped for NAT2 by sequencing the entire coding exon. Results: Twenty-three previously described NAT2 single nucleotide polymorphisms (SNPs) were identified, including the seven most common ones globally (c.191G>A, c.282C>T, c.341T>C, c.481C>T, c.590G>A, c.803A>G and c.857G>A). The main allelic groups were NAT2*5 (36%) and NAT2*6 (18.2%), followed to the reference allele NAT2*4 (20.4%). Combined into genotypes, the most prevalent allelic groups were NAT2*5/*5 (14.6%), NAT2*5/*6 (11.9%) and NAT2*6/*6 (6.2%). The genotype deduced NAT2 slow acetylation phenotype was predominant but showed significant variability between geographical regions. The prevalence of slow acetylation phenotype was higher in the Northeast, North and Midwest (51.3%, 45.5% and 41.5%, respectively) of the country. In the Southeast, the intermediate acetylation phenotype was the most prevalent (40.3%) and, in the South, the prevalence of rapid acetylation phenotype was significantly higher (36.7%), when compared to other Brazilian states (p < 0.0001). Comparison of the predicted acetylation profile among regions showed homogeneity among the North and Northeast but was significantly different when compared to the Southeast (p = 0.0396). The Southern region was significantly different from all other regions (p < 0.0001). Discussion: This study contributes not only to current knowledge of the NAT2 population genetic diversity in different geographical regions of Brazil, but also to the reconstruction of a more accurate phenotypic picture of NAT2 acetylator profiles in those regions.

Facial morphometric differences across face databases: influence of ethnicities and sex.

The scientific need for standardized, high-quality facial stimuli has driven the creation of several face image databases in recent years. These stimuli are particularly important in facial asymmetry research. However, previous studies have reported facial anthropometric differences across a variety of ethnicities. This highlights the need to investigate whether these differences can also impact the use of face image databases, particularly in facial asymmetry research. In this study, we investigated facial asymmetry-based morphometric differences between the multi-ethnic Chicago Face Database (CFD) and the LACOP Face Database, which is composed of Brazilian subjects. We found reliable differences in facial asymmetry between the two databases, which were related to ethnic groups. Specifically, differences in eye and mouth asymmetry seem to drive these differences. The asymmetry-based morphometric differences among databases and ethnicities found in this study reinforce the necessity of creating multi-ethnic face databases.

Association between HLA Class II Alleles/Haplotypes and Genomic Ancestry in Brazilian Patients with Type 1 Diabetes: A Nationwide Exploratory Study.

We aimed to identify HLA-DRB1, -DQA1, and -DQB1 alleles/haplotypes associated with European, African, or Native American genomic ancestry (GA) in admixed Brazilian patients with type 1 diabetes (T1D). This exploratory nationwide study enrolled 1599 participants. GA percentage was inferred using a panel of 46 ancestry informative marker-insertion/deletion. Receiver operating characteristic curve analysis (ROC) was applied to identify HLA class II alleles related to European, African, or Native American GA, and showed significant (p < 0.05) accuracy for identifying HLA risk alleles related to European GA: for DRB1*03:01, the area under the curve was (AUC) 0.533; for DRB1*04:01 AUC = 0.558, for DRB1*04:02 AUC = 0.545. A better accuracy for identifying African GA was observed for the risk allele DRB1*09:01AUC = 0.679 and for the protective alleles DRB1*03:02 AUC = 0.649, DRB1*11:02 AUC = 0.636, and DRB1*15:03 AUC = 0.690. Higher percentage of European GA was observed in patients with risk haplotypes (p < 0.05). African GA percentage was higher in patients with protective haplotypes (p < 0.05). Risk alleles and haplotypes were related to European GA and protective alleles/haplotypes to African GA. Future studies with other ancestry markers are warranted to fill the gap in knowledge regarding the genetic origin of T1D in highly admixed populations such as that found in Brazil.

Host genetic susceptibility underlying SARS-CoV-2-associated Multisystem Inflammatory Syndrome in Brazilian Children.

BACKGROUND: Multisystem Inflammatory Syndrome in Children (MIS-C) is a life-threatening complication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, which manifests as a hyper inflammatory process with multiorgan involvement in predominantly healthy children in the weeks following mild or asymptomatic coronavirus disease 2019 (COVID-19). However, host monogenic predisposing factors to MIS-C remain elusive. METHODS: Herein, we used whole exome sequencing (WES) on 16 MIS-C Brazilian patients to identify single nucleotide/InDels variants as predisposition factors associated with MIS-C. RESULTS: We identified ten very rare variants in eight genes (FREM1, MPO, POLG, C6, C9, ABCA4, ABCC6, and BSCL2) as the most promising candidates to be related to a higher risk of MIS-C development. These variants may propitiate a less effective immune response to infection or trigger the inflammatory response or yet a delayed hyperimmune response to SARS-CoV-2. Protein-Protein Interactions (PPIs) among the products of the mutated genes revealed an integrated network, enriched for immune and inflammatory response mechanisms with some of the direct partners representing gene products previously associated with MIS-C and Kawasaki disease (KD). In addition, the PPIs direct partners are also enriched for COVID-19-related gene sets. HLA alleles prediction from WES data allowed the identification of at least one risk allele in 100% of the MIS-C patients. CONCLUSIONS: This study is the first to explore host MIS-C-associated variants in a Latin American admixed population. Besides expanding the spectrum of MIS-C-associated variants, our findings highlight the relevance of using WES for characterising the genetic interindividual variability associated with COVID-19 complications and ratify the presence of overlapping/convergent mechanisms among MIS-C, KD and COVID-19, crucial for future therapeutic management.

Tierra Del Fuego: What Is Left from the Precolonial Male Lineages?

Similar to other South American regions, Tierra del Fuego has an admixed population characterized by distinct ancestors: Native Americans who first occupied the continent, European settlers who arrived from the late 15th century onwards, and Sub-Saharan Africans who were brought to the Americas for slave labor. To disclose the paternal lineages in the current population from Tierra del Fuego, 196 unrelated males were genotyped for 23 Y-STRs and 52 Y-SNPs. Haplotype and haplogroup diversities were high, indicating the absence of strong founder or drift events. A high frequency of Eurasian haplogroups was detected (94.4%), followed by Native American (5.1%) and African (0.5%) ones. The haplogroup R was the most abundant (48.5%), with the sub-haplogroup R-S116* taking up a quarter of the total dataset. Comparative analyses with other Latin American populations showed similarities with other admixed populations from Argentina. Regarding Eurasian populations, Tierra del Fuego presented similarities with Italian and Iberian populations. In an in-depth analysis of the haplogroup R-M269 and its subtypes, Tierra del Fuego displayed a close proximity to the Iberian Peninsula. The results from this study are in line with the historical records and reflect the severe demographic change led mainly by male newcomers with paternal European origin.

NAToRA, a relatedness-pruning method to minimize the loss of dataset size in genetic and omics analyses.

Genetic and omics analyses frequently require independent observations, which is not guaranteed in real datasets. When relatedness cannot be accounted for, solutions involve removing related individuals (or observations) and, consequently, a reduction of available data. We developed a network-based relatedness-pruning method that minimizes dataset reduction while removing unwanted relationships in a dataset. It uses node degree centrality metric to identify highly connected nodes (or individuals) and implements heuristics that approximate the minimal reduction of a dataset to allow its application to complex datasets. When compared with two other popular population genetics methodologies (PLINK and KING), NAToRA shows the best combination of removing all relatives while keeping the largest possible number of individuals in all datasets tested and also, with similar effects on the allele frequency spectrum and Principal Component Analysis than PLINK and KING. NAToRA is freely available, both as a standalone tool that can be easily incorporated as part of a pipeline, and as a graphical web tool that allows visualization of the relatedness networks. NAToRA also accepts a variety of relationship metrics as input, which facilitates its use. We also release a genealogies simulator software used for different tests performed in this study.

Impact of POR*28 Variant on Tacrolimus Pharmacokinetics in Kidney Transplant Patients with Different CYP3A5 Genotypes.

BACKGROUND: The introduction of tacrolimus (TAC) to clinical practice was essential to the establishment of transplantation as a therapy for patients with chronic renal disease. However, the higher interindividual variation of TAC metabolism has been an important limiting factor for its clinical use. Although the relationship between CYP3A5 polymorphisms and TAC pharmacokinetics (PK) is well established, the effects of other genetic variants on TAC metabolism, such as POR*28, still remain uncertain. OBJECTIVE: The study aimed to evaluate the impact of POR variants on TAC PK in renal transplant patients with different CYP3A5 genotypes (expressers and non-expressers). METHODS: A total of 115 patients were included in this study. Genomic DNA was isolated from peripheral blood, and the real-time PCR technique was used to analyze the polymorphism POR rs1057868; C>T. RESULTS: During the initial post-transplant period, variant allele carriers (*1/*28 and *28/*28) showed a lower TAC dose requirement than POR wild homozygotes (*1/*1). Regarding the influence of the different polymorphisms of POR within the CYP3A5 expresser and non-expresser groups, no differences were observed in any of the PK parameters analyzed during 12 months after transplantation. CONCLUSION: In the studied population, the variant allelic POR*28 was significantly associated with lower TAC dose requirements and higher Co/D ratio in the first-month post-transplant. However, the effects of this polymorphism on the CYP3A5 enzyme activity were not observed.

Indigenous Ancestry and Admixture in the Uruguayan Population.

The Amerindian group known as the Charrúas inhabited Uruguay at the timing of European colonial contact. Even though they were extinguished as an ethnic group as a result of a genocide, Charrúan heritage is part of the Uruguayan identity both culturally and genetically. While mitochondrial DNA studies have shown evidence of Amerindian ancestry in living Uruguayans, here we undertake whole-genome sequencing of 10 Uruguayan individuals with self-declared Charruan heritage. We detect chromosomal segments of Amerindian ancestry supporting the presence of indigenous genetic ancestry in living descendants. Specific haplotypes were found to be enriched in "Charrúas" and rare in the rest of the Amerindian groups studied. Some of these we interpret as the result of positive selection, as we identified selection signatures and they were located mostly within genes related to the infectivity of specific viruses. Historical records describe contacts of the Charrúas with other Amerindians, such as Guaraní, and patterns of genomic similarity observed here concur with genomic similarity between these groups. Less expected, we found a high genomic similarity of the Charrúas to Diaguita from Argentinian and Chile, which could be explained by geographically proximity. Finally, by fitting admixture models of Amerindian and European ancestry for the Uruguayan population, we were able to estimate the timing of the first pulse of admixture between European and Uruguayan indigenous peoples in approximately 1658 and the second migration pulse in 1683. Both dates roughly concurring with the Franciscan missions in 1662 and the foundation of the city of Colonia in 1680 by the Spanish.

Association Analysis of Candidate Variants in Admixed Brazilian Patients With Genetic Generalized Epilepsies.

Genetic generalized epilepsies (GGEs) include well-established epilepsy syndromes with generalized onset seizures: childhood absence epilepsy, juvenile myoclonic epilepsy (JME), juvenile absence epilepsy (JAE), myoclonic absence epilepsy, epilepsy with eyelid myoclonia (Jeavons syndrome), generalized tonic-clonic seizures, and generalized tonic-clonic seizures alone. Genome-wide association studies (GWASs) and exome sequencing have identified 48 single-nucleotide polymorphisms (SNPs) associated with GGE. However, these studies were mainly based on non-admixed, European, and Asian populations. Thus, it remains unclear whether these results apply to patients of other origins. This study aims to evaluate whether these previous results could be replicated in a cohort of admixed Brazilian patients with GGE. We obtained SNP-array data from 87 patients with GGE, compared with 340 controls from the BIPMed public dataset. We could directly access genotypes of 17 candidate SNPs, available in the SNP array, and the remaining 31 SNPs were imputed using the BEAGLE v5.1 software. We performed an association test by logistic regression analysis, including the first five principal components as covariates. Furthermore, to expand the analysis of the candidate regions, we also interrogated 14,047 SNPs that flank the candidate SNPs (1 Mb). The statistical power was evaluated in terms of odds ratio and minor allele frequency (MAF) by the genpwr package. Differences in SNP frequencies between Brazilian and Europeans, sub-Saharan African, and Native Americans were evaluated by a two-proportion Z-test. We identified nine flanking SNPs, located on eight candidate regions, which presented association signals that passed the Bonferroni correction (rs12726617; rs9428842; rs1915992; rs1464634; rs6459526; rs2510087; rs9551042; rs9888879; and rs8133217; p-values <3.55e-06). In addition, the two-proportion Z-test indicates that the lack of association of the remaining candidate SNPs could be due to different genomic backgrounds observed in admixed Brazilians. This is the first time that candidate SNPs for GGE are analyzed in an admixed Brazilian population, and we could successfully replicate the association signals in eight candidate regions. In addition, our results provide new insights on how we can account for population structure to improve risk stratification estimation in admixed individuals.

Evaluation of Promega PowerSeq™ Auto/Y systems prototype on an admixed sample of Rio de Janeiro, Brazil: Population data, sensitivity, stutter and mixture studies.

Forensic DNA typing typically relies on the length-based (LB) separation of PCR products containing short tandem repeat loci (STRs). Massively parallel sequencing (MPS) elucidates an additional level of STR motif and flanking region variation. Also, MPS enables simultaneous analysis of different marker-types - autosomal STRs, SNPs for lineage and identification purposes, reducing both the amount of sample used and the turn-around-time of analysis. Therefore, MPS methodologies are being considered as an additional tool in forensic genetic casework. The PowerSeq™ Auto/Y System (Promega Corp), a multiplex forensic kit for MPS, enables analysis of the 22 autosomal STR markers (plus Amelogenin) from the PowerPlex® Fusion 6C kit and 23 Y-STR markers from the PowerPlex® Y23 kit. Population data were generated from 140 individuals from an admixed sample from Rio de Janeiro, Brazil. All samples were processed according to the manufacturers' recommended protocols. Raw data (FastQ) were generated for each indexed sample and analyzed using STRait Razor v2s and PowerSeqv2.config file. The subsequent population data showed the largest increase in expected heterozygosity (23%), from LB to sequence-based (SB) analyses at the D5S818 locus. Unreported allele was found at the D21S11 locus. The random match probability across all loci decreased from 5.9 × 10-28 to 7.6 × 10-33. Sensitivity studies using 1, 0.25, 0.062 and 0.016 ng of DNA input were analyzed in triplicate. Full Y-STR profiles were detected in all samples, and no autosomal allele drop-out was observed with 62 pg of input DNA. For mixture studies, 1 ng of genomic DNA from a male and female sample at 1:1, 1:4, 1:9, 1:19 and 1:49 proportions were analyzed in triplicate. Clearly resolvable alleles (i.e., no stacking or shared alleles) were obtained at a 1:19 male to female contributor ratio. The minus one stutter (-1) increased with the longest uninterrupted stretch (LUS) allele size reads and according to simple or compound/complex repeats. The haplotype-specific stutter rates add more information for mixed samples interpretation. These data support the use of the PowerSeqTM Auto/Y systems prototype kit (22 autosomal STR loci, 23 Y-STR loci and Amelogenin) for forensic genetics applications.

Population structure and relatedness estimates in a Mexican sample.

Population stratification (PS) is a confounding factor in genome-wide association studies (GWASs) and also an interesting process itself. Latin American populations have mixed genetic ancestry, which may account for PS. We have analyzed the relatedness, by means of the identity-by-descent (IBD) estimations, in a sample of 1805 individuals and 1.006.703 autosomal mutations from a case-control study of colorectal cancer in Mexico. When using the recommended protocol for quality control assessment, 402 should have been removed due to relatedness. Our purpose was to analyze this value in the context of an admixed population. For that aim, we reanalyzed the sample using two software designed for admixed populations, obtaining estimates of 110 and 70 related individuals to remove. The results showed that the first estimation of relatedness was an effect of the higher Native American contribution in part of the data samples, being a confounding factor for IBD estimations. We conclude in the importance of considering PS and genetic ancestry in order to avoid spurious results, not only in GWAS but also in relatedness analysis.