Improving population scale statistical phasing with whole-genome sequencing data.

Haplotype estimation, or phasing, has gained significant traction in large-scale projects due to its valuable contributions to population genetics, variant analysis, and the creation of reference panels for imputation and phasing of new samples. To scale with the growing number of samples, haplotype estimation methods designed for population scale rely on highly optimized statistical models to phase genotype data, and usually ignore read-level information. Statistical methods excel in resolving common variants, however, they still struggle at rare variants due to the lack of statistical information. In this study we introduce SAPPHIRE, a new method that leverages whole-genome sequencing data to enhance the precision of haplotype calls produced by statistical phasing. SAPPHIRE achieves this by refining haplotype estimates through the realignment of sequencing reads, particularly targeting low-confidence phase calls. Our findings demonstrate that SAPPHIRE significantly enhances the accuracy of haplotypes obtained from state of the art methods and also provides the subset of phase calls that are validated by sequencing reads. Finally, we show that our method scales to large data sets by its successful application to the extensive 3.6 Petabytes of sequencing data of the last UK Biobank 200,031 sample release.

Risk of anthracycline-induced cardiac dysfunction in adolescent and young adult (AYA) cancer survivors: role of genetic susceptibility loci.

There is a known genetic susceptibility to anthracycline-induced cardiac dysfunction in childhood cancer survivors, but this has not been adequately shown in adolescent and young adult (AYA) patients. Our aim was to determine if the previously identified variants associated with cardiac dysfunction in childhood cancer patients affect AYA cancer patients similarly. Forty-five variants were selected for analysis in 253 AYAs previously treated with anthracyclines. We identified four variants that were associated with cardiac dysfunction: SLC10A2:rs7319981 (p = 0.017), SLC22A17:rs4982753 (p = 0.019), HAS3:rs2232228 (p = 0.023), and RARG:rs2229774 (p = 0.050). HAS3:rs2232228 and SLC10A2:rs7319981 displayed significant effects in our AYA cancer survivor population that were in the opposite direction than that reported in childhood cancer survivors. Genetic variants in the host genes were further analyzed for additional associations with cardiotoxicity in AYA cancer survivors. The host genes were then evaluated in a panel of induced pluripotent stem cell-derived cardiomyocytes to assess changes in levels of expression when treated with doxorubicin. Significant upregulation of HAS3 and SLC22A17 expression was observed (p < 0.05), with non-significant anthracycline-responsivity observed for RARG. Our study demonstrates that there is a genetic influence on cardiac dysfunction in AYA cancer patients, but there may be a difference in the role of genetics between childhood and AYA cancer survivors.

Transcription factor 7 like 2 gene polymorphism and advanced glycation end products as risk factors for diabetic nephropathy.

Hepatocellular carcinoma (HCC) is a multifactorial disease with both genetic and environmental factors contributing to its pathogenesis. ACYP2 is a gene that is related to cell differentiation, apoptosis and prevention of malignant tumors. The ACYP2 gene also affects telomere length. The aim of this study was to evaluate the association between ACYP2 single nucleotide polymorphisms (SNPs) (rs843711), and (rs843706) and incidence of HCC in Egyptian HCC patients. The study included 30 patients with HCC and 30 normal controls. Detection of ACYP2 gene SNPs rs843711, and rs843706 in all study participants was done using real time polymerase chain reaction (RT-PCR). The results showed that all participants including HCC patients and controls carried the heterozygous CA (100%) of the rs843706 SNP (p> 0.05). As for the rs843711, 3.3% of HCC patients had the homozygous TT genotype, 46.7% had the heterozygous CT genotype and 50% had the wild CC genotype, while in the control group, 60% had the heterozygous CT genotype and 40% had the wild CC genotype with no significant difference between both groups (p>0.05). We concluded that there was no association between SNPs ACYP2 rs843706 and rs843711 and occurrence of HCC.

Association of Body Mass Index and Parkinson Disease: A Bidirectional Mendelian Randomization Study.

BACKGROUND AND OBJECTIVES: The role of body mass index (BMI) in Parkinson disease (PD) is unclear. Based on the Comprehensive Unbiased Risk Factor Assessment for Genetics and Environment in PD (Courage-PD) consortium, we used 2-sample Mendelian randomization (MR) to replicate a previously reported inverse association of genetically predicted BMI with PD and investigated whether findings were robust in analyses addressing the potential for survival and incidence-prevalence biases. We also examined whether the BMI-PD relation is bidirectional by performing a reverse MR. METHODS: We used summary statistics from a genome-wide association study (GWAS) to extract the association of 501 single-nucleotide polymorphisms (SNPs) with BMI and from the Courage-PD and international Parkinson Disease Genomics Consortium (iPDGC) to estimate their association with PD. Analyses are based on participants of European ancestry. We used the inverse-weighted method to compute odds ratios (ORIVW per 4.8 kg/m2 [95% CI]) of PD and additional pleiotropy robust methods. We performed analyses stratified by age, disease duration, and sex. For reverse MR, we used SNPs associated with PD from 2 iPDGC GWAS to assess the effect of genetic liability toward PD on BMI. RESULTS: Summary statistics for BMI are based on 806,834 participants (54% women). Summary statistics for PD are based on 8,919 (40% women) cases and 7,600 (55% women) controls from Courage-PD, and 19,438 (38% women) cases and 24,388 (51% women) controls from iPDGC. In Courage-PD, we found an inverse association between genetically predicted BMI and PD (ORIVW 0.82 [0.70-0.97], p = 0.012) without evidence for pleiotropy. This association tended to be stronger in younger participants (≤67 years, ORIVW 0.71 [0.55-0.92]) and cases with shorter disease duration (≤7 years, ORIVW 0.75 [0.62-0.91]). In pooled Courage-PD + iPDGC analyses, the association was stronger in women (ORIVW 0.85 [0.74-0.99], p = 0.032) than men (ORIVW 0.92 [0.80-1.04], p = 0.18), but the interaction was not statistically significant (p-interaction = 0.48). In reverse MR, there was evidence for pleiotropy, but pleiotropy robust methods showed a significant inverse association. DISCUSSION: Using an independent data set (Courage-PD), we replicate an inverse association of genetically predicted BMI with PD, not explained by survival or incidence-prevalence biases. Moreover, reverse MR analyses support an inverse association between genetic liability toward PD and BMI, in favor of a bidirectional relation.

Impact of Toll-like receptor 2 (rs5743708) gene polymorphism in pediatric pneumonia: risk and severity.

Pediatric pneumonia is a common respiratory infection that affects children and is thought to be a major source of mortality and morbidity worldwide, particularly in low- and middle-income nations. Toll-like receptor2 (TLR2) is an important receptor involved in the recognition of bacterial pathogens and the activation of the immune response. Genetic variability in TLR2 may partially explain individual differences in susceptibility to infections. The purpose of this study was to investigate the possible contribution of the TLR2 (rs5743708) variant to the risk and severity of pediatric pneumonia infection. The study included 100 pediatric patients diagnosed with pneumonia and 100 normal controls who were age and gender matched. Real-time polymerase chain reaction (PCR) was used to genotype participants for the TLR2 (rs5743708) variant. The analysis revealed that children with the TLR2 (rs5743708) (G/A) genotype showed a 2.52-fold greater risk of having pneumonia (OR: 2.52; 95% CI: 1.32-4.79; p = 0.005) in comparison with patients who have wild homozygous genotypes. Furthermore, we observed that the TLR2 (rs5743708) (A) allele is connected to a greater risk of pneumonia infection in children (OR: 1.612; 95% CI: 1.07-2.43; p = 0.023) but did not significantly influence infection severity. In conclusion, children with the TLR2 (rs5743708) mutant (G/A) genotype and (A) allele had a significantly higher risk of having pneumonia, but they were not at high risk for the severity of the infection.

Contribution of plasma levels of VEGF-A and angiopoietin-2 in addition to a genetic variant in KCNAB1 to predict the risk of bevacizumab-induced hypertension.

Bevacizumab-induced hypertension poses a therapeutic challenge and identifying biomarkers for hypertension can enhance therapy safety. Lower plasma levels of VEGF-A, angiopoietin-2, and rs6770663 in KCNAB1 were previously associated with increased risk of bevacizumab-induced hypertension. This study investigated whether these factors independently contribute to grade 2-3 bevacizumab-induced hypertension risk in 277 cancer patients (CALGB/Alliance 90401). Multivariable analyses assessed the independent association of each factor and hypertension. Likelihood ratio test (LRT) evaluated the explanatory significance of combining protein levels and rs6770663 in predicting hypertension. Boostrap was employed to assess the mediation effect of protein levels on the rs6770663 association with hypertension. Lower protein levels and rs6770663 were independently associated with increased hypertension risk. Adding rs6770663 to protein levels improved the prediction of hypertension (LRT p = 0.0002), with no mediation effect observed. Protein levels of VEGF-A, angiopoietin-2 and rs6770663 in KCNAB1 are independent risk factors and, when combined, may improve prediction of bevacizumab-induced hypertension. ClinicalTrials.gov Identifier: NCT00110214.

Identification of consensus homozygous regions and their associations with growth and feed efficiency traits in American mink.

The recent chromosome-based genome assembly and the newly developed 70K single nucleotide polymorphism (SNP) array for American mink (Neogale vison) facilitate the identification of genetic variants underlying complex traits in this species. The objective of this study was to evaluate the association between consensus runs of homozygosity (ROH) with growth and feed efficiency traits in American mink. A subsample of two mink populations (n = 2,986) were genotyped using the Affymetrix Mink 70K SNP array. The identified ROH segments were included simultaneously, concatenated into consensus regions, and the ROH-based association studies were carried out with linear mixed models considering a genomic relationship matrix for 11 growth and feed efficiency traits implemented in ASReml-R version 4. In total, 298,313 ROH were identified across all individuals, with an average length and coverage of 4.16 Mb and 414.8 Mb, respectively. After merging ROH segments, 196 consensus ROH regions were detected and used for genome-wide ROH-based association analysis. Thirteen consensus ROH regions were significantly (P < 0.01) associated with growth and feed efficiency traits. Several candidate genes within the significant regions are known for their involvement in growth and body size development, including MEF2A, ADAMTS17, POU3F2, and TYRO3. In addition, we found ten consensus ROH regions, defined as ROH islands, with frequencies over 80% of the population. These islands harbored 12 annotated genes, some of which were related to immune system processes such as DTX3L, PARP9, PARP14, CD86, and HCLS1. This is the first study to explore the associations between homozygous regions with growth and feed efficiency traits in American mink. Our findings shed the light on the effects of homozygosity in the mink genome on growth and feed efficiency traits, that can be utilized in developing a sustainable breeding program for mink.

Combined Genome-Wide Association Study and Transcriptome Analysis Reveal Candidate Genes for Resistance to Rust (Puccinia graminis) in Dactylis glomerata.

Rust disease is a common plant disease that can cause wilting, slow growth of plant leaves, and even affect the growth and development of plants. Orchardgrass (Dactylis glomerata L.) is native to temperate regions of Europe, which has been introduced as a superior forage grass in temperate regions worldwide. Orchardgrass has rich genetic diversity and is widely distributed in the world, which may contain rust resistance genes not found in other crops. Therefore, we collected a total of 333 orchardgrass accessions from different regions around the world. Through a genome-wide association study (GWAS) analysis conducted in four different environments, 91 genes that overlap or are adjacent to significant single nucleotide polymorphisms (SNPs) were identified as potential rust disease resistance genes. Combining transcriptome data from susceptible (PI292589) and resistant (PI251814) accessions, the GWAS candidate gene DG5C04160.1 encoding glutathione S-transferase (GST) was found to be important for orchardgrass rust (Puccinia graminis) resistance. Interestingly, by comparing the number of GST gene family members in seven species, it was found that orchardgrass has the most GST gene family members, containing 119 GST genes. Among them, 23 GST genes showed significant differential expression after inoculation with the rust pathogen in resistant and susceptible accessions; 82% of the genes still showed significantly increased expression 14 days after inoculation in resistant accessions, while the expression level significantly decreased in susceptible accessions. These results indicate that GST genes play an important role in orchardgrass resistance to rust (P. graminis) stress by encoding GST to reduce its oxidative stress response.

Identification of candidate gene associated with maize northern leaf blight resistance in a multi-parent population.

KEY MESSAGE: QTL mapping combined with genome-wide association studies, revealed a potential candidate gene for  resistance to northern leaf blight in the tropical CATETO-related maize line YML226, providing a basis for marker-assisted selection of maize varieties Northern leaf blight (NLB) is a foliar disease that can cause severe yield losses in maize. Identifying and utilizing NLB-resistant genes is the most effective way to prevent and control this disease. In this study, five important inbred lines of maize were used as parental lines to construct a multi-parent population for the identification of NLB-resistant loci. QTL mapping and GWAS analysis revealed that QTL qtl_YML226_1, which had the largest phenotypic variance explanation (PVE) of 9.28%, and SNP 5-49,193,921 were co-located in the CATETO-related line YML226. This locus was associated with the candidate gene Zm00001d014471, which encodes a pentatricopeptide repeat (PPR) protein. In the coding region of Zm00001d014471, YML226 had more specific SNPs than the other parental lines. qRT-PCR showed that the relative expressions of Zm00001d014471 in inoculated and uninoculated leaves of YML226 were significantly higher, indicating that the expression of the candidate gene was correlated with NLB resistance. The analysis showed that the higher expression level in YML226 might be caused by SNP mutations. This study identified NLB resistance candidate loci and genes in the tropical maize inbred line YML226 derived from the CATETO germplasm, thereby providing a theoretical basis for using modern marker-assisted breeding techniques to select genetic resources resistant to NLB.

Association of interleukin6 rs1800796 gene polymorphism and serum level with bladder cancer in Egyptian population.

BACKGROUND: Urinary bladder cancer (UBC)is a common tumor of the urinary tract. OBJECTIVES: To assess the diagnostic significance of IL6 rs1800796 gene polymorphism and IL6 serum level among Egyptian patients with UBC. DESIGN AND METHODS: One hundred patients with UBC were selected from the Mansoura Urology and Nephrology Center, in addition to 100 healthy control subjects; using PCR and ELISA techniques for IL6 detection. RESULTS: The rs1800796 GC, CC genotypes, and C allele were significantly more prevalent in the cases with bladder cancer compared to the healthy group (p < 0.001, = 0.021, < 0.001 respectively). There was a clear association between elevated levels of IL6 and bladder cancer versus the control group (median = 4.2, 0.89 respectively, p < 0.001). Serum IL6 levels showed significantly higher levels in patients carrying CC, followed by GC then GG genotypes. No significant association was found between IL6 rs1800796 gene polymorphism or serum level with demographic or laboratory data. CONCLUSION: It is suggested that there is a clear link between elevated IL6 levels as well as IL6 rs1800796 gene polymorphism with bladder cancer, suggesting their potential utility as biomarkers for the disease.

Modifiable risk factors, oxidative stress markers, and SOD2 rs4880 SNP in coronary artery disease: an association study.

BACKGROUND: Coronary artery disease (CAD) has been linked to single nucleotide polymorphism (SNP) in superoxide dismutase 2 (SOD 2) gene. Additionally, several modifiable risk factors are also known to influence the CAD risk. AIM: To investigate the association between selected modifiable risk factors and oxidative stress markers with the SOD2 rs4880 SNP in CAD patients. METHODS: A cohort of 150 angiographically confirmed CAD patients, and 100 control subjects in the same geographic area were enrolled. SOD levels and lipid peroxidation were assessed in the blood samples using standard protocols. The genotyping of the SOD2 gene was conducted through the PCR-sequencing method. RESULTS: This study indicated that CAD patients with the rs4880 SNP having heterozygous AG and mutated homozygous GG genotypes have increased oxidative stress, decreased SOD activity, and a positive association with CAD risk (OR 2.85) in comparison with control individuals. The investigation among CAD patients was then carried out based on modifiable risk factors. The risk factors selected were clinical characteristics, physical habits, nutritional status, and body mass index. In all the cases, MDA levels showed a positive association, and SOD activity showed a negative association with the selected polymorphism. CONCLUSIONS: The study suggests that the selected modifiable risk factors have an important role in the higher oxidative stress found in patients, which may lead to SOD2 polymorphism. It also suggests that the SOD2 locus can be identified as a marker gene for CAD susceptibility.

Regulation of MicroRNA-4697-3p by Parkinson's disease-associated SNP rs329648 and its impact on SNCA112 mRNA.

BACKGROUND: Parkinson's disease (PD) is a common neurodegenerative disorder characterized by a multifaceted genetic foundation. Genome-Wide Association Studies (GWAS) have played a crucial role in pinpointing genetic variants linked to PD susceptibility. Current study aims to delve into the mechanistic aspects through which the PD-associated Single Nucleotide Polymorphism (SNP) rs329648, identified in prior GWAS, influences the pathogenesis of PD. METHODS AND RESULTS: Employing the CRISPR/Cas9-mediated genome editing mechanism, we demonstrated the association of the disease-associated allele of rs329648 with increased expression of miR-4697-3p in differentiated SH-SY5Y cells. We revealed that miR-4697-3p contributes to the formation of high molecular weight complexes of α-Synuclein (α-Syn), indicative of α-Syn aggregate formation, as evidenced by Western blot analysis. Furthermore, our study unveiled that miR-4697-3p elevates SNCA112 mRNA levels. The resultant protein product, α-Syn 112, a variant of α-Syn with 112 amino acids, is recognized for augmenting α-Syn aggregation. Notably, this regulatory effect minimally impacts the levels of full-length SNCA140 mRNA, as evidenced by qRT-PCR. Additionally, we observed a correlation between the disease-associated allele and miR-4697-3p with increased cell death, substantiated by assessments including cell viability assays, alterations in cell morphology, and TUNEL assays. CONCLUSION: Our research reveals that the disease-associated allele of rs329648 is linked to higher levels of miR-4697-3p. This increase in miR-4697-3p leads to elevated SNCA112 mRNA levels, consequently promoting the formation of α-Syn aggregates. Furthermore, miR-4697-3p appears to play a role in increased cell death, potentially contributing to the pathogenesis of PD.

Association of histidine-rich glycoprotein C633T single nucleotide polymorphism and recurrent miscarriage in Iranian women.

BACKGROUND: Recurrent miscarriage (RM) is defined as the occurrence of at least two or three subsequent miscarriages within the 20th -24th weeks of pregnancy. The primary objective of this study was to investigate whether histidine-rich glycoprotein C633T single nucleotide polymorphism (HRG C633T SNP) statistically correlates with the occurrence of RM among Iranian women. METHODS AND RESULTS: Blood samples from 200 women were taken at the outset of the study. Then, the blood samples of 100 women who had a record of RM (case group) were compared with the other 100 women's blood samples who had no record of RM (control group). Following DNA extraction, the polymorphism of histidine-rich glycoprotein C633T (HRG C633T) for every case was specified and all women were genotyped by the amplification-refractory mutation system (ARMS) method. The results of the study revealed that there was a statistically significant difference between T/T genotype (OR = 3.5, CI (1.39-8.77), p = 0.007), and C/T genotype (OR = 1.83, CI (0.99-3.37), p = 0.05) in the case and control groups. Also, a statistically significant association was observed in T allelic frequency in the RM participants compared to the control group (OR = 2.01, CI (1.31-3.09), p = 0.01). CONCLUSIONS: The present study determined that there was a statistically significant relationship between HRG C633T SNP and increased RM regarding allelic and genotypical aspects. Moreover, it became apparent that women with homozygous T/T genotype were more susceptible to the risk of RM.

Enhancing SNV identification in whole-genome sequencing data through the incorporation of known genetic variants into the minimap2 index.

MOTIVATION: Alignment of reads to a reference genome sequence is one of the key steps in the analysis of human whole-genome sequencing data obtained through Next-generation sequencing (NGS) technologies. The quality of the subsequent steps of the analysis, such as the results of clinical interpretation of genetic variants or the results of a genome-wide association study, depends on the correct identification of the position of the read as a result of its alignment. The amount of human NGS whole-genome sequencing data is constantly growing. There are a number of human genome sequencing projects worldwide that have resulted in the creation of large-scale databases of genetic variants of sequenced human genomes. Such information about known genetic variants can be used to improve the quality of alignment at the read alignment stage when analysing sequencing data obtained for a new individual, for example, by creating a genomic graph. While existing methods for aligning reads to a linear reference genome have high alignment speed, methods for aligning reads to a genomic graph have greater accuracy in variable regions of the genome. The development of a read alignment method that takes into account known genetic variants in the linear reference sequence index allows combining the advantages of both sets of methods. RESULTS: In this paper, we present the minimap2_index_modifier tool, which enables the construction of a modified index of a reference genome using known single nucleotide variants and insertions/deletions (indels) specific to a given human population. The use of the modified minimap2 index improves variant calling quality without modifying the bioinformatics pipeline and without significant additional computational overhead. Using the PrecisionFDA Truth Challenge V2 benchmark data (for HG002 short-read data aligned to the GRCh38 linear reference (GCA_000001405.15) with parameters k = 27 and w = 14) it was demonstrated that the number of false negative genetic variants decreased by more than 9500, and the number of false positives decreased by more than 7000 when modifying the index with genetic variants from the Human Pangenome Reference Consortium.

Mapping of Aegilops speltoides derived leaf rust and stripe rust resistance genes using 35K SNP array.

Wheat is an essential food commodity cultivated throughout the world. However, this crop faces continuous threats from fungal pathogens, leaf rust (LR) and stripe rust (YR). To continue feeding the growing population, these major destructors of wheat must be effectively countered by enhancing the genetic diversity of cultivated germplasm. In this study, an introgression line with hexaploid background (ILsp3603) carrying resistance against Pt pathotypes 77-5 (121R63-1), 77-9 (121R60-1) and Pst pathotypes 46S119 (46E159), 110S119 (110E159), 238S119 (238E159) was developed from donor wheat wild progenitor, Aegilops speltoides acc pau 3603. To understand the genetic basis of resistance and map these genes (named Lrsp3603 and Yrsp3603), inheritance studies were carried out in F6 and F7 mapping population, developed by crossing ILsp3603 with LR and YR susceptible cultivar WL711, which revealed a monogenic (single gene) inheritance pattern for each of these traits. Bulk segregant analysis combined with 35 K Axiom SNP array genotyping mapped both genes as separate entities on the short arm of chromosome 6B. A genetic linkage map, comprising five markers, 1 SNP, 1 PLUG and three gene based SSRs, covered a genetic distance of 12.65 cM. Lrsp3603 was flanked by markers Tag-SSR14 (located proximally at 2.42 cM) and SNP AX-94542331 (at 3.28 cM) while Yrsp3603 was mapped at one end closest to AX-94542331 at 6.62 cM distance. Functional annotation of Lrsp3603 target region (∼ 1 Mbp) revealed 10 gene IDs associated with disease resistance mechanisms including three encoding typical R gene domains.

Post-implantation analysis of genomic variations in the progeny from developing fetus to birth.

The analysis of genomic variations in offspring after implantation has been infrequently studied. In this study, we aim to investigate the extent of de novo mutations in humans from developing fetus to birth. Using high-depth whole-genome sequencing, 443 parent-offspring trios were studied to compare the results of de novo mutations (DNMs) between different groups. The focus was on fetuses and newborns, with DNA samples obtained from the families' blood and the aspirated embryonic tissues subjected to deep sequencing. It was observed that the average number of total DNMs in the newborns group was 56.26 (54.17-58.35), which appeared to be lower than that the multifetal reduction group, which was 76.05 (69.70-82.40) (F = 2.42, P = 0.12). However, after adjusting for parental age and maternal pre-pregnancy body mass index (BMI), significant differences were found between the two groups. The analysis was further divided into single nucleotide variants (SNVs) and insertion/deletion of a small number of bases (indels), and it was discovered that the average number of de novo SNVs associated with the multifetal reduction group and the newborn group was 49.89 (45.59-54.20) and 51.09 (49.22-52.96), respectively. No significant differences were noted between the groups (F = 1.01, P = 0.32). However, a significant difference was observed for de novo indels, with a higher average number found in the multifetal reduction group compared to the newborn group (F = 194.17, P < 0.001). The average number of de novo indels among the multifetal reduction group and the newborn group was 26.26 (23.27-29.05) and 5.17 (4.82-5.52), respectively. To conclude, it has been observed that the quantity of de novo indels in the newborns experiences a significant decrease when compared to that in the aspirated embryonic tissues (7-9 weeks). This phenomenon is evident across all genomic regions, highlighting the adverse effects of de novo indels on the fetus and emphasizing the significance of embryonic implantation and intrauterine growth in human genetic selection mechanisms.

Genome-wide association and polygenic risk score estimation of type 2 diabetes mellitus in Kinh Vietnamese-A pilot study.

A genome-wide association study (GWAS) is a powerful tool in investigating genetic contribution, which is a crucial factor in the development of complex multifactorial diseases, such as type 2 diabetes mellitus. Type 2 diabetes mellitus is a major healthcare burden in the Western Pacific region; however, there is limited availability of genetic-associated data for type 2 diabetes in Southeast Asia, especially among the Kinh Vietnamese population. This lack of information exacerbates global healthcare disparities. In this study, 997 Kinh Vietnamese individuals (503 with type 2 diabetes and 494 controls) were prospectively recruited and their clinical and paraclinical information was recorded. DNA samples were collected and whole genome genotyping was performed. Standard quality control and genetic imputation using the 1000 Genomes database were executed. A polygenic risk score for type 2 diabetes was generated in different models using East Asian, European, and mix ancestry GWAS summary statistics as training datasets. After quality control and genetic imputation, 107 polymorphisms reached suggestive statistical significance for GWAS (≤5 × 10-6) and rs11079784 was one of the potential markers strongly associated with type 2 diabetes in the studied population. The best polygenic risk score model predicting type 2 diabetes mellitus had AUC = 0.70 (95% confidence interval = 0.62-0.77) based on a mix of ancestral GWAS summary statistics. These data show promising results for genetic association with a polygenic risk score estimation in the Kinh Vietnamese population; the results also highlight the essential role of population diversity in a GWAS of type 2 diabetes mellitus.

Identification of RESP18 Gene Mutations Linked to Hereditary Non-Syndromic Cleft Lip and Palate in a Southern Chinese Family.

BACKGROUND Non-syndromic cleft lip with cleft palate (NSCLP) is one of the most common congenital birth defects worldwide; it causes lifelong problems and imposes burdens on patients and their families. This study aimed to describe the genomic analysis and identification of de novo regulated endocrine-specific protein 18 (RESP18) rs2385404 and rs2385405 gene polymorphisms associated with NSCLP in a southern Chinese family and to improve prevention, treatment, and prognosis of NSCLP. MATERIAL AND METHODS We performed a genome-wide association study (GWAS) to investigate the association of NSCLP phenotype with gene mutation. We investigated a 5-persons NSCLP family to screen the genetic variation of Han nationality in southern Chinese. Whole-genome sequencing (WGS) was used to detect all candidate genetic variants, and whole-exome sequencing (WES) was implemented to further verify mutations. The Clinical Variation Data Base (ClinVar) was employed for screening gene mutations. Finally, Sanger sequencing was applied to verify gene variations. RESULTS The combined analysis of WGS, WES, and ClinVar showed that a total of 9 variation positions overlapped among the 3 study cohorts. Sanger sequencing verified Glu amino acid variation in 2 mutation sites (rs2385404, rs2385405) from the RESP18 gene, which caused abnormal RESP18 function and was associated with hereditary NSCLP. CONCLUSIONS The combined genomic results showed that 2 mutations (rs2385404 and rs2385405) of the RESP18 gene were related to NSCLP in the family. The RESP18 gene may play an important role in the etiology and pathogenesis of cleft lip and palate.

Unveiling novel genetic variants in 370 challenging medically relevant genes using the long read sequencing data of 41 samples from 19 global populations.

BACKGROUND: A large number of challenging medically relevant genes (CMRGs) are situated in complex or highly repetitive regions of the human genome, hindering comprehensive characterization of genetic variants using next-generation sequencing technologies. In this study, we employed long-read sequencing technology, extensively utilized in studying complex genomic regions, to characterize genetic alterations, including short variants (single nucleotide variants and short insertions and deletions) and copy number variations, in 370 CMRGs across 41 individuals from 19 global populations. RESULTS: Our analysis revealed high levels of genetic variants in CMRGs, with 68.73% exhibiting copy number variations and 65.20% containing short variants that may disrupt protein function across individuals. Such variants can influence pharmacogenomics, genetic disease susceptibility, and other clinical outcomes. We observed significant differences in CMRG variation across populations, with individuals of African ancestry harboring the highest number of copy number variants and short variants compared to samples from other continents. Notably, 15.79% to 33.96% of short variants were exclusively detectable through long-read sequencing. While the T2T-CHM13 reference genome significantly improved the assembly of CMRG regions, thereby facilitating variant detection in these regions, some regions still lacked resolution. CONCLUSION: Our results provide an important reference for future clinical and pharmacogenetic studies, highlighting the need for a comprehensive representation of global genetic diversity in the reference genome and improved variant calling techniques to fully resolve medically relevant genes.