Fine-mapping Studies Research Articles

Exploring the genetic architecture of powdery mildew resistance in wheat through QTL meta-analysis.

Powdery mildew (PM), caused by Blumeria graminis f. sp. tritici, poses a significant threat to wheat production, necessitating the development of genetically resistant varieties for long-term control. Therefore, exploring genetic architecture of PM in wheat to uncover important genomic regions is an important area of wheat research. In recent years, the utilization of meta-QTL (MQTL) analysis has gained prominence as an essential tool for unraveling the complex genetic architecture underlying complex quantitative traits. The aim of this research was to conduct a QTL meta-analysis to pinpoint the specific genomic regions in wheat responsible for governing PM resistance. This study integrated 222 QTLs from 33 linkage-based studies using a consensus map with 54,672 markers. The analysis revealed 39 MQTLs, refined to 9 high-confidence MQTLs (hcMQTLs) with confidence intervals of 0.49 to 12.94 cM. The MQTLs had an average physical interval of 41.00 Mb, ranging from 0.000048 Mb to 380.71 Mb per MQTL. Importantly, 18 MQTLs co-localized with known resistance genes like Pm2, Pm3, Pm8, Pm21, Pm38, and Pm41. The study identified 256 gene models within hcMQTLs, providing potential targets for marker-assisted breeding and genomic prediction programs to enhance PM resistance. These MQTLs would serve as a foundation for fine mapping, gene isolation, and functional genomics studies, facilitating a deeper understanding of molecular mechanisms. The identification of candidate genes opens up exciting possibilities for the development of PM-resistant wheat varieties after validation.

Genetic association and transferability for urinary albumin-creatinine ratio as a marker of kidney disease in four Sub-Saharan African populations and non-continental individuals of African ancestry.

Genome-wide association studies (GWAS) have predominantly focused on populations of European and Asian ancestry, limiting our understanding of genetic factors influencing kidney disease in Sub-Saharan African (SSA) populations. This study presents the largest GWAS for urinary albumin-to-creatinine ratio (UACR) in SSA individuals, including 8,970 participants living in different African regions and an additional 9,705 non-resident individuals of African ancestry from the UK Biobank and African American cohorts. Urine biomarkers and genotype data were obtained from two SSA cohorts (AWI-Gen and ARK), and two non-resident African-ancestry studies (UK Biobank and CKD-Gen Consortium). Association testing and meta-analyses were conducted, with subsequent fine-mapping, conditional analyses, and replication studies. Polygenic scores (PGS) were assessed for transferability across populations. Two genome-wide significant (P < 5 × 10-8) UACR-associated loci were identified, one in the BMP6 region on chromosome 6, in the meta-analysis of resident African individuals, and another in the HBB region on chromosome 11 in the meta-analysis of non-resident SSA individuals, as well as the combined meta-analysis of all studies. Replication of previous significant results confirmed associations in known UACR-associated regions, including THB53, GATM, and ARL15. PGS estimated using previous studies from European ancestry, African ancestry, and multi-ancestry cohorts exhibited limited transferability of PGS across populations, with less than 1% of observed variance explained. This study contributes novel insights into the genetic architecture of kidney disease in SSA populations, emphasizing the need for conducting genetic research in diverse cohorts. The identified loci provide a foundation for future investigations into the genetic susceptibility to chronic kidney disease in underrepresented African populations Additionally, there is a need to develop integrated scores using multi-omics data and risk factors specific to the African context to improve the accuracy of predicting disease outcomes.

Genome-Wide Association Study Identifies IFIH1 and HLA-DQB1*05:02 Loci Associated With Anti-NMDAR Encephalitis.

Anti-N-methyl-d-aspartate receptor (NMDAR) encephalitis is a rare autoimmune neurologic disorder, the genetic etiology of which remains poorly understood. Our study aims to investigate the genetic basis of this disease in the Chinese Han population. We performed a genome-wide association study and fine-mapping study within the major histocompatibility complex (MHC) region of 413 Chinese patients with anti-NMDAR encephalitis recruited from 6 large tertiary hospitals and 7,127 healthy controls. Our genome-wide association analysis identified a strong association at the IFIH1 locus on chromosome 2q24.2 (rs3747517, p = 1.06 × 10-8, OR = 1.55, 95% CI, 1.34-1.80), outside of the human leukocyte antigen (HLA) region. Furthermore, through a fine-mapping study of the MHC region, we discovered associations for 3 specific HLA class I and II alleles. Notably, HLA-DQB1*05:02 (p = 1.43 × 10-12; OR, 2.10; 95% CI 1.70-2.59) demonstrates the strongest association among classical HLA alleles, closely followed by HLA-A*11:01 (p = 4.36 × 10-7; OR, 1.52; 95% CI 1.29-1.79) and HLA-A*02:07 (p = 1.28 × 10-8; OR, 1.87; 95% CI 1.50-2.31). In addition, we uncovered 2 main HLA amino acid variation associated with anti-NMDAR encephalitis including HLA-DQβ1-126H (p = 1.43 × 10-12; OR, 2.10; 95% CI 1.70-2.59), exhibiting a predisposing effect, and HLA-B-97R (p = 3.40 × 10-8; OR, 0.63; 95% CI 0.53-0.74), conferring a protective effect. Computational docking analysis suggested a close relationship between the NR1 subunit of NMDAR and DQB1*05:02. Our findings indicate that genetic variation in IFIH1, involved in the type I interferon signaling pathway and innate immunity, along with variations in the HLA class I and class II genes, has substantial implications for the susceptibility to anti-NMDAR encephalitis in the Chinese Han population.

Genomic regions associated with coat color in Gir cattle.

Indicine cattle breeds are adapted to the tropical climate, and their coat plays an important role in this process. Coat color influences thermoregulation and the adhesion of ectoparasites and may be associated with productive and reproductive traits. Furthermore, coat color is used for breed qualification, with breeders preferring certain colors. The Gir cattle is characterized by a wide variety of coat colors. Therefore, we performed genome-wide association studies to identify candidate genes for coat color in Gir cattle. Different phenotype scenarios were considered in the analyses and regions were identified on eight chromosomes. Some regions and many candidate genes are influencing coat color in the Gir cattle, which was found to be a polygenic trait. The candidate genes identified have been associated with white spotting patterns and base coat color in cattle and other species. In addition, a possible epistatic effect on coat color determination in the Gir cattle was suggested. This is the first published study that identified genomic regions and listed candidate genes associated with coat color in Gir cattle. The findings provided a better understanding of the genetic architecture of the trait in the breed and will allow to guide future fine-mapping studies for the development of genetic markers for selection.

Non-coding autoimmune risk variant defines role for ICOS in T peripheral helper cell development

Fine-mapping and functional studies implicate rs117701653, a non-coding single nucleotide polymorphism in the CD28/CTLA4/ICOS locus, as a risk variant for rheumatoid arthritis and type 1 diabetes. Here, using DNA pulldown, mass spectrometry, genome editing and eQTL analysis, we establish that the disease-associated risk allele is functional, reducing affinity for the inhibitory chromosomal regulator SMCHD1 to enhance expression of inducible T-cell costimulator (ICOS) in memory CD4+ T cells from healthy donors. Higher ICOS expression is paralleled by an increase in circulating T peripheral helper (Tph) cells and, in rheumatoid arthritis patients, of blood and joint fluid Tph cells as well as circulating plasmablasts. Correspondingly, ICOS ligation and carriage of the rs117701653 risk allele accelerate T cell differentiation into CXCR5-PD-1high Tph cells producing IL-21 and CXCL13. Thus, mechanistic dissection of a functional non-coding variant in human autoimmunity discloses a previously undefined pathway through which ICOS regulates Tph development and abundance.

A multi-ancestry genome-wide association study in type 1 diabetes.

Type 1 diabetes (T1D) is an autoimmune disease caused by destruction of the pancreatic β-cells. Genome-wide association (GWAS) and fine mapping studies have been conducted mainly in European ancestry (EUR) populations. We performed a multi-ancestry GWAS to identify SNPs and HLA alleles associated with T1D risk and age at onset. EUR families (N = 3223), and unrelated individuals of African (AFR, N = 891) and admixed (Hispanic/Latino) ancestry (AMR, N = 308) were genotyped using the Illumina HumanCoreExome BeadArray, with imputation to the TOPMed reference panel. The Multi-Ethnic HLA reference panel was utilized to impute HLA alleles and amino acid residues. Logistic mixed models (T1D risk) and frailty models (age at onset) were used for analysis. In GWAS meta-analysis, seven loci were associated with T1D risk at genome-wide significance: PTPN22, HLA-DQA1, IL2RA, RNLS, INS, IKZF4-RPS26-ERBB3, and SH2B3, with four associated with T1D age at onset (PTPN22, HLA-DQB1, INS, and ERBB3). AFR and AMR meta-analysis revealed NRP1 as associated with T1D risk and age at onset, although NRP1 variants were not associated in EUR ancestry. In contrast, the PTPN22 variant was significantly associated with risk only in EUR ancestry. HLA alleles and haplotypes most significantly associated with T1D risk in AFR and AMR ancestry differed from that seen in EUR ancestry; in addition, the HLA-DRB1*08:02-DQA1*04:01-DQB1*04:02 haplotype was 'protective' in AMR while HLA-DRB1*08:01-DQA1*04:01-DQB1*04:02 haplotype was 'risk' in EUR ancestry, differing only at HLA-DRB1*08. These results suggest that much larger sample sizes in non-EUR populations are required to capture novel loci associated with T1D risk.

PostGWAS: A web server for deciphering the causality post the genome-wide association studies

While genome-wide association studies (GWAS) have unequivocally identified vast disease susceptibility variants, a majority of them are situated in non-coding regions and are in high linkage disequilibrium (LD). To pave the way of translating GWAS signals to clinical drug targets, it is essential to identify the underlying causal variants and further causal genes. To this end, a myriad of post-GWAS methods have been devised, each grounded in distinct principles including fine-mapping, co-localization, and transcriptome-wide association study (TWAS) techniques. Yet, no platform currently exists that seamlessly integrates these diverse post-GWAS methodologies. In this work, we present a user-friendly web server for post-GWAS analysis, that seamlessly integrates 9 distinct methods with 12 models, categorized by fine-mapping, colocalization, and TWAS. The server mainly helps users decipher the causality hindered by complex GWAS signals, including casual variants and casual genes, without the burden of computational skills and complex environment configuration, and provides a convenient platform for post-GWAS analysis, result visualization, facilitating the understanding and interpretation of the genome-wide association studies. The postGWAS server is available at http://g2g.biographml.com/.

An Update on the Genetics of IgA Nephropathy.

Immunoglobulin A (IgA) nephropathy (IgAN), the most common form of glomerulonephritis, is one of the leading causes of end-stage kidney disease (ESKD). It is widely believed that genetic factors play a significant role in the development of IgAN. Previous studies of IgAN have provided important insights to unravel the genetic architecture of IgAN and its potential pathogenic mechanisms. The genome-wide association studies (GWASs) together have identified over 30 risk loci for IgAN, which emphasizes the importance of IgA production and regulation in the pathogenesis of IgAN. Follow-up fine-mapping studies help to elucidate the candidate causal variant and the potential pathogenic molecular pathway and provide new potential therapeutic targets. With the rapid development of next-generation sequencing technologies, linkage studies based on whole-genome sequencing (WGS)/whole-exome sequencing (WES) also identify rare variants associated with IgAN, accounting for some of the missing heritability. The complexity of pathogenesis and phenotypic variability may be better understood by integrating genetics, epigenetics, and environment. We have compiled a review summarizing the latest advancements in genetic studies on IgAN. We similarly summarized relevant studies examining the involvement of epigenetics in the pathogenesis of IgAN. Future directions and challenges in this field are also proposed.

Estimating gene-level false discovery probability improves eQTL statistical fine-mapping precision.

Statistical fine-mapping prioritizes putative causal variants from a large number of candidate variants, and is widely used in expression quantitative loci (eQTLs) studies. In eQTL fine-mapping, the existence of causal variants for gene expression is not guaranteed, since the genetic heritability of gene expression explained by nearby (cis-) variants is limited. Here we introduce a refined fine-mapping algorithm, named Knockoff-Finemap combination (KFc). KFc estimates the probability that the causal variant(s) exist in the cis-window of a gene through construction of knockoff genotypes (i.e. a set of synthetic genotypes that resembles the original genotypes), and uses it to adjust the posterior inclusion probabilities (PIPs). Utilizing simulated gene expression data, we show that KFc results in calibrated PIP distribution with improved precision. When applied to gene expression data of 465 genotyped samples from the Japan COVID-19 Task Force (JCTF), KFc resulted in significant enrichment of a functional score as well as reporter assay hits in the top PIP bins. When combined with functional priors derived from an external fine-mapping study (GTEx), KFc resulted in a significantly higher proportion of hematopoietic trait putative causal variants in the top PIP bins. Our work presents improvements in the precision of a major fine-mapping algorithm.

Deciphering Genomic Regions and Putative Candidate Genes for Grain Size and Shape Traits in Durum Wheat through GWAS

Durum wheat is an economically and nutritionally important cereal. The increase in durum wheat yield is mostly associated with improving grain traits, and the grain size- and shape-related traits are directly related to wheat yield. In addition, grain size influences the seed germination rate and seedling vigor, which play key roles in stand establishment and yield. Thus, it is important to investigate grain traits both agro-morphologically and genetically. In this study, a panel of durum wheat, consisting of 146 genotypes, was evaluated for grain traits agro-morphologically and a genome-wide association study (GWAS) was conducted to dissect the genomic regions associated with these traits. As a result of the GWAS, a total of 41 marker-trait associations (MTAs) were identified on different chromosomes of durum wheat. Of these MTAs, only 11 were stable across environments. A BLAST search for the flanking sequences of every stable MTA in the Svevo genome identified 18 putative candidate genes directly associated with the seed traits of different plants, particularly wheat seeds. In conclusion, the annotation results and literature information provide strong evidence that the identified stable MTAs and their candidate genes may have important functions in the formation of wheat grain traits. After the validation of these MTAs with different fine-mapping and functional characterization studies, these loci may provide valuable information for geneticists and breeders to improve wheat yield.

Functional Enrichment Analysis Identifying Regulatory Information Associated with Human Fracture.

Dozens of loci associated with fracture have been identified by genome-wide association studies (GWASs). However, most of these variants are located in the noncoding regions including introns, long terminal repeats, and intergenic regions. Although combining regulation information helps to identify the causal SNPs and interpret the involvement of these variants in the etiology of human fracture, regulation information which was truly associated with fracture was unknown. A novel functional enrichment method GARFIELD (GWAS Analysis of Regulatory of Functional Information Enrichment with LD correction) was applied to identify fracture-associated regulation information, including transcript factor binding sites, expression quantitative trait loci (eQTLs), chromatin states, enhancer, promoter, dyadic, super enhancer and Epigenome marks. Fracture SNPs were significantly enriched in exon (Bonferroni correction, p value < 7.14 × 10-3) at two GWAS p value thresholds through GARFIELD. High level of fold-enrichment was observed in super enhancer of monocyte and the enhancer of chondrocyte (Bonferroni correction, p value < 4.45 × 10-3). eQTLs of 44 tissues/cells and 10 transcription factors (TFs) were identified to be associated with human fracture. These results provide new insight into the etiology of human fracture, which might increase the identification of the causal SNPs through the fine-mapping study combined with functional annotation, as well as polygenic risk score.

Evaluating approaches for constructing polygenic risk scores for prostate cancer in men of African and European ancestry

Genome-wide polygenic risk scores (GW-PRSs) have been reported to have better predictive ability than PRSs based on genome-wide significance thresholds across numerous traits. We compared the predictive ability of several GW-PRS approaches to a recently developed PRS of 269 established prostate cancer-risk variants from multi-ancestry GWASs and fine-mapping studies (PRS269). GW-PRS models were trained with a large and diverse prostate cancer GWAS of 107,247 cases and 127,006 controls that we previously used to develop the multi-ancestry PRS269. Resulting models were independently tested in 1,586 cases and 1,047 controls of African ancestry from the California Uganda Study and 8,046 cases and 191,825 controls of European ancestry from the UK Biobank and further validated in 13,643 cases and 210,214 controls of European ancestry and 6,353 cases and 53,362 controls of African ancestry from the Million Veteran Program. Inthe testing data, the best performing GW-PRS approach had AUCs of 0.656 (95% CI= 0.635-0.677) in African and 0.844 (95%CI=0.840-0.848) in European ancestry men and corresponding prostate cancer ORs of 1.83 (95% CI= 1.67-2.00) and 2.19 (95% CI= 2.14-2.25), respectively, for each SD unit increase in the GW-PRS. Compared to the GW-PRS, in African and European ancestry men, the PRS269 had larger or similar AUCs (AUC= 0.679, 95% CI= 0.659-0.700 and AUC= 0.845, 95% CI= 0.841-0.849, respectively) and comparable prostate cancer ORs (OR= 2.05, 95% CI= 1.87-2.26 and OR= 2.21, 95% CI= 2.16-2.26, respectively). Findings were similar in the validation studies. This investigation suggests that current GW-PRS approaches may not improve the ability to predict prostate cancer risk compared to the PRS269 developed from multi-ancestry GWASs and fine-mapping.

The eQTL colocalization and transcriptome-wide association study identify potentially causal genes responsible for economic traits in Simmental beef cattle

BackgroundA detailed understanding of genetic variants that affect beef merit helps maximize the efficiency of breeding for improved production merit in beef cattle. To prioritize the putative variants and genes, we ran a comprehensive genome-wide association studies (GWAS) analysis for 21 agronomic traits using imputed whole-genome variants in Simmental beef cattle. Then, we applied expression quantitative trait loci (eQTL) mapping between the genotype variants and transcriptome of three tissues (longissimus dorsi muscle, backfat, and liver) in 120 cattle.ResultsWe identified 1,580 association signals for 21 beef agronomic traits using GWAS. We then illuminated 854,498 cis-eQTLs for 6,017 genes and 46,970 trans-eQTLs for 1,903 genes in three tissues and built a synergistic network by integrating transcriptomics with agronomic traits. These cis-eQTLs were preferentially close to the transcription start site and enriched in functional regulatory regions. We observed an average of 43.5% improvement in cis-eQTL discovery using multi-tissue eQTL mapping. Fine-mapping analysis revealed that 111, 192, and 194 variants were most likely to be causative to regulate gene expression in backfat, liver, and muscle, respectively. The transcriptome-wide association studies identified 722 genes significantly associated with 11 agronomic traits. Via the colocalization and Mendelian randomization analyses, we found that eQTLs of several genes were associated with the GWAS signals of agronomic traits in three tissues, which included genes, such as NADSYN1, NDUFS3, LTF and KIFC2 in liver, GRAMD1C, TMTC2 and ZNF613 in backfat, as well as TIGAR, NDUFS3 and L3HYPDH in muscle that could serve as the candidate genes for economic traits.ConclusionsThe extensive atlas of GWAS, eQTL, fine-mapping, and transcriptome-wide association studies aid in the suggestion of potentially functional variants and genes in cattle agronomic traits and will be an invaluable source for genomics and breeding in beef cattle.

Structural and immunological differences in Plasmodium falciparum sexual stage transmission-blocking vaccines comprised of Pfs25-EPA nanoparticles

Development of a malaria vaccine that blocks transmission of different parasite stages to humans and mosquitoes is considered critical for elimination efforts. A vaccine using Pfs25, a protein on the surface of zygotes and ookinetes, is under investigation as a transmission-blocking vaccine (TBV) that would interrupt parasite passage from mosquitoes to humans. The most extensively studied Pfs25 TBVs use Pichia pastoris-produced recombinant forms of Pfs25, chemically conjugated to a recombinant carrier protein, ExoProtein A (EPA). The recombinant form of Pfs25 first used in humans was identified as Pfs25H, which contained a total of 14 heterologous amino acid residues located at the amino- and carboxyl-termini including a His6 affinity tag. A second recombinant Pfs25, identified as Pfs25M, was produced to remove the heterologous amino acid residues and conjugated to EPA (Pfs25M-EPA). Here, monomeric Pfs25M was characterized biochemically and biophysically for identity, purity, and integrity including protein structure to assess its comparability with Pfs25H. Although the biological activities of Pfs25H and Pfs25M, whether generated by monomeric forms or conjugated nanoparticles, appeared similar, fine-mapping studies with two transmission-blocking monoclonal antibodies detected structural and immunological differences. In addition, evaluation of antisera generated against conjugated Pfs25H or Pfs25M nanoparticles in nonhuman primates identified polyclonal IgG that recognized these structural differences.

Systematic fine-mapping and functional studies of prostate cancer risk variants

To date, genome-wide association studies (GWAS) have revealed over 200 genetic risk loci associated with prostate cancer; yet, true disease-causing variants remain elusive. Identification of causal variants and their targets from association signals is complicated by high linkage disequilibrium and limited availability of functional genomics data for specific tissue/cell types. Here, we integrated statistical fine-mapping and functional annotation from prostate-specific epigenomic profiles, 3D genome features, and quantitative trait loci data to distinguish causal variants from associations and identify target genes. Our fine-mapping analysis yielded 3,395 likely causal variants, and multiscale functional annotation linked them to 487 target genes. We prioritized rs10486567 as a genome-wide top-ranked SNP and predicted HOTTIP as its target. Deletion of the rs10486567-associated enhancer in prostate cancer cells decreased their capacity for invasive migration. HOTTIP overexpression in enhancer-KO cell lines rescued defective invasive migration. Furthermore, we found that rs10486567 regulates HOTTIP through allele-specific long-range chromatin interaction.

TMEM232 promotes the inflammatory response in atopic dermatitis via the nuclear factor-κB and signal transducer and activator of transcription 3 signalling pathways.

Our group previously found that the transmembrane protein 232 (TMEM232) gene was associated with atopic dermatitis (AD) by genome-wide association study and fine mapping study. However, its function is unclear so far. To investigate the roles and mechanisms of TMEM232 in AD. The expression of TMEM232 was investigated in skin lesions of patients with AD, the MC903-induced AD mouse model, human primary keratinocytes and immortalized human keratinocyte cell line (HaCaT) cells stimulated with different inflammatory factors. The role of TMEM232 in AD was analysed in HaCaT cells and Tmem232 knockout (Tmem232-/-) mice. Tmem232-specific small interfering RNA (siRNA) was used to evaluate its therapeutic potential in the AD mouse model. The expression of TMEM232 was significantly increased in skin lesions of patients with AD, the MC903-induced AD mouse model and human primary keratinocytes and HaCaT cells stimulated with different inflammatory factors compared with controls. In the presence of MC903, Tmem232-/- mice exhibited significantly reduced dermatitis severity, mast-cell infiltration in the back, and expression of T-helper (Th)1 and Th2-related inflammatory factors in skin tissue compared with wild-type mice. In vitro and in vivo experiments further showed that upregulation of TMEM232 in AD exacerbated the inflammation response through activating the pathway of nuclear factor-κB and signal transducer and activator of transcription (STAT) 3, and was regulated by the interleukin-4/STAT6 axis, which formed a self-amplifying loop. Finally, topical application of Tmem232 siRNA markedly ameliorated AD-like lesions in the AD model. This study is the first to outline the function of TMEM232. It is involved in regulating inflammation in AD and may be a potential target for AD treatment.

Bayesian multivariant fine mapping using the Laplace prior.

Currently, the only effect size prior that is routinely implemented in a Bayesian fine-mapping multi-single-nucleotide polymorphism (SNP) analysis is the Gaussian prior. Here, we show how the Laplace prior can be deployed in Bayesian multi-SNP fine mapping studies. We compare the ranking performance of the posterior inclusion probability (PIP) using a Laplace prior with the ranking performance of the corresponding Gaussian prior and FINEMAP. Our results indicate that, for the simulation scenarios we consider here, the Laplace prior can lead to higher PIPs than either the Gaussian prior or FINEMAP, particularly for moderately sized fine-mapping studies. The Laplace prior also appears to have better worst-case scenario properties. We reanalyse the iCOGS case-control data from the CASP8 region on Chromosome 2. Even though this study has a total sample size of nearly 90,000individuals, there are still some differences in the top few ranked SNPs if the Laplace prior is used rather than the Gaussian prior. R code to implement the Laplace (and Gaussian) prior is available at https://github.com/Kevin-walters/lapmapr.

Identification of genetic variants associated with diabetic kidney disease in multiple Korean cohorts via a genome-wide association study mega-analysis

BackgroundThe pathogenesis of diabetic kidney disease (DKD) is complex, involving metabolic and hemodynamic factors. Although DKD has been established as a heritable disorder and several genetic studies have been conducted, the identification of unique genetic variants for DKD is limited by its multiplex classification based on the phenotypes of diabetes mellitus (DM) and chronic kidney disease (CKD). Thus, we aimed to identify the genetic variants related to DKD that differentiate it from type 2 DM and CKD.MethodsWe conducted a large-scale genome-wide association study mega-analysis, combining Korean multi-cohorts using multinomial logistic regression. A total of 33,879 patients were classified into four groups—normal, DM without CKD, CKD without DM, and DKD—and were further analyzed to identify novel single-nucleotide polymorphisms (SNPs) associated with DKD. Additionally, fine-mapping analysis was conducted to investigate whether the variants of interest contribute to a trait. Conditional analyses adjusting for the effect of type 1 DM (T1D)-associated HLA variants were also performed to remove confounding factors of genetic association with T1D. Moreover, analysis of expression quantitative trait loci (eQTL) was performed using the Genotype-Tissue Expression project. Differentially expressed genes (DEGs) were analyzed using the Gene Expression Omnibus database (GSE30529). The significant eQTL DEGs were used to explore the predicted interaction networks using search tools for the retrieval of interacting genes and proteins.ResultsWe identified three novel SNPs [rs3128852 (P = 8.21×10−25), rs117744700 (P = 8.28×10−10), and rs28366355 (P = 2.04×10−8)] associated with DKD. Moreover, the fine-mapping study validated the causal relationship between rs3128852 and DKD. rs3128852 is an eQTL for TRIM27 in whole blood tissues and HLA-A in adipose-subcutaneous tissues. rs28366355 is an eQTL for HLA-group genes present in most tissues.ConclusionsWe successfully identified SNPs (rs3128852, rs117744700, and rs28366355) associated with DKD and verified the causal association between rs3128852 and DKD. According to the in silico analysis, TRIM27 and HLA-A can define DKD pathophysiology and are associated with immune response and autophagy. However, further research is necessary to understand the mechanism of immunity and autophagy in the pathophysiology of DKD and to prevent and treat DKD.

Molecular analysis of resistance gene locus to bacterial canker and wilting disease in tomato mutants

The tomato plant is one of the most widely produced vegetables in the world. However, there are several disease factors which limit tomato production. The Clavibacter michiganensis subsp. michiganensis (Cmm) pathogen causes wilting and canker symptoms on the leaves of plants and bird’s eye symptoms on the tomato and then the plants completely die. Development of resistant tomato varieties is prerequisite due to absence of an effective control methods against the bacterial disease. The resistant M3-9 and M3-15 tomato plants have been developed because of mutation from susceptible NCEBR3 tomato seeds with ethyl methanosulfonate (EMS). For each chromosome of the tomato genome, 24 SSR markers were selected from each end of the haploid, 12 tomato chromosomes, and polymorphic differences between susceptible and resistant tomato plants were studied. Polymorphisms were found with SSR13 and SSRB18031 markers located on chromosome 5 with resistant mutants, M3-9, M3-15 and susceptible original NCEBR3 plants. It is envisaged that a resistance gene is located on the 5th chromosome of resistant M3-9 and M3-15 plants. Further fine mapping studies will reveal the location of the resistance gene(s) for controlling bacterial canker and wilting pathogen.

A Pleiotropic Flowering Time QTL Exhibits Gene-by-Environment Interaction for Fitness in a Perennial Grass.

Appropriate flowering time is a crucial adaptation impacting fitness in natural plant populations. Although the genetic basis of flowering variation has been extensively studied, its mechanisms in nonmodel organisms and its adaptive value in the field are still poorly understood. Here, we report new insights into the genetic basis of flowering time and its effect on fitness in Panicum hallii, a native perennial grass. Genetic mapping in populations derived from inland and coastal ecotypes identified flowering time quantitative trait loci (QTL) and many exhibited extensive QTL-by-environment interactions. Patterns of segregation within recombinant hybrids provide strong support for directional selection driving ecotypic divergence in flowering time. A major QTL on chromosome 5 (q-FT5) was detected in all experiments. Fine-mapping and expression studies identified a gene with orthology to a rice FLOWERING LOCUS T-like 9 (PhFTL9) as the candidate underlying q-FT5. We used a reciprocal transplant experiment to test for local adaptation and the specific impact of q-FT5 on performance. We did not observe local adaptation in terms of fitness tradeoffs when contrasting ecotypes in home versus away habitats. However, we observed that the coastal allele of q-FT5 conferred a fitness advantage only in its local habitat but not at the inland site. Sequence analyses identified an excess of low-frequency polymorphisms at the PhFTL9 promoter in the inland lineage, suggesting a role for either selection or population expansion on promoter evolution. Together, our findings demonstrate the genetic basis of flowering variation in a perennial grass and provide evidence for conditional neutrality underlying flowering time divergence.