Exomic variants of an elderly cohort of Brazilians in the ABraOM database

The Impact of TBXA2R Gene Variants on the Risk of Aspirin-Induced Upper Gastrointestinal Bleeding: A Case-Control Study

Objective: Upper gastrointestinal bleeding (UGIB) has been identified as a potential adverse drug reaction associated with the use of low-dose aspirin (LDA). This study aimed to investigate the relationship between variants in the TBXA2R gene, which is involved in platelet aggregation, and the risk of UGIB in patients with cardiovascular diseases treated with LDA. Methods: A case-control study was conducted at a Brazilian hospital complex. Three groups were defined: (1) case group (n = 50): patients with cardiovascular disease who used LDA and were diagnosed with UGIB of non-variceal etiology, (2) LDA control group (n = 50): patients with cardiovascular disease who used LDA without developing UGIB, and (3) healthy control group (n = 189). Data were collected through face-to-face interviews, and blood samples were collected for the analysis of Helicobacter pylori infection and genotyping of 3 genetic variants [rs2238631 (C > T), rs4807491 (A > G), and rs1131882 (A > G)]. Results: The case group had a significantly higher frequency of carriers of the rs4807491.G allele compared to the control group of LDA users (P-value = .004). No significant difference was observed in the proportion of carriers of the rs2238631.T and 1131882.G variants between the studied groups. Carriers of rs2238631.T (OR: 4.515, 95% CI: 1.37-14.89) and rs4807491.G allele (OR: 3.232, 95% CI: 1.12-9.37) exhibited a higher risk of UGIB. Conclusion: These findings suggest that the presence of the rs2238631 and rs4807491 variant alleles is associates with a 3- to 4-fold increased risk of UGIB in patients with cardiovascular diseases treated with LDA. Future studies with larger sample sizes should confirm these results and to better identify individuals who may benefit from chronic LDA use.

Molecular investigation in individuals with orofacial clefts and microphthalmia-anophthalmia-coloboma spectrum

This study describes genomic findings among individuals with both orofacial clefts (OC) and microphthalmia/anophthalmia/coloboma (MAC) recorded in the Brazilian Database on Craniofacial Anomalies (BDCA). Chromosomal microarray analysis (CMA) and Whole Exome Sequencing (WES) were performed in 17 individuals with OC-MAC. Clinical interpretation of molecular findings was based on data available at the BDCA and on re-examination. No copy number variants (CNVs) classified as likely pathogenic or pathogenic were detected by CMA. WES allowed a conclusive diagnosis in six individuals (35.29%), two of them with variants in the CHD7 gene, and the others with variants in the TFAP2A, POMT1, PTPN11, and TP63 genes with the following syndromes: CHARGE, CHD7-spectrum, Branchiooculofacial, POMT1-spectrum, LEOPARD, and ADULT. Variants of uncertain significance (VUS) possibly associated to the phenotypes were found in six other individuals. Among the individuals with VUSes, three individuals presented variants in genes associated to defects of cilia structure and/or function, including DYNC2H1, KIAA0586, WDR34, INTU, RPGRIP1L, KIF7, and LMNA. These results show that WES was the most effective molecular approach for OC-MAC in this cohort. This study also reinforces the genetic heterogeneity of OC-MAC, and the importance of genes related to ciliopathies in this phenotype.

Germline DNA Damage Repair Gene Alterations in Patients with Metachronous Breast and Colorectal Cancer

A hereditary component of breast (BC) and colorectal cancer (CRC) has been described in approximately one-third of these tumor types. BC patients have an increased risk of developing CRC as a second primary tumor and vice versa. Germline genomic variants (NextSeq550, Illumina) were investigated in 24 unrelated BC and/or CRC patients and 7 relatives from 3 index patients. Fifty-six pathogenic or likely pathogenic variants were identified in 19 of 24 patients. We detected single-nucleotide variants (SNVs) in CRC predisposition genes (MLH1 and MUTYH) and other promising candidates (CDK5RAP3, MAD1L1, NOS3, and POLM). Eighteen patients presented SNVs or copy number variants (CNVs) in DNA damage repair genes. We also identified SNVs recently associated with BC or CRC predisposition (PABPC1, TYRO3, MAP3K1, SLC15A4, and LAMA1). The PABPC1c.1255C>T variant was detected in nine unrelated patients. Each patient presented at least one SNV/CNV in a candidate gene, and most had alterations in more than one gene, reinforcing a polygenic model for BC/CRC predisposition. A significant fraction of BC/CRC patients with a family history of these tumors harbored deleterious germline variants in DNA repair genes. Our findings can lead to strategies to improve the diagnosis, genetic counseling, and treatment of patients and their relatives.

Language development in children from a public cochlear implant program

OBJECTIVE

To evaluate the rate of loss to follow-up in a cochlear implant program from the public health system in Southern Brazil as well as the characteristics of hearing loss, sociodemographic, sociocultural and the development of oral language in children with prelingual deafness.

METHODS

Retrospective cohort study with children who underwent CI surgery between 2010 and 2020. Data was collected through of interviews and review of medical records. The language development assessment was performed using the MUSS, MAIS and IT-MAIS scales. For the classification of language development, we used as parameters the values (mean ± SD) found in a previous national study. From those values, the Z-score for each patient at each hearing age (time of experience with the cochlear implant) was calculated.

RESULTS

Of the 225 children implanted between 2010-2020, 129 were included in this study. The rate of loss to follow-up in the program was 42.6%. The mean age at first surgery was 40.5 (±16.9) months, with 77.5% of patients having received a unilateral implant. Language results below the expected for hearing age (

CONCLUSIONS

Most patients had an elevated mean age at cochlear implantation and there was a high rate of loss to follow-up and low attendance to speech and programming sessions. An overall poor language performance was found for this pediatric cochlear implant program from the public health system in Southern Brazil.

LEVEL OF EVIDENCE

Level 3 (Non-randomized cohort study).

Collapse

Key Words

• Cochlear implant
• Language development disorders
• Prelingual deafness

Collapse

MESH Headings

• Humans
• Retrospective Studies
• Child, Preschool
• Male
• Female
• Brazil
• Language Development
• Cochlear Implants
• Cochlear Implantation
• Infant
• Deafness/surgery
• Child
• Socioeconomic Factors

Collapse

Grants Collapse

Affiliation(s)

Alice Lang Silva Universidade Federal do Rio Grande do Sul (UFRGS), Programa de Pós-Graduação em Saúde da Criança e do Adolescente, Porto Alegre, RS, Brazil.
Isadora Martins da Silva Stumpf Hospital de Clínicas de Porto Alegre (HCPA), Serviço de Otorrinolaringologia, Porto Alegre, RS, Brazil
Laura Prolla Lacroix Universidade Federal do Rio Grande do Sul (UFRGS), Faculdade de Medicina, Porto Alegre, RS, Brazil
Debora Milena Ferreira Alves Universidade Federal do Rio Grande do Sul (UFRGS), Faculdade de Medicina, Porto Alegre, RS, Brazil
Adriana Laybauer da Silveira Hospital de Clínicas de Porto Alegre (HCPA), Serviço de Fonoaudiologia, Porto Alegre, RS, Brazil
Sady Selaimen da Costa Universidade Federal do Rio Grande do Sul (UFRGS), Programa de Pós-Graduação em Saúde da Criança e do Adolescente, Porto Alegre, RS, Brazil; Hospital de Clínicas de Porto Alegre (HCPA), Serviço de Otorrinolaringologia, Porto Alegre, RS, Brazil
Letícia Petersen Schmidt Rosito Hospital de Clínicas de Porto Alegre (HCPA), Serviço de Otorrinolaringologia, Porto Alegre, RS, Brazil; Universidade Federal do Rio Grande do Sul (UFRGS), Faculdade de Medicina, Porto Alegre, RS, Brazil

Collapse

Maple syrup urine disease diagnosis in Brazilian patients by massive parallel sequencing

Biallelic pathogenic variants cause maple syrup urine disease (MSUD) in one of the branched-chain α-keto acid dehydrogenase (BCKDH) complex genes (BCKDHA, BCKDHB, DBT, DLD, and PPM1K) leading to the accumulation of leucine, isoleucine, and valine. This study aimed to perform a molecular diagnosis of Brazilian patients with MSUD using gene panels and massive parallel sequencing. Eighteen Brazilian patients with a biochemical diagnosis of MSUD were analyzed by massive parallel sequencing in the Ion PGM Torrent Server using a gene panel with the BCKDHA, BCKDHB, and DBT genes. The American College of Medical Genetics and Genomics guidelines were used to determine variant pathogenicity. Thirteen patients had both variants found by massive parallel sequencing, whereas 3 patients had only one variant found. In 2 patients, the variants were not found by this analysis. These 5 patients required additional Sanger sequencing to confirm their genotype. Twenty-five pathogenic variants were identified in the 3 MSUD-related genes (BCKDHA, BCKDHB, and DBT). Most variants were present in the BCKDHB gene, and no common variants were found. Nine novel variants were observed: c.922 A > G, c.964C > A, and c.1237 T > C in the BCKDHA gene; and c.80_90dup, c.384delA, c.478 A > T, c.528C > G, c.977 T > C, and c.1039-2 A > G in the BCKDHB gene. All novel variants were classified as pathogenic. Molecular modeling of the novel variants indicated that the binding of monomers was affected in the BCKDH complex tetramer, which could lead to a change in the stability and activity of the enzyme. Massive parallel sequencing with targeted gene panels seems to be a cost-effective method that can provide a molecular diagnosis of MSUD.

Clinical Characterization and Underlying Genetic Findings in Brazilian Patients with Syndromic Microcephaly Associated with Neurodevelopmental Disorders

Microcephaly is characterized by an occipitofrontal circumference at least two standard deviations below the mean for age and sex. Neurodevelopmental disorders (NDD) are commonly associated with microcephaly, due to perturbations in brain development and functioning. Given the extensive genetic heterogeneity of microcephaly, managing patients is hindered by the broad spectrum of diagnostic possibilities that exist before conducting molecular testing. We investigated the genetic basis of syndromic microcephaly accompanied by NDD in a Brazilian cohort of 45 individuals and characterized associated clinical features, as well as evaluated the effectiveness of whole-exome sequencing (WES) as a diagnostic tool for this condition. Patients previously negative for pathogenic copy number variants underwent WES, which was performed using a trio approach for isolated index cases (n = 31), only the index in isolated cases with parental consanguinity (n = 8) or affected siblings in familial cases (n = 3). Pathogenic/likely pathogenic variants were identified in 19 families (18 genes) with a diagnostic yield of approximately 45%. Nearly 86% of the individuals had global developmental delay/intellectual disability and 51% presented with behavioral disturbances. Additional frequent clinical features included facial dysmorphisms (80%), brain malformations (67%), musculoskeletal (71%) or cardiovascular (47%) defects, and short stature (54%). Our findings unraveled the underlying genetic basis of microcephaly in half of the patients, demonstrating a high diagnostic yield of WES for microcephaly and reinforcing its genetic heterogeneity. We expanded the phenotypic spectrum associated with the condition and identified a potentially novel gene (CCDC17) for congenital microcephaly.

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.

Variant effect predictors: a systematic review and practical guide

Large-scale association analyses using whole-genome sequence data have become feasible, but understanding the functional impacts of these associations remains challenging. Although many tools are available to predict the functional impacts of genetic variants, it is unclear which tool should be used in practice. This work provides a practical guide to assist in selecting appropriate tools for variant annotation. We conducted a MEDLINE search up to November 10, 2023, and included tools that are applicable to a broad range of phenotypes, can be used locally, and have been recently updated. Tools were categorized based on the types of variants they accept and the functional impacts they predict. Sequence Ontology terms were used for standardization. We identified 118 databases and software packages, encompassing 36 variant types and 161 functional impacts. Combining only three tools, namely SnpEff, FAVOR, and SparkINFERNO, allows predicting 99 (61%) distinct functional impacts. Thirty-seven tools predict 89 functional impacts that are not supported by any other tool, while 75 tools predict pathogenicity and can be used within the ACMG/AMP guidelines in a clinical context. We launched a website allowing researchers to select tools based on desired variants and impacts. In summary, more than 100 tools are already available to predict approximately 160 functional impacts. About 60% of the functional impacts can be predicted by the combination of three tools. Unexpectedly, recent tools do not predict more impacts than older ones. Future research should allow predicting the functionality of so far unsupported variant types, such as gene fusions.URL: https://cardio-care.shinyapps.io/VEP_Finder/ .Registration: OSF Registries on November 10, 2023, https://osf.io/s2gct .

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.

Mapping the APOE structurally on missense variants in elderly Brazilians

Cardiovascular diseases (CVDs) pose a significant global health threat, with familial hypercholesterolemia (FH) being a key genetic contributor. The apolipoprotein E (APOE) gene plays a vital role in lipid metabolism, and its variants are associated with CVD risk. This study explores prevalent APOE variants (p.R163C, p.R176C, p.R246C and p.V254E) using genetic and structural analyses. The research, initiated by identifying high-frequency APOE variants through the ABraOM database, utilizes homology modeling and molecular dynamics (MD) simulations to understand the structural consequences. The major lipid-binding region, a critical domain for lipid metabolism, was a focal point. Structural dynamics, including principal component analyses and domain movement analyses, highlighted distinct patterns in APOE variants compared to the wild type (WT). Results revealed significant differences in the structural behavior of variants, particularly in the Major lipid-binding region. The identification of an 'elbow' structure with two states (Elbow I and Elbow II) provided insights into conformational changes. Notably, variants exhibited unique patterns in hydrogen bonding (hb) and hydrophobic interactions, indicating potential functional consequences. The study further associated APOE variants with clinical outcomes, including cognitive impairment and cholesterol levels. Specific variants demonstrated correlations with cognitive decline and variations in lipid profiles, emphasizing their relevance to cardiovascular and neurobiological health. In conclusion, this integrated approach enhances our understanding of APOE variants, shedding light on their role in lipid metabolism and cardiovascular health. The identified structural 'elbows' and their association with clinical outcomes offer a nuanced perspective, guiding future research toward targeted interventions for diseases linked to lipid metabolism and neurobiology.Communicated by Ramaswamy H. Sarma.

Functional pri-miR-34b/c rs4938723 and KRAS 3'UTR rs61764370 SNPs: Novel phenotype modifiers in Li-Fraumeni Syndrome?

PURPOSE

Li-Fraumeni Syndrome (LFS) is a rare cancer predisposing condition caused by germline pathogenic TP53 variants, in which core tumors comprise sarcomas, breast, brain and adrenocortical neoplasms. Clinical manifestations are highly variable in carriers of the Brazilian germline founder variant TP53 p.R337H, possibly due to the influence of modifier genes such as miRNA genes involved in the regulation of the p53 pathway. Herein, we investigated the potential phenotypic effects of two miRNA-related functional SNPs, pri-miR-34b/c rs4938723 and 3'UTR KRAS rs61764370, in a cohort of 273 LFS patients from Southern and Southeastern Brazil.

METHODS

The genotyping of selected SNPs was performed by TaqMan® allelic discrimination and subsequently custom TaqMan® genotyping results were confirmed by Sanger sequencing in all SNP-positive LFS patients.

RESULTS

Although the KRAS SNP showed no effect as a phenotype modulator, the rs4938723 CC genotype was significantly associated with development of LFS non-core tumors (first tumor diagnosis) in p.R337H carriers (p = 0.039). Non-core tumors were also more frequently diagnosed in carriers of germline TP53 DNA binding domain variants harboring the rs4938723 C variant allele. Previous studies described pri-miR-34b/c rs4938723 C as a risk allele for sporadic occurrence of thyroid and prostate cancers (non-core tumors of the LFS spectrum).

CONCLUSION

With this study, we presented additional evidence about the importance of analyzing miRNA genes that could indirectly regulate p53 expression, and, therefore, may modulate the LFS phenotype, such as those of the miR-34 family.

Genomic disparity impacts variant classification of cancer susceptibility genes in Turkish breast cancer patients

OBJECTIVE

Turkish genome is underrepresented in large genomic databases. This study aims to evaluate the effect of allele frequency in the Turkish population in determining the clinical utility of germline findings in breast cancer, including invasive lobular carcinoma (ILC), mixed invasive ductal and lobular carcinoma (IDC-L), and ductal carcinoma (DC).

METHODS

Two clinic-based cohorts from the Umraniye Research and Training Hospital (URTH) were used in this study: a cohort consisting of 132 women with breast cancer and a non-cancer cohort consisting of 492 participants. The evaluation of the germline landscape was performed by analysis of 27 cancer genes. The frequency and type of variants in the breast cancer cohort were compared to those in the non-cancer cohort to investigate the effect of population genetics. The variant allele frequencies in Turkish Variome and gnomAD were statistically evaluated.

RESULTS

The genetic analysis identified 121 variants in the breast cancer cohort (actionable = 32, VUS = 89) and 223 variants in the non-cancer cohort (actionable = 25, VUS = 188). The occurrence of 21 variants in both suggested a possible genetic population effect. Evaluation of allele frequency of 121 variants from the breast cancer cohort showed 22% had a significantly higher value in Turkish Variome compared to gnomAD (p < 0.0001, 95% CI) with a mean difference of 60 times (ranging from 1.37-354.4). After adjusting for variant allele frequency using the ancestry-appropriate database, 6.7% (5/75) of VUS was reclassified to likely benign.

CONCLUSION

To our knowledge, this is the first study of population genetic effects in breast cancer subtypes in Turkish women. Our findings underscore the need for a large genomic database representing Turkish population-specific variants. It further highlights the significance of the ancestry-appropriate population database for accurate variant assessment in clinical settings.

Exploring the frequency of a TP53 polyadenylation signal variant in tumor DNA from patients diagnosed with lung adenocarcinomas, sarcomas and uterine leiomyomas

The TP53 3'UTR variant rs78378222 A>C has been detected in different tumor types as a somatic alteration that reduces p53 expression through modification of polyadenylation and miRNA regulation. Its prevalence is not yet known in all tumors. Herein, we examine tumor tissue prevalence of rs7837822 in Brazilian cohorts of patients from south and southeast regions diagnosed with lung adenocarcinoma (LUAD, n=586), sarcoma (SARC, n=188) and uterine leiomyoma (ULM, n=41). The minor allele (C) was identified in heterozygosity in 6/586 LUAD tumors (prevalence = 1.02 %) and none of the SARC and ULM samples. Additionally, next generation sequencing analysis revealed that all variant-positive tumors (n=4) with sample availability had additional pathogenic or likely pathogenic somatic variants in the TP53 coding regions. Among them, 3/4 (75 %) had the same pathogenic or likely pathogenic sequence variant (allele frequency <0.05 in tumor DNA) namely c.751A>C (p.Ile251Leu). Our results indicate a low somatic prevalence of rs78378222 in LUAD, ULM and SARC tumors from Brazilian patients, which suggests that no further analysis of this variant in the specific studied regions of Brazil is warranted. However, these findings should not exclude tumor molecular testing of this TP53 3'UTR functional variant for different populations.

Waardenburg Syndrome: The Contribution of Next-Generation Sequencing to the Identification of Novel Causative Variants

Waardenburg syndrome (WS) is characterized by hearing loss and pigmentary abnormalities of the eyes, hair, and skin. The condition is genetically heterogeneous, and is classified into four clinical types differentiated by the presence of dystopia canthorum in type 1 and its absence in type 2. Additionally, limb musculoskeletal abnormalities and Hirschsprung disease differentiate types 3 and 4, respectively. Genes PAX3, MITF, SOX10, KITLG, EDNRB, and EDN3 are already known to be associated with WS. In WS, a certain degree of molecularly undetected patients remains, especially in type 2. This study aims to pinpoint causative variants using different NGS approaches in a cohort of 26 Brazilian probands with possible/probable diagnosis of WS1 (8) or WS2 (18). DNA from the patients was first analyzed by exome sequencing. Seven of these families were submitted to trio analysis. For inconclusive cases, we applied a targeted NGS panel targeting WS/neurocristopathies genes. Causative variants were detected in 20 of the 26 probands analyzed, these being five in PAX3, eight in MITF, two in SOX10, four in EDNRB, and one in ACTG1 (type 2 Baraitser-Winter syndrome, BWS2). In conclusion, in our cohort of patients, the detection rate of the causative variant was 77%, confirming the superior detection power of NGS in genetically heterogeneous diseases.

Unveiling the clinical and electrophysiological profile of CMTX6: Insights from two Brazilian families

BACKGROUND AND AIMS

X-linked Charcot-Marie-Tooth disease type 6 (CMTX6) is an extremely rare condition associated with mutations in the PDK3 gene. To date, only three families from different countries have been reported (Australia, South Korea, and Germany). In this study, we sought to provide a comprehensive clinical and electrophysiological characterization of two Brazilian families.

METHODS

We conducted comprehensive clinical assessments, extensive electrophysiological evaluations, and performed whole-exome sequencing in the probands to investigate the genetic basis of the disease.

RESULTS

Males in the family carrying the Arg162His mutation displayed early-onset motor and/or sensory axonal neuropathy, absence of tendon jerks, pes cavus, and frequently reported pain. Females in the same family exhibited a milder phenotype of the disease with later onset and some remained asymptomatic into their 50s. In the unrelated family with a single affected male, the clinical presentation was characterized by severe progressive sensorimotor polyneuropathy accompanied by neuropathic pain.

INTERPRETATION

We report two Brazilian families with CMTX6 including one harboring a previously unpublished variant in the PDK3 gene, which co-segregates with the disease as expected in a X-linked disease. Notably, the clinical presentations across the five families with available descriptions, including our study, share striking similarities. Furthermore, the proximity of the three reported mutations suggests potential functional similarities and common underlying mechanisms. This study contributes to the growing knowledge of CMTX6 and underscores the importance of international collaborations in studying rare genetic disorders.

Characterization of Potential Melanoma Predisposition Genes in High-Risk Brazilian Patients

Increased genetic risk for melanoma can occur in the context of germline pathogenic variants in high-penetrance genes, such as CDKN2A and CDK4, risk variants in low- to moderate-penetrance genes (MC1R and MITF), and possibly due to variants in emerging genes, such as ACD, TERF2IP, and TERT. We aimed to identify germline variants in high- and low- to moderate-penetrance melanoma risk genes in Brazilian patients with clinical criteria for familial melanoma syndrome. We selected patients with three or more melanomas or melanoma patients from families with three tumors (melanoma and pancreatic cancer) in first- or second-degree relatives. Genetic testing was performed with a nine-gene panel (ACD, BAP1, CDK4, CDKN2A, POT1, TERT, TERF2IP, MC1R, and MITF). In 36 patients, we identified 2 (5.6%) with germline pathogenic variants in CDKN2A and BAP1 and 4 (11.1%) with variants of uncertain significance in the high-penetrance genes. MC1R variants were found in 86.5%, and both red hair color variants and unknown risk variants were enriched in patients compared to a control group. The low frequency of germline pathogenic variants in the high-penetrance genes and the high prevalence of MC1R variants found in our cohort show the importance of the MC1R genotype in determining the risk of melanoma in the Brazilian melanoma-prone families.

Identification of a second genetic alteration in patients with SHOX deficiency individuals: a potential explanation for phenotype variability

OBJECTIVE

Our study aimed to assess the impact of genetic modifiers on the significant variation in phenotype that is observed in individuals with SHOX deficiency, which is the most prevalent monogenic cause of short stature.

DESIGN AND METHODS

We performed a genetic analysis in 98 individuals from 48 families with SHOX deficiency with a target panel designed to capture the entire SHOX genomic region and 114 other genes that modulate growth and/or SHOX action. We prioritized rare potentially deleterious variants.

RESULTS

We did not identify potential deleterious variants in the promoter or intronic regions of the SHOX genomic locus. In contrast, we found eight heterozygous variants in 11 individuals from nine families in genes with a potential role as genetic modifiers. In addition to a previously described likely pathogenic (LP) variant in CYP26C1 observed in two families, we identified LP variants in PTHLH and ACAN, and variants of uncertain significance in NPR2, RUNX2, and TP53 in more affected individuals from families with SHOX deficiency. Families with a SHOX alteration restricted to the regulatory region had a higher prevalence of a second likely pathogenic variant (27%) than families with an alteration compromising the SHOX coding region (2.9%, P = .04).

CONCLUSION

In conclusion, variants in genes related to the growth plate have a potential role as genetic modifiers of the phenotype in individuals with SHOX deficiency. In individuals with SHOX alterations restricted to the regulatory region, a second alteration could be critical to determine the penetrance and expression of the phenotype.

Mutational signatures and increased retrotransposon insertions in xeroderma pigmentosum variant skin tumors

Xeroderma pigmentosum variant (XP-V) is an autosomal recessive disease with an increased risk of developing cutaneous neoplasms in sunlight-exposed regions. These cells are deficient in the translesion synthesis (TLS) DNA polymerase eta, responsible for bypassing different types of DNA lesions. From the exome sequencing of 11 skin tumors of a genetic XP-V patients' cluster, classical mutational signatures related to sunlight exposure, such as C>T transitions targeted to pyrimidine dimers, were identified. However, basal cell carcinomas also showed distinct C>A mutation spectra reflecting a mutational signature possibly related to sunlight-induced oxidative stress. Moreover, four samples carry different mutational signatures, with C>A mutations associated with tobacco chewing or smoking usage. Thus, XP-V patients should be warned of the risk of these habits. Surprisingly, higher levels of retrotransposon somatic insertions were also detected when the tumors were compared with non-XP skin tumors, revealing other possible causes for XP-V tumors and novel functions for the TLS polymerase eta in suppressing retrotransposition. Finally, the expected high mutation burden found in most of these tumors renders these XP patients good candidates for checkpoint blockade immunotherapy.

Co-Occurrence of Germline Genomic Variants and Copy Number Variations in Hereditary Breast and Colorectal Cancer Patients

Hereditary Breast and Ovarian Cancer (HBOC) syndrome is an autosomal dominant disease associated with a high risk of developing breast, ovarian, and other malignancies. Lynch syndrome is caused by mutations in mismatch repair genes predisposing to colorectal and endometrial cancers, among others. A rare phenotype overlapping hereditary colorectal and breast cancer syndromes is poorly characterized. Three breast and colorectal cancer unrelated patients fulfilling clinical criteria for HBOC were tested by whole exome sequencing. A family history of colorectal cancer was reported in two patients (cases 2 and 3). Several variants and copy number variations were identified, which potentially contribute to the cancer risk or prognosis. All patients presented copy number imbalances encompassing PMS2 (two deletions and one duplication), a known gene involved in the DNA mismatch repair pathway. Two patients showed gains covering the POLE2 (cases 1 and 3), which is associated with DNA replication. Germline potentially damaging variants were found in PTCH1 (patient 3), MAT1A, and WRN (patient 2). Overall, concurrent genomic alterations were described that may increase the risk of cancer appearance in HBOC patients with breast and colorectal cancers.

Analysis of the Mutational Landscape of Osteosarcomas Identifies Genes Related to Metastasis and Prognosis and Disrupted Biological Pathways of Immune Response and Bone Development

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.

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 Novel Variant in VPS13B Underlying Cohen Syndrome

Pathogenic variants in vacuolar protein sorting 13 homolog B (VPS13B) cause Cohen syndrome (CS), a clinically diverse neurodevelopmental disorder. We used whole exome and Sanger sequencing to identify disease-causing variants in a Pakistani family with intellectual disability, microcephaly, facial dysmorphism, neutropenia, truncal obesity, speech delay, motor delay, and insomnia. We identified a novel homozygous nonsense variant c.8841G > A: p.(W2947∗) in VPS13B (NM_017890.5) which segregated with the disease. Sleep disturbances are commonly seen in neurodevelopmental disorders and can exacerbate medical issues if left untreated. We demonstrate that individuals with Cohen syndrome may also be affected by sleep disturbances. In conclusion, we expand the genetic and phenotypic features of Cohen syndrome in the Pakistani population.

DNA methylation patterns suggest the involvement of DNMT3B and TET1 in osteosarcoma development

DNA methylation may be involved in the development of osteosarcomas. Osteosarcomas commonly arise during the bone growth and remodeling in puberty, making it plausible to infer the involvement of epigenetic alterations in their development. As a highly studied epigenetic mechanism, we investigated DNA methylation and related genetic variants in 28 primary osteosarcomas aiming to identify deregulated driver alterations. Methylation and genomic data were obtained using the Illumina HM450K beadchips and the TruSight One sequencing panel, respectively. Aberrant DNA methylation was spread throughout the osteosarcomas genomes. We identified 3146 differentially methylated CpGs comparing osteosarcomas and bone tissue samples, with high methylation heterogeneity, global hypomethylation and focal hypermethylation at CpG islands. Differentially methylated regions (DMR) were detected in 585 loci (319 hypomethylated and 266 hypermethylated), mapped to the promoter regions of 350 genes. These DMR genes were enriched for biological processes related to skeletal system morphogenesis, proliferation, inflammatory response, and signal transduction. Both methylation and expression data were validated in independent groups of cases. Six tumor suppressor genes harbored deletions or promoter hypermethylation (DLEC1, GJB2, HIC1, MIR149, PAX6, and WNT5A), and four oncogenes presented gains or hypomethylation (ASPSCR1, NOTCH4, PRDM16, and RUNX3). Our analysis also revealed hypomethylation at 6p22, a region that contains several histone genes. Copy-number changes in DNMT3B (gain) and TET1 (loss), as well as overexpression of DNMT3B in osteosarcomas provide a possible explanation for the observed phenotype of CpG island hypermethylation. While the detected open-sea hypomethylation likely contributes to the well-known osteosarcoma genomic instability, enriched CpG island hypermethylation suggests an underlying mechanism possibly driven by overexpression of DNMT3B likely resulting in silencing of tumor suppressors and DNA repair genes.

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.

A systematic evaluation of normalization methods and probe replicability using infinium EPIC methylation data

BACKGROUND

The Infinium EPIC array measures the methylation status of > 850,000 CpG sites. The EPIC BeadChip uses a two-array design: Infinium Type I and Type II probes. These probe types exhibit different technical characteristics which may confound analyses. Numerous normalization and pre-processing methods have been developed to reduce probe type bias as well as other issues such as background and dye bias.

METHODS

This study evaluates the performance of various normalization methods using 16 replicated samples and three metrics: absolute beta-value difference, overlap of non-replicated CpGs between replicate pairs, and effect on beta-value distributions. Additionally, we carried out Pearson's correlation and intraclass correlation coefficient (ICC) analyses using both raw and SeSAMe 2 normalized data.

RESULTS

The method we define as SeSAMe 2, which consists of the application of the regular SeSAMe pipeline with an additional round of QC, pOOBAH masking, was found to be the best performing normalization method, while quantile-based methods were found to be the worst performing methods. Whole-array Pearson's correlations were found to be high. However, in agreement with previous studies, a substantial proportion of the probes on the EPIC array showed poor reproducibility (ICC < 0.50). The majority of poor performing probes have beta values close to either 0 or 1, and relatively low standard deviations. These results suggest that probe reliability is largely the result of limited biological variation rather than technical measurement variation. Importantly, normalizing the data with SeSAMe 2 dramatically improved ICC estimates, with the proportion of probes with ICC values > 0.50 increasing from 45.18% (raw data) to 61.35% (SeSAMe 2).

Protein domains provide a new layer of information for classifying human variations in rare diseases

Introduction: Using the ACMG-AMP guidelines for the interpretation of sequence variants, it remains difficult to meet the criterion associated with the protein domain, PM1, which is assigned in only about 10% of cases, whereas the criteria related to variant frequency, PM2/BA1/BS1, is reported in 50% of cases. To improve the classification of human missense variants using protein domains information, we developed the DOLPHIN system (https://dolphin.mmg-gbit.eu). Methods: We used Pfam alignments of eukaryotes to define DOLPHIN scores to identify protein domain residues and variants that have a significant impact. In parallel, we enriched gnomAD variants frequencies for each domains' residue. These were validated using ClinVar data. Results: We applied this method to all potential human transcripts' variants, resulting in 30.0% being assigned a PM1 label, whereas 33.2% were eligible for a new benign support criterion, BP8. We also showed that DOLPHIN provides an extrapolated frequency for 31.8% of the variants, compared to the original frequency available in gnomAD for 7.6% of them. Discussion: Overall, DOLPHIN allows a simplified use of the PM1 criterion, an expanded application of the PM2/BS1 criteria and the creation of a new BP8 criterion. DOLPHIN could facilitate the classification of amino acid substitutions in protein domains that cover nearly 40% of proteins and represent the sites of most pathogenic variants.

A truncating variant of RAD51B associated with primary ovarian insufficiency provides insights into its meiotic and somatic functions

Primary ovarian insufficiency (POI) causes female infertility by abolishing normal ovarian function. Although its genetic etiology has been extensively investigated, most POI cases remain unexplained. Using whole-exome sequencing, we identified a homozygous variant in RAD51B -(c.92delT) in two sisters with POI. In vitro studies revealed that this variant leads to translation reinitiation at methionine 64. Here, we show that this is a pathogenic hypomorphic variant in a mouse model. Rad51bc.92delT/c.92delT mice exhibited meiotic DNA repair defects due to RAD51 and HSF2BP/BMRE1 accumulation in the chromosome axes leading to a reduction in the number of crossovers. Interestingly, the interaction of RAD51B-c.92delT with RAD51C and with its newly identified interactors RAD51 and HELQ was abrogated or diminished. Repair of mitomycin-C-induced chromosomal aberrations was impaired in RAD51B/Rad51b-c.92delT human and mouse somatic cells in vitro and in explanted mouse bone marrow cells. Accordingly, Rad51b-c.92delT variant reduced replication fork progression of patient-derived lymphoblastoid cell lines and pluripotent reprogramming efficiency of primary mouse embryonic fibroblasts. Finally, Rad51bc.92delT/c.92delT mice displayed increased incidence of pituitary gland hyperplasia. These results provide new mechanistic insights into the role of RAD51B not only in meiosis but in the maintenance of somatic genome stability.

Diagnostic yield of a multigene sequencing approach in children classified as idiopathic short stature

Objective

Most children with short stature remain without an etiologic diagnosis after extensive clinical and laboratory evaluation and are classified as idiopathic short stature (ISS). This study aimed to determine the diagnostic yield of a multigene analysis in children classified as ISS.

Design and methods

We selected 102 children with ISS and performed the genetic analysis as part of the initial investigation. We developed customized targeted panel sequencing, including all genes already implicated in the isolated short-stature phenotype. Rare and deleterious single nucleotide or copy number variants were assessed by bioinformatic tools.

Results

We identified 20 heterozygous pathogenic (P) or likely pathogenic (LP) genetic variants in 17 of 102 patients (diagnostic yield = 16.7%). Three patients had more than one P/LP genetic alteration. Most of the findings were in genes associated with the growth plate differentiation: IHH (n = 4), SHOX (n = 3), FGFR3 (n = 2), NPR2 (n = 2), ACAN (n = 2), and COL2A1 (n = 1) or involved in the RAS/MAPK pathway: NF1 (n = 2), PTPN11 (n = 1), CBL (n = 1), and BRAF (n = 1). None of these patients had clinical findings to guide a candidate gene approach. The diagnostic yield was higher among children with severe short stature (35% vs 12.2% for height SDS ≤ or > -3; P = 0.034). The genetic diagnosis had an impact on clinical management for four children.

Conclusion

A multigene sequencing approach can determine the genetic etiology of short stature in up to one in six children with ISS, removing the term idiopathic from their clinical classification.

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.

Diagnostic reasoning in neurogenetics: a general approach

Establishing the definitive diagnosis of a neurogenetic disease is usually a complex task. However, like any type of clinical diagnostic reasoning, an organized process of development and consideration of diagnostic hypotheses may guide neurologists and medical geneticists to solve this difficult task. The aim of the present review is to propose a general method for diagnostic reasoning in neurogenetics, with the definition of the main neurological syndrome and its associated topographical diagnosis, followed by the identification of major and secondary neurological syndromes, extraneurological findings, and inheritance pattern. We also discuss general rules and knowledge requirements of the ordering physician to request genetic testing and information on how to interpret genetic variants in a genetic report. By guiding the requests for genetic testing according to an organized model of diagnostic reasoning and with the availability of specific treatments, clinicians may find greater resoluteness and efficacy in the diagnostic investigation, shortening the struggle of patients for a definitive diagnosis.

Nemaline Myopathy in Brazilian Patients: Molecular and Clinical Characterization

Nemaline myopathy (NM), a structural congenital myopathy, presents a significant clinical and genetic heterogeneity. Here, we compiled molecular and clinical data of 30 Brazilian patients from 25 unrelated families. Next-generation sequencing was able to genetically classify all patients: sixteen families (64%) with mutation in NEB, five (20%) in ACTA1, two (8%) in KLHL40, and one in TPM2 (4%) and TPM3 (4%). In the NEB-related families, 25 different variants, 11 of them novel, were identified; splice site (10/25) and frame shift (9/25) mutations were the most common. Mutation c.24579 G>C was recurrent in three unrelated patients from the same region, suggesting a common ancestor. Clinically, the “typical” form was the more frequent and caused by mutations in the different NM genes. Phenotypic heterogeneity was observed among patients with mutations in the same gene. Respiratory involvement was very common and often out of proportion with limb weakness. Muscle MRI patterns showed variability within the forms and genes, which was related to the severity of the weakness. Considering the high frequency of NEB mutations and the complexity of this gene, NGS tools should be combined with CNV identification, especially in patients with a likely non-identified second mutation.

Chagas Disease Megaesophagus Patients Carrying Variant MRPS18B P260A Display Nitro-Oxidative Stress and Mitochondrial Dysfunction in Response to IFN-γ Stimulus

Chagas disease (CD), caused by the protozoan parasite Trypanosoma cruzi, affects 8 million people, and around 1/3 develop chronic cardiac (CCC) or digestive disease (megaesophagus/megacolon), while the majority remain asymptomatic, in the indeterminate form of Chagas disease (ASY). Most CCC cases in families with multiple Chagas disease patients carry damaging mutations in mitochondrial genes. We searched for exonic mutations associated to chagasic megaesophagus (CME) in genes essential to mitochondrial processes. We performed whole exome sequencing of 13 CME and 45 ASY patients. We found the damaging variant MRPS18B 688C > G P230A, in five out of the 13 CME patients (one of them being homozygous; 38.4%), while the variant appeared in one out of 45 ASY patients (2.2%). We analyzed the interferon (IFN)-γ-induced nitro-oxidative stress and mitochondrial function of EBV-transformed lymphoblastoid cell lines. We found the CME carriers of the mutation displayed increased levels of nitrite and nitrated proteins; in addition, the homozygous (G/G) CME patient also showed increased mitochondrial superoxide and reduced levels of ATP production. The results suggest that pathogenic mitochondrial mutations may contribute to cytokine-induced nitro-oxidative stress and mitochondrial dysfunction. We hypothesize that, in mutation carriers, IFN-γ produced in the esophageal myenteric plexus might cause nitro-oxidative stress and mitochondrial dysfunction in neurons, contributing to megaesophagus.

Genomic study of nonsyndromic hearing loss in unaffected individuals: Frequency of pathogenic and likely pathogenic variants in a Brazilian cohort of 2,097 genomes

Hearing loss (HL) is a common sensory deficit in humans and represents an important clinical and social burden. We studied whole-genome sequencing data of a cohort of 2,097 individuals from the Brazilian Rare Genomes Project who were unaffected by hearing loss to investigate pathogenic and likely pathogenic variants associated with nonsyndromic hearing loss (NSHL). We found relevant frequencies of individuals harboring these alterations: 222 heterozygotes (10.59%) for sequence variants, 54 heterozygotes (2.58%) for copy-number variants (CNV), and four homozygotes (0.19%) for sequence variants. The top five most frequent genes and their corresponding combined allelic frequencies (AF) were GJB2 (AF = 1.57%), STRC (AF = 1%), OTOA (AF = 0.69%), TMPRSS3 (AF = 0.41%), and OTOF (AF = 0.29%). The most frequent sequence variant was GJB2:c.35del (AF = 0.72%), followed by OTOA:p. (Glu787Ter) (AF = 0.61%), while the most recurrent CNV was a microdeletion of 57.9 kb involving the STRC gene (AF = 0.91%). An important fraction of these individuals (n = 104; 4.96%) presented variants associated with autosomal dominant forms of NSHL, which may imply the development of some hearing impairment in the future. Using data from the heterozygous individuals for recessive forms and the Hardy–Weinberg equation, we estimated the population frequency of affected individuals with autosomal recessive NSHL to be 1:2,222. Considering that the overall prevalence of HL in adults ranges from 4–15% worldwide, our data indicate that an important fraction of this condition may be associated with a monogenic origin and dominant inheritance.

Comparative Analyses of Turkish Variome and Widely Used Genomic Variation Databases for the Evaluation of Rare Sequence Variants in Turkish Individuals: Idiopathic Hypogonadotropic Hypogonadism as a Disease Model

Objective

With the increasing use of whole-exome sequencing, one of the challenges in identifying the causal allele for a Mendelian disease is the lack of availability of population-specific human genetic variation reference databases. The people of Turkey were not represented in GnomAD or other publicly available large databases until recently, when the first comprehensive genomic variation database, Turkish Variome (TRV), was published. The aim of this study was to evaluate whether TRV or other publicly available large genomic variation databases can reliably be used for rare disease variant evaluation in Turkish individuals.

Methods

Sixty non-disease-causing, non-synonymous variants (minor allele frequencies >1%) were identified in 58 genes that are known to be associated with idiopathic hypogonadotropic hypogonadism from a large Turkish patient cohort. The allelic frequencies of these variants were then compared with those in various public genomic variation databases, including TRV.

Results

Our cohort variants showed the highest correlations with those in the TRV, Iranome, and The Greater Middle East Variome, in decreasing order.

Conclusion

These results suggest that the TRV is the appropriate database to use for rare genomic variant evaluations in the Turkish population. Our data also suggest that variomes from geographic neighborhoods may serve as substitute references for populations devoid of their own genomic variation databases.

High Throughput Sequencing to Identify Monogenic Etiologies in a Preselected Polycystic Ovary Syndrome Cohort

Context

Polycystic ovary syndrome (PCOS) etiology remains to be elucidated, but familial clustering and twin studies have shown a strong heritable component.

Objective

The purpose of this study was to identify rare genetic variants that are associated with the etiology of PCOS in a preselected cohort.

Methods

This prospective study was conducted among a selected group of women with PCOS. The study’s inclusion criteria were patients with PCOS diagnosed by the Rotterdam criteria with the following phenotypes: severe insulin resistance (IR), normoandrogenic–normometabolic phenotype, adrenal hyperandrogenism, primary amenorrhea, and familial PCOS. Forty-five patients were studied by target sequencing, while 8 familial cases were studied by whole exome sequencing.

Results

Patients were grouped according to the inclusion criteria with the following distribution: 22 (41.5%) with severe IR, 13 (24.5%) with adrenal hyperandrogenism, 7 (13.2%) with normoandrogenic phenotype, 3 (5.7%) with primary amenorrhea, and 8 (15.1%) familial cases. DNA sequencing analysis identified 1 pathogenic variant in LMNA, 3 likely pathogenic variants in INSR, PIK3R1, and DLK1, and 6 variants of uncertain significance level with interesting biologic rationale in 5 genes (LMNA, GATA4, NR5A1, BMP15, and FSHR). LMNA was the most prevalent affected gene in this cohort (3 variants).

Conclusion

Several rare variants in genes related to IR were identified in women with PCOS. Although IR is a common feature of PCOS, patients with extreme or atypical phenotype should be carefully evaluated to rule out monogenic conditions.

Environmental Contaminants Modulate Breast Cancer Development and Outcome in TP53 p.R337H Carriers and Noncarriers

Two major concerns associated with cancer development in Paraná state, South Brazil, are environmental pollution and the germline TP53 p.R337H variant found in 0.27−0.30% of the population. We assessed breast cancer (BC) risk in rural (C1 and C2) and industrialized (C3) subregions, previously classified by geochemistry, agricultural productivity, and population density. C2 presents lower organochloride levels in rivers and lower agricultural outputs than C1, and lower levels of chlorine anions in rivers and lower industrial activities than C3. TP53 p.R337H status was assessed in 4658 women aged >30 years from C1, C2, and C3, subsequent to a genetic screening (Group 1, longitudinal study). BC risk in this group was 4.58 times higher among TP53 p.R337H carriers. BC prevalence and risk were significantly lower in C2 compared to that in C3. Mortality rate and risk associated with BC in women aged >30 years (n = 8181 deceased women; Group 2) were also lower in C2 than those in C3 and C1. These results suggest that environmental factors modulate BC risk and outcome in carriers and noncarriers.

Variant-level matching for diagnosis and discovery: Challenges and opportunities

Here we describe MyGene2, Geno2MP, VariantMatcher, and Franklin; databases that provide variant-level information and phenotypic features to researchers, clinicians, healthcare providers and patients. Following the footsteps of the Matchmaker Exchange project that connects exome, genome, and phenotype databases at the gene level, these databases have as one goal to facilitate connection to one another using Data Connect, a standard for discovery and search of biomedical data from the Global Alliance for Genomics and Health (GA4GH).

Exome and Tissue-Associated Microbiota as Predictive Markers of Response to Neoadjuvant Treatment in Locally Advanced Rectal Cancer

The clinical and pathological responses to multimodal neoadjuvant therapy in locally advanced rectal cancers (LARCs) remain unpredictable, and robust biomarkers are still lacking. Recent studies have shown that tumors present somatic molecular alterations related to better treatment response, and it is also clear that tumor-associated bacteria are modulators of chemotherapy and immunotherapy efficacy, therefore having implications for long-term survivorship and a good potential as the biomarkers of outcome. Here, we performed whole exome sequencing and 16S ribosomal RNA (rRNA) amplicon sequencing from 44 pre-treatment LARC biopsies from Argentinian and Brazilian patients, treated with neoadjuvant chemoradiotherapy or total neoadjuvant treatment, searching for predictive biomarkers of response (responders, n = 17; non-responders, n = 27). In general, the somatic landscape of LARC was not capable to predict a response; however, a significant enrichment in mutational signature SBS5 was observed in non-responders (p = 0.0021), as well as the co-occurrence of APC and FAT4 mutations (p < 0.05). Microbiota studies revealed a similar alpha and beta diversity of bacteria between response groups. Yet, the linear discriminant analysis (LDA) of effect size indicated an enrichment of Hungatella, Flavonifractor, and Methanosphaera (LDA score ≥3) in the pre-treatment biopsies of responders, while non-responders had a higher abundance of Enhydrobacter, Paraprevotella (LDA score ≥3) and Finegoldia (LDA score ≥4). Altogether, the evaluation of these biomarkers in pre-treatment biopsies could eventually predict a neoadjuvant treatment response, while in post-treatment samples, it could help in guiding non-operative treatment strategies.

Causative Variants for Inherited Cardiac Conditions in a Southeast Asian Population Cohort

Background:

Variable penetrance and late-onset phenotypes are key challenges for classifying causal as well as incidental findings in inherited cardiac conditions. Allele frequencies of variants in ancestry-specific populations, along with clinical variant analysis and interpretation, are critical to determine their true significance.

Methods:

Here, we carefully reviewed and classified variants in genes associated with inherited cardiac conditions based on a population whole-genome sequencing cohort of 4810 Singaporeans representing Southeast Asian ancestries.

Results:

Eighty-nine (1.85%) individuals carried either pathogenic or likely pathogenic variants across 25 genes. Forty-six (51.7%) had variants in causal genes for familial hyperlipidemia, but there were also recurrent variants in SCN5A and MYBPC3 , causal genes for inherited arrhythmia and cardiomyopathy, which, despite previous reports, we determined to lack criteria for pathogenicity.

Conclusions:

Our findings highlight the incidence of disease-related variants in inherited cardiac conditions and emphasize the value of large-scale sequencing in specific ancestries. Follow-up detailed phenotyping and analysis of pedigrees are crucial because assigning pathogenicity will significantly affect clinical management for individuals and their family members.

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.

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.

Whole genome sequencing (WGS) data on non-European and admixed individuals remains scarce. Here, the authors analyse WGS data from 1,171 admixed elderly Brazilians from a census cohort, characterising population-specific genetic variation and exploring the clinical utility of this expanded dataset.

New insights on familial colorectal cancer type X syndrome

Familial colorectal cancer type X (FCCTX) is a heterogeneous colorectal cancer predisposition syndrome that, although displays a cancer pattern similar to Lynch syndrome, is mismatch repair proficient and does not exhibit microsatellite instability. Besides, its genetic etiology remains to be elucidated. In this study we performed germline exome sequencing of 39 cancer-affected patients from 34 families at risk for FCCTX. Variant classification followed the American College of Medical Genetics and Genomics (ACMG) guidelines. Pathogenic/likely pathogenic variants were identified in 17.65% of the families. Rare and potentially pathogenic alterations were identified in known hereditary cancer genes (CHEK2), in putative FCCTX candidate genes (OGG1 and FAN1) and in other cancer-related genes such as ATR, ASXL1, PARK2, SLX4 and TREX1. This study provides novel important clues that can contribute to the understanding of FCCTX genetic basis.

Characterization of the 3'UTR of the BTD gene and identification of regulatory elements and microRNAs

Reduced biotinidase activity is associated with a spectrum of deficiency ranging from total deficiency to heterozygous levels, a finding that is not always explained by the pathogenic variants observed in the BTD gene. The investigation of miRNAs, regulatory elements and variants in the 3’UTR region may present relevance in understanding the genotype-phenotype association. The aims of the study were to characterize the regulatory elements of the 3’UTR of the BTD gene and identify variants and miRNAs which may explain the discrepancies observed between genotype and biochemical phenotype. We evaluated 92 individuals with reduced biotinidase activity (level of heterozygotes = 33, borderline = 35, partial DB = 20 or total DB= 4) with previously determined BTD genotype. The 3’UTR of the BTD gene was Sanger sequenced. In silico analysis was performed to identify miRNAs and regulatory elements. No variants were found in the 3’UTR. We found 97 possible miRNAs associated with the BTD gene, 49 predicted miRNAs involved in the alanine, biotin, citrate and pyruvate metabolic pathways and 5 genes involved in biotin metabolism. Six AU-rich elements were found. Our data suggest variants in the 3'UTR of BTD do not explain the genotype-phenotype discrepancies found in Brazilian individuals with reduced biotinidase.

Haplotypes of single cancer driver genes and their local ancestry in a highly admixed long-lived population of Northeast Brazil

Admixed populations have not been examined in detail in cancer genetic studies. Here, we inferred the local ancestry of cancer-associated single nucleotide polymorphisms (SNPs) and haplotypes of a highly admixed Brazilian population. SNP array was used to genotype 73 unrelated individuals aged 80-102 years. Local ancestry inference was performed by merging genotyped regions with phase three data from the 1000 Genomes Project Consortium using RFmix. The average ancestry tract length was 9.12-81.71 megabases. Strong linkage disequilibrium was detected in 48 haplotypes containing 35 SNPs in 10 cancer driver genes. All together, 19 risk and eight protective alleles were identified in 23 out of 48 haplotypes. Homozygous individuals were mainly of European ancestry, whereas heterozygotes had at least one Native American and one African ancestry tract. Native-American ancestry for homozygous individuals with risk alleles for HNF1B, CDH1, and BRCA1 was inferred for the first time. Results indicated that analysis of SNP polymorphism in the present admixed population has a high potential to identify new ancestry-associated alleles and haplotypes that modify cancer susceptibility differentially in distinct human populations. Future case-control studies with populations with a complex history of admixture could help elucidate ancestry-associated biological differences in cancer incidence and therapeutic outcomes.

Polymorphisms at CYP enzymes, NR1I2 and NR1I3 in association with virologic response to antiretroviral therapy in Brazilian HIV-positive individuals

Virologic failure of antiretroviral therapy (ART) may be explained by single nucleotide polymorphisms (SNPs) in drug absorption and metabolism genes. Here, we characterized the associations between polymorphisms in cytochrome P450 enzymes' genes CYP2B6 and CYP3A4/A5, nuclear receptor genes NR1I2/3, and initial ART efficacy among 203 HIV-positive individuals from Rio de Janeiro. Association between SNPs and virologic control was evaluated after 6 and 12 months of follow-up using Cox regression models. The SNP rs2307424 (NR1I3) was associated with increased virologic response after 12 months of treatment, while rs1523127 (NR1I2), rs3003596, and rs2502815 (NR1I3) were associated with decreased response. Increased virologic response after 12 months (_adjHR = 1.54; p = 0.02) was also observed among carriers of the NR1I3 haplotype rs2502815G-rs3003596A-rs2307424A versus the reference haplotype G-A-G. Our results suggest that NR1I2 and NR1I3 variants are associated with virologic responses to ART among Brazilians.

In silico analysis of potential off-target sites to gene editing for Mucopolysaccharidosis type I using the CRISPR/Cas9 system: Implications for population-specific treatments

Mucopolysaccharidosis type I (MPS I) is caused by alpha-L-iduronidase deficiency encoded by the IDUA gene. Therapy with CRISPR/Cas9 is being developed for treatment, however a detailed investigation of off-target effects must be performed. This study aims to evaluate possible off-targets for a sgRNA aiming to correct the most common variant found in MPS I patients (p.Trp402*). A total of 272 potential off-target sequences was obtained and 84 polymorphic sites were identified in these sequences with a frequency equal to or greater than 1% in at least one of the populations. In the majority of cases, polymorphic sites decrease the chance of off-target cleavage and a new PAM was created, which indicates the importance of such analysis. This study highlights the importance of screening off-targets in a population-specific context using Mucopolysaccharidosis type I as an example of a problem that concerns all therapeutic treatments. Our results can have broader applications for other targets already clinically in use, as they could affect CRISPR/Cas9 safety and efficiency.

Pharmacogenetic differentiation across Latin America

Aim: To explore the pharmacogenetic differentiation across Latin American populations, using the fixation index statistics (F_ST). Materials & methods: F_ST analyses were applied to 1519 pharmacogenetic markers in the 1000 Genomes Admixed American superpopulation (1KG_AMR) and an admixed Brazilian sample. Results: Allele-specific F_ST values for the overall cohort point to little overall pharmacogenetic differentiation (average F_ST = 0.017); however, moderate differentiation (F_ST = 0.05-0.15) was observed for 83 markers, while large differentiation (F_ST = 0.15-0.25) was restricted to three markers. Pairwise F_ST analysis identified three markers with very large differentiation (F_ST >0.25). Conclusion: The present study verifies and extends previous reports of little overall pharmacogenetic divergence across Latin America, although a number of markers display substantial differentiation.

Neuroprogenitor Cells From Patients With TBCK Encephalopathy Suggest Deregulation of Early Secretory Vesicle Transport

Biallelic pathogenic variants in TBCK cause encephaloneuropathy, infantile hypotonia with psychomotor retardation, and characteristic facies 3 (IHPRF3). The molecular mechanisms underlying its neuronal phenotype are largely unexplored. In this study, we reported two sisters, who harbored biallelic variants in TBCK and met diagnostic criteria for IHPRF3. We provided evidence that TBCK may play an important role in the early secretory pathway in neuroprogenitor cells (iNPC) differentiated from induced pluripotent stem cells (iPSC). Lack of functional TBCK protein in iNPC is associated with impaired endoplasmic reticulum-to-Golgi vesicle transport and autophagosome biogenesis, as well as altered cell cycle progression and severe impairment in the capacity of migration. Alteration in these processes, which are crucial for neurogenesis, neuronal migration, and cytoarchitecture organization, may represent an important causative mechanism of both neurodevelopmental and neurodegenerative phenotypes observed in IHPRF3. Whether reduced mechanistic target of rapamycin (mTOR) signaling is secondary to impaired TBCK function over other secretory transport regulators still needs further investigation.

The recurrent homozygous translation start site variant in CCDC134 in an individual with severe osteogenesis imperfecta of non-Morrocan ancestry

Osteogenesis imperfecta (OI) is a rare low-bone mass skeletal Mendelian disorder characterized by bone fragility leading to bone fractures, with deformities and stunted growth in the more severe phenotypes. Other common, nonskeletal findings include blue sclerae and dentinogenesis imperfecta. It is caused mainly by quantitative or structural defects in type I collagen, although dysregulation of different signaling pathways that play a role in bone morphogenesis has been described to be associated with a small fraction of individuals with OI. Recently, a homozygous variant in the translation start site of CCDC134, showing increased activation of the RAS/MAPK signaling pathway, has been reported in three families of Moroccan origin with a severe, deforming form of OI. We report on a 9-year-old Brazilian boy, harboring the same homozygous variant in CCDC134, also presenting severe bone involvement. This report contributes to the phenotypic delineation of this novel autosomal recessive form of OI, which presents with high prevalence of nonunion fractures considered rare events in OI in general. In addition, it expands the phenotype to include base skull anomalies, potentially leading to serious complications, as seen in severe forms of OI. A poor response to bisphosphonate therapy was observed in these individuals. As the variant in CCDC134 leads to dysregulation of the RAS/MAPK signaling pathway, drugs targeted to this pathway could be an alternative to achieve a better management of these individuals.

Identification of Variants Responsible for Monogenic Forms of Diabetes in Brazil

Monogenic forms of diabetes mellitus may affect a significant number of patients of this disease, and it is an important molecular cause to be investigated. However, studies of the genetic causes of monogenic diabetes, especially in populations with mixed ethnic backgrounds, such as the one in Brazil, are scarce. The aim of this study was to screen several genes associated with monogenic diabetes in fifty-seven Brazilian patients with recurrence of the disease in their families and thirty-four relatives. Inclusion criteria were: Age of onset ≤ 40 years old, BMI < 30 kg/m², at least two affected generations and negative anti-GAD and anti-IA2 antibodies. MODY genes HNF4A, GCK, HNF1A, HNF1B, NEUROD1, KLF11, PAX4, INS, KCNJ11, and MT-TL1 were sequenced by Sanger sequencing. We identified a total of 20 patients with variants, 13 GCK-MODY, four HNF1A-MODY, and one variant in each of the following genes, HNF4A, HNF1B and MT-TL1. Segregation analysis was performed in 13 families. Four variants were novel, two in GCK (p.(Met115Val) [c.343A>G] and p.(Asp365GlufsTer95) [c.1094_1095insGCGA]) and two in HNF1A (p.(Tyr163Ter) [c.489C>G] and p.(Val380CysfsTer39) [c.1136_1137insC]). Here we highlight the importance of screening for monogenic diabetes in admixed populations.