The GenomeAsia 100K Project enables genetic discoveries across Asia

Including diverse and admixed populations in genetic epidemiology research

The inclusion of ancestrally diverse participants in genetic studies can lead to new discoveries and is important to ensure equitable health care benefit from research advances. Here, members of the Ethical, Legal, Social, Implications (ELSI) committee of the International Genetic Epidemiology Society (IGES) offer perspectives on methods and analysis tools for the conduct of inclusive genetic epidemiology research, with a focus on admixed and ancestrally diverse populations in support of reproducible research practices. We emphasize the importance of distinguishing socially defined population categorizations from genetic ancestry in the design, analysis, reporting, and interpretation of genetic epidemiology research findings. Finally, we discuss the current state of genomic resources used in genetic association studies, functional interpretation, and clinical and public health translation of genomic findings with respect to diverse populations.

Updated Genome-Wide Association Study of Intracranial Aneurysms by Genotype Correction and Imputation in Koreans

OBJECTIVE

Compared to European, Japanese, and Chinese populations, genetic studies on intracranial aneurysms (IAs) in Koreans are lacking. We conducted an updated genome-wide association study (GWAS) to more accurately identify candidate variations predicting IA by genotype correction and imputation than in the first Korean GWAS.

METHODS

We performed a high-throughput imputation of single-nucleotide polymorphisms (SNPs) and genotype missing values for 250 IA and 296 controls. Out of a total of 7,333,746 sites with an imputation R² score of ≥0.5, 6,105,212 SNPs were analyzed. A high-throughput GWAS was performed after adjusting for clinical variables and 4 principal component analysis values.

RESULTS

A total of 39 SNPs reached a significant genome-wide threshold (P < 5 × 10^-8). After pruning by pairwise linkage disequilibrium (r² < 0.8), 11 SNPs were consistently associated with IA. Six tagging SNPs, including rs3120004, rs1851347, rs1522095, rs7779989, rs12935558, rs3826442, and rs2440154, showed strong linkage disequilibrium tower tagging haplotype structures. Among them, rs3120004 tagged a large and strong haplotype structure between LOC440704 and RGS18 genes in 1q31.2 (odds ratio, 2.34; 95% confidence interval, 1.74-3.14; P = 1.4 × 10^-8). The rs2440154 (SLC47A1, 17p11.2) SNP increased the risk of IA most significantly (odds ratio, 2.90; 95% confidence interval, 2.07-4.08; P = 8.2 × 10^-10). The region encompassing rs3826442 (MYH13, 17p13.1) showed a high recombination rate of approximately 70 cM/Mbp.

CONCLUSIONS

Our updated GWAS using high-throughput imputation approaches can be an informative milestone in understanding IA formation via susceptibility loci in this stage before large-scale genome-wide association meta-analysis.

Genome interpretation using in silico predictors of variant impact

Estimating the effects of variants found in disease driver genes opens the door to personalized therapeutic opportunities. Clinical associations and laboratory experiments can only characterize a tiny fraction of all the available variants, leaving the majority as variants of unknown significance (VUS). In silico methods bridge this gap by providing instant estimates on a large scale, most often based on the numerous genetic differences between species. Despite concerns that these methods may lack reliability in individual subjects, their numerous practical applications over cohorts suggest they are already helpful and have a role to play in genome interpretation when used at the proper scale and context. In this review, we aim to gain insights into the training and validation of these variant effect predicting methods and illustrate representative types of experimental and clinical applications. Objective performance assessments using various datasets that are not yet published indicate the strengths and limitations of each method. These show that cautious use of in silico variant impact predictors is essential for addressing genome interpretation challenges.

Use of race, ethnicity, and ancestry data in health research

Race, ethnicity, and ancestry are common classification variables used in health research. However, there has been no formal agreement on the definitions of these terms, resulting in misuse, confusion, and a lack of clarity surrounding these concepts for researchers and their readers. This article examines past and current understandings of race, ethnicity, and ancestry in research, identifies the distinctions between these terms, examines the reliability of these terms, and provides researchers with guidance on how to use these terms. Although race, ethnicity, and ancestry are often treated synonymously, they should be considered as distinct terms in the context of health research. Researchers should carefully consider which term is most appropriate for their study, define and use the terms consistently, and consider how their classification may be used in future research by others. The classification should be self-reported rather than assigned by an observer wherever possible.

Genetics of Chronic Kidney Disease in Low-Resource Settings

Advances in kidney genomics in the past 20 years has opened the door for more precise diagnosis of kidney disease and identification of new and specific therapeutic agents. Despite these advances, an imbalance exists between low-resource and affluent regions of the world. Individuals of European ancestry from the United States, United Kingdom, and Iceland account for 16% of the world's population, but represent more than 80% of all genome-wide association studies. South Asia, Southeast Asia, Latin America, and Africa together account for 57% of the world population but less than 5% of genome-wide association studies. Implications of this difference include limitations in new variant discovery, inaccurate interpretation of the effect of genetic variants in non-European populations, and unequal access to genomic testing and novel therapies in resource-poor regions. It also further introduces ethical, legal, and social pitfalls, and ultimately may propagate global health inequities. Ongoing efforts to reduce the imbalance in low-resource regions include funding and capacity building, population-based genome sequencing, population-based genome registries, and genetic research networks. More funding, training, and capacity building for infrastructure and expertise is needed in resource-poor regions. Focusing on this will ensure multiple-fold returns on investments in genomic research and technology.

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.

Autosomal recessive inheritance of a novel missense mutation of ITGB4 for Epidermolysis-Bullosa pyloric-atresia: a case report

Epidermolysis-Bullosa (EB), a rare Mendelian disorder, exhibits complex phenotypic and locus-heterogeneity. We identified a nuclear family of clinically unaffected parents with two offsprings manifesting EB-Pyloric-Atresia (EB-PA), with a variable clinical severity. We generated whole exome sequence data on all four individuals to (1) identify the causal mutation behind EB-PA (2) understand the background genetic variation for phenotype variability of the siblings. We assumed an autosomal recessive mode of inheritance and used suites of bioinformatic and computational tools to collate information through global databases to identify the causal genetic variant for the disease. We also investigated variations in key genes that are likely to impact phenotype severity. We identified a novel missense mutation in the ITGB4 gene (p.Ala1227Asp), for which the parents were heterozygous and the children homozygous. The mutation in ITGB4 gene, predicted to reduce the stability of the primary alpha6beta4-plectin complex compared to all previously studied mutations on ITGB4 reported to cause EB.

GAWMerge expands GWAS sample size and diversity by combining array-based genotyping and whole-genome sequencing

Genome-wide association studies (GWAS) have made impactful discoveries for complex diseases, often by amassing very large sample sizes. Yet, GWAS of many diseases remain underpowered, especially for non-European ancestries. One cost-effective approach to increase sample size is to combine existing cohorts, which may have limited sample size or be case-only, with public controls, but this approach is limited by the need for a large overlap in variants across genotyping arrays and the scarcity of non-European controls. We developed and validated a protocol, Genotyping Array-WGS Merge (GAWMerge), for combining genotypes from arrays and whole-genome sequencing, ensuring complete variant overlap, and allowing for diverse samples like Trans-Omics for Precision Medicine to be used. Our protocol involves phasing, imputation, and filtering. We illustrated its ability to control technology driven artifacts and type-I error, as well as recover known disease-associated signals across technologies, independent datasets, and ancestries in smoking-related cohorts. GAWMerge enables genetic studies to leverage existing cohorts to validly increase sample size and enhance discovery for understudied traits and ancestries.

Transferability of genetic loci and polygenic scores for cardiometabolic traits in British Pakistani and Bangladeshi individuals

Individuals with South Asian ancestry have a higher risk of heart disease than other groups but have been largely excluded from genetic research. Using data from 22,000 British Pakistani and Bangladeshi individuals with linked electronic health records from the Genes & Health cohort, we conducted genome-wide association studies of coronary artery disease and its key risk factors. Using power-adjusted transferability ratios, we found evidence for transferability for the majority of cardiometabolic loci powered to replicate. The performance of polygenic scores was high for lipids and blood pressure, but lower for BMI and coronary artery disease. Adding a polygenic score for coronary artery disease to clinical risk factors showed significant improvement in reclassification. In Mendelian randomisation using transferable loci as instruments, our findings were consistent with results in European-ancestry individuals. Taken together, trait-specific transferability of trait loci between populations is an important consideration with implications for risk prediction and causal inference.

1029 genomes of self-declared healthy individuals from India reveal prevalent and clinically relevant cardiac ion channelopathy variants

BACKGROUND

The prevalence and genetic spectrum of cardiac channelopathies exhibit population-specific differences. We aimed to understand the spectrum of cardiac channelopathy-associated variations in India, which is characterised by a genetically diverse population and is largely understudied in the context of these disorders.

RESULTS

We utilised the IndiGenomes dataset comprising 1029 whole genomes from self-declared healthy individuals as a template to filter variants in 36 genes known to cause cardiac channelopathies. Our analysis revealed 186,782 variants, of which we filtered 470 variants that were identified as possibly pathogenic (440 nonsynonymous, 30 high-confidence predicted loss of function ). About 26% (124 out of 470) of these variants were unique to the Indian population as they were not reported in the global population datasets and published literature. Classification of 470 variants by ACMG/AMP guidelines unveiled 13 pathogenic/likely pathogenic (P/LP) variants mapping to 19 out of the 1029 individuals. Further query of 53 probands in an independent cohort of cardiac channelopathy, using exome sequencing, revealed the presence of 3 out of the 13 P/LP variants. The identification of p.G179Sfs*62, p.R823W and c.420 + 2 T > C variants in KCNQ1, KCNH2 and CASQ2 genes, respectively, validate the significance of the P/LP variants in the context of clinical applicability as well as for large-scale population analysis.

CONCLUSION

A compendium of ACMG/AMP classified cardiac channelopathy variants in 1029 self-declared healthy Indian population was created. A conservative genotypic prevalence was estimated to be 0.9-1.8% which poses a huge public health burden for a country with large population size like India. In the majority of cases, these disorders are manageable and the risk of sudden cardiac death can be alleviated by appropriate lifestyle modifications as well as treatment regimens/clinical interventions. Clinical utility of the obtained variants was demonstrated using a cardiac channelopathy patient cohort. Our study emphasises the need for large-scale population screening to identify at-risk individuals and take preventive measures. However, we suggest cautious clinical interpretation to be exercised by taking other cardiac channelopathy risk factors into account.

CFTR genotype analysis of Asians in international registries highlights disparities in the diagnosis and treatment of Asian patients with cystic fibrosis

PURPOSE

Cystic fibrosis (CF) is not well-characterized in Asians, potentially resulting in delayed diagnosis and poor prognosis. We characterized CF in Asian subgroups to address these disparities.

METHODS

De-identified ethnicity and CFTR variant data were obtained from the United States, United Kingdom, and Canadian CF registries. We measured the prevalence of CF, CFTR variant allele frequencies, effectiveness of screening panels, and eligibility for modulator therapies.

RESULTS

The prevalence of CF was 1 in 74,982 people (Canada) to 1 in 13,340 people (United Kingdom) for South Asians and 1 in 256,541 (Canada) to 1 in 52,563 (United Kingdom) for other Asians, suggesting 26,000 to 146,000 patients with CF in South Asia. p.(F508del) variant was markedly less frequent in Asians than in non-Hispanic Whites. Splicing and nonsense variants occurred at high allelic frequencies in Asians, resulting in 41% to 49% of South Asians and 21% to 39% of other Asians being ineligible for CFTR modulator therapies. Hologic/EU2v1 panels failed to identify 37% to 47% of South Asian and 23% to 46% of other Asian patients with CF.

CONCLUSIONS

Among Asians, CF appears to be more common in South Asians. A significant CF population may exist in South Asia. CFTR variants in South and other Asians markedly differ from non-Hispanic Whites causing inequities in newborn screening, diagnosis, and treatment. New strategies are necessary to mitigate these health care disparities.

Advancing precision medicines for ocular disorders: Diagnostic genomics to tailored therapies

Successful sequencing of the human genome and evolving functional knowledge of gene products has taken genomic medicine to the forefront, soon combining broadly with traditional diagnostics, therapeutics, and prognostics in patients. Recent years have witnessed an extraordinary leap in our understanding of ocular diseases and their respective genetic underpinnings. As we are entering the age of genomic medicine, rapid advances in genome sequencing, gene delivery, genome surgery, and computational genomics enable an ever-increasing capacity to provide a precise and robust diagnosis of diseases and the development of targeted treatment strategies. Inherited retinal diseases are a major source of blindness around the world where a large number of causative genes have been identified, paving the way for personalized diagnostics in the clinic. Developments in functional genetics and gene transfer techniques has also led to the first FDA approval of gene therapy for LCA, a childhood blindness. Many such retinal diseases are the focus of various clinical trials, making clinical diagnoses of retinal diseases, their underlying genetics and the studies of natural history important. Here, we review methodologies for identifying new genes and variants associated with various ocular disorders and the complexities associated with them. Thereafter we discuss briefly, various retinal diseases and the application of genomic technologies in their diagnosis. We also discuss the strategies, challenges, and potential of gene therapy for the treatment of inherited and acquired retinal diseases. Additionally, we discuss the translational aspects of gene therapy, the important vector types and considerations for human trials that may help advance personalized therapeutics in ophthalmology. Retinal disease research has led the application of precision diagnostics and precision therapies; therefore, this review provides a general understanding of the current status of precision medicine in ophthalmology.

Novel Missense CNTNAP2 Variant Identified in Two Consanguineous Pakistani Families With Developmental Delay, Epilepsy, Intellectual Disability, and Aggressive Behavior

We report the genetic analysis of two consanguineous pedigrees of Pakistani ancestry in which two siblings in each family exhibited developmental delay, epilepsy, intellectual disability and aggressive behavior. Whole-genome sequencing was performed in Family 1, and we identified ~80,000 variants located in regions of homozygosity. Of these, 615 variants had a minor allele frequency ≤ 0.001, and 21 variants had CADD scores ≥ 15. Four homozygous exonic variants were identified in both affected siblings: PDZD7 (c.1348_1350delGAG, p.Glu450del), ALG6 (c.1033G&gt;C, p.Glu345Gln), RBM20 (c.1587C&gt;G, p.Ser529Arg), and CNTNAP2 (c.785G&gt;A, p.Gly228Arg). Sanger sequencing revealed co-segregation of the PDZD7, RBM20, and CNTNAP2 variants with disease in Family 1. Pathogenic variants in PDZD7 and RBM20 are associated with autosomal recessive non-syndromic hearing loss and autosomal dominant dilated cardiomyopathy, respectively, suggesting that these variants are unlikely likely to contribute to the clinical presentation. Gene panel analysis was performed on the two affected siblings in Family 2, and they were found to also be homozygous for the p.Gly228Arg CNTNAP2 variant. Together these families provide a LOD score 2.9 toward p.Gly228Arg CNTNAP2 being a completely penetrant recessive cause of this disease. The clinical presentation of the affected siblings in both families is also consistent with previous reports from individuals with homozygous CNTNAP2 variants where at least one allele was a nonsense variant, frameshift or small deletion. Our data suggests that homozygous CNTNAP2 missense variants can also contribute to disease, thereby expanding the genetic landscape of CNTNAP2 dysfunction.

Patterns and distribution of de novo mutations in multiplex Middle Eastern families

While de novo mutations (DNMs) are key to genetic diversity, they are also responsible for a high number of rare disorders. To date, no study has systematically examined the rate and distribution of DNMs in multiplex families in highly consanguineous populations. Leveraging WGS profiles of 645 individuals in 146 families, we implemented a combinatorial approach using 3 complementary tools for DNM discovery in 353 unique trio combinations. We found a total of 27,168 DNMs (median: 70 single-nucleotide and 6 insertion-deletions per individual). Phasing revealed around 80% of DNMs were paternal in origin. Notably, using whole-genome methylation data of spermatogonial stem cells, these DNMs were significantly more likely to occur at highly methylated CpGs (OR: 2.03; p value = 6.62 × 10⁻¹¹). We then examined the effects of consanguinity and ethnicity on DNMs, and found that consanguinity does not seem to correlate with DNM rate, and special attention has to be considered while measuring such a correlation. Additionally, we found that Middle-Eastern families with Arab ancestry had fewer DNMs than African families, although not significant (p value = 0.16). Finally, for families with diseased probands, we examined the difference in DNM counts and putative impact across affected and unaffected siblings, but did not find significant differences between disease groups, likely owing to the enrichment for recessive disorders in this part of the world, or the small sample size per clinical condition. This study serves as a reference for DNM discovery in multiplex families from the globally under-represented populations of the Middle-East.

A de novo start-loss in EFTUD2 associated with mandibulofacial dysostosis with microcephaly: case report

Mandibulofacial dysostosis with microcephaly (MFDM) is a rare genetic disorder inherited in an autosomal dominant pattern. Major characteristics include developmental delay, craniofacial malformations such as malar and mandibular hypoplasia, and ear anomalies. Here, we report a 4.5-yr-old female patient with symptoms fitting MFDM. Using whole-genome sequencing, we identified a de novo start-codon loss (c.3G > T) in the EFTUD2. We examined EFTUD2 expression in the patient by RNA sequencing and observed a notable functional consequence of the variant on gene expression in the patient. We identified a novel variant for the development of MFDM in humans. To the best of our knowledge, this is the first report of a start-codon loss in EFTUD2 associated with MFDM.

Cross-continental admixture in the Kho population from northwest Pakistan

Northern Pakistan is home to many diverse ethnicities and languages. The region acted as a prime corridor for ancient invasions and population migrations between Western Eurasia and South Asia. Kho, one of the major ethnic groups living in this region, resides in the remote and isolated mountainous region in the Chitral Valley of the Hindu Kush Mountain range. They are culturally and linguistically distinct from the rest of the Pakistani population groups and their genetic ancestry is still unknown. In this study, we generated genome-wide genotype data of ~1 M loci (Illumina WeGene array) for 116 unrelated Kho individuals and carried out comprehensive analyses in the context of worldwide extant and ancient anatomically modern human populations across Eurasia. The results inferred that the Kho can trace a large proportion of their ancestry to the population who migrated south from the Southern Siberian steppes during the second millennium BCE ~110 generations ago. An additional wave of gene flow from a population carrying East Asian ancestry was also identified in the Kho that occurred ~60 generations ago and may possibly be linked to the expansion of the Tibetan Empire during 7th to 9th centuries CE (current era) in the northwestern regions of the Indian sub-continent. We identified several candidate regions suggestive of positive selection in the Kho, that included genes mainly involved in pigmentation, immune responses, muscular development, DNA repair, and tumor suppression.

Population-based incidence and time to classification of systemic lupus erythematosus by three different classification criteria: a Lupus Midwest Network (LUMEN) study

OBJECTIVE

To estimate the incidence and time-to-classification of SLE by the 1997 ACR (ACR97) criteria, the SLICC criteria, and the European Alliance of Associations for Rheumatology/ACR (EULAR/ACR) criteria.

METHODS

We identified all incident SLE cases from 2000-2018 in the well-defined Olmsted County population. Clinical data included in the ACR97, SLICC and EULAR/ACR criteria were manually abstracted from medical records. All incident cases met at least one of the three classification criteria. Time-to-classification was estimated from the first documented lupus-attributable disease manifestation to the time of criteria fulfilment by each of the three definitions. Annual incidence rates were age or age/sex adjusted to the 2000 US population.

RESULTS

Of 139 incident cases there were 126 cases by the EULAR/ACR criteria, corresponding to an age/sex-adjusted incidence of 4.5 per 100 000 population (95% CI: 3.7, 5.2). The age/sex-incidence was higher than that of the SLICC criteria (113 cases; 4.0 per 100 000 [95% CI: 3.3, 4.7], P = 0.020) and the ACR97 (92 cases; 3.3 per 100 000 [95% CI: 2.6, 3.9], P < 0.001). The median time from first disease manifestation to criteria fulfilment was shorter for the EULAR/ACR criteria (29.4 months) than the ACR97 criteria (47.0 months, P < 0.001) and similar to the SLICC criteria (30.6 months, P = 0.83).

CONCLUSION

The incidence of SLE was higher by the EULAR/ACR criteria compared with the ACR97 and the SLICC criteria, and the EULAR/ACR criteria classified patients earlier that the ACR97 criteria but similar to the SLICC criteria.

Screening of the TMEM151A Gene in Patients With Paroxysmal Kinesigenic Dyskinesia and Other Movement Disorders

Background

Paroxysmal kinesigenic dyskinesia (PKD) is a rare neurological disorder characterized by recurrent involuntary movements usually triggered by sudden movements. Mutations in the TMEM151A gene were found to be the causative factor of PKD in recent studies. It has also been revealed that loss-of-function is the mechanism by which TMEM151A mutations cause PKD.

Methods

To investigate the genetic basis of PKD and broaden the clinical spectrum of the TMEM151A mutations, we recruited 181 patients of Chinese origin with movement disorders (MDs), including 39 PRRT2-negative PKD, 3 paroxysmal exercise-induced dyskinesia (PED), 2 paroxysmal non-kinesigenic dyskinesia (PNKD), 127 isolated dystonia, 8 choreas, and 2 myoclonus-dystonia syndromes. Whole-exome sequencing was applied to identify their possible disease-causing mutations. Then, Sanger sequencing was performed for validation and co-segregation analysis. Genetic analysis was also performed on additional family members of patients with TMEM151A mutations. Clinical manifestations of all PKD cases with mutations in TMEM151A reported, so far, were reviewed.

Results

Two novel variants of the TMEM151A gene (NM_153266.4, NP_694998.1), c.627_643dup (p.A215Gfs^*53) and c.627delG (p.L210Wfs^*52), were identified in 2 patients with PKD by whole-exome sequencing and further Sanger sequencing. Both variants were inherited by the patients from their respective mothers. No mutation of the TMEM151A gene was found in the other type of movement disorders. In reviewing the clinical presentation of TMEM151A-related PKD, no statistically significant difference in the age of onset, family history, duration of attacks, laterality, and phenotype was found between genders. More male patients received treatment and had a good response. A higher proportion of female patients did not receive any treatment, possibly because they had a milder condition of the disease.

Conclusions

This study further validated the role of TMEM151A in PKD. Future studies on protein function will be needed to ascertain the pathogenesis of TMEM151A in PKD.

Genomic analyses of 10,376 individuals in the Westlake BioBank for Chinese (WBBC) pilot project

We initiate the Westlake BioBank for Chinese (WBBC) pilot project with 4,535 whole-genome sequencing (WGS) individuals and 5,841 high-density genotyping individuals, and identify 81.5 million SNPs and INDELs, of which 38.5% are absent in dbSNP Build 151. We provide a population-specific reference panel and an online imputation server ( https://wbbc.westlake.edu.cn/ ) which could yield substantial improvement of imputation performance in Chinese population, especially for low-frequency and rare variants. By analyzing the singleton density of the WGS data, we find selection signatures in SNX29, DNAH1 and WDR1 genes, and the derived alleles of the alcohol metabolism genes (ADH1A and ADH1B) emerge around 7,000 years ago and tend to be more common from 4,000 years ago in East Asia. Genetic evidence supports the corresponding geographical boundaries of the Qinling-Huaihe Line and Nanling Mountains, which separate the Han Chinese into subgroups, and we reveal that North Han was more homogeneous than South Han.

Genetic Determinants of Survival in Parkinson's Disease in the Asian Population

BACKGROUND

Patients with Parkinson's disease (PD) have reduced life expectancy compared to the general population. Genetic variation was shown to play a role in the heterogeneity of survival for patients with PD, although the underlying genetic background remains poorly studied.

OBJECTIVE

The aim was to explore the genetic determinants influencing the survival of PD.

METHODS

We performed a genome-wide association analysis using a Cox proportional hazards model in a longitudinal cohort of 1080 Chinese patients with PD. Furthermore, we built a clinical-genetic model to predict the survival of patients using clinical variables combined with polygenic risk score (PRS) of survival of PD.

RESULTS

The cohort was followed up for an average of 7.13 years, with 85 incidents of death. One locus rs12628329 (RPL3) was significantly associated with reduced survival time by ~10.8 months (P = 2.72E-08, β = 1.79, standard error = 0.32). Functional exploration suggested this variant could upregulate the expression of RPL3 and induce apoptosis and cell death. In addition, adding PRS of survival in the prediction model substantially improved survival predictability (concordance index [Cindex]: 0.936) compared with the clinical model (Cindex: 0.860).

CONCLUSIONS

These findings improve the current understanding of the genetic cause of survival of PD and provide a novel target RPL3 for further research on PD pathogenesis and potential therapeutic options. Our results also demonstrate the potential utility of PRS of survival in identifying patients with shorter survival and providing personalized clinical monitoring and treatment. © 2022 International Parkinson and Movement Disorder Society.

TMC-SNPdb 2.0: an ethnic-specific database of Indian germline variants

Cancer is a somatic disease. The lack of Indian-specific reference germline variation resources limits the ability to identify true cancer-associated somatic variants among Indian cancer patients. We integrate two recent studies, the GenomeAsia 100K and the Genomics for Public Health in India (IndiGen) program, describing genome sequence variations across 598 and 1029 healthy individuals of Indian origin, respectively, along with the unique variants generated from our in-house 173 normal germline samples derived from cancer patients to generate the Tata Memorial Centre-SNP database (TMC-SNPdb) 2.0. To show its utility, GATK/Mutect2-based somatic variant calling was performed on 224 in-house tumor samples to demonstrate a reduction in false-positive somatic variants. In addition to the ethnic-specific variants from GenomeAsia 100K and IndiGenomes databases, 305 132 unique variants generated from 173 in-house normal germline samples derived from cancer patients of Indian origin constitute the Indian specific, TMC-SNPdb 2.0. Of 305 132 unique variants, 11.13% were found in the coding region with missense variants (31.3%) as the most predominant category. Among the non-coding variations, intronic variants (49%) were the highest contributors. The non-synonymous to synonymous SNP ratio was observed to be 1.9, consistent with the previous version of TMC-SNPdb and literature. Using TMC SNPdb 2.0, we analyzed a whole-exome sequence from 224 in-house tumor samples (180 paired and 44 orphans). We show an average depletion of 3.44% variants per paired tumor and significantly higher depletion (P-value < 0.001) for orphan tumors (4.21%), demonstrating the utility of the rare, unique variants found in the ethnic-specific variant datasets in reducing the false-positive somatic mutations. TMC-SNPdb 2.0 is the most exhaustive open-source reference database of germline variants occurring across 1800 Indian individuals to analyze cancer genomes and other genetic disorders. The database and toolkit package is available for download at the following: Database URL  http://www.actrec.gov.in/pi-webpages/AmitDutt/TMCSNPdb2/TMCSNPdb2.html.

North Asian population relationships in a global context

Population genetic studies of North Asian ethnic groups have focused on genetic variation of sex chromosomes and mitochondria. Studies of the extensive variation available from autosomal variation have appeared infrequently. We focus on relationships among population samples using new North Asia microhaplotype data. We combined genotypes from our laboratory on 58 microhaplotypes, distributed across 18 autosomes, on 3945 individuals from 75 populations with corresponding data extracted for 26 populations from the Thousand Genomes consortium and for 22 populations from the GenomeAsia 100 K project. A total of 7107 individuals in 122 total populations are analyzed using STRUCTURE, Principal Component Analysis, and phylogenetic tree analyses. North Asia populations sampled in Mongolia include: Buryats, Mongolians, Altai Kazakhs, and Tsaatans. Available Siberians include samples of Yakut, Khanty, and Komi Zyriane. Analyses of all 122 populations confirm many known relationships and show that most populations from North Asia form a cluster distinct from all other groups. Refinement of analyses on smaller subsets of populations reinforces the distinctiveness of North Asia and shows that the North Asia cluster identifies a region that is ancestral to Native Americans.

Genetic insight into Birt-Hogg-Dubé syndrome in Indian patients reveals novel mutations at FLCN

Background

Birt-Hogg-Dubé syndrome (BHDS) is a rare monogenic condition mostly associated with germline mutations at FLCN. It is characterized by either one or more manifestations of primary spontaneous pneumothorax (PSP), skin fibrofolliculomas and renal carcinoma (chromophobe). Here, we comprehensively studied the mutational background of 31 clinically diagnosed BHDS patients and their 74 asymptomatic related members from 15 Indian families.

Results

Targeted amplicon next-generation sequencing (NGS) and Sanger sequencing of FLCN in patients and asymptomatic members revealed a total of 76 variants. Among these variants, six different types of pathogenic FLCN mutations were detected in 26 patients and some asymptomatic family members. Two of the variants were novel mutations: an 11-nucleotide deletion (c.1150_1160delGTCCAGTCAGC) and a splice acceptor mutation (c.1301-1G > A). Two variants were Clinvar reported pathogenic mutations: a stop-gain (c.634C > T) and a 4-nucleotide duplication (c.1329_1332dupAGCC). Two known variants were: hotspot deletion (c.1285delC) and a splice donor mutation (c.1300 + 1G > A). FLCN mutations could not be detected in patients and asymptomatic members from 5 families. All these mutations greatly affected the protein stability and FLCN-FNIP2 interaction as observed by molecular docking method. Family-based association study inferred pathogenic FLCN mutations are significantly associated with BHDS.

Conclusion

Six pathogenic FLCN mutations were detected in patients from 10 families out of 15 families in the cohort. Therefore, genetic screening is necessary to validate the clinical diagnosis. The pathogenic mutations at FLCN affects the protein–protein interaction, which plays key roles in various metabolic pathways. Since, pathogenic mutations could not be detected in exonic regions of FLCN in 5 families, whole genome sequencing is necessary to detect all mutations at FLCN and/or any undescribed gene/s that may also be implicated in BHDS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-022-02326-5.

Integrative analysis of genomic and transcriptomic data of normal, tumour and co-occurring leukoplakia tissue triads drawn from patients with gingivobuccal oral cancer identifies signatures of tumour initiation and progression

A thickened, white patch — leukoplakia — in the oral cavity is usually benign, but sometimes (in ~9% of individuals) it progresses to malignant tumour. Because the genomic basis of this progression is poorly understood, we undertook this study and collected samples of four tissues — leukoplakia, tumour, adjacent normal, and blood — from each of 28 patients suffering from gingivobuccal oral cancer. We performed multiomics analysis of the 112 collected tissues (four tissues per patient from 28 patients) and integrated information on progressive changes in the mutational and transcriptional profiles of each patient to create this genomic narrative. Additionally, we generated and analysed whole‐exome sequence data from leukoplakia tissues collected from 11 individuals not suffering from oral cancer. Nonsynonymous somatic mutations in the CASP8 gene were identified as the likely events to initiate malignant transformation, since these were frequently shared between tumour and co‐occurring leukoplakia. CASP8 alterations were also shown to enhance expressions of genes that favour lateral spread of mutant cells. During malignant transformation, additional pathogenic mutations are acquired in key genes (TP53, NOTCH1, HRAS) (41% of patients); chromosomal‐instability (arm‐level deletions of 19p and q, focal‐deletion of DNA‐repair pathway genes and NOTCH1, amplification of EGFR) (77%), and increased APOBEC‐activity (23%) are also observed. These additional alterations were present singly (18% of patients) or in combination (68%). Some of these alterations likely impact immune‐dynamics of the evolving transformed tissue; progression to malignancy is associated with immune suppression through infiltration of regulatory T‐cells (56%), depletion of cytotoxic T‐cells (68%), and antigen‐presenting dendritic cells (72%), with a concomitant increase in inflammation (92%). Patients can be grouped into three clusters by the estimated time to development of cancer from precancer by acquiring additional mutations (range: 4–10 years). Our findings provide deep molecular insights into the evolutionary processes and trajectories of oral cancer initiation and progression. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Variant landscape of the RYR1 gene based on whole genome sequencing of the Singaporean population

The RYR1 gene codes for a ryanodine receptor which is a calcium release channel in the skeletal muscle sarcoplasmic reticulum. It is associated with Malignant Hyperthermia (MH) and congenital myopathies including Central Core Disease (CCD), Multiminicore Disease (MMD) and Congenital Fibre-Type Disproportion (CFTD). There is currently little information on the epidemiology of RYR1 variants in Asians. Our study aims to describe the RYR1 variant landscape in a Singapore cohort unselected for RYR1-associated conditions. Data was retrieved from the SG10K pilot project, where whole genome sequencing was performed on volunteers unselected and undetermined for RYR1-associated conditions. Variants were classified based on pathogenicity using databases ClinVar and InterVar. Allele frequencies of pathogenic variants were compared between Chinese, Indians and Malays. Using databases ExAC, GnomAD and GenomeAsia 100k study, we further compared local allele frequencies to those in Europe, America and Asia. Data was analysed using R Commander. Significant P value was set at p < 0.05. Majority of the RYR1 variants were missense mutations. We identified four pathogenic and four likely pathogenic RYR1 variants. All were related to the aforementioned RYR1-associated conditions. There were 6 carriers of RYR1 pathogenic variants amongst 4810 individuals, corresponding to an allele frequency of 0.06%. The prevalence of pathogenic variants was the highest amongst Indians (4 in 1127 individuals) (p = 0.030). Majority of pathogenic and likely pathogenic mutations were missense and located in mutational hotspots. These variants also occurred at higher frequencies in Asians than globally. This study describes the variant landscape of the RYR1 gene in Singapore. This knowledge will facilitate genetic screening for RYR1-related conditions.

Understanding signatures of positive natural selection in human zinc transporter genes

Zinc is an essential micronutrient with a tightly regulated systemic and cellular homeostasis. In humans, some zinc transporter genes (ZTGs) have been previously reported as candidates for strong geographically restricted selective sweeps. However, since zinc homeostasis is maintained by the joint action of 24 ZTGs, other more subtle modes of selection could have also facilitated human adaptation to zinc availability. Here, we studied whether the complete set of ZTGs are enriched for signals of positive selection in worldwide populations and population groups from South Asia. ZTGs showed higher levels of genetic differentiation between African and non-African populations than would be randomly expected, as well as other signals of polygenic selection outside Africa. Moreover, in several South Asian population groups, ZTGs were significantly enriched for SNPs with unusually extended haplotypes and displayed SNP genotype-environmental correlations when considering zinc deficiency levels in soil in that geographical area. Our study replicated some well-characterized targets for positive selection in East Asia and sub-Saharan Africa, and proposes new candidates for follow-up in South Asia (SLC39A5) and Africa (SLC39A7). Finally, we identified candidate variants for adaptation in ZTGs that could contribute to different disease susceptibilities and zinc-related human health traits.

A homozygous hypomorphic BNIP1 variant causes an increase in autophagosomes and reduced autophagic flux and results in a spondylo-epiphyseal dysplasia

BNIP1 (BCL2 interacting protein 1) is a soluble N-ethylmaleimide-sensitive factor-attachment protein receptor involved in ER membrane fusion. We identified the homozygous BNIP1 intronic variant c.84+3A>T in the apparently unrelated patients 1 and 2 with disproportionate short stature. Radiographs showed abnormalities affecting both the axial and appendicular skeleton and spondylo-epiphyseal dysplasia. We detected ~80% aberrantly spliced BNIP1 pre-mRNAs, reduced BNIP1 mRNA level to ~80%, and BNIP1 protein level reduction by ~50% in patient 1 compared to control fibroblasts. The BNIP1 ortholog in drosophila, Sec. 20, regulates autophagy and lysosomal degradation. We assessed lysosome positioning and identified a decrease in lysosomes in the perinuclear region and an increase in the cell periphery in patient 1 cells. Immunofluorescence microscopy and immunoblotting demonstrated an increase in LC3B-positive structures and LC3B-II levels, respectively, in patient 1 fibroblasts under steady-state condition. Treatment of serum-starved fibroblasts with or without bafilomycin A1 identified significantly decreased autophagic flux in patient 1 cells. Our data suggest a block at the terminal stage of autolysosome formation and/or clearance in patient fibroblasts. BNIP1 together with RAB33B and VPS16, disease genes for Smith-McCort dysplasia 2 and a multisystem disorder with short stature, respectively, highlight the importance of autophagy in skeletal development. This article is protected by copyright. All rights reserved.

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.

HKG: an open genetic variant database of 205 Hong Kong cantonese exomes

HKG is the first fully accessible variant database for Hong Kong Cantonese, constructed from 205 novel whole-exome sequencing data. There has long been a research gap in the understanding of the genetic architecture of southern Chinese subgroups, including Hong Kong Cantonese. HKG detected 196 325 high-quality variants with 5.93% being novel, and 25 472 variants were found to be unique in HKG compared to three Chinese populations sampled from 1000 Genomes (CHN). PCA illustrates the uniqueness of HKG in CHN, and the admixture study estimated the ancestral composition of HKG and CHN, with a gradient change from north to south, consistent with their geological distribution. ClinVar, CIViC and PharmGKB annotated 599 clinically significant variants and 360 putative loss-of-function variants, substantiating our understanding of population characteristics for future medical development. Among the novel variants, 96.57% were singleton and 6.85% were of high impact. With a good representation of Hong Kong Cantonese, we demonstrated better variant imputation using reference with the addition of HKG data, thus successfully filling the data gap in southern Chinese to facilitate the regional and global development of population genetics.

The current understanding and future directions for sputum microbiome profiling in chronic obstructive pulmonary disease

PURPOSE OF REVIEW

Next-generation sequencing (NGS) has deepened our understanding of the respiratory microbiome in health and disease. The number of microbiome studies employing sputum as an airway surrogate has continued to increase over the past decade to include multiple large multicentre and longitudinal studies of the microbiome in chronic obstructive pulmonary disease (COPD). In this review, we summarize the recent advances to our understanding of the bacteriome, virome and mycobiome in COPD.

RECENT FINDINGS

Diverse microbiome profiles are reported in COPD. The neutrophilic Haemophilus-predominant bacteriome remains a prominent COPD phenotype, relatively stable over time and during exacerbations. Studies of the virome remain limited but reveal a potential involvement of viruses and bacteriophages particularly during COPD exacerbations and advancing disease severity. Mycobiome signatures, even in stable COPD are associated with poorer clinical outcomes including mortality.

SUMMARY

The sputum microbiome in COPD is being increasingly recognized for its clinical relevance, even in the stable state. Future studies integrating microbial kingdoms holistically (i.e. bacterial, viral and fungal) will provide deeper insight into its functionality including the relevance of microbial interactions and effect of treatment on microbiome-associated clinical outcomes.

A guide for the diagnosis of rare and undiagnosed disease: beyond the exome

Rare diseases affect 30 million people in the USA and more than 300-400 million worldwide, often causing chronic illness, disability, and premature death. Traditional diagnostic techniques rely heavily on heuristic approaches, coupling clinical experience from prior rare disease presentations with the medical literature. A large number of rare disease patients remain undiagnosed for years and many even die without an accurate diagnosis. In recent years, gene panels, microarrays, and exome sequencing have helped to identify the molecular cause of such rare and undiagnosed diseases. These technologies have allowed diagnoses for a sizable proportion (25-35%) of undiagnosed patients, often with actionable findings. However, a large proportion of these patients remain undiagnosed. In this review, we focus on technologies that can be adopted if exome sequencing is unrevealing. We discuss the benefits of sequencing the whole genome and the additional benefit that may be offered by long-read technology, pan-genome reference, transcriptomics, metabolomics, proteomics, and methyl profiling. We highlight computational methods to help identify regionally distant patients with similar phenotypes or similar genetic mutations. Finally, we describe approaches to automate and accelerate genomic analysis. The strategies discussed here are intended to serve as a guide for clinicians and researchers in the next steps when encountering patients with non-diagnostic exomes.

Leveraging deep phenotyping from health check-up cohort with 10,000 Korean individuals for phenome-wide association study of 136 traits

The expanding use of the phenome-wide association study (PheWAS) faces challenges in the context of using International Classification of Diseases billing codes for phenotype definition, imbalanced study population ethnicity, and constrained application of the results in research. We performed a PheWAS utilizing 136 deep phenotypes corroborated by comprehensive health check-ups in a Korean population, along with trans-ethnic comparisons through using the UK Biobank and Biobank Japan Project. Meta-analysis with Korean and Japanese population was done. The PheWAS associated 65 phenotypes with 14,101 significant variants (P < 4.92 × 10-10). Network analysis, visualization of cross-phenotype mapping, and causal inference mapping with Mendelian randomization were conducted. Among phenotype pairs from the genotype-driven cross-phenotype associations, we evaluated penetrance in correlation analysis using a clinical database. We focused on the application of PheWAS in order to make it robust and to aid the derivation of biological meaning post-PheWAS. This comprehensive analysis of PheWAS results based on a health check-up database will provide researchers and clinicians with a panoramic overview of the networks among multiple phenotypes and genetic variants, laying groundwork for the practical application of precision medicine.

Genetic Connections and Convergent Evolution of Tropical Indigenous Peoples in Asia

Tropical indigenous peoples in Asia (TIA) attract much attention for their unique appearance, whereas their genetic history and adaptive evolution remain mysteries. We conducted a comprehensive study to characterize the genetic distinction and connection of broad geographical TIAs. Despite the diverse genetic makeup and large interarea genetic differentiation between the TIA groups, we identified a basal Asian ancestry (bASN) specifically shared by these populations. The bASN ancestry was relatively enriched in ancient Asian human genomes dated as early as ∼50,000 years before the present and diminished in more recent history. Notably, the bASN ancestry is unlikely to be derived from archaic hominins. Instead, we suggest it may be better modeled as a survived lineage of the initial peopling of Asia. Shared adaptations inherited from the ancient Asian ancestry were detected among the TIA groups (e.g., LIMS1 for hair morphology, and COL24A1 for bone formation), and they are enriched in neurological functions either at an identical locus (e.g., NKAIN3), or different loci in an identical gene (e.g., TENM4). The bASN ancestry could also have formed the substrate of the genetic architecture of the dark pigmentation observed in the TIA peoples. We hypothesize that phenotypic convergence of the dark pigmentation in TIAs could have resulted from parallel (e.g., DDB1/DAK) or genetic convergence driven by admixture (e.g., MTHFD1 and RAD18), new mutations (e.g., STK11), or notably purifying selection (e.g., MC1R). Our results provide new insights into the initial peopling of Asia and an advanced understanding of the phenotypic convergence of the TIA peoples.

Etiological roles of core promoter variation in triple-negative breast cancer

Abnormal gene expression plays key role in cancer development. A core promoter is located around the transcriptional start site. Through interaction between core promoter sequences and transcriptional factors, core promoter controls transcriptional initiation. We hypothesized that in cancer, core promoter sequences could be mutated to interfere the interaction with transcriptional factors, resulting in altered transcriptional initiation and abnormal gene expression and cancer development. We used triple-negative breast cancer (TNBC) as a model to test our hypothesis. We collected genome-wide core promoter variants from 279 TNBC genomes. After extensive filtering of normal genomic polymorphism, we identified 19,427 recurrent somatic variants in 1,238 core promoters of 1,274 genes and 1,694 recurrent germline variants in 272 core promoters of 294 genes. Many of the affected genes were oncogenes and tumor suppressors. Analysis of RNA-seq data from the same patient cohort identified increased or decreased gene expression in 439 somatic and 85 germline variants-affected genes, and the results were validated by luciferase reporter assay. By comparing with the core promoter variation data from 610 unclassified breast cancer, we observed that core promoter variants in TNBC were highly TNBC-specific. We further identified the drugs targeting the genes with core promoter variation. Our study demonstrates that core promoter is highly mutable in cancer, and can play etiological roles in TNBC and other types of cancer through influencing transcriptional initiation.

Who's afraid of genetic tests?: An assessment of Singapore's public attitudes and changes in attitudes after taking a genetic test

BACKGROUND

As a consequence of precision medicine initiatives, genomic technologies have rapidly spread around the world, raising questions about genetic privacy and the ethics of data sharing. Previous scholarship in bioethics and science and technology studies has made clear that different nations have varying expectations about trust, transparency, and public reason in relation to emerging technologies and their governance. The key aims of this article are to assess genetic literacy, perceptions of genetic testing, privacy concerns, and governing norms amongst the Singapore population by collecting surveys.

METHODS

This study investigated genetic literacy and broad public attitudes toward genetic tests in Singapore with an online public survey (n = 560). To assess potential changes in attitudes following receipt of results from a genetic test, we also surveyed undergraduate students who underwent a genetic screen as part of a university class before and after they received their test results (n = 25).

RESULTS

Public participants showed broad support for the use of genetic tests; scored an average of 48.9% in genetic literacy; and expressed privacy concerns over data sharing and a desire for control over their genetic data. After taking a genetic test and receiving genetic test results, students reported less fear of genetic tests while other attitudes did not change significantly.

CONCLUSION

These findings highlight the potential of genetic education and active engagement with genetic testing to increase support and participation in genomic projects, PM, and biobanking initiatives; and they suggest that data privacy protections could potentially reduce discrimination by giving participants control over who can access their data. More specifically, these findings and the dataset we provide may be helpful in formulating culturally sensitive education programs and regulations concerning genomic technologies and data privacy.

Imputation Performance in Latin American Populations: Improving Rare Variants Representation With the Inclusion of Native American Genomes

Current Genome-Wide Association Studies (GWAS) rely on genotype imputation to increase statistical power, improve fine-mapping of association signals, and facilitate meta-analyses. Due to the complex demographic history of Latin America and the lack of balanced representation of Native American genomes in current imputation panels, the discovery of locally relevant disease variants is likely to be missed, limiting the scope and impact of biomedical research in these populations. Therefore, the necessity of better diversity representation in genomic databases is a scientific imperative. Here, we expand the 1,000 Genomes reference panel (1KGP) with 134 Native American genomes (1KGP + NAT) to assess imputation performance in Latin American individuals of mixed ancestry. Our panel increased the number of SNPs above the GWAS quality threshold, thus improving statistical power for association studies in the region. It also increased imputation accuracy, particularly in low-frequency variants segregating in Native American ancestry tracts. The improvement is subtle but consistent across countries and proportional to the number of genomes added from local source populations. To project the potential improvement with a higher number of reference genomes, we performed simulations and found that at least 3,000 Native American genomes are needed to equal the imputation performance of variants in European ancestry tracts. This reflects the concerning imbalance of diversity in current references and highlights the contribution of our work to reducing it while complementing efforts to improve global equity in genomic research.

Haplotype-based inference of the distribution of fitness effects

Recent genome sequencing studies with large sample sizes in humans have discovered a vast quantity of low-frequency variants, providing an important source of information to analyze how selection is acting on human genetic variation. In order to estimate the strength of natural selection acting on low-frequency variants, we have developed a likelihood-based method that uses the lengths of pairwise identity-by-state between haplotypes carrying low-frequency variants. We show that in some non-equilibrium populations (such as those that have had recent population expansions) it is possible to distinguish between positive or negative selection acting on a set of variants. With our new framework, one can infer a fixed selection intensity acting on a set of variants at a particular frequency, or a distribution of selection coefficients for standing variants and new mutations. We show an application of our method to the UK10K phased haplotype dataset of individuals.

Analyzing the Korean reference genome with meta-imputation increased the imputation accuracy and spectrum of rare variants in the Korean population

Genotype imputation is essential for enhancing the power of association-mapping and discovering rare and indels that are missed by most genotyping arrays. Imputation analysis can be more accurate with a population-specific reference panel or a multi-ethnic reference panel with numerous samples. The National Institute of Health, Republic of Korea, initiated the Korean Reference Genome (KRG) project to identify variants in whole-genome sequences of ∼20,000 Korean participants. In the pilot phase, we analyzed the data from 1,490 participants. The genetic characteristics and imputation performance of the KRG were compared with those of the 1,000 Genomes Project Phase 3, GenomeAsia 100K Project, ChinaMAP, NARD, and TOPMed reference panels. For comparison analysis, genotype panels were artificially generated using whole-genome sequencing data from combinations of four different ancestries (Korean, Japanese, Chinese, and European) and two population-specific optimized microarrays (Korea Biobank Array and UK Biobank Array). The KRG reference panel performed best for the Korean population (R ² = 0.78-0.84, percentage of well-imputed is 91.9% for allele frequency >5%), although the other reference panels comprised a larger number of samples with genetically different background. By comparing multiple reference panels and multi-ethnic genotype panels, optimal imputation was obtained using reference panels from genetically related populations and a population-optimized microarray. Indeed, the reference panels of KRG and TOPMed showed the best performance when applied to the genotype panels of KBA (R ² = 0.84) and UKB (R ² = 0.87), respectively. Using a meta-imputation approach to merge imputation results from different reference panels increased the imputation accuracy for rare variants (∼7%) and provided additional well-imputed variants (∼20%) with comparable imputation accuracy to that of the KRG. Our results demonstrate the importance of using a population-specific reference panel and meta-imputation to assess a substantial number of accurately imputed rare variants.

COVID-19: Impact on linguistic and genetic isolates of India

The rapid expansion of coronavirus SARS-CoV-2 has impacted various ethnic groups all over the world. The burden of infectious diseases including COVID-19 are generally reported to be higher for the Indigenous people. The historical knowledge have also suggested that the indigenous populations suffer more than the general populations in the pandemic. Recently, it has been reported that the indigenous groups of Brazil have been massively affected by COVID-19. Series of studies have shown that many of the indigenous communities reached at the verge of extinction due to this pandemic. Importantly, South Asia also has several indigenous and smaller communities, that are living in isolation. Till date, despite the two consecutive waves in India, there is no report on the impact of COVID-19 for indigenous tribes. Since smaller populations experiencing drift may have greater risk of such pandemic, we have analysed Runs of Homozygosity (ROH) among South Asian populations and identified several populations with longer homozygous segments. The longer runs of homozygosity at certain genomic regions may increases the susceptibility for COVID-19. Thus, we suggest extreme careful management of this pandemic among isolated populations of South Asia.

BGvar: A comprehensive resource for blood group immunogenetics

BACKGROUND

Blood groups form the basis of effective and safe blood transfusion. There are about 43 well-recognised human blood group systems presently known. Blood groups are molecularly determined by the presence of specific antigens on the red blood cells and are genetically determined and inherited following Mendelian principles. The lack of a comprehensive, relevant, manually compiled and genome-ready dataset of red cell antigens limited the widespread application of genomic technologies to characterise and interpret the blood group complement of an individual from genomic datasets.

MATERIALS AND METHODS

A range of public datasets was used to systematically annotate the variation compendium for its functionality and allele frequencies across global populations. Details on phenotype or relevant clinical importance were collated from reported literature evidence.

RESULTS

We have compiled the Blood Group Associated Genomic Variant Resource (BGvar), a manually curated online resource comprising all known human blood group related allelic variants including a total of 1700 International Society of Blood Transfusion approved alleles and 1706 alleles predicted and curated from literature reports. This repository includes 1682 single nucleotide variations (SNVs), 310 Insertions, Deletions (InDels) and Duplications (Copy Number Variations) and about 1360 combination mutations corresponding to 43 human blood group systems and 2 transcription factors. This compendium also encompasses gene fusion and rearrangement events occurring in human blood group genes.

CONCLUSION

To the best of our knowledge, BGvar is a comprehensive and a user-friendly resource with most relevant collation of blood group alleles in humans. BGvar is accessible online at URL: http://clingen.igib.res.in/bgvar/.

Loss of Y chromosome: An emerging next-generation biomarker for disease prediction and early detection?

As human life expectancy increases substantially and aging is the primary risk factor for most chronic diseases, there is an urgent need for advancing the development of post-genomic era biomarkers that can be used for disease prediction and early detection (DPED). Mosaic loss of Y chromosome (LOY) is the state of nullisomy Y in sub-groups of somatic cells acquired from different post-zygotic development stages and onwards throughout the lifespan. Multiple large-cohort based epidemiology studies have found that LOY in blood cells is a significant risk factor for future mortality and various diseases in males. Many features intrinsic to LOY analysis may be leveraged to enhance its use as a non-invasive, sensitive, reliable, high throughput-biomarker for DPED. Here, we review the emerging literatures in LOY studies and highlight ten strengths for using LOY as a novel biomarker for genomics-driven DPED diagnostics. Meanwhile, the current limitations in this area are also discussed. We conclude by identifying some important knowledge gaps regarding the consequences of malsegregation of the Y chromosome and propose further steps that are required before clinical implementation of LOY. Taken together, we think that LOY has substantial potential as a biomarker for DPED, despite some hurdles that still need to be addressed before its integration into healthcare becomes acceptable.

Fine-scale population structure and demographic history of British Pakistanis

Previous genetic and public health research in the Pakistani population has focused on the role of consanguinity in increasing recessive disease risk, but little is known about its recent population history or the effects of endogamy. Here, we investigate fine-scale population structure, history and consanguinity patterns using genotype chip data from 2,200 British Pakistanis. We reveal strong recent population structure driven by the biraderi social stratification system. We find that all subgroups have had low recent effective population sizes (Ne), with some showing a decrease 15‒20 generations ago that has resulted in extensive identity-by-descent sharing and homozygosity, increasing the risk of recessive disorders. Our results from two orthogonal methods (one using machine learning and the other coalescent-based) suggest that the detailed reporting of parental relatedness for mothers in the cohort under-represents the true levels of consanguinity. These results demonstrate the impact of cultural practices on population structure and genomic diversity in Pakistanis, and have important implications for medical genetic studies.

[CASE REPORT] Homozygous N-terminal missense variant in PLEKHG5 associated with intermediate CMT: a case report

Mutations in PLEKHG5, a pleckstrin homology domain containing member of the GEF family, are associated with distal spinal muscular atrophy and intermediate Charcot-Marie-Tooth disease. Here, we describe an isolated case with distal intermediate neuropathy with scapular winging. By whole exome sequencing, we identified the homozygous PLEKHG5 Arg97Gln missense mutation, located in the N-terminal region of the protein. This mutation resides between a zinc-finger motif and a RBD domain, involved in binding rnd3, a RhoA effector protein. We conclude that based on the characteristic phenotype presented by the patient and the supportive genetic findings, the PLEKHG5 mutation is the causative variant.

The perils of generalization: Rethinking breast cancer screening guidelines for young women of color

Women-of-Color are more likely to be diagnosed with breast cancer under the age of 40 and more likely to die (vs. non-Hispanic Whites). These finding support the concern that the United States Preventive Services Task Force breast cancer adversely impacts Women-of-Color.

The global carrier frequency and genetic prevalence of Upshaw-Schulman syndrome

Background

Upshaw–Schulman syndrome (USS) is an autosomal recessive disease characterized by thrombotic microangiopathies caused by pathogenic variants in ADAMTS13. We aimed to (1) curate the ADAMTS13 gene pathogenic variant dataset and (2) estimate the carrier frequency and genetic prevalence of USS using Genome Aggregation Database (gnomAD) data.

Methods

Studies were comprehensively retrieved. All previously reported pathogenic ADAMTS13 variants were compiled and annotated with gnomAD allele frequencies. The pooled global and population-specific carrier frequencies and genetic prevalence of USS were calculated using the Hardy-Weinberg equation.

Results

We mined reported disease-causing variants that were present in the gnomAD v2.1.1, filtered by allele frequency. The pathogenicity of variants was classified according to the American College of Medical Genetics and Genomics criteria. The genetic prevalence and carrier frequency of USS were 0.43 per 1 million (95% CI: [0.36, 0.55]) and 1.31 per 1 thousand population, respectively. When the novel pathogenic/likely pathogenic variants were included, the genetic prevalence and carrier frequency were 1.1 per 1 million (95% CI: [0.89, 1.37]) and 2.1 per 1 thousand population, respectively.

Conclusions

The genetic prevalence and carrier frequency of USS were within the ranges of previous estimates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-01010-0.

NyuWa Genome resource: A deep whole-genome sequencing-based variation profile and reference panel for the Chinese population

The lack of haplotype reference panels and whole-genome sequencing resources specific to the Chinese population has greatly hindered genetic studies in the world's largest population. Here, we present the NyuWa genome resource, based on deep (26.2×) sequencing of 2,999 Chinese individuals, and construct a NyuWa reference panel of 5,804 haplotypes and 19.3 million variants, which is a high-quality publicly available Chinese population-specific reference panel with thousands of samples. Compared with other panels, the NyuWa reference panel reduces the Han Chinese imputation error rate by a margin ranging from 30% to 51%. Population structure and imputation simulation tests support the applicability of one integrated reference panel for northern and southern Chinese. In addition, a total of 22,504 loss-of-function variants in coding and noncoding genes are identified, including 11,493 novel variants. These results highlight the value of the NyuWa genome resource in facilitating genetic research in Chinese and Asian populations.

MicroRNA Variants and HLA-miRNA Interactions are Novel Rheumatoid Arthritis Susceptibility Factors

Genome-wide association studies have identified >100 genetic risk factors for rheumatoid arthritis. However, the reported genetic variants could only explain less than 40% heritability of rheumatoid arthritis. The majority of the heritability is still missing and needs to be identified with more studies with different approaches and populations. In order to identify novel function SNPs to explain missing heritability and reveal novel mechanism pathogenesis of rheumatoid arthritis, 4 HLA SNPs (HLA-DRB1, HLA-DRB9, HLA-DQB1, and TNFAIP3) and 225 common SNPs located in miRNA, which might influence the miRNA target binding or pre-miRNA stability, were genotyped in 1,607 rheumatoid arthritis and 1,580 matched normal individuals. We identified 2 novel SNPs as significantly associated with rheumatoid arthritis including rs1414273 (miR-548ac, OR = 0.84, p = 8.26 × 10^-4) and rs2620381 (miR-627, OR = 0.77, p = 2.55 × 10^-3). We also identified that rs5997893 (miR-3928) showed significant epistasis effect with rs4947332 (HLA-DRB1, OR = 4.23, p = 0.04) and rs2967897 (miR-5695) with rs7752903 (TNFAIP3, OR = 4.43, p = 0.03). In addition, we found that individuals who carried 8 risk alleles showed 15.38 (95%CI: 4.69-50.49, p < 1.0 × 10^-6) times more risk of being affected by RA. Finally, we demonstrated that the targets of the significant miRNAs showed enrichment in immune related genes (p = 2.0 × 10^-5) and FDA approved drug target genes (p = 0.014). Overall, 6 novel miRNA SNPs including rs1414273 (miR-548ac, p = 8.26 × 10^-4), rs2620381 (miR-627, p = 2.55 × 10^-3), rs4285314 (miR-3135b, p = 1.10 × 10^-13), rs28477407 (miR-4308, p = 3.44 × 10^-5), rs5997893 (miR-3928, p = 5.9 × 10^-3) and rs45596840 (miR-4482, p = 6.6 × 10^-3) were confirmed to be significantly associated with RA in a Chinese population. Our study suggests that miRNAs might be interesting targets to accelerate understanding of the pathogenesis and drug development for rheumatoid arthritis.