Large-scale cross-ancestry genome-wide meta-analysis of serum urate

Hyperuricemia is an essential causal risk factor for gout and is associated with cardiometabolic diseases. Given the limited contribution of East Asian ancestry to genome-wide association studies of serum urate, the genetic architecture of serum urate requires exploration. A large-scale cross-ancestry genome-wide association meta-analysis of 1,029,323 individuals and ancestry-specific meta-analysis identifies a total of 351 loci, including 17 previously unreported loci. The genetic architecture of serum urate control is similar between European and East Asian populations. A transcriptome-wide association study, enrichment analysis, and colocalization analysis in relevant tissues identify candidate serum urate-associated genes, including CTBP1, SKIV2L, and WWP2. A phenome-wide association study using polygenic risk scores identifies serum urate-correlated diseases including heart failure and hypertension. Mendelian randomization and mediation analyses show that serum urate-associated genes might have a causal relationship with serum urate-correlated diseases via mediation effects. This study elucidates our understanding of the genetic architecture of serum urate control.

Genetic drivers of heterogeneity in type 2 diabetes pathophysiology

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10-8) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care.

Multi-ancestry genome-wide study in >2.5 million individuals reveals heterogeneity in mechanistic pathways of type 2 diabetes and complications

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes. To characterise the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study (GWAS) data from 2,535,601 individuals (39.7% non-European ancestry), including 428,452 T2D cases. We identify 1,289 independent association signals at genome-wide significance (P<5×10-8) that map to 611 loci, of which 145 loci are previously unreported. We define eight non-overlapping clusters of T2D signals characterised by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial, and enteroendocrine cells. We build cluster-specific partitioned genetic risk scores (GRS) in an additional 137,559 individuals of diverse ancestry, including 10,159 T2D cases, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned GRS are more strongly associated with coronary artery disease and end-stage diabetic nephropathy than an overall T2D GRS across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings demonstrate the value of integrating multi-ancestry GWAS with single-cell epigenomics to disentangle the aetiological heterogeneity driving the development and progression of T2D, which may offer a route to optimise global access to genetically-informed diabetes care.

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

A saturated map of common genetic variants associated with human height

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

Biological insights into systemic lupus erythematosus through an immune cell-specific transcriptome-wide association study

OBJECTIVE

Genome-wide association studies (GWAS) have identified >100 risk loci for systemic lupus erythematosus (SLE), but the disease genes at most loci remain unclear, hampering translation of these genetic discoveries. We aimed to prioritise genes underlying the 110 SLE loci that were identified in the latest East Asian GWAS meta-analysis.

METHODS

We built gene expression predictive models in blood B cells, CD4+ and CD8+ T cells, monocytes, natural killer cells and peripheral blood cells of 105 Japanese individuals. We performed a transcriptome-wide association study (TWAS) using data from the latest genome-wide association meta-analysis of 208 370 East Asians and searched for candidate genes using TWAS and three data-driven computational approaches.

RESULTS

TWAS identified 171 genes for SLE (p<1.0×10-5); 114 (66.7%) showed significance only in a single cell type; 127 (74.3%) were in SLE GWAS loci. TWAS identified a strong association between CD83 and SLE (p<7.7×10-8). Meta-analysis of genetic associations in the existing 208 370 East Asian and additional 1498 cases and 3330 controls found a novel single-variant association at rs72836542 (OR=1.11, p=4.5×10-9) around CD83. For the 110 SLE loci, we identified 276 gene candidates, including 104 genes at recently-identified SLE novel loci. We demonstrated in vitro that putative causal variant rs61759532 exhibited an allele-specific regulatory effect on ACAP1, and that presence of the SLE risk allele decreased ACAP1 expression.

CONCLUSIONS

Cell-level TWAS in six types of immune cells complemented SLE gene discovery and guided the identification of novel genetic associations. The gene findings shed biological insights into SLE genetic associations.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation

We assembled an ancestrally diverse collection of genome-wide association studies (GWAS) of type 2 diabetes (T2D) in 180,834 affected individuals and 1,159,055 controls (48.9% non-European descent) through the Diabetes Meta-Analysis of Trans-Ethnic association studies (DIAMANTE) Consortium. Multi-ancestry GWAS meta-analysis identified 237 loci attaining stringent genome-wide significance (P < 5 × 10-9), which were delineated to 338 distinct association signals. Fine-mapping of these signals was enhanced by the increased sample size and expanded population diversity of the multi-ancestry meta-analysis, which localized 54.4% of T2D associations to a single variant with >50% posterior probability. This improved fine-mapping enabled systematic assessment of candidate causal genes and molecular mechanisms through which T2D associations are mediated, laying the foundations for functional investigations. Multi-ancestry genetic risk scores enhanced transferability of T2D prediction across diverse populations. Our study provides a step toward more effective clinical translation of T2D GWAS to improve global health for all, irrespective of genetic background.

Phenome-wide association study of the major histocompatibility complex region in the Korean population identifies novel association signals

Human leukocyte antigen (HLA) gene variants in the major histocompatibility complex (MHC) region are associated with numerous complex human diseases and quantitative traits. Previous phenome-wide association studies (PheWAS) for this region demonstrated that HLA association patterns to the phenome have both population-specific and population-shared components. We performed MHC PheWAS in the Korean population by analyzing associations between phenotypes and genetic variants in the MHC region using the Korea Biobank Array project data samples from the Korean Genome and Epidemiology Study (KoGES) cohorts. Using this single-population dataset, we curated and analyzed 82 phenotypes for 125 673 Korean individuals after imputing HLA using CookHLA, a recently developed imputation framework. More than one-third of these phenotypes showed significant associations, confirming 56 known associations and discovering 13 novel association signals that were not reported previously. In addition, we analyzed heritability explained by the variants in the MHC region and genetic correlations among phenotypes based on the MHC variants.

Validation and genetic heritability estimation of known type 2 diabetes related variants in the Korean population

Genome-wide association studies (GWASs) facilitated the discovery of countless disease-associated variants. However, GWASs have mostly been conducted in European ancestry samples. Recent studies have reported that these European-based association results may reduce disease prediction accuracy when applied in non-Europeans. Therefore, previously reported variants should be validated in non-European populations to establish reliable scientific evidence for precision medicine. In this study, we validated known associations with type 2 diabetes (T2D) and related metabolic traits in 125,850 samples from a Korean population genotyped by the Korea Biobank Array (KBA). At the end of December 2020, there were 8,823 variants associated with glycemic traits, lipids, liver enzymes, and T2D in the GWAS catalog. Considering the availability of imputed datasets in the KBA genome data, publicly available East-Asian T2D summary statistics, and the linkage disequilibrium among the variants (r2 < 0.2), 2,900 independent variants were selected for further analysis. Among these, 1,837 variants (63.3%) were statistically significant (p ≤ 0.05). Most of the non-replicated variants (n = 1,063) showed insufficient statistical power and decreased minor allele frequencies compared with the replicated variants. Moreover, most of known variants showed <10% genetic heritability. These results could provide valuable scientific evidence for future study designs, the current power of GWASs, and future applications in precision medicine in the Korean population.

Rare coding variants in 35 genes associate with circulating lipid levels-A multi-ancestry analysis of 170,000 exomes

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

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.

The power of genetic diversity in genome-wide association studies of lipids

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

Determinants of penetrance and variable expressivity in monogenic metabolic conditions across 77,184 exomes

Hundreds of thousands of genetic variants have been reported to cause severe monogenic diseases, but the probability that a variant carrier develops the disease (termed penetrance) is unknown for virtually all of them. Additionally, the clinical utility of common polygenetic variation remains uncertain. Using exome sequencing from 77,184 adult individuals (38,618 multi-ancestral individuals from a type 2 diabetes case-control study and 38,566 participants from the UK Biobank, for whom genotype array data were also available), we apply clinical standard-of-care gene variant curation for eight monogenic metabolic conditions. Rare variants causing monogenic diabetes and dyslipidemias display effect sizes significantly larger than the top 1% of the corresponding polygenic scores. Nevertheless, penetrance estimates for monogenic variant carriers average 60% or lower for most conditions. We assess epidemiologic and genetic factors contributing to risk prediction in monogenic variant carriers, demonstrating that inclusion of polygenic variation significantly improves biomarker estimation for two monogenic dyslipidemias.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Meta-analysis of 208370 East Asians identifies 113 susceptibility loci for systemic lupus erythematosus

OBJECTIVE

Systemic lupus erythematosus (SLE), an autoimmune disorder, has been associated with nearly 100 susceptibility loci. Nevertheless, these loci only partially explain SLE heritability and their putative causal variants are rarely prioritised, which make challenging to elucidate disease biology. To detect new SLE loci and causal variants, we performed the largest genome-wide meta-analysis for SLE in East Asian populations.

METHODS

We newly genotyped 10 029 SLE cases and 180 167 controls and subsequently meta-analysed them jointly with 3348 SLE cases and 14 826 controls from published studies in East Asians. We further applied a Bayesian statistical approach to localise the putative causal variants for SLE associations.

RESULTS

We identified 113 genetic regions including 46 novel loci at genome-wide significance (p<5×10-8). Conditional analysis detected 233 association signals within these loci, which suggest widespread allelic heterogeneity. We detected genome-wide associations at six new missense variants. Bayesian statistical fine-mapping analysis prioritised the putative causal variants to a small set of variants (95% credible set size ≤10) for 28 association signals. We identified 110 putative causal variants with posterior probabilities ≥0.1 for 57 SLE loci, among which we prioritised 10 most likely putative causal variants (posterior probability ≥0.8). Linkage disequilibrium score regression detected genetic correlations for SLE with albumin/globulin ratio (rg=-0.242) and non-albumin protein (rg=0.238).

CONCLUSION

This study reiterates the power of large-scale genome-wide meta-analysis for novel genetic discovery. These findings shed light on genetic and biological understandings of SLE.

Genome-wide association study in a Korean population identifies six novel susceptibility loci for rheumatoid arthritis

OBJECTIVE

Genome-wide association studies (GWAS) in rheumatoid arthritis (RA) have discovered over 100 RA loci, explaining patient-relevant RA pathogenesis but showing a large fraction of missing heritability. As a continuous effort, we conducted GWAS in a large Korean RA case-control population.

METHODS

We newly generated genome-wide variant data in two independent Korean cohorts comprising 4068 RA cases and 36 487 controls, followed by a whole-genome imputation and a meta-analysis of the disease association results in the two cohorts. By integrating publicly available omics data with the GWAS results, a series of bioinformatic analyses were conducted to prioritise the RA-risk genes in RA loci and to dissect biological mechanisms underlying disease associations.

RESULTS

We identified six new RA-risk loci (SLAMF6, CXCL13, SWAP70, NFKBIA, ZFP36L1 and LINC00158) with pmeta<5×10-8 and consistent disease effect sizes in the two cohorts. A total of 122 genes were prioritised from the 6 novel and 13 replicated RA loci based on physical distance, regulatory variants and chromatin interaction. Bioinformatics analyses highlighted potentially RA-relevant tissues (including immune tissues, lung and small intestine) with tissue-specific expression of RA-associated genes and suggested the immune-related gene sets (such as CD40 pathway, IL-21-mediated pathway and citrullination) and the risk-allele sharing with other diseases.

CONCLUSION

This study identified six new RA-associated loci that contributed to better understanding of the genetic aetiology and biology in RA.

Identification of type 2 diabetes loci in 433,540 East Asian individuals

Meta-analyses of genome-wide association studies (GWAS) have identified more than 240 loci that are associated with type 2 diabetes (T2D)1,2; however, most of these loci have been identified in analyses of individuals with European ancestry. Here, to examine T2D risk in East Asian individuals, we carried out a meta-analysis of GWAS data from 77,418 individuals with T2D and 356,122 healthy control individuals. In the main analysis, we identified 301 distinct association signals at 183 loci, and across T2D association models with and without consideration of body mass index and sex, we identified 61 loci that are newly implicated in predisposition to T2D. Common variants associated with T2D in both East Asian and European populations exhibited strongly correlated effect sizes. Previously undescribed associations include signals in or near GDAP1, PTF1A, SIX3, ALDH2, a microRNA cluster, and genes that affect the differentiation of muscle and adipose cells3. At another locus, expression quantitative trait loci at two overlapping T2D signals affect two genes-NKX6-3 and ANK1-in different tissues4-6. Association studies in diverse populations identify additional loci and elucidate disease-associated genes, biology, and pathways.

Recapitulation of previously reported associations for type 2 diabetes and metabolic traits in the 126K East Asians

Over the last decade, genome-wide association studies (GWASs) have provided an unprecedented amount of genetic variations that are associated with various phenotypes. However, previous GWAS were mostly conducted in European populations, and these biased results for non-Europeans may result in a significant reduction in risk prediction for non-Europeans. An issue with the early GWAS was the winner’s curse problem, which led to misleading results when constructing the polygenic risk scores (PRS). Therefore, more non-European population-based studies are needed to validate reported variants and improve genetic risk assessment across diverse populations. In this study, we validated 422 variants independently associated with glycemic indexes, liver enzymes, and type 2 diabetes in 125,872 samples from a Korean population, and further validated the results by assessing publicly available summary statistics from European GWAS (n = 898,130). Among the 422 independently associated variants, 284, 320, and 361 variants were replicated in Koreans, Europeans, and either one of the two populations. In addition, the effect sizes for Koreans and Europeans were moderately correlated (r = 0.33–0.68). However, 61 variants were not replicated in both Koreans and Europeans. Our findings provide valuable information on effect sizes and statistical significance, which is essential to improve the assessment of disease risk using PRS analysis.

Programmable Nuclease-Based Integration into Novel Extragenic Genomic Safe Harbor Identified from Korean Population-Based CNV Analysis

Here, we found two genomic safe harbor (GSH) candidates from chromosomes 3 and 8, based on large-scale population-based cohort data from 4,694 Koreans by CNV analysis. Furthermore, estimated genotype of these CNVRs was validated by quantitative real-time PCR, and epidemiological data examined no significant genetic association between diseases or traits and two CNVRs. After screening the GSH candidates by in silico approaches, we designed TALEN pairs to integrate EGFP expression cassette into human cell lines in order to confirm the functionality of GSH candidates in an in vitro setting. As a result, transgene insertion into one of the two loci using TALEN showed robust transgene expression comparable to that with an AAVS1 site without significantly perturbing neighboring genes. Changing the promoter or cell type did not noticeably disturb this trend. Thus, we could validate two CNVRs as a site for effective and safe transgene insertion in human cells.

Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls

Protein-coding genetic variants that strongly affect disease risk can yield relevant clues to disease pathogenesis. Here we report exome-sequencing analyses of 20,791 individuals with type 2 diabetes (T2D) and 24,440 non-diabetic control participants from 5 ancestries. We identify gene-level associations of rare variants (with minor allele frequencies of less than 0.5%) in 4 genes at exome-wide significance, including a series of more than 30 SLC30A8 alleles that conveys protection against T2D, and in 12 gene sets, including those corresponding to T2D drug targets (P = 6.1 × 10-3) and candidate genes from knockout mice (P = 5.2 × 10-3). Within our study, the strongest T2D gene-level signals for rare variants explain at most 25% of the heritability of the strongest common single-variant signals, and the gene-level effect sizes of the rare variants that we observed in established T2D drug targets will require 75,000-185,000 sequenced cases to achieve exome-wide significance. We propose a method to interpret these modest rare-variant associations and to incorporate these associations into future target or gene prioritization efforts.

Exome chip-driven association study of lipidemia in >14,000 Koreans and evaluation of genetic effect on identified variants between different ethnic groups

Lipid levels in blood are widely used to diagnose and monitor chronic diseases. It is essential to identify the genetic traits involved in lipid metabolism for understanding chronic diseases. However, the influence of genetic traits varies depending on race, sex, age, and ethnicity. Therefore, research focusing on populations of individual countries is required, and the results can be used as a basis for comparison of results of other studies at the cross-racial and cross-country levels. In the present study, we selected lipid-related variants and evaluated their effects on lipid-related diseases in more than 14,000 subjects of three cohorts using the Illumina Human Exome Beadchip. A genome-wide association study was conducted using EPACTs after adjusting for age, sex, and recruitment area. A genome-wide significance cutoff was defined as p < 5E-08 in all the three cohorts. Sixteen variants represented the lipid traits and were classified as vulnerable to borderline hypertriglyceridemia, hyper-LDL-cholesterolemia, or hypo-HDL-cholesterolemia. Moreover, we compared the genetic effects of the 16 variants between ethnic groups and identified the missense variants in apolipoprotein A-V, cholesterol ester transfer protein, and apolipoprotein E as Asian-specific. Our study provides candidate genes as markers for chronic diseases through the evaluation of genetic effects.

The Korea Biobank Array: Design and Identification of Coding Variants Associated with Blood Biochemical Traits

We introduce the design and implementation of a new array, the Korea Biobank Array (referred to as KoreanChip), optimized for the Korean population and demonstrate findings from GWAS of blood biochemical traits. KoreanChip comprised >833,000 markers including >247,000 rare-frequency or functional variants estimated from >2,500 sequencing data in Koreans. Of the 833 K markers, 208 K functional markers were directly genotyped. Particularly, >89 K markers were presented in East Asians. KoreanChip achieved higher imputation performance owing to the excellent genomic coverage of 95.38% for common and 73.65% for low-frequency variants. From GWAS (Genome-wide association study) using 6,949 individuals, 28 associations were successfully recapitulated. Moreover, 9 missense variants were newly identified, of which we identified new associations between a common population-specific missense variant, rs671 (p.Glu457Lys) of ALDH2, and two traits including aspartate aminotransferase (P = 5.20 × 10⁻¹³) and alanine aminotransferase (P = 4.98 × 10⁻⁸). Furthermore, two novel missense variants of GPT with rare frequency in East Asians but extreme rarity in other populations were associated with alanine aminotransferase (rs200088103; p.Arg133Trp, P = 2.02 × 10⁻⁹ and rs748547625; p.Arg143Cys, P = 1.41 × 10⁻⁶). These variants were successfully replicated in 6,000 individuals (P = 5.30 × 10⁻⁸ and P = 1.24 × 10⁻⁶). GWAS results suggest the promising utility of KoreanChip with a substantial number of damaging variants to identify new population-specific disease-associated rare/functional variants.

Nonsynonymous Variants in PAX4 and GLP1R Are Associated With Type 2 Diabetes in an East Asian Population

We investigated ethnicity-specific exonic variants of type 2 diabetes (T2D) and its related clinical phenotypes in an East Asian population. We performed whole-exome sequencing in 917 T2D case and control subjects, and the findings were validated by exome array genotyping in 3,026 participants. In silico replication was conducted for seven nonsynonymous variants in an additional 13,122 participants. Single-variant and gene-based association tests for T2D were analyzed. A total of 728,838 variants were identified by whole-exome sequencing. Among nonsynonymous variants, PAX4 Arg192His increased risk of T2D and GLP1R Arg131Gln decreased risk of T2D in genome-wide significance (odds ratio [OR] 1.48, P = 4.47 × 10^-16 and OR 0.84, P = 3.55 × 10^-8, respectively). Another variant at PAX4 192 codon Arg192Ser was nominally associated with T2D (OR 1.62, P = 5.18 × 10^-4). In T2D patients, PAX4 Arg192His was associated with earlier age at diagnosis, and GLP1R Arg131Gln was associated with decreased risk of cardiovascular disease. In control subjects without diabetes, the PAX4 Arg192His was associated with higher fasting glucose and GLP1R Arg131Gln was associated with lower fasting glucose and HbA_1c level. Gene-based analysis revealed that SLC30A8 was most significantly associated with decreased risk of T2D (P = 1.0 × 10^-4). In summary, we have identified nonsynonymous variants associated with risk of T2D and related phenotypes in Koreans.

Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels

Large-scale meta-analyses of genome-wide association studies (GWAS) have identified >175 loci associated with fasting cholesterol levels, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). With differences in linkage disequilibrium (LD) structure and allele frequencies between ancestry groups, studies in additional large samples may detect new associations. We conducted staged GWAS meta-analyses in up to 69,414 East Asian individuals from 24 studies with participants from Japan, the Philippines, Korea, China, Singapore, and Taiwan. These meta-analyses identified (P < 5 × 10-8) three novel loci associated with HDL-C near CD163-APOBEC1 (P = 7.4 × 10-9), NCOA2 (P = 1.6 × 10-8), and NID2-PTGDR (P = 4.2 × 10-8), and one novel locus associated with TG near WDR11-FGFR2 (P = 2.7 × 10-10). Conditional analyses identified a second signal near CD163-APOBEC1. We then combined results from the East Asian meta-analysis with association results from up to 187,365 European individuals from the Global Lipids Genetics Consortium in a trans-ancestry meta-analysis. This analysis identified (log10Bayes Factor ≥6.1) eight additional novel lipid loci. Among the twelve total loci identified, the index variants at eight loci have demonstrated at least nominal significance with other metabolic traits in prior studies, and two loci exhibited coincident eQTLs (P < 1 × 10-5) in subcutaneous adipose tissue for BPTF and PDGFC. Taken together, these analyses identified multiple novel lipid loci, providing new potential therapeutic targets.

Positive correlation of cg16672562 methylation with obesity-related traits in childhood obesity, and its independence with underlying HIF3A (hypoxia-inducible factor 3a) genetic background

Differential methylations of the HIF3A (hypoxia-inducible factor 3a) gene have been linked to body mass index (BMI). To explore the association of these methylations to childhood obesity, we measured 5 CpG methylation sites (cg27146050, cg46801562, cg22891070, cg16672562 and cg46801675) in intron 1 of the HIF3A gene by pyrosequencing, in the Korean population (mean age: 13.9 yrs, 305 obese cases and 387 controls). Two CpG methylations, cg46801562 and cg16672562, had statistically significant association with childhood obesity (P = 2.09E-9 and 1.66E-7, respectively). Notably, in the case of cg16672562, all correlations were significantly positive with BMI (beta = 0.285, P = 1.652E-13), waist-hip ratio (beta = 0.0028, P = 1.42E-15) and fasting plasma glucose level (beta = 0.0645, P = 2.61E-4), when analyzed by linear regression, with age and sex as covariates. We investigated any genetic effect of cg16672562 methylation by using 14 single nucleotide polymorphisms (SNP) identified by exome sequencing of the HIF3A locus cg16672562 methylation showed no statistically significant changes due to the 14 polymorphisms. In this study, we show that cg16672562 is the most significant blood DNA methylation marker for childhood obesity in the Korean population, and might be independent of any underlying HIF3A genetic background.

On the association analysis of CNV data: a fast and robust family-based association method

BACKGROUND

Copy number variation (CNV) is known to play an important role in the genetics of complex diseases and several methods have been proposed to detect association of CNV with phenotypes of interest. Statistical methods for CNV association analysis can be categorized into two different strategies. First, the copy number is estimated by maximum likelihood and association of the expected copy number with the phenotype is tested. Second, the observed probe intensity measurements can be directly used to detect association of CNV with the phenotypes of interest.

RESULTS

For each strategy we provide a statistic that can be applied to extended families. The computational efficiency of the proposed methods enables genome-wide association analysis and we show with simulation studies that the proposed methods outperform other existing approaches. In particular, we found that the first strategy is always more efficient than the second strategy no matter whether copy numbers for each individual are well identified or not. With the proposed methods, we performed genome-wide CNV association analyses of hematological trait, hematocrit, on 521 Korean family samples.

CONCLUSIONS

We found that statistical analysis with the expected copy number is more powerful than the statistic with the probe intensity measurements regardless of the accuracy of the estimation of copy numbers.

Identification of a Copy Number Variation on Chromosome 20q13.12 Associated with Osteoporotic Fractures in the Korean Population

Osteoporotic fractures (OFs) are critical hard outcomes of osteoporosis and are characterized by decreased bone strength induced by low bone density and microarchitectural deterioration in bone tissue. Most OFs cause acute pain, hospitalization, immobilization, and slow recovery in patients and are associated with increased mortality. A variety of genetic studies have suggested associations of genetic variants with the risk of OF. Genome-wide association studies have reported various single-nucleotide polymorphisms and copy number variations (CNVs) in European and Asian populations. To identify CNV regions associated with OF risk, we conducted a genome-wide CNV study in a Korean population. We performed logistic regression analyses in 1,537 Korean subjects (299 OF cases and 1,238 healthy controls) and identified a total of 8 CNV regions significantly associated with OF (p < 0.05). Then, one CNV region located on chromosome 20q13.12 was selected for experimental validation. The selected CNV region was experimentally validated by quantitative polymerase chain reaction. The CNV region of chromosome 20q13.12 is positioned upstream of a family of long non-coding RNAs, LINC01260. Our findings could provide new information on the genetic factors associated with the risk of OF.

Evaluation of pleiotropic effects among common genetic loci identified for cardio-metabolic traits in a Korean population

Background

The genetic contribution to complex diseases or traits, including cardio-metabolic traits, has been elucidated recently by large-scale genome-wide association studies. These genome-wide association studies have indicated that most pleiotropic loci contain genes associated with lipids. Clinically, lipid related abnormalities are strongly associated with other diseases such as type 2 diabetes, coronary artery disease and hypertension. The aim of this study was to evaluate the shared genetic background of lipids and other cardio-metabolic traits.

Methods

We conducted meta-analyses of the association between 157 published lipid-associated loci and 10 cardio-metabolic traits in 14,028 Korean individuals genotyped using the Exome chip (Illumina HumanExome BeadChip). We also examined whether the pleiotropic effects of such loci constituted independent (i.e., biological) pleiotropy or mediated pleiotropy in these metabolic pathways.

Results

Eighteen lipid-associated loci were significantly associated with one of six cardio-metabolic traits after correction for multiple testing (P < 3.70 × 10⁻⁴). Region 12q24.12 had pleiotropic effects on fasting plasma glucose, blood pressure and obesity-related traits (body mass index and waist-hip ratio) independent of its effects on the lipid profile. Lipid risk scores, calculated according to whether or not subjects carried the risk allele for lipid traits, were significantly associated with fasting plasma glucose, blood pressure and obesity-related traits.

Conclusions

The 12q24.12 region showed ethnic-specific genetic pleiotropy among cardio-metabolic traits in this study. Our findings may help to account for molecular mechanisms based on shared genetic background underlying not only dyslipidemia, but also cardiovascular disease and type 2 diabetes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12933-016-0337-1) contains supplementary material, which is available to authorized users.

A genome-wide association study of copy-number variation identifies putative loci associated with osteoarthritis in Koreans

Background

OA is a complex disease caused by environmental and genetic risk factors. The purpose of this study is to identify candidate copy number variations (CNVs) associated with OA.

Methods

We performed a genome-wide association study of CNV to identify potential loci that confer susceptibility to or protection from OA. CNV genotyping was conducted using NimbleGen HD2 3 × 720K comparative hybridization array and included samples from 371 OA patients and 467 healthy controls. The putative CNV regions identified were confirmed with a TaqMan assay.

Results

We identified six genomic regions associated with OA encompassing CNV loci. None of six loci had previously been reported in genome-wide association studies with OA, although a genetic analysis suggested that they have functional effects. The protein product of a candidate risk gene for obesity, TNKS, targets Wnt inhibition, and this gene was significantly associated with hand and knee OA. Copy number deletion on TNKS was associated with a 1.37-fold decreased risk for OA. In addition, CA10, which shows a strong association with osteoporosis, was also significant in our study. Copy number deletion on this gene was associated with a 1.69-fold decreased risk for OA.

Conclusion

We identified several CNV loci that may contribute to OA susceptibility in Koreans. Further functional investigations of these genes are warranted to fully characterize their putative association.

Electronic supplementary material

The online version of this article (doi:10.1186/s12891-015-0531-4) contains supplementary material, which is available to authorized users.

Combinatorial approach to estimate copy number genotype using whole-exome sequencing data

Copy number variations (CNVs) are known risk factors in complex diseases. Array-based approaches have been widely used to detect CNVs, but limitations of array-based CNV detection methods, such as noisy signal and low resolution, have hindered detection of small CNVs. Recently, the development of next-generation sequencing techniques has increased rapidly owing to declines in cost. Particularly, whole-exome sequencing has proved useful for finding causal genes and variants in complex diseases. Because gene copy number may affect expression, CNV genotyping can be very valuable in disease association studies. However, almost all current CNV detection tools consider only two types of CNV genotypes. In this study, we propose a CNV genotype estimation approach using a combination of existing methods. Our approach was comprehensively compared with the customized Agilent array-comparative genomic hybridization. We found that our genotyping approach proved to be accurate, and reproducible, suggesting that it can complement existing CNV genotyping methods.

Identification of a genetic variant at 2q12.1 associated with blood pressure in East Asians by genome-wide scan including gene-environment interactions

Background

Genome-wide association studies have identified many genetic loci associated with blood pressure (BP). Genetic effects on BP can be altered by environmental exposures via multiple biological pathways. Especially, obesity is one of important environmental risk factors that can have considerable effect on BP and it may interact with genetic factors. Given that, we aimed to test whether genetic factors and obesity may jointly influence BP.

Methods

We performed meta-analyses of genome-wide association data for systolic blood pressure (SBP) and diastolic blood pressure (DBP) that included analyses of interaction between single nucleotide polymorphisms (SNPs) and the obesity-related anthropometric measures, body mass index (BMI), height, weight, and waist/hip ratio (WHR) in East-Asians (n = 12,030).

Results

We identified that rs13390641 on 2q12.1 demonstrated significant association with SBP when the interaction between SNPs and BMI was considered (P < 5 × 10 ^-8). The gene located nearest to rs13390641, TMEM182, encodes transmembrane protein 182. In stratified analyses, the effect of rs13390641 on BP was much stronger in obese individuals (BMI ≥ 30) than non-obese individuals and the effect of BMI on BP was strongest in individuals with the homozygous A allele of rs13390641.

Conclusions

Our analyses that included interactions between SNPs and environmental factors identified a genetic variant associated with BP that was overlooked in standard analyses in which only genetic factors were included. This result also revealed a potential mechanism that integrates genetic factors and obesity related traits in the development of high BP.

KGVDB: a population-based genomic map of CNVs tagged by SNPs in Koreans

SUMMARY

Despite a growing interest in a correlation between copy number variations (CNVs) and flanking single nucleotide polymorphisms, few databases provide such information. In particular, most information on CNV available so far was obtained in Caucasian and Yoruba populations, and little is known about CNV in Asian populations. This article presents a database that provides CNV regions tagged by single nucleotide polymorphisms in about 4700 Koreans, which were detected under strict quality control, manually curated and experimentally validated.

AVAILABILITY

KGVDB is freely available for non-commercial use at http://biomi.cdc.go.kr/KGVDB.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Gene-based copy number variation study reveals a microdeletion at 12q24 that influences height in the Korean population

Height is a classic polygenic trait with high heritability (h(2)=0.8). Recent genome-wide association studies have revealed many independent loci associated with human height. In addition, although many studies have reported an association between copy number variation (CNV) and complex diseases, few have explored the relationship between CNV and height. Recent studies reported that single nucleotide polymorphisms (SNPs) are highly correlated with common CNVs, suggesting that it is warranted to survey CNVs to identify additional genetic factors affecting heritable traits such as height. This study tested the hypothesis that there would be CNV regions (CNVRs) associated with height nearby genes from the GWASs known to affect height. We identified regions containing >1% copy number deletion frequency from 3667 population-based cohort samples using the Illumina HumanOmni1-Quad BeadChip. Among the identified CNVRs, we selected 15 candidate regions that were located within 1Mb of 283 previously reported genes. To assess the effect of these CNVRs on height, statistical analyses were conducted with samples from a case group of 370 taller (upper 10%) individuals and a control group of 1828 individuals (lower 50%). We found that a newly identified 17.7 kb deletion at chromosomal position 12q24.33, approximately 171.6 kb downstream of GPR133, significantly correlated with height; this finding was validated using quantitative PCR. These results suggest that CNVs are potentially important in determining height and may contribute to height variation in human populations.

MR findings of the spinal paraganglioma : report of three cases

Extraadrenal paragangliomas involving the spine is less common and usually takes the form of intradural compression of the cauda equina. The authors report three cases of spinal paragangliomas resulting in extradural spinal cord compression and their MR findings. The MR imaging revealed a well-demarcated extradural mass with low to intermediate signal intensity on T1-weighted images and intermediate to high signal intensity on T2-weighted images compared to paravertebral muscles. After Gd-DTPA administration, heterogeneous and intense enhancement was found. Multiple punctate and serpiginous structures of signal void due to high-velocity flow were noted around and within the tumors on all sequences. In one case, the signal void structures were well corresponded with feeding arteries on angiography. These may be the characteristic findings of the extraadrenal paraganglioma involving the spine.

Expressed sequence tags (ESTs) analysis of Acanthamoeba healyi

Randomly selected 435 clones from Acanthamoeba healyi cDNA library were sequenced and a total of 387 expressed sequence tags (ESTs) had been generated. Based on the results of BLAST search, 130 clones (34.4%) were identified as the genes encoding surface proteins, enzymes for DNA, energy production or other metabolism, kinases and phosphatases, protease, proteins for signal transduction, structural and cytoskeletal proteins, cell cycle related proteins, transcription factors, transcription and translational machineries, and transporter proteins. Most of the genes (88.5%) are newly identified in the genus Acanthamoeba. Although 15 clones matched the genes of Acanthamoeba located in the public databases, twelve clones were actin gene which was the most frequently expressed gene in this study. These ESTs of Acanthamoeba would give valuable information to study the organism as a model system for biological investigations such as cytoskeleton or cell movement, signal transduction, transcriptional and translational regulations. These results would also provide clues to elucidate factors for pathogenesis in human granulomatous amoebic encephalitis or keratitis by Acanthamoeba.

Clinical and neuropsychological profiles of obsessive-compulsive schizophrenia: a pilot study

This pilot study compared characteristics of obsessive-compulsive (OC) schizophrenic patients and a matched non-OC schizophrenic control group. The OC-schizophrenic group required more intensive clinical interventions and had a poorer clinical course, lower levels of functioning, and longer periods of hospitalization. They showed greater negative symptoms and more impaired executive functioning. These findings suggest OC-schizophrenic patients may have an atypical set of clinical and neuropsychiatric characteristics, perhaps constituting a subgroup within the schizophrenia spectrum. Pathophysiology and possible treatment implications require further study.

Effect of arachidonic acid metabolic inhibitors on hypoxia/reoxygenation-induced renal cell injury

The present study was undertaken to examine the role of arachidonic acid (AA) metabolites in hypoxia/reoxygenation (H/R)-induced renal cell injury in rabbit renal cortical slices using AA metabolic inhibitors. Inhibitors of cyclooxygenase (indomethacin and diclofenac sodium) and lipoxygenase pathways (nordihydroguaiaretic acid, caffeic acid, and eicosapentaenoic acid) reduced H/R-induced LDH release in a dose-dependent manner, whereas an inhibitor of cytochrome P-450 monooxygenase pathway ethoxyresorufin was not effective. AA increased LDH release in control slices, and the effect was not altered by indomethacin and nordihydroguaiaretic acid. The protective effect of indomethacin was not affected by addition of PGE2, a main product of cyclooxygenase pathway in the kidney. H2O2-induced LDH release was prevented by inhibitors of lipoxygenase but not by inhibitors of cyclooxygenase and cytochrome P-450 monooxygenase H/R-induced LDH release was not altered by iron chelators, phenanthroline and deferoxamine, and a potent antioxidant, N,N'-diphenyl-p-phenylenediamine, suggesting that the H/R-induced cell injury is not attributed to a generation of reactive oxygen species. Morphological studies showed that H/R-induced structural changes including cell necrosis were significantly prevented by indomethacin. These results suggest that inhibitors of cyclooxygenase and lipoxygenase pathways exert a direct protective effect against the H/R-induced cell injury in renal tubules. Whether these effects are mediated by alterations of AA metabolic pathways is not certain.

Phylogenetic relationships among Acanthamoeba spp. based on PCR-RFLP analyses of mitochondrial small subunit rRNA gene

We investigated the value of mitochondrial small subunit rRNA gene (mt SSU rDNA) PCR-RFLP as a taxonomic tool for Acanthamoeba isolates with close inter-relationships. Twenty-five isolates representing 20 species were included in the analysis. As in nuclear 18S rDNA analysis, two type strains (A. astronyxis and A. tubiashi) of morphological group 1 diverged earliest from the other strains, but the divergence between them was less than in 18S riboprinting. Acanthamoeba griffini of morphological group 2 branched between pathogenic (A. culbertsoni A-1 and A. healyi OC-3A) and nonpathogenic (A. palestinensis Reich, A. pustulosa GE-3a, A. royreba Oak Ridge, and A lenticulata PD2S) strains of morphological group 3. Among the remaining isolates of morphological group 2, the Chang strain had the identical mitochondrial riboprints as the type strain of A. hatchetti. AA2 and AA1, the type strains of A. divionensis and A. paradivionensis, respectively, had the identical riboprints as A. quina Vil3 and A. castellanii Ma. Although the branching orders of A. castellanii Neff, A. polyphaga P23, A. triangularis SH621, and A. lugdunensis L3a were different from those in 18S riboprinting analysis, the results obtained from this study generally coincided well with those from 18S riboprinting. Mitochondrial riboprinting may have an advantage over nuclear 18S rDNA riboprinting because the mt SSU rDNAs do not seem to have introns that are found in the 18S genes of Acanthamoeba and that distort phylogenetic analyses.

Generation and chromosome mapping of expressed sequence tags (ESTs) from a human infant thymus

In an effort to identify novel genes that are expressed differentially in an infant thymus, we constructed an oligo-d(T) primed cDNA library from a human infant thymus followed by single-run partial sequencing to generate expressed sequence tags (ESTs). Characterization of more than 1400 sequences enabled us to convert human thymus transcripts into 1223 useful ESTs. These ESTs consisted of 613 (50.1%) showing homology to known human genes, 51 (4.2%) matching to genes from other species, 289 (23.6%) matching ESTs of unknown functions, and 182 (14.9%) being novel transcripts. The expression profile of an infant thymus features a high number of genes related to cell division-DNA synthesis and gene-protein expression, indicating the active growth stage of an infant thymus. To identify the chromosomal localization of 43 thymus ESTs, PCR-based mapping was performed using a human-rodent somatic cell hybrid or radiation hybrid mapping panel. The results indicated that several novel genes were determined to be located in the vicinity of previously mapped disease loci; histidinemia loci, plasminogen Tochigi disease loci, Ehlers-Danlos syndrome, hypertriglyceridemia, thyroid resistance locus, ocular albinism, galactosemia, porphyria variegata, Charcot-Marie-tooth disease, FEOM (fibrosis of extraocular muscles), Prader-Willi syndrome.

Subgenus classification of Acanthamoeba by riboprinting

Subgenus classification of Acanthamoeba remains uncertain. Twenty-three reference strains of Acanthamoeba including 18 (neo)type-strains were subjected for classification at the subgenus level by riboprinting. PCR/RFLP analysis of 18S rRNA gene (rDNA). On the dendrogram reconstructed on the basis of riboprint analyses, two type-strains (A. astronyxis and A. tubiashi) of morphological group 1 diverged early from the other strains and were quite distinct from each other. Four type-strains of morphological group 3, A. culbertsoni, A. palestinensis, A. healyi were considered taxonomically valid, but A. pustulosa was regarded as an invalid synonym of A. palestinensis. Strains of morphological group 2 were classified into 6 subgroups. Among them, A. griffini which has an intron in its 18S rDNA was the most divergent from the remaining strains. Acanthamoeba castellanii Castellani, A. quina Vil3, A. lugdunensis L3a, A. polyphaga Jones, A. triangularis SH621, and A. castellanii Ma strains belonged to a subgroup, A. castellanii complex. However, A. quina and A. lugdunensis were regarded as synonyms of A. castellanii. The Chang strain could be regarded as A. hatchetti. Acanthamoeba mauritaniensis, A. divionensis, A. paradivionensis could be considered as synonyms of A. rhysodes. Neff strain was regarded as A. polyphaga rather than as A. castellanii. It is likely that riboprinting can be applied for rapid identification of Acanthamoeba isolated from the clinical specimens and environments.

[Bacterial endosymbiosis within the cytoplasm of Acanthamoeba lugdunensis isolated from a contact lens storage case]

Transmission electron microscopy of an Acanthamoeba isolate (KA/L5) from a contact lens case revealed bacterial endosymbionts within cytoplasm of the amoebae. The Acanthamoeba isolate belonged to the morphological group II. Based on the polymerase chain reaction (PCR)-restriction fragment length polymorphism (RFLP) of 18S ribosomal RNA coding DNA (rDNA), the isolate was identified as A. lugdunensis. Strain typing by isoenzyme analysis using isoelectric focusing (IEF) and mitochondrial (Mt) DNA RFLP revealed that the isolate was closely related with KA/L1, the most predominant type of isolates from contact lens storage cases, KA/E2, a clinical isolate, KA/W4, previously reported to host endosymbionts, and L3a strains of A. lugdunensis. The endosymbionts were similar to those of KA/W4 in aspects that they were randomly distributed in both trophozoites and cysts, and were rod-shaped bacteria measuring approximately 1.38 x 0.50 microns. But the number of endosymbionts per amoeba was significantly lower than that of KA/W4. They were neither limited by phagosomal membranes nor included in lacunaelike structure.

Close relatedness of Acanthamoeba pustulosa with Acanthamoeba palestinensis based on isoenzyme profiles and rDNA PCR-RFLP patterns

The taxonomic validity of morphological group III Acanthamoeba spp. is uncertain. In the present study, six type strains of group III Acanthamoeba spp., A. culbertsoni, A. healyi, A. pustulosa, A. palestinensis, A. royreba and A. lenticulata were subjected for the evaluation of their taxonomic validity by comparison of the isoenzyme patterns by isoelectic focusing on polyacrylamide gels, mitochondrial DNA (Mt DNA) restriction fragment length polymorphism (RFLP), and small subunit ribosomal DNA (ssu rDNA) PCR-RFLP patterns. The Mt DNA RFLP patterns were heterogeneous between the species. The type strains of A. palestinensis and A. pustulosa showed almost identical patterns of isoenzymes and rDNA PCR RFLP with an estimated sequence divergence of 2.6%. The other species showed heterogeneous patterns of isoenzymes and rDNA PCR-RFLP. It is likely that A. pustulosa is closely related with A. palestinensis and that the former may be regarded as a junior synonym of the latter.

Effect of high-velocity missile injury on plasma cortisol concentration in dogs

Plasma cortisol concentrations were measured in nine dogs after high-velocity missile injury induced by spherical steel bullets. All nine dogs sustained soft-tissue wounds in the thighs without damage to large vessels, femur, or nerves. When the dogs were considered as one group the mean plasma cortisol concentration increased promptly 10 min after injury (p less than 0.05) and reached its peak at 1 hour and then decreased. The mean concentration at 6 hours after injury was still higher than that before injury (p less than 0.05). Even though the severities of the injuries varied rather much there was a significant difference (p less than 0.01) between the plasma cortisol concentrations at 1 hour of the two groups into which the dogs were divided according to the severity of their injuries. The changes of the plasma cortisol concentrations may be parallel with the extent of the injuries and there was a positive correlation with the dimensions of the entrance wounds.