Valeriana sambucifolia f. dageletiana (Nakai. ex Maekawa) Hara is a broad-leaved valerian endemic to Ulleung Island, a noted hot spot of endemism in Korea. However, despite its widespread pharmacological use, this plant remains comparatively understudied. Plant cells generally contain two types of organellar genomes (the plastome and the mitogenome) that have undergone independent evolution, which accordingly can provide valuable information for elucidating the phylogenetic relationships and evolutionary histories of terrestrial plants. Moreover, the extensive mega-data available for plant genomes, particularly those of plastomes, can enable researchers to gain an in-depth understanding of the transfer of genes between different types of genomes. In this study, we analyzed two organellar genomes (the 155,179 bp plastome and the 1,187,459 bp mitogenome) of V. sambucifolia f. dageletiana and detected extensive changes throughout the plastome sequence, including rapid structural mutations associated with inverted repeat (IR) contraction and genetic variation. We also described features characterizing the first reported mitogenome sequence obtained for a plant in the order Dipsacales and confirmed frequent gene transfer in this mitogenome. We identified eight non-plastome-originated regions (NPRs) distributed within the plastome of this endemic plant, for six of which there were no corresponding sequences in the current nucleotide sequence databases. Indeed, one of these unidentified NPRs unexpectedly showed certain similarities to sequences from bony fish. Although this is ostensibly difficult to explain, we suggest that this surprising association may conceivably reflect the occurrence of gene transfer from a bony fish to the plastome of an ancestor of V. sambucifolia f. dageletiana mediated by either fungi or bacteria.

Complete and accurate identification of genetic variants associated with specific phenotypes can be challenging when there is a high level of genomic divergence between individuals in a study and the corresponding reference genome. We have applied the Cas9-mediated enrichment coupled to nanopore sequencing to perform a targeted de novo assembly and accurately reconstruct a genomic region of interest. This approach was used to reconstruct a 250-kbp target region on chromosome 5 of the common bean genome (Phaseolus vulgaris) associated with the shattering phenotype. Comparing a non-shattering cultivar (Midas) with the reference genome revealed many single-nucleotide variants and structural variants in this region. We cut five 50-kbp tiled sub-regions of Midas genomic DNA using Cas9, followed by sequencing on a MinION device and de novo assembly, generating a single contig spanning the whole 250-kbp region. This assembly increased the number of Illumina reads mapping to genes in the region, improving their genotypability for downstream analysis. The Cas9 tiling approach for target enrichment and sequencing is a valuable alternative to whole-genome sequencing for the assembly of ultra-long regions of interest, improving the accuracy of downstream genotype-phenotype association analysis.

The phenotype of the 5α-reductase type 2 deficiency (5αRD2) by the SRD5A2 gene mutation varies, and although there have been many attempts, the genotype-phenotype correlation still has not yet been adequately evaluated. Recently, the crystal structure of the 5α-reductase type 2 isozyme (SRD5A2) has been determined. Therefore, the present study retrospectively evaluated the genotype-phenotype correlation from a structural perspective in 19 Korean patients with 5αRD2. Additionally, variants were classified according to structural categories, and phenotypic severity was compared with previously published data. The p.R227Q variant, which belongs to the NADPH-binding residue mutation category, exhibited a more masculine phenotype (higher external masculinization score) than other variants. Furthermore, compound heterozygous mutations with p.R227Q mitigated phenotypic severity. Similarly, other mutations in this category showed mild to moderate phenotypes. Conversely, the variants categorized as structure-destabilizing and small to bulky residue mutations showed moderate to severe phenotypes, and those categorized as catalytic site and helix-breaking mutations exhibited severe phenotypes. Therefore, the SRD5A2 structural approach suggested that a genotype-phenotype correlation does exist in 5αRD2. Furthermore, the categorization of SRD5A2 gene variants according to the SRD5A2 structure facilitates the prediction of the severity of 5αRD2 and the management and genetic counseling of patients affected by it.

As a major noncoding fragment, the control region (CR) of mtDNA is responsible for the initiation of mitogenome transcription and replication. Several structural features of CR sequences have been reported in many insects. However, comprehensive analyses on the structural organization and phylogenetic utility, as well as the role of tandem replications (TRs) on length variation, high A+T content, and shift of base skew of CR sequences are poorly investigated in hemipteran insects. In this study, we conducted a series of comparative analyses, using 116 samples covering all 11 infraorders of the five currently recognized monophyletic groups in the Hemiptera. Several structural elements (mononucleotide stretches containing conserved sequence blocks (CSBs), TRs, and GA-rich region) were identified in the mitochondrial control region in hemipteran insects, without showing a consistent location. The presence and absence of certain specific structural elements in CR sequences show the various structural organizations of that segment among the five monophyletic groups, which indicates the diversification of the control region’s structural organization in Hemiptera. Among the many groups within Hemiptera, eight monophyletic groups and three consistent phylogenetic trees were recovered, using CSBs datasets by maximum likelihood and Bayesian methods, which suggests the possible utility of CR sequences for phylogenetic reconstruction in certain groups of Hemiptera. Statistical analyses showed that TRs may contribute to the length variation, high AT content, and the shift of base skewing of CR sequences toward high AT content in the Hemiptera. Our findings enrich the knowledge of structural organization, phylogenetic utility, and roles of tandem replication of hemipteran CR, and provide a possible framework for mitochondrial control region analyses in hemimetabolous insects.

Calcineurin, a Ca(2+)/calmodulin-dependent protein phosphatase, plays a critical role in controlling skeletal muscle fiber type. However, little information is available concerning the expression of calcineurin in goat. Therefore, protein phosphatase 3 catalytic subunit alpha isoform (PPP3CA) gene, also called calcineurin Aα, was cloned and its expression characterized in Tianfu goat muscle. Real time quantitative polymerase chain reaction (RT-qPCR) analyses revealed that Tianfu goat PPP3CA was detected in cardiac muscle, biceps femoris muscle, abdominal muscle, longissimus dors muscle, and soleus muscle. High expression levels were found in biceps femoris muscle, longissimus muscle and abdominal muscle (p < 0.01), and low expression levels were seen in cardiac muscle and soleus muscle (p > 0.05). In addition, the spatial-temporal mRNA expression levels showed different variation trends in different muscles with the age of the goats. Western blotting further revealed that PPP3CA protein was expressed in the above-mentioned tissues, with the highest level in biceps femoris muscle, and the lowest level in soleus muscle. In this study, we isolated the full-length coding sequence of Tianfu goat PPP3CA gene, analyzed its structure, and investigated its expression in different muscle tissues from different age stages. These results provide a foundation for understanding the function of the PPP3CA gene in goats.

In plants, the Nuclear Factor-Y (NF-Y) transcription factors (TFs), which include three distinct types of NF-YA, NF-YB, and NF-YC TFs, have been identified to play key roles in the regulation of various plant growth and developmental processes under both normal and environmental stress conditions. In this work, a total of 40 CaNF-Y-encoding genes, including eight CaNF-YAs, 21 CaNF-YBs, and 11 CaNF-YCs, were identified in chickpea, and their major gene and protein characteristics were subsequently obtained using various web-based tools. Of our interest, a phylogenetically-based analysis predicted 18 CaNF-Ys (eight CaNF-YAs, seven CaNF-YBs, and three CaNF-YCs) that potentially play roles in chickpea responses to dehydration according to their close relationship with the well-characterized GmNF-Ys in soybean. These results were in good agreement with the enrichment of drought-responsive cis-regulatory motifs and expression patterns obtained from in silico analyses using publically available transcriptome data. Most of the phylogenetically predicted drought-responsive CaNF-Y genes (15 of 18) were quantitatively validated to significantly respond to dehydration treatment in leaves and/or roots, further supporting the results of in silico analyses. Among these CaNF-Y genes, the transcript levels of CaNF-YA01 and CaNF-YC10 were the most highly accumulated in leaves (by approximately eight-fold) and roots (by approximately 18-fold), respectively, by dehydration. Furthermore, 12 of the 18 CaNF-Y genes were found to be responsive to the most well-known stress hormone, namely abscisic acid (ABA), in leaves and/or roots, suggesting that these genes may act in chickpea response to dehydration in ABA-dependent manner. Taken together, our study has provided a comprehensive and fundamental information for further functional analyses of selected CaNF-Y candidate genes, ultimately leading to the improvement of chickpea growth under water-limited conditions.

DNA repair processes are involved in both the onset and treatment efficacy of colorectal cancer (CRC). A change of a single nucleotide causing an amino acid substitution in the corresponding protein may alter the efficiency of DNA repair, thus modifying the CRC susceptibility and clinical outcome. We performed a candidate gene approach in order to analyze the association of non-synonymous single nucleotide polymorphisms (nsSNPs) in the genes covering the main DNA repair pathways with CRC risk and clinical outcome modifications. Our candidate polymorphisms were selected according to the foremost genomic and functional prediction databases. Sixteen nsSNPs in 12 DNA repair genes were evaluated in cohorts from the Czech Republic and Austria. Apart from the tumor-node-metastasis (TNM) stage, which occurred as the main prognostic factor in all of the performed analyses, we observed several significant associations of different nsSNPs with survival and clinical outcomes in both cohorts. However, only some of the genes (REV3L, POLQ, and NEIL3) were prominently defined as prediction factors in the classification and regression tree analysis; therefore, the study suggests their association for patient survival. In summary, we provide observational and bioinformatics evidence that even subtle alterations in specific proteins of the DNA repair pathways may contribute to CRC susceptibility and clinical outcome.

Cerebral amyloid angiopathy (CAA) is a cerebrovascular disorder caused by the deposition of amyloid beta-peptide (Aβ) aggregates. Aβ aggregates lead to vessel rupture and intracerebral hemorrhages, detected by magnetic resonance imaging (MRI). Presenile CAA is usually genetically determined by mutations in the amyloid precursor protein (APP) gene. However, mutations after codon 200 in the presenilin 1 (PSEN1) gene have been reported to facilitate CAA onset. Here, we analyzed the genetic bases in a patient of 55 years old affected by CAA and cognitive decline. DNA was isolated and genetic analysis was performed by Next-Generation Sequencing (NGS). RNA was extracted and retro-transcribed to perform segregation analysis by TOPO-TA cloning. WB analysis was carried out to check the impact of the mutations on protein. Two compound heterozygous mutations in PSEN1 exon 10, such as a novel stop-gain mutation (c.1070C > G) and a pathogenic splice variant (c.1129A > T), were found by NGS. Both mutations altered the presenilin 1 protein, truncating its C-terminal portion. This is the first case of CAA and cognitive decline caused by two compound mutations in PSEN1. With this report, we suggest extending the genetic analysis to PSEN1 when cerebral microbleeds are observed by MRI investigation in a patient affected by presenile cognitive decline.

Amyloid precursor protein (APP), presenilin 1 (PSEN1), and presenilin 2 (PSEN2) are associated with autosomal-dominant early-onset Alzheimer's disease (AD). Most mutations have been identified in the PSEN1 gene. We discovered a PSEN1 mutation (Tyr389His) in a Korean patient with early-onset AD who presented memory decline at 41 years of age followed by language, memory, and visuospatial dysfunctions. As this is the third such patient identified in Korea, this mutation may be involved in AD pathogenesis, suggesting that routine screening is necessary in this population. Altered intra-molecular interactions with the mutated amino acid may result in the destabilization of γ-secretase. In the future, a panel incorporating genes with relatively high-frequency rare variants, along with the APOE4 gene, may predict the onset of AD and facilitate customized treatment.

Chromosomal instability (CIN) induces a high rate of losses or gains of whole chromosomes or parts of chromosomes. It is a hallmark of most human cancers and one of the causes of aneuploidy and intra-tumor heterogeneity. The present study aimed to evaluate the potential prognostic role of CIN in NB patients at diagnosis. We performed array comparative genomic hybridization analyses on 451 primary NB patients at the onset of the disease. To assess global chromosomal instability with high precision, we focused on the total number of DNA breakpoints of gains or losses of chromosome arms. For each tumor, an array-CGH-based breakpoint instability index (BPI) was assigned which defined the total number of chromosomal breakpoints per genome. This approach allowed us to quantify CIN related to whole genome disruption in all NB cases analyzed. We found differences in chromosomal breakages among the NB clinical risk groups. High BPI values are negatively associated with survival of NB patients. This association remains significant when correcting for stage, age, and MYCN status in the Cox model. Stratified analysis confirms the prognostic effect of BPI index in low-risk NB patients with non-amplified MYCN and with segmental chromosome aberrations.

Growth traits represent a main goal in aquaculture breeding programs and may be related to adaptive variation in wild fisheries. Integrating quantitative trait loci (QTL) mapping and next generation sequencing can greatly help to identify variation in candidate genes, which can result in marker-assisted selection and better genetic structure information. Turbot is a commercially important flatfish in Europe and China, with available genomic information on QTLs and genome mapping. Muscle and liver RNA-seq from 18 individuals was carried out to obtain gene sequences and markers functionally related to growth, resulting in a total of 20,447 genes and 85,344 single nucleotide polymorphisms (SNPs). Many growth-related genes and SNPs were identified and placed in the turbot genome and genetic map to explore their co-localization with growth-QTL markers. Forty-five SNPs on growth-related genes were selected based on QTL co-localization and relevant function for growth traits. Forty-three SNPs were technically feasible and validated in a wild Atlantic population, where 91% were polymorphic. The integration of functional and structural genomic resources in turbot provides a practical approach for QTL mining in this species. Validated SNPs represent a useful set of growth-related gene markers for future association, functional and population studies in this flatfish species.

Aceruloplasminemia is a rare autosomal recessive genetic disease characterized by mild microcytic anemia, diabetes, retinopathy, liver disease, and progressive neurological symptoms due to iron accumulation in pancreas, retina, liver, and brain. The disease is caused by mutations in the Ceruloplasmin (CP) gene that produce a strong reduction or absence of ceruloplasmin ferroxidase activity, leading to an impairment of iron metabolism. Most patients described so far are from Japan. Prompt diagnosis and therapy are crucial to prevent neurological complications since, once established, they are usually irreversible. Here, we describe the largest series of non-Japanese patients with aceruloplasminemia published so far, including 13 individuals from 11 families carrying 13 mutations in the CP gene (7 missense, 3 frameshifts, and 3 splicing mutations), 10 of which are novel. All missense mutations were studied by computational modeling. Clinical manifestations were heterogeneous, but anemia, often but not necessarily microcytic, was frequently the earliest one. This study confirms the clinical and genetic heterogeneity of aceruloplasminemia, a disease expected to be increasingly diagnosed in the Next-Generation Sequencing (NGS) era. Unexplained anemia with low transferrin saturation and high ferritin levels without inflammation should prompt the suspicion of aceruloplasminemia, which can be easily confirmed by low serum ceruloplasmin levels. Collaborative joint efforts are needed to better understand the pathophysiology of this potentially disabling disease.

Designer effectors based on the DNA binding domain (DBD) of Xanthomonas transcription activator-like effectors (TALEs) are powerful sequence-specific tools with an excellent reputation for their specificity in editing the genome, transcriptome, and more recently the epigenome in multiple cellular systems. However, the repetitive structure of the TALE arrays composing the DBD impedes their generation as gene synthesis product and prevents the delivery of TALE-based genes using lentiviral vectors (LVs), a widely used system for human gene therapy. To overcome these limitations, we aimed at chimerizing the DNA sequence encoding for the TALE-DBDs by introducing sufficient diversity to facilitate both their gene synthesis and enable their lentiviral delivery. To this end, we replaced three out of 17 Xanthomonas TALE repeats with TALE-like units from the bacterium Burkholderia rhizoxinica. This was combined with extensive codon variation and specific amino acid substitutions throughout the DBD in order to maximize intra- and inter-repeat sequence variability. We demonstrate that chimerized TALEs can be easily generated using conventional Golden Gate cloning strategy or gene synthesis. Moreover, chimerization enabled the delivery of TALE-based designer nucleases, transcriptome and epigenome editors using lentiviral vectors. When delivered as plasmid DNA, chimerized TALEs targeting the CCR5 and CXCR4 loci showed comparable activities in human cells. However, lentiviral delivery of TALE-based transcriptional activators was only successful in the chimerized form. Similarly, delivery of a chimerized CXCR4-specific epigenome editor resulted in rapid silencing of endogenous CXCR4 expression. In conclusion, extensive codon variation and chimerization of TALE-based DBDs enables both the simplified generation and the lentiviral delivery of designer TALEs, and therefore facilitates the clinical application of these tools to precisely edit the genome, transcriptome and epigenome.

Long non-coding, tandem-repetitive regions in mitochondrial (mt) genomes of many metazoans have been notoriously difficult to characterise accurately using conventional sequencing methods. Here, we show how the use of a third-generation (long-read) sequencing and informatic approach can overcome this problem. We employed Oxford Nanopore technology to sequence genomic DNAs from a pool of adult worms of the carcinogenic parasite, Schistosoma haematobium, and used an informatic workflow to define the complete mt non-coding region(s). Using long-read data of high coverage, we defined six dominant mt genomes of 33.4 kb to 22.6 kb. Although no variation was detected in the order or lengths of the protein-coding genes, there was marked length (18.5 kb to 7.6 kb) and structural variation in the non-coding region, raising questions about the evolution and function of what might be a control region that regulates mt transcription and/or replication. The discovery here of the largest tandem-repetitive, non-coding region (18.5 kb) in a metazoan organism also raises a question about the completeness of some of the mt genomes of animals reported to date, and stimulates further explorations using a Nanopore-informatic workflow.

The red or purple color of radish (Raphanus sativus L.) taproots is due to anthocyanins, which have nutritional and aesthetic value, as well as antioxidant properties. Moreover, the varied patterns and levels of anthocyanin accumulation in radish roots make them an interesting system for studying the transcriptional regulation of anthocyanin biosynthesis. The R2R3 MYB transcription factor RsMYB1 is a key positive regulator of anthocyanin biosynthesis in radish. Here, we isolated an allele of RsMYB1, named RsMYB1Short, in radish cultivars with white taproots. The RsMYB1Short allele carried a 4 bp insertion in the first exon causing a frame-shift mutation of RsMYB1, generating a truncated protein with only a partial R2 domain at the N-terminus. Unlike RsMYB1Full, RsMYB1Short was localized to the nucleus and the cytoplasm and failed to interact with their cognate partner RsTT8. Transient expression of genomic or cDNA sequences for RsMYB1Short in radish cotyledons failed to induce anthocyanin accumulation, but that for RsMYB1Full activated it. Additionally, RsMYB1Short showed the lost ability to induce pigment accumulation and to enhance the transcript level of anthocyanin biosynthetic genes, while RsMYB1Full promoted both processes when co-expressed with RsTT8 in tobacco leaves. As the result of the transient assay, co-expressing RsTT8 and RsMYB1Full, but not RsMYB1Short, also enhanced the promoter activity of RsCHS and RsDFR. We designed a molecular marker for RsMYB1 genotyping, and revealed that the RsMYB1Short allele is common in white radish cultivars, underscoring the importance of variation at the RsMYB1 locus in anthocyanin biosynthesis in the radish taproot. Together, these results indicate that the nonsense mutation of RsMYB1 generated the truncated protein, RsMYB1Short, that had the loss of ability to regulate anthocyanin biosynthesis. Our findings highlight that the frame shift mutation of RsMYB1 plays a key role in anthocyanin biosynthesis in the radish taproot.

Three Apiaceae species Ledebouriella seseloides, Peucedanum japonicum, and Glehnia littoralis are used as Asian herbal medicines, with the confusingly similar common name "Bang-poong". We characterized the complete chloroplast (cp) genomes and 45S nuclear ribosomal DNA (45S nrDNA) sequences of two accessions for each species. The complete cp genomes of G. littoralis, L. seseloides, and P. japonicum were 147,467, 147,830, and 164,633 bp, respectively. Compared to the other species, the P. japonicum cp genome had a huge inverted repeat expansion and a segmental inversion. The 45S nrDNA cistron sequences of the three species were almost identical in size and structure. Despite the structural variation in the P. japonicum cp genome, phylogenetic analysis revealed that G. littoralis diverged 5-6 million years ago (Mya), while P. japonicum diverged from L. seseloides only 2-3 Mya. Abundant copy number variations including tandem repeats, insertion/deletions, and single nucleotide polymorphisms, were found at the interspecies level. Intraspecies-level polymorphism was also found for L. seseloides and G. littoralis. We developed nine PCR barcode markers to authenticate all three species. This study characterizes the genomic differences between L. seseloides, P. japonicum, and G. littoralis; provides a method of species identification; and sheds light on the evolutionary history of these three species.

Severe bioprosthetic mitral valve calcification is a significant problem in cardiovascular surgery. Unfortunately, clinical markers did not demonstrate efficacy in prediction of severe bioprosthetic mitral valve calcification. Here, we examined whether a genomics-based approach is efficient in predicting the risk of severe bioprosthetic mitral valve calcification. A total of 124 consecutive Russian patients who underwent mitral valve replacement surgery were recruited. We investigated the associations of the inherited variation in innate immunity, lipid metabolism and calcium metabolism genes with severe bioprosthetic mitral valve calcification. Genotyping was conducted utilizing the TaqMan assay. Eight gene polymorphisms were significantly associated with severe bioprosthetic mitral valve calcification and were therefore included into stepwise logistic regression which identified male gender, the T/T genotype of the rs3775073 polymorphism within the TLR6 gene, the C/T genotype of the rs2229238 polymorphism within the IL6R gene, and the A/A genotype of the rs10455872 polymorphism within the LPA gene as independent predictors of severe bioprosthetic mitral valve calcification. The developed genomics-based model had fair predictive value with area under the receiver operating characteristic (ROC) curve of 0.73. In conclusion, our genomics-based approach is efficient for the prediction of severe bioprosthetic mitral valve calcification.

The KMT2A/AFF1 rearrangement is associated with an unfavorable prognosis in infant acute lymphocytic leukemia (ALL). Discordant ALL in monozygotic twins is uncommon and represents an attractive resource to evaluate intrauterine environment-genetic interplay in ALL. Mutational and epigenetic profiles were characterized for a discordant KMT2A/AFF1-rearranged infant monozygotic twin pair and their parents, and they were compared to three independent KMT2A/AFF1-positive ALL infants, in which the DNA methylation and gene expression profiles were investigated. A de novo Q61H NRAS mutation was detected in the affected twin at diagnosis and backtracked in both twins at birth. The KMT2A/AFF1 rearrangement was absent at birth in both twins. Genetic analyses conducted at birth gave more insights into the timing of the mutation hit. We identified correlations between DNA methylation and gene expression changes for 32 genes in the three independent affected versus remitted patients. The strongest correlations were observed for the RAB32, PDK4, CXCL3, RANBP17, and MACROD2 genes. This epigenetic signature could be a putative target for the development of novel epigenetic-based therapies and could help in explaining the molecular mechanisms characterizing ALL infants with KMT2A/AFF1 fusions.

Airway remodeling in asthma is characterized by reticular basement membrane (RBM) thickening, likely related to epithelial structural and functional changes. Gene expression profiling of the airway epithelium might identify genes involved in bronchial structural alterations. We analyzed bronchial wall geometry (computed tomography (CT)), RBM thickness (histology), and the bronchial epithelium transcriptome profile (gene expression array) in moderate to severe persistent (n = 21) vs. no persistent (n = 19) airflow limitation asthmatics. RBM thickness was similar in the two studied subgroups. Among the genes associated with increased RBM thickness, the most essential were those engaged in cell activation, proliferation, and growth (e.g., CDK20, TACC2, ORC5, and NEK5) and inhibiting apoptosis (e.g., higher mRNA expression of RFN34, BIRC3, NAA16, and lower of RNF13, MRPL37, CACNA1G). Additionally, RBM thickness correlated with the expression of genes encoding extracellular matrix (ECM) components (LAMA3, USH2A), involved in ECM remodeling (LTBP1), neovascularization (FGD5, HPRT1), nerve functioning (TPH1, PCDHGC4), oxidative stress adaptation (RIT1, HSP90AB1), epigenetic modifications (OLMALINC, DNMT3A), and the innate immune response (STAP1, OAS2). Cluster analysis revealed that genes linked with RBM thickness were also related to thicker bronchial walls in CT. Our study suggests that the pro-fibrotic profile in the airway epithelial cell transcriptome is associated with a thicker RBM, and thus, may contribute to asthma airway remodeling.

The phenotypic and genetic links between body fat phenotypes and primary open-angle glaucoma (POAG) are unclear. We conducted a meta-analysis of relevant longitudinal epidemiological studies to evaluate the phenotypic link. To identify genetic links, we performed genetic correlation analysis and pleiotropy analysis of genome-wide association study summary statistics datasets of POAG, intraocular pressure (IOP), vertical cup-to-disc ratio, obesity, body mass index (BMI), and waist-to-hip ratio. In the meta-analysis, we first established that obese and underweight populations have a significantly higher risk of POAG using longitudinal data. We also discovered positive genetic correlations between POAG and BMI and obesity phenotypes. Finally, we identified over 20 genomic loci jointly associated with POAG/IOP and BMI. Among them, the genes loci CADM2, RP3-335N17.2, RP11-793K1.1, RPS17P5, and CASC20 showed the lowest false discovery rate. These findings support the connection between body fat phenotypes and POAG. The newly identified genomic loci and genes render further functional investigation.