Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

The Polygenic and Monogenic Basis of Blood Traits and Diseases.

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.

Genetic diversity fuels gene discovery for tobacco and alcohol use.

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

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.

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.

Chromosomal alterations among age-related haematopoietic clones in Japan.

The extent to which the biology of oncogenesis and ageing are shaped by factors that distinguish human populations is unknown. Haematopoietic clones with acquired mutations become common with advancing age and can lead to blood cancers1-10. Here we describe shared and population-specific patterns of genomic mutations and clonal selection in haematopoietic cells on the basis of 33,250 autosomal mosaic chromosomal alterations that we detected in 179,417 Japanese participants in the BioBank Japan cohort and compared with analogous data from the UK Biobank. In this long-lived Japanese population, mosaic chromosomal alterations were detected in more than 35.0% (s.e.m., 1.4%) of individuals older than 90 years, which suggests that such clones trend towards inevitability with advancing age. Japanese and European individuals exhibited key differences in the genomic locations of mutations in their respective haematopoietic clones; these differences predicted the relative rates of chronic lymphocytic leukaemia (which is more common among European individuals) and T cell leukaemia (which is more common among Japanese individuals) in these populations. Three different mutational precursors of chronic lymphocytic leukaemia (including trisomy 12, loss of chromosomes 13q and 13q, and copy-neutral loss of heterozygosity) were between two and six times less common among Japanese individuals, which suggests that the Japanese and European populations differ in selective pressures on clones long before the development of clinically apparent chronic lymphocytic leukaemia. Japanese and British populations also exhibited very different rates of clones that arose from B and T cell lineages, which predicted the relative rates of B and T cell cancers in these populations. We identified six previously undescribed loci at which inherited variants predispose to mosaic chromosomal alterations that duplicate or remove the inherited risk alleles, including large-effect rare variants at NBN, MRE11 and CTU2 (odds ratio, 28-91). We suggest that selective pressures on clones are modulated by factors that are specific to human populations. Further genomic characterization of clonal selection and cancer in populations from around the world is therefore warranted.

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.

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.

Genotype imputation accuracy and the quality metrics of the minor ancestry in multi-ancestry reference panels.

Large-scale imputation reference panels are currently available and have contributed to efficient genome-wide association studies through genotype imputation. However, whether large-size multi-ancestry or small-size population-specific reference panels are the optimal choices for under-represented populations continues to be debated. We imputed genotypes of East Asian (180k Japanese) subjects using the Trans-Omics for Precision Medicine reference panel and found that the standard imputation quality metric (Rsq) overestimated dosage r2 (squared correlation between imputed dosage and true genotype) particularly in marginal-quality bins. Variance component analysis of Rsq revealed that the increased imputed-genotype certainty (dosages closer to 0, 1 or 2) caused upward bias, indicating some systemic bias in the imputation. Through systematic simulations using different template switching rates (θ value) in the hidden Markov model, we revealed that the lower θ value increased the imputed-genotype certainty and Rsq; however, dosage r2 was insensitive to the θ value, thereby causing a deviation. In simulated reference panels with different sizes and ancestral diversities, the θ value estimates from Minimac decreased with the size of a single ancestry and increased with the ancestral diversity. Thus, Rsq could be deviated from dosage r2 for a subpopulation in the multi-ancestry panel, and the deviation represents different imputed-dosage distributions. Finally, despite the impact of the θ value, distant ancestries in the reference panel contributed only a few additional variants passing a predefined Rsq threshold. We conclude that the θ value substantially impacts the imputed dosage and the imputation quality metric value.

Disease-specific variant interpretation highlighted the genetic findings in 2325 Japanese patients with retinitis pigmentosa and allied diseases.

BACKGROUND: As gene-specific therapy for inherited retinal dystrophy (IRD) advances, unified variant interpretation across institutes is becoming increasingly important. This study aims to update the genetic findings of 86 retinitis pigmentosa (RP)-related genes in a large number of Japanese patients with RP by applying the standardised variant interpretation guidelines for Japanese patients with IRD (J-IRD-VI guidelines) built upon the American College of Medical Genetics and Genomics and the Association for Molecular Pathology rules, and assess the contribution of these genes in RP-allied diseases. METHODS: We assessed 2325 probands with RP (n=2155, including n=1204 sequenced previously with the same sequencing panel) and allied diseases (n=170, newly analysed), including Usher syndrome, Leber congenital amaurosis and cone-rod dystrophy (CRD). Target sequencing using a panel of 86 genes was performed. The variants were interpreted according to the J-IRD-VI guidelines. RESULTS: A total of 3564 variants were detected, of which 524 variants were interpreted as pathogenic or likely pathogenic. Among these 524 variants, 280 (53.4%) had been either undetected or interpreted as variants of unknown significance or benign variants in our earlier study of 1204 patients with RP. This led to a genetic diagnostic rate in 38.6% of patients with RP, with EYS accounting for 46.7% of the genetically solved patients, showing a 9% increase in diagnostic rate from our earlier study. The genetic diagnostic rate for patients with CRD was 28.2%, with RP-related genes significantly contributing over other allied diseases. CONCLUSION: A large-scale genetic analysis using the J-IRD-VI guidelines highlighted the population-specific genetic findings for Japanese patients with IRD; these findings serve as a foundation for the clinical application of gene-specific therapies.

TNFRSF10A downregulation induces retinal pigment epithelium degeneration during the pathogenesis of age-related macular degeneration and central serous chorioretinopathy.

Age-related macular degeneration (AMD) and central serous chorioretinopathy (CSC) are common diseases that can cause vision loss in older and younger populations. These diseases share pathophysiological conditions derived from retinal pigment epithelium (RPE) dysfunction. Tumor necrosis factor receptor superfamily 10A (TNFRSF10A)-LOC389641 with the same lead single-nucleotide polymorphism (SNP) (rs13278062) is the only overlapped susceptibility locus found in both AMD and CSC through genome-wide association studies. This lead SNP has been reported to alter the transcriptional activity of TNFRSF10A. This study aimed to elucidate the function of TNFRSF10A in RPE degeneration using human primary RPE cells and Tnfrsf10 knockout (Tnfrsf10-/-) mice. TNFRSF10A was found to be localized in human RPE. In vitro assays revealed that a T allele of rs13278062, the risk allele for AMD and CSC, downregulated TNFRSF10A transcription in RPE, leading to decreased cell viability and increased apoptosis through protein kinase C-α (PKCA) downregulation. Treatment with phorbol 12-myristate 13-acetate, a PKC activator, rescued the cell viability. Morphological RPE abnormality was found in the retina of Tnfrsf10-/- mice. Our data suggest that downregulation of TNFRSF10A expression inactivates PKCA signaling and causes cellular vulnerability of the RPE, which may contribute to the pathogenesis of AMD and CSC.

Arteriolar hyalinization at 0-hour biopsy predicts long-term graft function in deceased kidney transplantation.

OBJECTIVES: Regarding the relationship between donor kidney quality and renal graft function after deceased kidney transplantation (KTx) following donation after cardiac death (DCD), the evaluation timing varies depending on the study. Evaluation of histology and changes in long-term renal graft function is limited. METHODS: A retrospective single-center study included 71 recipients who underwent 0-hour biopsy for KTx from DCD. The recipients were divided into two groups to evaluate factors related to renal graft function (study1). The two groups were categorized as stable graft function and poor graft function with the change of estimated glomerular filtration rate (eGFR) after KTx. The recipients were then divided into four groups to assess whether the factors identified in study1 were related to the change in long-term renal graft function (study2). They were categorized as follows: Improved, Stable, Deteriorated, and Primary non-function with the change of eGFR after KTx. RESULTS: In study1, donor age ≥ 50 years (29.5% vs. 65.2%; p = 0.09), banff arteriolar hyalinosis (ah) score (0.66 ± 0.78 vs. 1.2 ± 1.0; p = 0.018), and presence of glomerulosclerosis (43.2% vs. 76.2%; p = 0.017) were significant risk factors for poor long-term graft function. When the recipients were divided into four groups, the severity of ah correlated well with changes in long-term renal function. CONCLUSIONS: We can predict the shift in long-term renal graft function after KTx from DCD according to the severity of ah by 0-hour biopsy.

Evidence of Polygenic Adaptation in Sardinia at Height-Associated Loci Ascertained from the Biobank Japan.

Adult height is one of the earliest putative examples of polygenic adaptation in humans. However, this conclusion was recently challenged because residual uncorrected stratification from large-scale consortium studies was considered responsible for the previously noted genetic difference. It thus remains an open question whether height loci exhibit signals of polygenic adaptation in any human population. We re-examined this question, focusing on one of the shortest European populations, the Sardinians, in addition to mainland European populations. We utilized height-associated loci from the Biobank Japan (BBJ) dataset to further alleviate concerns of biased ascertainment of GWAS loci and showed that the Sardinians remain significantly shorter than expected under neutrality (∼0.22 standard deviation shorter than Utah residents with ancestry from northern and western Europe [CEU] on the basis of polygenic height scores, p = 3.89 × 10-4). We also found the trajectory of polygenic height scores between the Sardinian and the British populations diverged over at least the last 10,000 years (p = 0.0082), consistent with a signature of polygenic adaptation driven primarily by the Sardinian population. Although the polygenic score-based analysis showed a much subtler signature in mainland European populations, we found a clear and robust adaptive signature in the UK population by using a haplotype-based statistic, the trait singleton density score (tSDS), driven by the height-increasing alleles (p = 9.1 × 10-4). In summary, by ascertaining height loci in a distant East Asian population, we further supported the evidence of polygenic adaptation at height-associated loci among the Sardinians. In mainland Europeans, the adaptive signature was detected in haplotype-based analysis but not in polygenic score-based analysis.

Genome-Wide Association Study of Age-Related Macular Degeneration Reveals 2 New Loci Implying Shared Genetic Components with Central Serous Chorioretinopathy.

PURPOSE: To investigate the genetic architecture of age-related macular degeneration (AMD) in a Japanese population. DESIGN: Genome-wide association study (GWAS). PARTICIPANTS: Three thousand seven hundred seventy-two patients with AMD and 16 770 control participants from the Japanese population were enrolled in the association analyses. METHODS: We conducted a meta-analysis of 2 independent GWASs that included a total of 2663 patients with AMD and 9471 control participants using the imputation reference panel for genotype imputation specified for the Japanese population (n = 3541). A replication study was performed using an independent set of 1109 patients with AMD and 7299 control participants. MAIN OUTCOME MEASURES: Associations of genetic variants with AMD. RESULTS: A meta-analysis of the 2 GWASs identified 6 loci significantly associated with AMD (P < 5.0 × 10-8). Of these loci, 4 were known to be associated with AMD (CFH, C2/FB, TNFRSF10A, and ARMS2), and 2 were novel (rs4147157 near WBP1L and rs76228488 near GATA5). The newly identified associations were confirmed in a replication study (P < 0.01). After the meta-analysis of all datasets, we observed strong associations in these loci (P = 1.88 × 10-12 and P = 1.35 × 10-9 for meta-analysis for rs4147157 and rs76228488, respectively). When we looked up the associations in the reported central serous chorioretinopathy (CSC) GWAS conducted in the Japanese population, both loci were associated significantly with CSC (P = 4.86 × 10-3 and P = 4.28 × 10-3 for rs4147157 and rs76228488, respectively). We performed a genetic colocalization analysis for these loci and estimated that the posterior probabilities of shared causal variants between AMD and CSC were 0.39 and 0.60 for WBP1L and GATA5, respectively. Genetic correlation analysis focusing on the epidemiologically suggested clinical risk factors implicated shared polygenic architecture between AMD and smoking cessation (rg [the measure of genetic correlation] = -0.33; P = 0.01; false discovery rate, 0.099). CONCLUSIONS: Our findings imply shared genetic components conferring the risk of both AMD and CSC. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found after the references.

Mapping the genetic architecture of suicide attempt and suicide death using polygenic risk scores for clinically-related psychiatric disorders and traits.

BACKGROUND: Suicidal behavior is moderately heritable and a consequence of a combination of the diathesis traits for suicidal behavior and suicide-related major psychiatric disorders. Here, we sought to examine shared polygenic effects between various psychiatric disorders/traits and suicidal behavior and to compare the shared polygenic effects of various psychiatric disorders/traits on non-fatal suicide attempt and suicide death. METHODS: We used our genotyped European ancestry sample of 260 non-fatal suicide attempters, 317 suicide decedents and 874 non-psychiatric controls to test whether polygenic risk scores (PRSs) obtained from large GWASs for 22 suicide-related psychiatric disorders/traits were associated with suicidal behavior. Results were compared between non-fatal suicide attempt and suicide death in a sensitivity analysis. RESULTS: PRSs for major depressive disorder, bipolar disorder, schizophrenia, ADHD, alcohol dependence, sensitivity to environmental stress and adversity, educational attainment, cognitive performance, and IQ were associated with suicidal behavior (Bonferroni-corrected p < 2.5 × 10-4). The polygenic effects of all 22 psychiatric disorders/traits had the same direction (p for binomial tests = 4.8 × 10-7) and were correlated (Spearman's ρ = 0.85) between non-fatal suicide attempters and suicide decedents. CONCLUSIONS: We found that polygenic effects for major psychiatric disorders and diathesis-related traits including stress responsiveness and intellect/cognitive function contributed to suicidal behavior. While we found comparable polygenic architecture between non-fatal suicide attempters and suicide decedents based on correlations with PRSs of suicide-related psychiatric disorders/traits, our analyses are limited by small sample size resulting in low statistical power to detect difference between non-fatal suicide attempt and suicide death.

Likely pathogenic structural variants in genetically unsolved patients with retinitis pigmentosa revealed by long-read sequencing.

Despite the successful identification of causative genes and genetic variants of retinitis pigmentosa (RP), many patients have not been molecularly diagnosed. Our recent study using targeted short-read sequencing showed that the proportion of carriers of pathogenic variants in EYS, the cause of autosomal recessive RP, was unexpectedly high in Japanese patients with unsolved RP. This result suggested that causative genetic variants, which are difficult to detect by short-read sequencing, exist in such patients. Using long-read sequencing technology (Oxford Nanopore), we analysed the whole genomes of 15 patients with RP with one heterozygous pathogenic variant in EYS detected in our previous study along with structural variants (SVs) in EYS and another 88 RP-associated genes. Two large exon-overlapping deletions involving six exons were identified in EYS in two patients with unsolved RP. An analysis of an independent patient set (n=1189) suggested that these two deletions are not founder mutations. Our results suggest that searching for SVs by long-read sequencing in genetically unsolved cases benefits the molecular diagnosis of RP.

One-year efficacy of "rescue photodynamic therapy" for patients with typical age-related macular degeneration, polypoidal choroidal vasculopathy, and pachychoroid neovasculopathy refractory to anti-vascular endothelial growth factor therapy.

PURPOSE: This study aimed to evaluate the one-year outcomes of photodynamic therapy (PDT) as a rescue treatment for age-related macular degeneration (AMD) refractory to anti-vascular endothelial growth factor (VEGF) therapy. METHODS: Patients with AMD refractory to anti-VEGF therapy, treated with "rescue-PDT" were retrospectively investigated. The time of PDT was defined as the baseline value. Baseline characteristics including sex, age, best-corrected visual acuity (BCVA), central macular thickness (CMT), and foveal choroidal thickness (FCT) were examined. The changes in BCVA, CMT, and recurrence were also assessed at the 1-year follow-up. The logMAR VA change of 0.3 or more was defined as "improved" or "declined." RESULTS: Twenty-three consecutive eyes (typical AMD: 10 eyes, polypoidal choroidal vasculopathy: 10 eyes, and pachychoroid neovasculopathy: 3 eyes), which underwent "rescue-PDT," were analyzed in this study. The BCVA was improved in three patients and maintained in 20 patients at 12 months after PDT (mean BCVA change: 0.11 ± 0.19). The CMT improved in 19 patients (82.6%), and the mean CMT changed from 318.5 ± 93.7 µm to 225.9 ± 51.6 µm (p < 0.01) 12 months after PDT. "Retreatment" of anti-VEGF drug injections was considered if the retinal fluid or retinal hemorrhage recurred after PDT. The baseline FCT of the "retreatment group (15 eyes)" was significantly lower than that of the "no retreatment group (8 eyes)" (206.3 ± 50.7 µm vs 293.9 ± 85.7 µm: p = 0.033). CONCLUSIONS: PDT could be an effective treatment option for anti-VEGF refractory AMD to maintain visual acuity and control retinal fluid for up to 12 months.

RECURRENCE RATE OF CYSTOID MACULAR EDEMA WITH TOPICAL DORZOLAMIDE TREATMENT AND ITS RISK FACTORS IN RETINITIS PIGMENTOSA.

PURPOSE: To investigate the rate of the recurrence of cystoid macular edema (CME) secondary to retinitis pigmentosa (RP) after the initiation of topical dorzolamide and the recurrence risk factors. METHODS: We retrospectively analyzed the data of RP patients at Kyushu University Hospital. We included patients who showed a treatment response to 1.0% topical dorzolamide. The day of treatment initiation was set as the baseline. Topical dorzolamide treatment was continued during the follow-up. The recurrence of CME (defined as a >20% increase in central subfield thickness compared to previous visit, or a central subfield thickness value that exceed baseline value) was evaluated at each follow-up visit. Risk factors for RP-CME recurrence were analyzed by Cox proportional hazards modeling. A Kaplan-Meier survival analysis was used to evaluate the time to recurrent RP-CME. RESULTS: Forty RP-CME patients showed a treatment response to topical dorzolamide. During the mean 3.9-year follow-up, 14 patients exhibited recurrence; its rate was 15.6%, 34.7%, and 48.7% at 1, 3, and 5 years, respectively. A high baseline central subfield thickness was significantly associated with recurrent (hazard ratio 1.11, 95% CI: 1.05-1.18, P = 0.0004). CONCLUSION: The recurrence rate of RP-CME increased with time. A high baseline central subfield thickness value was a risk factor for recurrence.

Multi-omics study for interpretation of genome-wide association study.

Genome-wide association studies (GWASs) have identified thousands of genetic loci associated with complex traits, including a wide variety of diseases. Despite the successful identification of associated loci, interpreting GWAS findings remains challenging and requires other biological resources. Omics, including genomics, transcriptomics, proteomics, metabolomics, and epigenomics, are increasingly used in a broad range of research fields. Integrative analyses applying GWAS with these omics data are expected to expand our knowledge of complex traits and provide insight into the pathogenesis of complex diseases and their causative factors. Recently, associations between genetic variants and omics data have been comprehensively evaluated, providing new information on the influence of genetic variants on omics. Furthermore, recent advances in analytic methods, including single-cell technologies, have revealed previously unknown disease mechanisms. To advance our understanding of complex traits, integrative analysis using GWAS with multi-omics data is needed. In this review, I describe successful examples of integrative analyses based on omics and GWAS, discuss the limitations of current multi-omics analyses, and provide a perspective on future integrative studies.