Marília Oliveira Scliar

Osteosarcoma (OS) is the most prevalent type of bone tumor, but slow progress has been achieved in disentangling the full set of genomic events involved in its initiation and progression. We assessed by NGS the mutational spectrum of 28 primary OSs from Brazilian patients, and identified 445 potentially deleterious SNVs/indels and 1176 copy number alterations (CNAs). TP53 was the most recurrently mutated gene, with an overall rate of ~60%, considering SNVs/indels and CNAs. The most frequent CNAs (~60%) were gains at 1q21.2q21.3, 6p21.1, and 8q13.3q24.22, and losses at 10q26 and 13q14.3q21.1. Seven cases presented CNA patterns reminiscent of complex events (chromothripsis and chromoanasynthesis). Putative RB1 and TP53 germline variants were found in five samples associated with metastasis at diagnosis along with complex genomic patterns of CNAs. PTPRQ, KNL1, ZFHX4, and DMD alterations were prevalent in metastatic or deceased patients, being potentially indicative of poor prognosis. TNFRSF11B, involved in skeletal system development and maintenance, emerged as a candidate for osteosarcomagenesis due to its biological function and a high frequency of copy number gains. A protein–protein network enrichment highlighted biological pathways involved in immunity and bone development. Our findings reinforced the high genomic OS instability and heterogeneity, and led to the identification of novel disrupted genes deserving further evaluation as biomarkers due to their association with poor outcomes.

Human genomics has quickly evolved, powering genome-wide association studies (GWASs). SNP-based GWASs cannot capture the intense polymorphism of HLA genes, highly associated with disease susceptibility. There are methods to statistically impute HLA genotypes from SNP-genotypes data, but lack of diversity in reference panels hinders their performance. We evaluated the accuracy of the 1,000 Genomes data as a reference panel for imputing HLA from admixed individuals of African and European ancestries, focusing on (a) the full dataset, (b) 10 replications from 6 populations, (c) 19 conditions for the custom reference panels. The full dataset outperformed smaller models, with a good F1-score of 0.66 for HLA-B . However, custom models outperformed the multiethnic or population models of similar size (F1-scores up to 0.53, against up to 0.42). We demonstrated the importance of using genetically specific models for imputing admixed populations, which are currently underrepresented in public datasets, opening the door to HLA imputation for every genetic population.

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.

Dementia is more prevalent in Blacks than in Whites, likely due to a combination of environmental and biological factors. Paradoxically, clinical studies suggest an attenuation of APOE ε4 risk of dementia in African ancestry (AFR), but lack of neuropathological data preclude the interpretation of the biological factors underlying these findings, including the association between APOE ε4 risk and Alzheimer’s disease (AD) pathology, the most frequent cause of dementia. We investigated the interaction between African ancestry, AD-related neuropathology, APOE genotype , and functional cognition in a postmortem sample of 400 individuals with a range of AD pathology severity and lack of comorbid neuropathology from a cohort of community-dwelling, admixed Brazilians. Increasing proportions of African ancestry (AFR) correlated with a lower burden of neuritic plaques (NP). However, for individuals with high levels of NP and neurofibrillary tangles (NFT), AFR proportion was associated with worse Clinical Dementia Rating sum of boxes (CDR-SOB). Among APOE ε4 carriers, the association between AFR proportion and CDR-SOB disappeared. APOE local ancestry inference of a subset of 309 individuals revealed that, in APOE ε4 noncarriers, non-European APOE background associated with lower NP burden, but with worst cognitive outcomes compared to European APOE when adjusting by the similar NP burden. Finally, APOE ε4 was associated with worse AD neuropathological burden only in a European APOE background. APOE genotype and its association with AD neuropathology and clinical pattern are highly influenced by ancestry, with AFR associated with lower NP burden and attenuated APOE4 risk compared to European ancestry.

Background Despite the high number of individuals infected by SARS-CoV-2 who develop COVID-19 symptoms worldwide, many exposed individuals remain asymptomatic and/or stay uninfected. This could be explained by a combination of environmental (exposure, previous infection), epigenetic, and genetic factors. Aiming to identify genetic variants involved in SARS-CoV-2 resistance, we analyzed 86 discordant Brazilian couples where one was infected and symptomatic while the partner remained asymptomatic and seronegative despite sharing the same bedroom during the infection. The discordant partners had comparable ages, and genetic ancestry proportions. Methods Whole-exome sequencing followed by a state-of-the-art method to call genotypes and haplotypes across the highly polymorphic MHC and LRC. Results We observed a minor impact in antigen-presentation genes and KIR genes associated with resistance. Interestingly, genes related to immune modulation, mainly involved in NK cell killing activation/inhibition harbor variants potentially contributing to infection resistance. We hypothesize that individuals prone to produce higher amounts of MICA (possibly soluble), LILRB1, LILRB2, and low amounts of MICB, would be more susceptible to infection. Conclusion According to this hypothesis, quantitative differences in these NK activity-related molecules could contribute to resistance to COVID-19 down regulating NK cell cytotoxic activity in infected individuals but not in resistant partners.

Background: Clinical recurrence of COVID-19 in convalescent patients has been reported, which immune mechanisms have not been thoroughly investigated. Presence of neutralizing antibodies suggests other types of immune response are involved. Methods: We assessed the innate type I/III IFN response, T cell responses to SARS-CoV-2 with IFNγ ELISPOT, binding and neutralizing antibody assays, in two monozygotic twin pairs with one COVID-19 recurrence case. Results: In pair 1, four months after a first mild episode of infection for both siblings, one displayed severe clinical recurrence of COVID-19. Twin pair 2 of siblings underwent non-recurring asymptomatic infection. All fours individuals presented similar overall responses, except for remarkably difference found in specific cellular responses. Recurring sibling presented a reduced number of recognized T cell epitopes as compared to the other three including her non-recurring sibling. Conclusions: Our results suggest that an effective SARS-CoV-2-specific T cell immune response is key for complete viral control and avoidance of clinical recurrence of COVID-19. Besides, adaptive immunity can be distinct in MZ twins. Given the rising concern about SARS-CoV-2 variants that evade neutralizing antibodies elicited by vaccination or infection, our study stresses the importance of T cell responses in protection against recurrence/reinfection. Key points: Immune parameters leading to COVID-19 recurrence/reinfection are incompletely understood. A COVID-19 recurrence case in a monozygotic twin pair is described with an intact antibody and innate type I/III Interferon response and drastically reduced number of recognized SARS-CoV-2 T cell epitopes.

As whole-genome sequencing (WGS) becomes the gold standard tool for studying population genomics and medical applications, data on diverse non-European and admixed individuals are still scarce. Here, we present a high-coverage WGS dataset of 1,171 highly admixed elderly Brazilians from a census-based cohort, providing over 76 million variants, of which ~2 million are absent from large public databases. WGS enabled identifying ~2,000 novel mobile element insertions, nearly 5Mb of genomic segments absent from human genome reference, and over 140 novel alleles from HLA genes. We reclassified and curated nearly four hundred variant's pathogenicity assertions in genes associated with dominantly inherited Mendelian disorders and calculated the incidence for selected recessive disorders, demonstrating the clinical usefulness of the present study. Finally, we observed that whole-genome and HLA imputation could be significantly improved compared to available datasets since rare variation represents the largest proportion of input from WGS. These results demonstrate that even smaller sample sizes of underrepresented populations bring relevant data for genomic studies, especially when exploring analyses allowed only by WGS.

Western South America was one of the worldwide cradles of civilization. The well known Inca Empire was thetip of the icebergof a cultural and biological evolutionary process that started 14-11 thousand years ago. Genetic data from 18 Peruvian populations reveal that: (1) The between-population homogenization of the central-southern Andes and its differentiation with respect to Amazonian populations of similar latitudes do not extend northward. Instead, longitudinal gene flow between the northern coast of Peru, Andes and Amazonia accompanied cultural and socioeconomic interactions revealed by archeological studies. This pattern recapitulates the environmental and cultural differentiation between the fertile north, where altitudes are lower; and the arid south, where the Andes are higher, acting as a genetic barrier between the sharply different environments of the Andes and Amazonia (2). The genetic homogenization between the populations of the arid Andes is not only due to migration during the Inca Empire or the subsequent colonial period. It started at least during the earlier expansion of the pre-Inca Wari Empire (600-1000 YBP) (3) This demographic history allowed for cases of positive natural selection in the high and arid Andes vs. the low Amazon tropical forest: in the Andes,HAND2-AS1(heart and neural crest derivatives expressed 2 antisense RNA1, related with cardiovascular function) andDUOX2(dual oxidase 2, related to thyroid function and innate immunity) genes; in the Amazon, the gene encoding for the CD45 protein, essential for antigen recognition by T/B lymphocytes in viral-host interaction, consistent with thehost-virus arms racehypothesis.

Admixed populations are a resource to study the global genetic architecture of complex phenotypes, which is critical, considering that non-European populations are severely under-represented in genomic studies. Leveraging admixture in Brazilians, whose chromosomes are mosaics of fragments of Native American, European and African origins, we used genome-wide data to perform admixture mapping/fine-mapping of Body Mass Index (BMI) in three population-based cohorts from Northeast (Salvador), Southeast (Bambuí) and South (Pelotas) of the country. We found significant associations with African-associated alleles in children from Salvador (PALD1 and ZMIZ1 genes), and in young adults from Pelotas (NOD2 and MTUS2 genes). More importantly, in Pelotas, rs114066381, mapped in a potential regulatory region, is significantly associated only in females (p= 2.76 e-06). This variant is very rare in Europeans but with frequencies of ~3% in West Africa, and has a strong female-specific effect (95%CI: 2.32-5.65 kg/m2 per each A allele). We confirmed this sex-specific association and replicated its strong effect for an adjusted fat-mass index in the same Pelotas cohort, and for BMI in another Brazilian cohort from São Paulo (Southeast Brazil). A meta-analysis confirmed the significant association. Remarkably, we observed that while the frequency of rs114066381-A allele ranges from 0.8 to 2.1% in the studied populations, it attains ~9% among morbidly obese women from Pelotas, São Paulo, and Bambuí. The effect size of rs114066381 is at least five-times the effect size of the FTO SNPs rs9939609 and rs1558902, already emblematic for their high effects, and for which we replicated associations in Pelotas. We demonstrate how, after a decade of GWAS mostly performed in European-ancestry populations, non-European and admixed populations are a source of new relevant phenotype-associated genetic variants.