Whole-genome sequencing of 1,171 elderly admixed individuals from São Paulo, Brazil

Assessing the risk stratification of breast cancer polygenic risk scores in a Brazilian cohort

Polygenic risk scores (PRS) for breast cancer (BC) have a clear clinical utility in risk prediction. PRS transferability across populations and ancestry groups is hampered by population-specific factors, ultimately leading to differences in variant effects, such as linkage disequilibrium (LD) and differences in variant frequency (AF-diff). Thus, locally-sourced population-based phenotypic and genomic datasets are essential to assess the validity of PRS derived from signals detected across populations. Here, we assess the transferability of a BC PRS composed of 313 risk variants (313-PRS) in a Brazilian tri-hybrid admixed ancestries (European, African and Native American) whole-genome sequenced cohort, GRAR. We computed 313-PRS in GRAR (n=853) using the UK Biobank (UKBB, n=264,307) as reference. We show that although the Brazilian cohorts have a high European (EA) component, with AF-diff and to a lesser extent LD patterns similar to those found in EA populations, the 313-PRS distribution is inflated when compared to that of the UKBB, leading to potential overestimation of PRS-based risk if EA is taken as a standard. Interestingly, we find that case-controls lead to equivalent predictive power when compared to UKBB-EA samples with AUROC values of 0.66-0.62 compared to 0.63 for UKBB.

Misclassification of a frequent variant from PMS2CL pseudogene as a PMS2 loss of function variant in Brazilian patients

PMS2, a Lynch Syndrome gene, presents challenges in genetic testing due to the existence of multiple pseudogenes. This study aims to describe a series of cases harboring a variant in the PMS2CL pseudogene that has been incorrectly assigned to PMS2 with different nomenclatures. We reviewed data from 647 Brazilian patients who underwent multigene genetic testing at a single center to identify those harboring the PMS2 V1:c.2186_2187delTC or V2:c.2182_2184delACTinsG variants, allegedly located at PMS2 exon 13. Gene-specific PCR and transcript sequencing was performed. Among the 647 individuals, 1.8% (12) carried the investigated variants, with variant allele frequencies ranging from 15 to 34%. By visually inspecting the alignments, we confirmed that both V1 and V2 represented the same variant and through gene-specific PCR and PMS2 transcript analysis, we demonstrated that V1/V2 is actually located in the PMS2CL pseudogene. Genomic databases (ExAC and gnomAD) report an incidence of 2.5 - 5.3% of this variant in the African population. Currently, V1 is classified as "uncertain significance" and V2 as "conflicting" in ClinVar, with several laboratories classifying them as "pathogenic". We identified a frequent African PMS2CL variant in the Brazilian population that is misclassified as a PMS2 variant. It is likely that V1/V2 have been erroneously assigned to PMS2 in several manuscripts and by clinical laboratories, underscoring a disparity-induced matter. Considering the limitations of short-read NGS differentiating between certain regions of PMS2 and PMS2CL, using complementary methodologies is imperative to provide an accurate diagnosis.

Single-Base Gene Variants in MIR-146A and SCN1A Genes Related to the Epileptogenic Process in Drug-Responsive and Drug-Resistant Temporal Lobe Epilepsy-A Preliminary Study in a Brazilian Cohort Sample

The drug-resistant temporal lobe epilepsy (TLE) has recently been associated with single nucleotide variants (SNVs) in microRNA(miR)-146a (MIR-146A) (rs2910164) and Sodium Voltage-Gated Channel Alpha Subunit 1 (SCN1A) (rs2298771 and rs3812718) genes. Moreover, no studies have shown an association between these SNVs and susceptibility to drug-resistant and drug-responsive TLE in Brazil. Thus, deoxyribonucleic acid (DNA) samples from 120 patients with TLE (55 drug-responsive and 65 drug-resistant) were evaluated by real-time polymerase chain reaction (RT-PCR). A total of 1171 healthy blood donor individuals from the Online Archive of Brazilian Mutations (ABraOM, from Portuguese Arquivo Brasileiro On-line de Mutações), a repository containing genomic variants of the Brazilian population, were added as a control population for the studied SNVs. MIR-146A and SCN1A relative expression was performed by quantitative RT-PCR (qRT-PCR). The statistical analysis protocol was performed using an alpha error of 0.05. TLE patient samples and ABraOM control samples were in Hardy-Weinberg equilibrium for all studied SNVs. For rs2910164, the frequencies of the homozygous genotype (CC) (15.00% vs. 9.65%) and C allele (37.80% vs. 29.97%) were superior in patients with TLE compared to controls with a higher risk for TLE disease [odds ratio (OR) = 1.89 (95% confidence interval (95%CI) = 1.06-3.37); OR = 1.38 (95%CI = 1.04-1.82), respectively]. Drug-responsive patients also presented higher frequencies of the CC genotype [21.81% vs. 9.65%; OR = 2.58 (95%CI = 1.25-5.30)] and C allele [39.09% vs. 29.97%; OR = 1.50 (95%CI = 1.01-2.22)] compared to controls. For rs2298771, the frequency of the heterozygous genotype (AG) (51.67% vs. 40.40%) was superior in patients with TLE compared to controls with a higher risk for TLE disease [OR = 2.42 (95%CI = 1.08-5.41)]. Drug-resistant patients presented a higher AG frequency [56.92% vs. 40.40%; OR = 3.36 (95%CI = 1.04-17.30)] compared to the control group. For rs3812718, the prevalence of genotypes and alleles were similar in both studied groups. The MIR-146A relative expression level was lower in drug-resistant compared to drug-responsive patients for GC (1.6 vs. 0.1, p-value = 0.049) and CC (1.8 vs. 0.6, p-value = 0.039). Also, the SCN1A relative expression levels in samples from TLE patients were significantly higher in AG [2.09 vs. 1.10, p-value = 0.038] and GG (3.19 vs. 1.10, p-value < 0.001) compared to the AA genotype. In conclusion, the rs2910164-CC and rs2298771-AG genotypes are exerting significant risk influence, respectively, on responsive disease and resistant disease, probably due to an upregulated nuclear factor kappa B (NF-kB) and SCN1A loss of function.

Assessing the Occurrence and Influence of Cancer Chemotherapy-Related Pharmacogenetic Alleles in the Chilean Population

BACKGROUND

Pharmacogenomic knowledge as a biomarker for cancer care has transformed clinical practice, however, as current guidelines are primarily derived from Eurocentric populations, this limits their application in Latin America, particularly among Hispanic or Latino groups. Despite advancements, systemic chemotherapy still poses challenges in drug toxicity and suboptimal response. This study explores pharmacogenetic markers related to anticancer drugs in a Chilean cohort, filling a gap in Latin American research. Notably, the influence of native South American Mapuche-Huilliche ancestry.

METHODS

To explore pharmacogenetic markers related to anticancer drugs, we utilized an ethnically Admixed Chilean genome-wide association studies (GWAS) dataset of 1095 unrelated individuals. Pharmacogenomic markers were selected from PharmGKB, totaling 36 level 1 and 2 evidence single nucleotide polymorphisms (SNPs) and 571 level 3 SNPs. Comparative analyses involved assessing SNP frequencies across diverse populations from the 1000 Genomes Project. Haplotypes were estimated, and linkage disequilibrium was examined. Ancestry-based association analyses explored relationships between SNPs and Mapuche-Huilliche and European ancestries. Chi-square distribution with p ≤ 0.05 and Bonferroni's multiple adjustment tests determined statistical differences between allele frequencies.

RESULTS

Our study reveals significant disparities in SNP frequency within the Chilean population. Notably, dihydropyrimidine dehydrogenase (DPYD) variants (rs75017182 and rs67376798), linked to an increased risk of severe fluoropyrimidine toxicity, exhibit an exceptionally low frequency (minor allele frequency (MAF) < 0.005). Nudix hydrolase 15 (NUDT15) rs116855232, associated with hematological mercaptopurine toxicity, is relatively common (MAF = 0.062), and is further linked to Mapuche-Huilliche ancestry. Thiopurine methyltransferase enzyme (TPMT), implicated in severe toxicity to mercaptopurines, SNPs rs1142345 and rs1800460 of TMPT gene demonstrate higher MAFs in Admixed Americans and the Chilean population (MAF range 0.031-0.057). Finally, the variant in the UDP-glucuronosyltransferase 1 gene (UGT1A1) rs4148323, correlated with irinotecan neutropenia, exhibits the highest MAF in East Asian (MAF = 0.136) and Chilean (MAF = 0.025) populations, distinguishing them from other investigated populations.

CONCLUSIONS

This study provides the first comprehensive pharmacogenetic characterization of cancer therapy-related SNPs and highlights significant disparities in SNP frequencies within the Chilean population. Our findings underscore the necessity for inclusive research and personalized therapeutic strategies to ensure the equitable and effective application of precision medicine across diverse global communities.

Comparing the effect of imputation reference panel composition in four distinct Latin American cohorts

Genome-wide association studies have been useful in identifying genetic risk factors for various phenotypes. These studies rely on imputation and many existing panels are largely composed of individuals of European ancestry, resulting in lower levels of imputation quality in underrepresented populations. We aim to analyze how the composition of imputation reference panels affects imputation quality in four target Latin American cohorts. We compared imputation quality for chromosomes 7 and X when altering the imputation reference panel by: 1) increasing the number of Latin American individuals; 2) excluding either Latin American, African, or European individuals, or 3) increasing the Indigenous American (IA) admixture proportions of included Latin Americans. We found that increasing the number of Latin Americans in the reference panel improved imputation quality in the four populations; however, there were differences between chromosomes 7 and X in some cohorts. Excluding Latin Americans from analysis resulted in worse imputation quality in every cohort, while differential effects were seen when excluding Europeans and Africans between and within cohorts and between chromosomes 7 and X. Finally, increasing IA-like admixture proportions in the reference panel increased imputation quality at different levels in different populations. The difference in results between populations and chromosomes suggests that existing and future reference panels containing Latin American individuals are likely to perform differently in different Latin American populations.

Genomic ancestry and cancer among Latin Americans

Latin American populations, characterized by intricate admixture patterns resulting from the intermingling of ancestries from European, Native American (NA) Asian, and African ancestries which result in a vast and complex genetic landscape, harboring unique combinations of novel variants. This genetic diversity not only poses challenges in traditional population genetics methods but also opens avenues for a deeper understanding of its implications in health. In cancer, the interplay between genetic ancestry, lifestyle factors, and healthcare disparities adds a layer of complexity to the varying incidence and mortality rates observed across different Latin American subpopulations. This complex interdependence has been unveiled through numerous studies, whether conducted on Latin American patients residing on the continent or abroad, revealing discernible differences in germline composition that influence divergent disease phenotypes such as higher incidence of Luminal B and Her2 breast tumors, EGFR and KRAS mutated lung adenocarcinomas in addition to an enrichment in BRCA1/2 pathogenic variants and a higher than expected prevalence of variants in colorectal cancer associated genes such as APC and MLH1. In prostate cancer novel risk variants have also been solely identified in Latin American populations. Due to the complexity of genetic divergence, inputs from each individual ancestry seem to carry independent contributions that interplay in the development of these complex disease phenotypes. By understanding these unique population characteristics, genomic ancestries hold a promising avenue for tailoring prognostic assessments and optimizing responses to oncological interventions.

Germline mutations in cancer predisposition genes among pediatric patients with cancer and congenital anomalies

BACKGROUND

Childhood cancer has a poorly known etiology, and investigating the underlying genetic background may provide novel insights. A recognized association exists between non-chromosomal birth defects and childhood cancer susceptibility.

METHODS

We performed whole-exome sequencing and chromosomal microarray analysis in a cohort of childhood cancer (22 individuals, 50% with congenital anomalies) to unravel deleterious germline variants.

RESULTS

A diagnostic yield of 14% was found, encompassing heterozygous variants in bona fide dominant Cancer Predisposition Genes (CPGs). Considering candidate and recessive CPGs harboring monoallelic variants, which were also deemed to play a role in the phenotype, the yield escalated to 45%. Most of the deleterious variants were mapped in genes not conventionally linked to the patient's tumor type. Relevant findings were detected in 55% of the syndromic individuals, mostly variants potentially underlying both phenotypes.

CONCLUSION

We uncovered a remarkable prevalence of germline deleterious CPG variants, highlighting the significance of a comprehensive genetic analysis in pediatric cancer, especially when coupled with additional clinical signs. Moreover, our findings emphasized the potential for oligogenic inheritance, wherein multiple genes synergistically increase cancer risk. Lastly, our investigation unveiled potentially novel genotype-phenotype associations, such as SETD5 in neuroblastoma, KAT6A in gliomas, JAG1 in hepatoblastomas, and TNFRSF13B in Langerhans cell histiocytosis.

IMPACT

Novel gene-phenotype associations and candidate genes for pediatric cancer were unraveled, such as KAT6A in gliomas, SETD5 in neuroblastoma, JAG1 in hepatoblastomas, and TNFRSF13B in Langerhans cell histiocytosis. Our analysis revealed a high frequency of deleterious germline variants, particularly in cases accompanied by additional clinical signs, highlighting the importance of a comprehensive genetic evaluation in childhood cancer. Our findings also underscored the potential for oligogenic inheritance in pediatric cancer risk. Understanding the cancer etiology is crucial for genetic counseling, often influencing therapeutic decisions and offering valuable insights into molecular targets for the development of oncological therapies.

Variants in Candidate Genes for Phenotype Heterogeneity in Patients with the 22q11.2 Deletion Syndrome

22q11.2 deletion syndrome (22q11.2DS) is a microdeletion syndrome with a broad and heterogeneous phenotype, even though most of the deletions present similar sizes, involving ∼3 Mb of DNA. In a relatively large population of a Brazilian 22q11.2DS cohort (60 patients), we investigated genetic variants that could act as genetic modifiers and contribute to the phenotypic heterogeneity, using a targeted NGS (Next Generation Sequencing) with a specific Ion AmpliSeq panel to sequence nine candidate genes (CRKL, MAPK1, HIRA, TANGO2, PI4KA, HDAC1, ZDHHC8, ZFPM2, and JAM3), mapped in and outside the 22q11.2 hemizygous deleted region. In silico prediction was performed, and the whole-genome sequencing annotation analysis package (WGSA) was used to predict the possible pathogenic effect of single nucleotide variants (SNVs). For the in silico prediction of the indels, we used the genomic variants filtered by a deep learning model in NGS (GARFIELD-NGS). We identified six variants, 4 SNVs and 2 indels, in MAPK1, JAM3, and ZFPM2 genes with possibly synergistic deleterious effects in the context of the 22q11.2 deletion. Our results provide the opportunity for the discovery of the co-occurrence of genetic variants with 22q11.2 deletions, which may influence the patients´ phenotype.

A comparative medical genomics approach may facilitate the interpretation of rare missense variation

PURPOSE

To determine the degree to which likely causal missense variants of single-locus traits in domesticated species have features suggestive of pathogenicity in a human genomic context.

METHODS

We extracted missense variants from the Online Mendelian Inheritance in Animals database for nine animals (cat, cattle, chicken, dog, goat, horse, pig, rabbit and sheep), mapped coordinates to the human reference genome and annotated variants using genome analysis tools. We also searched a private commercial laboratory database of genetic testing results from >400 000 individuals with suspected rare disorders.

RESULTS

Of 339 variants that were mappable to the same residue and gene in the human genome, 56 had been previously classified with respect to pathogenicity: 31 (55.4%) pathogenic/likely pathogenic, 1 (1.8%) benign/likely benign and 24 (42.9%) uncertain/other. The odds ratio for a pathogenic/likely pathogenic classification in ClinVar was 7.0 (95% CI 4.1 to 12.0, p<0.0001), compared with all other germline missense variants in these same 220 genes. The remaining 283 variants disproportionately had allele frequencies and REVEL scores that supported pathogenicity.

CONCLUSION

Cross-species comparisons could facilitate the interpretation of rare missense variation. These results provide further support for comparative medical genomics approaches that connect big data initiatives in human and veterinary genetics.

Hereditary breast cancer next-generation sequencing (NGS) panel evaluation in the south region of Brazil: a novel BRCA2 candidate pathogenic variant is reported

In this article, we delineate a loosely selected cohort comprising patients with a history of early-onset breast cancer and/or a familial occurrence of cancer. The aim of this study was to gain insights into the presence of breast cancer-related gene variants in a population from a micro-region in southern Brazil, specifically the Metropolitan Region of Curitiba. This area exhibits a highly genetically mixed population, mirroring the general characteristics of the Brazilian people. Comprehensive next-generation sequencing (NGS) multigene panel testing was conducted, involving the evaluation of twelve patients. Two pathogenic variants and one candidate pathogenic variant were identified: BRCA2:c.8878C>T, p.Gln2960Ter; CHEK2:c.1100delAG>A, p.Thr367Metfs*15 and BRCA2:c.3482dupG>GA, p.Asp1161Glufs*3, a novel variant, previously unpublished, is reported.

Optimal population-specific HLA imputation with dimension reduction

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 1000 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, and (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 populations, which are currently underrepresented in public datasets, opening the door to HLA imputation for every genetic population.

BRAF-Mutated Melanoma Journey in Latin America: Expert Recommendations From Diagnosis to Treatment

OBJECTIVES

• Gather a panel of Latin American experts in testing and treating BRAF-melanoma. • Describe the current landscape of BRAF-mutated melanoma in Latin America. • Outline the current gaps in testing and recommend improvements for testing and treating BRAF-mutated melanoma in the region.

INTRODUCTION

Melanoma prevalence in Latin America is lower than in high- and middle-income countries. However, recent data indicate that the region's incidence and mortality are rising, with more stage IV patients being diagnosed. According to international clinical practice guidelines, conducting BRAF-mutation testing in patients with stage III or stage IV melanoma and high-risk resected disease is imperative. Still, BRAF-mutation testing and targeted therapies are inconsistently available in the region.

METHODS

Americas Health Foundation convened a meeting of Latin American experts on BRAF-mutated melanoma to develop guidelines and recommendations for diagnosis through treatment.

RESULTS AND CONCLUSIONS

Some recommendations for improving diagnostics through improving access and reducing the cost of BRAF-mutation testing, enhancing efficiency in pathology laboratories, and creating country-specific local guidelines. The panel also gave treatment recommendations for neo-adjuvant therapy, adjuvant therapy, and therapy for patients with metastatic disease in Latin America.

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.

Whole genomes from Angola and Mozambique inform about the origins and dispersals of major African migrations

As the continent of origin for our species, Africa harbours the highest levels of diversity anywhere on Earth. However, many regions of Africa remain under-sampled genetically. Here we present 350 whole genomes from Angola and Mozambique belonging to ten Bantu ethnolinguistic groups, enabling the construction of a reference variation catalogue including 2.9 million novel SNPs. We investigate the emergence of Bantu speaker population structure, admixture involving migrations across sub-Saharan Africa and model the demographic histories of Angolan and Mozambican Bantu speakers. Our results bring together concordant views from genomics, archaeology, and linguistics to paint an updated view of the complexity of the Bantu Expansion. Moreover, we generate reference panels that better represents the diversity of African populations involved in the trans-Atlantic slave trade, improving imputation accuracy in African Americans and Brazilians. We anticipate that our collection of genomes will form the foundation for future African genomic healthcare initiatives.

Neuromuscular disease genetics in under-represented populations: increasing data diversity

Neuromuscular diseases (NMDs) affect ∼15 million people globally. In high income settings DNA-based diagnosis has transformed care pathways and led to gene-specific therapies. However, most affected families are in low-to-middle income countries (LMICs) with limited access to DNA-based diagnosis. Most (86%) published genetic data is derived from European ancestry. This marked genetic data inequality hampers understanding of genetic diversity and hinders accurate genetic diagnosis in all income settings. We developed a cloud-based transcontinental partnership to build diverse, deeply-phenotyped and genetically characterized cohorts to improve genetic architecture knowledge, and potentially advance diagnosis and clinical management. We connected 18 centres in Brazil, India, South Africa, Turkey, Zambia, Netherlands and the UK. We co-developed a cloud-based data solution and trained 17 international neurology fellows in clinical genomic data interpretation. Single gene and whole exome data were analysed via a bespoke bioinformatics pipeline and reviewed alongside clinical and phenotypic data in global webinars to inform genetic outcome decisions. We recruited 6001 participants in the first 43 months. Initial genetic analyses 'solved' or 'possibly solved' ∼56% probands overall. In-depth genetic data review of the four commonest clinical categories (limb girdle muscular dystrophy, inherited peripheral neuropathies, congenital myopathy/muscular dystrophies and Duchenne/Becker muscular dystrophy) delivered a ∼59% 'solved' and ∼13% 'possibly solved' outcome. Almost 29% of disease causing variants were novel, increasing diverse pathogenic variant knowledge. Unsolved participants represent a new discovery cohort. The dataset provides a large resource from under-represented populations for genetic and translational research. In conclusion, we established a remote transcontinental partnership to assess genetic architecture of NMDs across diverse populations. It supported DNA-based diagnosis, potentially enabling genetic counselling, care pathways and eligibility for gene-specific trials. Similar virtual partnerships could be adopted by other areas of global genomic neurological practice to reduce genetic data inequality and benefit patients globally.

The hospital Israelita Albert Einstein standards for constitutional sequence variants classification: version 2023

BACKGROUND

Next-generation sequencing has had a significant impact on genetic disease diagnosis, but the interpretation of the vast amount of genomic data it generates can be challenging. To address this, the American College of Medical Genetics and Genomics and the Association for Molecular Pathology have established guidelines for standardized variant interpretation. In this manuscript, we present the updated Hospital Israelita Albert Einstein Standards for Constitutional Sequence Variants Classification, incorporating modifications from leading genetics societies and the ClinGen initiative.

RESULTS

First, we standardized the scientific publications, documents, and other reliable sources for this document to ensure an evidence-based approach. Next, we defined the databases that would provide variant information for the classification process, established the terminology for molecular findings, set standards for disease-gene associations, and determined the nomenclature for classification criteria. Subsequently, we defined the general rules for variant classification and the Bayesian statistical reasoning principles to enhance this process. We also defined bioinformatics standards for automated classification. Our workgroup adhered to gene-specific rules and workflows curated by the ClinGen Variant Curation Expert Panels whenever available. Additionally, a distinct set of specifications for criteria modulation was created for cancer genes, recognizing their unique characteristics.

CONCLUSIONS

The development of an internal consensus and standards for constitutional sequence variant classification, specifically adapted to the Brazilian population, further contributes to the continuous refinement of variant classification practices. The aim of these efforts from the workgroup is to enhance the reliability and uniformity of variant classification.

Whole-Genome Sequencing of 502 Individuals from Latvia: The First Step towards a Population-Specific Reference of Genetic Variation

Despite rapid improvements in the accessibility of whole-genome sequencing (WGS), understanding the extent of human genetic variation is limited by the scarce availability of genome sequences from underrepresented populations. Developing the population-scale reference database of Latvian genetic variation may fill the gap in European genomes and improve human genomics research. In this study, we analysed a high-coverage WGS dataset comprising 502 individuals selected from the Genome Database of the Latvian Population. An assessment of variant type, location in the genome, function, medical relevance, and novelty was performed, and a population-specific imputation reference panel (IRP) was developed. We identified more than 18.2 million variants in total, of which 3.3% so far are not represented in gnomAD and dbSNP databases. Moreover, we observed a notable though distinct clustering of the Latvian cohort within the European subpopulations. Finally, our findings demonstrate the improved performance of imputation of variants using the Latvian population-specific reference panel in the Latvian population compared to established IRPs. In summary, our study provides the first WGS data for a regional reference genome that will serve as a resource for the development of precision medicine and complement the global genome dataset, improving the understanding of human genetic variation.

Osteopathia striata with cranial sclerosis causing a compressive optic neuropathy

BACKGROUND

Osteopathia striata combined with cranial sclerosis (OS-CS) is an inherited skeletal dysplasia that manifests with macrocephaly, orofacial abnormalities, thickened craniofacial bones, and vertically oriented radiodensities of the long bones.

CASE REPORT

Here, we present a severe case of OS-CS in a 4-year-old girl causing optic neuropathy as shown by radiographic evidence, ophthalmic findings, and histopathology. Previous genetic testing in this patient revealed a de novo heterozygous mutation in AMER1 (c.1057C>T, p.Arg353Ter). Although the patient had a pre-existing, appropriately functioning, ventriculoperitoneal (VP) shunt, a subsequent MRI of the brain and orbits showed narrowing of the bilateral optic nerve canals secondary to osseous thickening causing bilateral optic nerve atrophy, worse on the left. The patient underwent staged bilateral orbital osteotomies, optic canal decompression, and bilateral frontal craniotomy, and at 11 months postoperatively, her vision remained stable. Conclusions: While up to 50% of the patients with OS-CS may experience hearing loss due to cranial nerve compression, we present a case of severe visual loss secondary to OS-CS-associated optic nerve compression.

Glucocorticoid Receptor Gene (NR3C1) Polymorphisms and Metabolic Syndrome: Insights from the Mennonite Population

The regulation of the hypothalamic-pituitary-adrenal (HPA) axis is associated with polymorphisms and the methylation degree of the glucocorticoid receptor gene (NR3C1) and is potentially involved in the development of metabolic syndrome (MetS). In order to evaluate the association between MetS with the polymorphisms, methylation, and gene expression of the NR3C1 in the genetically isolated Brazilian Mennonite population, we genotyped 20 NR3C1 polymorphisms in 74 affected (MetS) and 138 unaffected individuals without affected first-degree relatives (Co), using exome sequencing, as well as five variants from non-exonic regions, in 70 MetS and 166 Co, using mass spectrometry. The methylation levels of 11 1F CpG sites were quantified using pyrosequencing (66 MetS and 141 Co), and the NR3C1 expression was evaluated via RT-qPCR (14 MetS and 25 Co). Age, physical activity, and family environment during childhood were associated with MetS. Susceptibility to MetS, independent of these factors, was associated with homozygosity for rs10482605*C (OR = 4.74, pcorr = 0.024) and the haplotype containing TTCGTTGATT (rs3806855*T_ rs3806854*T_rs10482605*C_rs10482614*G_rs6188*T_rs258813*T_rs33944801*G_rs34176759*A_rs17209258*T_rs6196*T, OR = 4.74, pcorr = 0.048), as well as for the CCT haplotype (rs41423247*C_ rs6877893*C_rs258763*T), OR = 6.02, pcorr = 0.030), but not to the differences in methylation or gene expression. Thus, NR3C1 polymorphisms seem to modulate the susceptibility to MetS in Mennonites, independently of lifestyle and early childhood events, and their role seems to be unrelated to DNA methylation and gene expression.

Clinical and Genetic Characterization of Familial Central Precocious Puberty

CONTEXT

Central precocious puberty (CPP) can have a familial form in approximately one-quarter of the children. The recognition of this inherited condition increased after the identification of autosomal dominant CPP with paternal transmission caused by mutations in the MKRN3 and DLK1 genes.

OBJECTIVE

We aimed to characterize the inheritance and estimate the prevalence of familial CPP in a large multiethnic cohort; to compare clinical and hormonal features, as well as treatment response to GnRH analogs (GnRHa), in children with distinct modes of transmission; and to investigate the genetic basis of familial CPP.

METHODS

We retrospectively studied 586 children with a diagnosis of CPP. Patients with familial CPP (n = 276) were selected for clinical and genetic analysis. Data from previous studies were grouped, encompassing sequencing of MKRN3 and DLK1 genes in 204 patients. Large-scale parallel sequencing was performed in 48 individuals from 34 families.

RESULTS

The prevalence of familial CPP was estimated at 22%, with a similar frequency of maternal and paternal transmission. Pedigree analyses of families with maternal transmission suggested an autosomal dominant inheritance. Clinical and hormonal features, as well as treatment response to GnRHa, were similar among patients with different forms of transmission of familial CPP. MKRN3 loss-of-function mutations were the most prevalent cause of familial CPP, followed by DLK1 loss-of-function mutations, affecting, respectively, 22% and 4% of the studied families; both affected exclusively families with paternal transmission. Rare variants of uncertain significance were identified in CPP families with maternal transmission.

CONCLUSION

We demonstrated a similar prevalence of familial CPP with maternal and paternal transmission. MKRN3 and DLK1 loss-of-function mutations were the major causes of familial CPP with paternal transmission.

A brazilian nationwide multicenter study on deficiency of deaminase-2 (DADA2)

INTRODUCTION

The deficiency of ADA2 (DADA2) is a rare autoinflammatory disease provoked by mutations in the ADA2 gene inherited in a recessive fashion. Up to this moment there is no consensus for the treatment of DADA2 and anti-TNF is the therapy of choice for chronic management whereas bone marrow transplantation is considered for refractory or severe phenotypes. Data from Brazil is scarce and this multicentric study reports 18 patients with DADA2 from Brazil.

PATIENTS AND METHODS

This is a multicentric study proposed by the Center for Rare and Immunological Disorders of the Hospital 9 de Julho - DASA, São Paulo - Brazil. Patients of any age with a confirmed diagnosis of DADA2 were eligible for this project and data on clinical, laboratory, genetics and treatment were collected.

RESULTS

Eighteen patients from 10 different centers are reported here. All patients had disease onset at the pediatric age (median of 5 years) and most of them from the state of São Paulo. Vasculopathy with recurrent stroke was the most common phenotype but atypical phenotypes compatible with ALPS-like and Common Variable Immunodeficiency (CVID) was also found. All patients carried pathogenic mutations in the ADA2 gene. Acute management of vasculitis was not satisfactory with steroids in many patients and all those who used anti-TNF had favorable responses.

CONCLUSION

The low number of patients diagnosed with DADA2 in Brazil reinforces the need for disease awareness for this condition. Moreover, the absence of guidelines for diagnosis and management is also necessary (t).

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.

A thermosensitive PCNA allele underlies an ataxia-telangiectasia-like disorder

Proliferating cell nuclear antigen (PCNA) is a sliding clamp protein that coordinates DNA replication with various DNA maintenance events that are critical for human health. Recently, a hypomorphic homozygous serine to isoleucine (S228I) substitution in PCNA was described to underlie a rare DNA repair disorder known as PCNA-associated DNA repair disorder (PARD). PARD symptoms range from UV sensitivity, neurodegeneration, telangiectasia, and premature aging. We, and others, previously showed that the S228I variant changes the protein-binding pocket of PCNA to a conformation that impairs interactions with specific partners. Here, we report a second PCNA substitution (C148S) that also causes PARD. Unlike PCNA-S228I, PCNA-C148S has WT-like structure and affinity toward partners. In contrast, both disease-associated variants possess a thermostability defect. Furthermore, patient-derived cells homozygous for the C148S allele exhibit low levels of chromatin-bound PCNA and display temperature-dependent phenotypes. The stability defect of both PARD variants indicates that PCNA levels are likely an important driver of PARD disease. These results significantly advance our understanding of PARD and will likely stimulate additional work focused on clinical, diagnostic, and therapeutic aspects of this severe disease.

Dynamics of the personalities of PCSK9 on missense variants (rs505151 and rs562556) from elderly cohort studies in Brazil

The Proprotein Convertase Subtilisin/Kexin Type 9 (PCSK9) promotes the degradation of the low-density lipoprotein receptors (LDLR). Gain-of-function (GOF) variants of PCSK9 significantly affects lipid metabolism leading to coronary artery disease (CAD), owing to the raising the plasma low-density lipoprotein (LDL). Considering the public health matter, large-scale genomic studies have been conducted worldwide to provide the genetic architecture of populations for the implementation of precision medicine actions. Nevertheless, despite the advances in genomic studies, non-European populations are still underrepresented in public genomic data banks. Despite this, we found two high-frequency variants (rs505151 and rs562556) in the ABraOM databank (Brazilian genomic variants) from a cohort SABE study conducted in the largest city of Brazil, São Paulo. Here, we assessed the structural and dynamical features of these variants against WT through a molecular dynamics study. We sought fundamental dynamical interdomain relations through Perturb Response Scanning (PRS) and we found an interesting change of dynamical relation between prodomain and Cysteine-Histidine-Rich-Domain (CHRD) in the variants. The results highlight the pivotal role of prodomain in the PCSK9 dynamic and the implications for the development of new drugs depending on patient group genotype.

Immunological evaluation of young unvaccinated patients with Turner syndrome after COVID-19

OBJECTIVES

The X-chromosome contains the largest number of immune-related genes, which play a major role in COVID-19 symptomatology and susceptibility. Here, we had a unique opportunity to investigate, for the first time, COVID-19 outcomes in six unvaccinated young Brazilian patients with Turner syndrome (TS; 45, X0), including one case of critical illness in a child aged 10 years, to evaluate their immune response according to their genetic profile.

METHODS

A serological analysis of humoral immune response against SARS-CoV-2, phenotypic characterization of antiviral responses in peripheral blood mononuclear cells after stimuli, and the production of cytotoxic cytokines of T lymphocytes and natural killer cells were performed in blood samples collected from the patients with TS during the convalescence period. Whole exome sequencing was also performed.

RESULTS

Our volunteers with TS showed a delayed or insufficient humoral immune response to SARS-CoV-2 (particularly immunoglobulin G) and a decrease in interferon-γ production by cluster of differentiation (CD)4+ and CD8+ T lymphocytes after stimulation with toll-like receptors 7/8 agonists. In contrast, we observed a higher cytotoxic activity in the volunteers with TS than the volunteers without TS after phorbol myristate acetate/ionomycin stimulation, particularly granzyme B and perforin by CD8+ and natural killer cells. Interestingly, two volunteers with TS carry rare genetic variants in genes that regulate type I and III interferon immunity.

CONCLUSION

Following previous reports in the literature for other conditions, our data showed that patients with TS may have an impaired immune response against SARS-CoV-2. Furthermore, other medical conditions associated with TS could make them more vulnerable to COVID-19.

Adaptive Savitzky–Golay Filters for Analysis of Copy Number Variation Peaks from Whole-Exome Sequencing Data

Copy number variation (CNV) is a form of structural variation in the human genome that provides medical insight into complex human diseases; while whole-genome sequencing is becoming more affordable, whole-exome sequencing (WES) remains an important tool in clinical diagnostics. Because of its discontinuous nature and unique characteristics of sparse target-enrichment-based WES data, the analysis and detection of CNV peaks remain difficult tasks. The Savitzky–Golay (SG) smoothing is well known as a fast and efficient smoothing method. However, no study has documented the use of this technique for CNV peak detection. It is well known that the effectiveness of the classical SG filter depends on the proper selection of the window length and polynomial degree, which should correspond with the scale of the peak because, in the case of peaks with a high rate of change, the effectiveness of the filter could be restricted. Based on the Savitzky–Golay algorithm, this paper introduces a novel adaptive method to smooth irregular peak distributions. The proposed method ensures high-precision noise reduction by dynamically modifying the results of the prior smoothing to automatically adjust parameters. Our method offers an additional feature extraction technique based on density and Euclidean distance. In comparison to classical Savitzky–Golay filtering and other peer filtering methods, the performance evaluation demonstrates that adaptive Savitzky–Golay filtering performs better. According to experimental results, our method effectively detects CNV peaks across all genomic segments for both short and long tags, with minimal peak height fidelity values (i.e., low estimation bias). As a result, we clearly demonstrate how well the adaptive Savitzky–Golay filtering method works and how its use in the detection of CNV peaks can complement the existing techniques used in CNV peak analysis.

SCAF4-related syndromic intellectual disability

The causal link between variants in the SCAF4 gene and a syndromic form of intellectual disability (ID) was established in 2020 by Fliedner et al. Since then, no additional cases have been reported. We performed exome sequencing in a 16-year-old Brazilian male presenting with ID, epilepsy, behavioral problems, speech impairment, facial dysmorphisms, heart malformations, and obesity. A de novo pathogenic variant [SCAF4(NM_020706.2):c.374_375dup(p.Glu126LeufsTer20)] was identified. This is the second study reporting the involvement of SCAF4 in syndromic ID, and the description of the patient's clinical features contributes to defining the phenotypic spectrum of this recently described Mendelian disorder.

A review of ancestrality and admixture in Latin America and the caribbean focusing on native American and African descendant populations

Genomics can reveal essential features about the demographic evolution of a population that may not be apparent from historical elements. In recent years, there has been a significant increase in the number of studies applying genomic epidemiological approaches to understand the genetic structure and diversity of human populations in the context of demographic history and for implementing precision medicine. These efforts have traditionally been applied predominantly to populations of European origin. More recently, initiatives in the United States and Africa are including more diverse populations, establishing new horizons for research in human populations with African and/or Native ancestries. Still, even in the most recent projects, the under-representation of genomic data from Latin America and the Caribbean (LAC) is remarkable. In addition, because the region presents the most recent global miscegenation, genomics data from LAC may add relevant information to understand population admixture better. Admixture in LAC started during the colonial period, in the 15th century, with intense miscegenation between European settlers, mainly from Portugal and Spain, with local indigenous and sub-Saharan Africans brought through the slave trade. Since, there are descendants of formerly enslaved and Native American populations in the LAC territory; they are considered vulnerable populations because of their history and current living conditions. In this context, studying LAC Native American and African descendant populations is important for several reasons. First, studying human populations from different origins makes it possible to understand the diversity of the human genome better. Second, it also has an immediate application to these populations, such as empowering communities with the knowledge of their ancestral origins. Furthermore, because knowledge of the population genomic structure is an essential requirement for implementing genomic medicine and precision health practices, population genomics studies may ensure that these communities have access to genomic information for risk assessment, prevention, and the delivery of optimized treatment; thus, helping to reduce inequalities in the Western Hemisphere. Hoping to set the stage for future studies, we review different aspects related to genetic and genomic research in vulnerable populations from LAC countries.

Dual Molecular Diagnoses of Recessive Disorders in a Child from Consanguineous Parents: Case Report and Literature Review

The widespread use of whole exome sequencing (WES) resulted in the discovery of multilocus pathogenic variations (MPV), defined as two or more distinct or overlapping Mendelian disorders occurring in a patient, leading to a blended phenotype. In this study, we report on a child with autosomal recessive primary microcephaly-5 (MCPH5) and nephropathic cystinosis. The proband is the first child of consanguineous parents, presenting a complex phenotype including neurodevelopmental delay, microcephaly, growth restriction, significant delay of bone maturation, lissencephaly, and abnormality of neuronal migration, photophobia, and renal tubular acidosis. WES revealed two pathogenic and homozygous variants: a c.4174C>T variant in the ASPM gene and a c.382C>T variant in the CTNS gene, explaining the complex phenotype. The literature review showed that most of the patients harboring two variants in recessive disease genes are born to consanguineous parents. To the best of our knowledge, the patient herein described is the first one harboring pathogenic variants in both the ASPM and CTNS genes. These findings highlight the importance of searching for MPV in patients with complex phenotypes investigated by genome-wide testing methods, especially for those patients born to consanguineous parents.

Challenges in selecting admixture models and marker sets to infer genetic ancestry in a Brazilian admixed population

The inference of genetic ancestry plays an increasingly prominent role in clinical, population, and forensic genetics studies. Several genotyping strategies and analytical methodologies have been developed over the last few decades to assign individuals to specific biogeographic regions. However, despite these efforts, ancestry inference in populations with a recent history of admixture, such as those in Brazil, remains a challenge. In admixed populations, proportion and components of genetic ancestry vary on different levels: (i) between populations; (ii) between individuals of the same population, and (iii) throughout the individual's genome. The present study evaluated 1171 admixed Brazilian samples to compare the genetic ancestry inferred by tri-/tetra-hybrid admixture models and evaluated different marker sets from those with small numbers of ancestry informative markers panels (AIMs), to high-density SNPs (HDSNP) and whole-genome-sequence (WGS) data. Analyses revealed greater variation in the correlation coefficient of ancestry components within and between admixed populations, especially for minority ancestral components. We also observed positive correlation between the number of markers in the AIMs panel and HDSNP/WGS. Furthermore, the greater the number of markers, the more accurate the tri-/tetra-hybrid admixture models.

A Polymorphism in the TMPRSS2 Gene Increases the Risk of Death in Older Patients Hospitalized with COVID-19

BACKGROUND

Transmembrane serine protease type 2 (TMPRSS2) and angiotensin-converting enzyme 2 (ACE2) are the main molecules involved in the entry of SARS-CoV-2 into host cells. Changes in TMPRSS2 expression levels caused by single nucleotide polymorphisms (SNPs) may contribute to the outcome of COVID-19. The aim was to investigate the association between TMPRSS2 gene polymorphisms and the risk of death in hospitalized patients with COVID-19.

METHODS

We included patients with confirmed COVID-19, recruited from two hospitals in northeastern Brazil from August 2020 to July 2021. Two functional polymorphisms (rs2070788 and rs12329760) in TMPRSS2 were evaluated by real-time PCR. The Kaplan-Meier method was used to estimate death. The Cox's proportional hazards model was used to adjust for potentially confounding factors.

RESULTS

A total of 402 patients were followed prospectively. Survival analysis demonstrated that older patients carrying the rs2070788 GG genotype had shorter survival times when compared to those with AG or AA genotypes (p = 0.009). In multivariable analysis, the GG genotype was a factor independently associated with the risk of death in older individuals (hazard ratio = 4.03, 95% confidence interval 1.49 to 10.84).

CONCLUSIONS

The rs2070788 polymorphism in TMPRSS2 increases risk of death four-fold in older patients hospitalized with COVID-19.

The oldest unvaccinated Covid-19 survivors in South America

BACKGROUND

Although older adults are at a high risk of severe or critical Covid-19, there are many cases of unvaccinated centenarians who had a silent infection or recovered from mild or moderate Covid-19. We studied three Brazilian supercentenarians, older than 110 years, who survived Covid-19 in 2020 before being vaccinated.

RESULTS

Despite their advanced age, humoral immune response analysis showed that these individuals displayed robust levels of IgG and neutralizing antibodies (NAbs) against SARS-CoV-2. Enrichment of plasma proteins and metabolites related to innate immune response and host defense was also observed. None presented autoantibodies (auto-Abs) to type I interferon (IFN). Furthermore, these supercentenarians do not carry rare variants in genes underlying the known inborn errors of immunity, including particular inborn errors of type I IFN.

CONCLUSION

These observations suggest that their Covid-19 resilience might be a combination of their genetic background and their innate and adaptive immunity.

MUC22, HLA-A, and HLA-DOB variants and COVID-19 in resilient super-agers from Brazil

BACKGROUND

Although aging correlates with a worse prognosis for Covid-19, super elderly still unvaccinated individuals presenting mild or no symptoms have been reported worldwide. Most of the reported genetic variants responsible for increased disease susceptibility are associated with immune response, involving type I IFN immunity and modulation; HLA cluster genes; inflammasome activation; genes of interleukins; and chemokines receptors. On the other hand, little is known about the resistance mechanisms against SARS-CoV-2 infection. Here, we addressed polymorphisms in the MHC region associated with Covid-19 outcome in super elderly resilient patients as compared to younger patients with a severe outcome.

METHODS

SARS-CoV-2 infection was confirmed by RT-PCR test. Aiming to identify candidate genes associated with host resistance, we investigated 87 individuals older than 90 years who recovered from Covid-19 with mild symptoms or who remained asymptomatic following positive test for SARS-CoV-2 as compared to 55 individuals younger than 60 years who had a severe disease or died due to Covid-19, as well as to the general elderly population from the same city. Whole-exome sequencing and an in-depth analysis of the MHC region was performed. All samples were collected in early 2020 and before the local vaccination programs started.

RESULTS

We found that the resilient super elderly group displayed a higher frequency of some missense variants in the MUC22 gene (a member of the mucins' family) as one of the strongest signals in the MHC region as compared to the severe Covid-19 group and the general elderly control population. For example, the missense variant rs62399430 at MUC22 is two times more frequent among the resilient super elderly (p = 0.00002, OR = 2.24).

CONCLUSION

Since the pro-inflammatory basal state in the elderly may enhance the susceptibility to severe Covid-19, we hypothesized that MUC22 might play an important protective role against severe Covid-19, by reducing overactive immune responses in the senior population.

Artificial intelligence and Big Data in neurology

Recent advances in technology have allowed us access to a multitude of datasets pertaining to various dimensions in neurology. Together with the enormous opportunities, we also face challenges related to data quality, ethics and intrinsic difficulties related to the application of data science in healthcare. In this article we will describe the main advances in the field of artificial intelligence and Big Data applied to neurology with a focus on neurosciences based on medical images. Real-World Data (RWD) and analytics related to large volumes of information will be described as well as some of the most relevant scientific initiatives at the time of this writing.