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.

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.

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.

Corrigendum: Follow-up of young adult monozygotic twins after simultaneous critical coronavirus disease 2019: a case report.

[This corrects the article DOI: 10.3389/fmed.2022.1008585.].

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.

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.

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.

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.

Follow-up of young adult monozygotic twins after simultaneous critical coronavirus disease 2019: A case report.

Background: The influence of the host genome on coronavirus disease 2019 (COVID-19) susceptibility and severity is supported by reports on monozygotic (MZ) twins where both were infected simultaneously with similar disease outcomes, including several who died due to the SARS-CoV-2 infection within days apart. However, successive exposures to pathogens throughout life along with other environmental factors make the immune response unique for each individual, even among MZ twins. Case presentation and methods: Here we report a case of a young adult monozygotic twin pair, who caught attention since both presented simultaneously severe COVID-19 with the need for oxygen support despite age and good health conditions. One of the twins, who spent more time hospitalized, reported symptoms of long-COVID even 7 months after infection. Immune cell profile and specific responses to SARS-CoV-2 were evaluated as well as whole exome sequencing. Conclusion: Although the MZ twin brothers shared the same genetic mutations which may be associated with their increased risk of developing severe COVID-19, their clinical progression was different, reinforcing the role of both immune response and genetics in the COVID-19 presentation and course. Besides, post-COVID syndrome was observed in one of them, corroborating an association between the duration of hospitalization and the occurrence of long-COVID symptoms.

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

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 enables identification of ~2,000 previously undescribed mobile element insertions without previous description, nearly 5 Mb of genomic segments absent from the human genome reference, and over 140 alleles from HLA genes absent from public resources. We reclassify and curate pathogenicity assertions for nearly four hundred variants in genes associated with dominantly-inherited Mendelian disorders and calculate the incidence for selected recessive disorders, demonstrating the clinical usefulness of the present study. Finally, we observe 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.

The 90 plus: longevity and COVID-19 survival.

The world population is getting older and studies aiming to enhance our comprehension of the underlying mechanisms responsible for health span are of utmost interest for longevity and as a measure for health care. In this review, we summarized previous genetic association studies (GWAS) and next-generation sequencing (NGS) of elderly cohorts. We also present the updated hypothesis for the aging process, together with the factors associated with healthy aging. We discuss the relevance of studying older individuals and build databanks to characterize the presence and resistance against late-onset disorders. The identification of about 2 million novel variants in our cohort of more than 1000 elderly Brazilians illustrates the importance of studying highly admixed populations of non-European ancestry. Finally, the ascertainment of nonagenarians and particularly of centenarians who were recovered from COVID-19 or remained asymptomatic opens new avenues of research aiming to enhance our comprehension of biological mechanisms associated with resistance against pathogens.

Recurrence of COVID-19 associated with reduced T-cell responses in a monozygotic twin pair.

Recurrence of COVID-19 in recovered patients has been increasingly reported. However, the immune mechanisms behind the recurrence have not been thoroughly investigated. The presence of neutralizing antibodies (nAbs) in recurrence/reinfection cases suggests that other types of immune response are involved in protection against recurrence. Here, we investigated the innate type I/III interferon (IFN) response, binding and nAb assays and T-cell responses to severe acute respiratory distress syndrome coronavirus 2 (SARS-CoV-2) with IFN gamma (IFNγ) enzyme-linked spot assay (ELISPOT) in three pairs of young adult monozygotic (MZ) twins with previous confirmed COVID-19, one of them presenting a severe recurrence four months after the initial infection. Twin studies have been of paramount importance to comprehend the immunogenetics of infectious diseases. Each MZ twin pair was previously exposed to SARS-CoV-2, as seen by clinical reports. The six individuals presented similar overall recovered immune responses except for the recurrence case, who presented a drastically reduced number of recognized SARS-CoV-2 T-cell epitopes on ELISPOT as compared to her twin sister and the other twin pairs. Our results suggest that the lack of a broad T-cell response to initial infection may have led to recurrence, emphasizing 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.

Men are the main COVID-19 transmitters: behavior or biology?

BACKGROUND: COVID-19 has affected millions of people worldwide. Clinical manifestations range from severe cases with lethal outcome to mild or asymptomatic cases. Although the proportion of infected individuals does not differ between sexes, men are more susceptible to severe COVID-19, with a higher risk of death than women. Also, men are pointed out as more lax regarding protective measures, mask wearing and vaccination. Thus, we questioned whether sex-bias may be explained by biological pathways and/or behavioral aspects or both. METHODS: Between July 2020 and July 2021, we performed an epidemiological survey including 1744 unvaccinated adult Brazilian couples, with there was at least one infected symptomatic member, who were living together during the COVID-19 infection without protective measures. Presence or absence of infection was confirmed by RT-PCR and/or serology results. Couples were divided into two groups: (1) both partners were infected (concordant couples) and (2) one partner was infected and the spouse remained asymptomatic despite the close contact with the COVID-19 symptomatic partner (discordant couples). Statistical analysis of the collected data was performed aiming to verify a differential transmission potential between genders in couples keeping contact without protective measures. RESULTS: The combination of our collected data showed that the man is the first (or the only) affected member in most cases when compared to women and that this difference may be explained by biological and behavioral factors. CONCLUSIONS: The present study confirmed the existence of gender differences not only for susceptibility to infection and resistance to COVID-19 but also in its transmission rate.

Longitudinal 16S rRNA gut microbiota data of infant triplets show partial susceptibility to host genetics.

The question of whether host genetics plays a role in the development of the infant gut microbiota does not, as yet, have a clear answer. In order to throw additional light on this question, we have analyzed 16S rRNA amplicon sequences from 99 valid fecal samples of five sets of dichorionic triplet babies born by C-section from 1 to 36 months of age. Beta diversity analysis showed that monozygotic twins were more similar to each other than their dizygotic siblings. Monozygotic twins also tended to share more amplicon sequence variants between them. Heritability analysis showed that the genera Bacteroides and Veillonella are particularly susceptible to host genetics. We conclude that infant gut microbiota development is influenced by host genetics, but this effect is subtle and may affect only certain bacterial taxa during a limited time period early in life.

HLA-G genetic diversity and evolutive aspects in worldwide populations.

HLA-G is a promiscuous immune checkpoint molecule. The HLA-G gene presents substantial nucleotide variability in its regulatory regions. However, it encodes a limited number of proteins compared to classical HLA class I genes. We characterized the HLA-G genetic variability in 4640 individuals from 88 different population samples across the globe by using a state-of-the-art method to characterize polymorphisms and haplotypes from high-coverage next-generation sequencing data. We also provide insights regarding the HLA-G genetic diversity and a resource for future studies evaluating HLA-G polymorphisms in different populations and association studies. Despite the great haplotype variability, we demonstrated that: (1) most of the HLA-G polymorphisms are in introns and regulatory sequences, and these are the sites with evidence of balancing selection, (2) linkage disequilibrium is high throughout the gene, extending up to HLA-A, (3) there are few proteins frequently observed in worldwide populations, with lack of variation in residues associated with major HLA-G biological properties (dimer formation, interaction with leukocyte receptors). These observations corroborate the role of HLA-G as an immune checkpoint molecule rather than as an antigen-presenting molecule. Understanding HLA-G variability across populations is relevant for disease association and functional studies.

MHC Variants Associated With Symptomatic Versus Asymptomatic SARS-CoV-2 Infection in Highly Exposed Individuals.

Despite the high number of individuals infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) who develop coronavirus disease 2019 (COVID-19) symptoms worldwide, many exposed individuals remain asymptomatic and/or uninfected and seronegative. This could be explained by a combination of environmental (exposure), immunological (previous infection), epigenetic, and genetic factors. Aiming to identify genetic factors involved in immune response in symptomatic COVID-19 as compared to asymptomatic exposed individuals, we analyzed 83 Brazilian couples where one individual was infected and symptomatic while the partner remained asymptomatic and serum-negative for at least 6 months despite sharing the same bedroom during the infection. We refer to these as "discordant couples". We performed whole-exome sequencing followed by a state-of-the-art method to call genotypes and haplotypes across the highly polymorphic major histocompatibility complex (MHC) region. The discordant partners had comparable ages and genetic ancestry, but women were overrepresented (65%) in the asymptomatic group. In the antigen-presentation pathway, we observed an association between HLA-DRB1 alleles encoding Lys at residue 71 (mostly DRB1*03:01 and DRB1*04:01) and DOB*01:02 with symptomatic infections and HLA-A alleles encoding 144Q/151R with asymptomatic seronegative women. Among the genes related to immune modulation, we detected variants in MICA and MICB associated with symptomatic infections. These variants are related to higher expression of soluble MICA and low expression of MICB. Thus, quantitative differences in these molecules that modulate natural killer (NK) activity could contribute to susceptibility to COVID-19 by downregulating NK cell cytotoxic activity in infected individuals but not in the asymptomatic partners.

Biased pathogenic assertions of loss of function variants challenge molecular diagnosis of admixed individuals.

Diagnosis of individuals affected by monogenic disorders was significantly improved by next-generation sequencing targeting clinically relevant genes. Whole exomes yield a large number of variants that require several filtering steps, prioritization, and pathogenicity classification. Among the criteria recommended by ACMG, those that rely on population databases critically affect analyses of individuals with underrepresented ancestries. Population-specific allelic frequencies need consideration when characterizing potential deleteriousness of variants. An orthogonal input for classification is annotation of variants previously classified as pathogenic as a criterion that provide supporting evidence widely sourced at ClinVar. We used a whole-genome dataset from a census-based cohort of 1,171 elderly individuals from São Paulo, Brazil, highly admixed, and unaffected by severe monogenic disorders, to investigate if pathogenic assertions in ClinVar are enriched with higher proportions of European ancestry, indicating bias. Potential loss of function (pLOF) variants were filtered from 4,250 genes associated with Mendelian disorders and annotated with ClinVar assertions. Over 1,800 single nucleotide pLOF variants were included, 381 had non-benign assertions. Among carriers (N = 463), average European ancestry was significantly higher than noncarriers (N = 708; p = .011). pLOFs in genomic contexts of non-European local ancestries were nearly three times less likely to have any ClinVar entry (OR = 0.353; p <.0001). Independent pathogenicity assertions are useful for variant classification in molecular diagnosis. However, European overrepresentation of assertions can promote distortions when classifying variants in non-European individuals, even in admixed samples with a relatively high proportion of European ancestry. The investigation and deposit of clinically relevant findings of diverse populations is fundamental improve this scenario.

Admixture/fine-mapping in Brazilians reveals a West African associated potential regulatory variant (rs114066381) with a strong female-specific effect on body mass and fat mass indexes.

BACKGROUND/OBJECTIVES: Admixed populations are a resource to study the global genetic architecture of complex phenotypes, which is critical, considering that non-European populations are severely underrepresented in genomic studies. Here, we study the genetic architecture of BMI in children, young adults, and elderly individuals from the admixed population of Brazil. SUBJECTS/METHODS: 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 Brazilian population-based cohorts from Northeast (Salvador), Southeast (Bambuí), and South (Pelotas). RESULTS: 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.76e-06). This variant is 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 women with morbid obesity from Pelotas, São Paulo, and Bambuí. The effect size of rs114066381 is at least five times higher than the FTO SNPs rs9939609 and rs1558902, already emblematic for their high effects. CONCLUSIONS: We identified six candidate SNPs associated with BMI. rs114066381 stands out for its high effect that was replicated and its high frequency in women with morbid obesity. We demonstrate how admixed populations are a source of new relevant phenotype-associated genetic variants.