Quantitative Trait Locus Association Research Articles

Concurrently mapping quantitative trait loci associations from multiple subspecies within hybrid populations

Many of the world’s agriculturally important plant and animal populations consist of hybrids of subspecies. Cattle in tropical and sub-tropical regions for example, originate from two subspecies, Bos taurus indicus (Bos indicus) and Bos taurus taurus (Bos taurus). Methods to derive the underlying genetic architecture for these two subspecies are essential to develop accurate genomic predictions in these hybrid populations. We propose a novel method to achieve this. First, we use haplotypes to assign SNP alleles to ancestral subspecies of origin in a multi-breed and multi-subspecies population. Then we use a BayesR framework to allow SNP alleles originating from the different subspecies differing effects. Applying this method in a composite population of B. indicus and B. taurus hybrids, our results show that there are underlying genomic differences between the two subspecies, and these effects are not identified in multi-breed genomic evaluations that do not account for subspecies of origin effects. The method slightly improved the accuracy of genomic prediction. More significantly, by allocating SNP alleles to ancestral subspecies of origin, we were able to identify four SNP with high posterior probabilities of inclusion that have not been previously associated with cattle fertility and were close to genes associated with fertility in other species. These results show that haplotypes can be used to trace subspecies of origin through the genome of this hybrid population and, in conjunction with our novel Bayesian analysis, subspecies SNP allele allocation can be used to increase the accuracy of QTL association mapping in genetically diverse populations.

APOA2: New Target for Molecular Hydrogen Therapy in Sepsis-Related Lung Injury Based on Proteomic and Genomic Analysis.

Target biomarkers for H2 at both the protein and genome levels are still unclear. In this study, quantitative proteomics acquired from a mouse model were first analyzed. At the same time, functional pathway analysis helped identify functional pathways at the protein level. Then, bioinformatics on mRNA sequencing data were conducted between sepsis and normal mouse models. Differential expressional genes with the closest relationship to disease status and development were identified through module correlation analysis. Then, common biomarkers in proteomics and transcriptomics were extracted as target biomarkers. Through analyzing expression quantitative trait locus (eQTL) and genome-wide association studies (GWAS), colocalization analysis on Apoa2 and sepsis phenotype was conducted by summary-data-based Mendelian randomization (SMR). Then, two-sample and drug-target, syndrome Mendelian randomization (MR) analyses were all conducted using the Twosample R package. For protein level, protein quantitative trait loci (pQTLs) of the target biomarker were also included in MR. Animal experiments helped validate these results. As a result, Apoa2 protein or mRNA was identified as a target biomarker for H2 with a protective, causal relationship with sepsis. HDL and type 2 diabetes were proven to possess causal relationships with sepsis. The agitation and inhibition of Apoa2 were indicated to influence sepsis and related syndromes. In conclusion, we first proposed Apoa2 as a target for H2 treatment.

Molecular and functional dissection of LIGULELESS1 (LG1) in plants.

Plant architecture is a culmination of the features necessary for capturing light energy and adapting to the environment. An ideal architecture can promote an increase in planting density, light penetration to the lower canopy, airflow as well as heat distribution to achieve an increase in crop yield. A number of plant architecture-related genes have been identified by map cloning, quantitative trait locus (QTL) and genome-wide association study (GWAS) analysis. LIGULELESS1 (LG1) belongs to the squamosa promoter-binding protein (SBP) family of transcription factors (TFs) that are key regulators for plant growth and development, especially leaf angle (LA) and flower development. The DRL1/2-LG1-RAVL pathway is involved in brassinosteroid (BR) signaling to regulate the LA in maize, which has facilitated the regulation of plant architecture. Therefore, exploring the gene regulatory functions of LG1, especially its relationship with LA genes, can help achieve the precise regulation of plant phenotypes adapted to varied environments, thereby increasing the yield. This review comprehensively summarizes the advances in LG1 research, including its effect on LA and flower development. Finally, we discuss the current challenges and future research goals associate with LG1.

Genome-Wide Association Study Points to Novel Locus for Gilles de la Tourette Syndrome

BackgroundTourette syndrome (TS) is a childhood-onset neurodevelopmental disorder of complex genetic architecture and is characterized by multiple motor tics and at least one vocal tic persisting for more than 1 year. MethodsWe performed a genome-wide meta-analysis integrating a novel TS cohort with previously published data, resulting in a sample size of 6133 individuals with TS and 13,565 ancestry-matched control participants. ResultsWe identified a genome-wide significant locus on chromosome 5q15. Integration of expression quantitative trait locus, Hi-C (high-throughput chromosome conformation capture), and genome-wide association study data implicated the NR2F1 gene and associated long noncoding RNAs within the 5q15 locus. Heritability partitioning identified statistically significant enrichment in brain tissue histone marks, while polygenic risk scoring of brain volume data identified statistically significant associations with right and left thalamus volumes and right putamen volume. ConclusionsOur work presents novel insights into the neurobiology of TS, thereby opening up new directions for future studies.

Combined linkage mapping and association analysis uncovers candidate genes for 25 leaf-related traits across three environments in maize.

Combined linkage and association analysis revealed five co-localized genetic loci across multiple environments. The key gene Zm00001d026491 was further verified to influence leaf length by candidate gene association analysis. Leaf morphology and number determine the canopy structure and thus affect crop yield. Herein, the genetic basis and key genes for 25 leaf-related traits, including leaf lengths (LL), leaf widths (LW), and leaf areas (LA) of eight continuous leaves under the tassel, and the number of leaves above the primary ear (LAE), were dissected by using an association panel and a biparental population. Using an intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, 290 quantitative trait loci (QTL) controlling these traits were detected across different locations, among which 115 QTL were individually repeatedly identified in at least two environments. Using the association panel, 165 unique significant single-nucleotide polymorphisms (SNPs) were associated with target traits (P < 2.15E-06), of which 35 were separately detected across multiple environments. In total, 42 pleiotropic QTL/SNPs (pQTL/SNPs) were responsible for at least two of the LL, LW, LA, and LAE traits across multiple environments. Combining the QTL mapping and association study, five unique SNPs were located within the confidence intervals of seven QTL, and 77 genes were identified based on the linkage disequilibrium regions of co-localized SNP loci. Gene-based association studies verified that the intragenic variants in the candidate gene Zm00001d026491 influenced LL of the third leaf counted from the top node. These findings will provide vital information to understanding the genetic basis of leaf-related traits and help to cultivate maize varieties with ideal plant architecture.

Epigenome-wide association study of serum folate in maternal peripheral blood leukocytes.

Aim: To perform an epigenome-wide association study (EWAS) of serum folate in maternal blood. Methods: Cross-ancestry (Europeans=302, South Asians=161) and ancestry-specific EWAS in the EPIPREG cohort were performed, followed by methyl quantitative trait loci analysis and association with cardiometabolic phenotypes. Replication was attempted using maternal folate intake and blood methylation data from theMoBa study and verified if the findings were significant in a previous EWAS of maternal serum folate in cord blood. Results & conclusion: cg19888088 (cross-ancestry) in EBF3, cg01952260 (Europeans) and cg07077240 (South Asians) in HERC3 were associated with serum folate. cg19888088 and cg01952260 were associated with diastolic blood pressure. cg07077240 was associated with variants in CASC15. The findings were not replicated and were not significant in cord blood.

Runs of Homozygosity and Quantitative Trait Locus/Association for Semen Parameters in Selected Chinese and South African Beef Cattle.

Simple SummaryGenome-wide runs of homozygosity (ROH) are excellent in understanding population history, estimating genomic inbreeding, and deciphering genetic architecture of complex traits and diseases, as well as identifying genes of agro-economic traits. ROH are defined as continuous region on the chromosome where an individual is homozygous across the genome. This study identified distribution of ROH in the six selected beef cattle breeds, Chinese Simentaler (CSI), Belgian Blue (BEL), and South African Angus (ANG), Nguni (NGU), Bonsmara (BON), and Simentaler (SIM), using Bovine BeadChip markers. Furthermore, nine candidate genes, CDF9, MARCH1, WDR19, SLOICI, ST7, DOP1B, CFAF9, INHBA, and ADAMTS1, were suggested to be associated with semen QTL traits and reported moderate inbreeding in some breeds with high to low correlation inbreeding between breeds. The study findings will allow proper guidelines for breeder’s societies.In this study, runs of homozygosity (ROH) and quantitative trait locus/association (QTL) for semen parameters in selected Chinese and South African beef cattle breed were estimated. The computed results showed 7516 ROH were observed between classes 0–5 Mb with no ROH observed in classes >40 Mb. Distribution of ROH showed high level of genomic coverage for ANG, NGU, CSI, and BEL breeds. Approximately 13 genomic regions with QTL were controlling sperm motility, sperm concentration, semen volume, sperm count, sperm head abnormalities, sperm tail abnormalities, sperm integrity, and percentage of abnormal sperm traits. Nine candidate genes, CDF9, MARCH1, WDR19, SLOICI, ST7, DOP1B, CFAF9, INHBA, and ADAMTS1, were suggested to be associated with above mentioned QTL traits. The results for inbreeding coefficient showed moderate correlation between FROH vs FHOM at 0.603 and high correlation between FROH 0–5 Mb 0.929, and lowest correlation for 0–>40 Mb 0.400. This study suggested recent inbreeding in CSI, BEL, ANG, BON, SIM, and NGU breeds. Furthermore, it highlighted varied inbreeding levels and identified QTL for semen traits and genes of association. These results can assist in implementation of genetic improvement strategies for bulls and provide awareness and proper guidelines in developing breeding programs.

Golgin A7 family member B (GOLGA7B) is a plausible novel gene associating high glycaemic index diet with acne vulgaris.

While the IGF1/FoxO1/mTORC1signalling pathway is a well-established nutrigenomic link between high glycaemic index (GI)/glycaemic load (GL) diet and acne vulgaris, other signalling pathways remain elusive. Therefore, we aimed to investigate other genes that are involved in the high GI/GL diet-acne link, using our Singapore/Malaysia population epidemiological, genomics and transcriptomics data. High GI/GL dietary habit of 3207 acne cases (1869 and 1341 further classified into severity and scarring grades, respectively) and 2521 controls were evaluated based on Quality of Diet based on Glycaemic Index Score (QDGIS). Overlapping concordant differentially expressed genes (DEGs) between acne case-controls and QDGIS poor-moderate/good classes were identified from whole-transcriptome sequencing data of PBMC of a subset of participants. Finally, we assessed the expression quantitative trait loci (eQTL) association of single nucleotide polymorphisms (SNPs) of the concordant DEGs. Daily intake of fruits significantly reduced the risk of acne presentation, severity and scarring by up to 48.5%. Those with good QDGIS had significantly lower risk of mild and moderate/severe acne, and grade 1/2 acne scarring. Sequential filtering identified four overlapping concordant DEGs that were significantly associated with acne and QDGIS, namely GOLGA7B, SNCB, LOC102723849 and LOC283683. Combining transcriptome and genetic association data, we identified intronic SNP rs1953947 in GOLGA7B as an eQTL for acne. In conclusion, we identified GOLGA7B as a plausible novel gene that links high GI/GL with acne, and hence propose a model for the involvement of Golga7b in high GI/GL diet-acne pathogenesis, which includes palmitoyl acyltransferase zDHHC5, fatty acid translocase CD36 and palmitic acid.

Metagenomic and metabolomic remodeling in nonagenarians and centenarians and its association with genetic and socioeconomic factors.

A better understanding of the biological and environmental variables that contribute to exceptional longevity has the potential to inform the treatment of geriatric diseases and help achieve healthy aging. Here, we compared the gut microbiome and blood metabolome of extremely long-lived individuals (94-105 years old) to that of their children (50-79 years old) in 116 Han Chinese families. We found extensive metagenomic and metabolomic remodeling in advanced age and observed a generational divergence in the correlations with socioeconomic factors. An analysis of quantitative trait loci revealed that genetic associations with metagenomic and metabolomic features were largely generation-specific, but we also found 131 plasma metabolic quantitative trait loci associations that were cross-generational with the genetic variants concentrated in six loci. These included associations between FADS1/2 and arachidonate, PTPA and succinylcarnitine and FLVCR1 and choline. Our characterization of the extensive metagenomic and metabolomic remodeling that occurs in people reaching extreme ages may offer new targets for aging-related interventions.

QtlXplorer: an online systems genetics browser in the Eucalyptus Genome Integrative Explorer (EucGenIE)

BackgroundAffordable high-throughput DNA and RNA sequencing technologies are allowing genomic analysis of plant and animal populations and as a result empowering new systems genetics approaches to study complex traits. The availability of intuitive tools to browse and analyze the resulting large-scale genetic and genomic datasets remain a significant challenge. Furthermore, these integrative genomics approaches require innovative methods to dissect the flow and interconnectedness of biological information underlying complex trait variation. The Plant Genome Integrative Explorer (PlantGenIE.org) is a multi-species database and domain that houses online tools for model and woody plant species including Eucalyptus. Since the Eucalyptus Genome Integrative Explorer (EucGenIE) is integrated within PlantGenIE, it shares genome and expression analysis tools previously implemented within the various subdomains (ConGenIE, PopGenIE and AtGenIE). Despite the success in setting up integrative genomics databases, online tools for systems genetics modelling and high-resolution dissection of complex trait variation in plant populations have been lacking.ResultsWe have developed qtlXplorer (https://eucgenie.org/QTLXplorer) for visualizing and exploring systems genetics data from genome-wide association studies including quantitative trait loci (QTLs) and expression-based QTL (eQTL) associations. This module allows users to, for example, find co-located QTLs and eQTLs using an interactive version of Circos, or explore underlying genes using JBrowse. It provides users with a means to build systems genetics models and generate hypotheses from large-scale population genomics data. We also substantially upgraded the EucGenIE resource and show how it enables users to combine genomics and systems genetics approaches to discover candidate genes involved in biotic stress responses and wood formation by focusing on two multigene families, laccases and peroxidases.ConclusionsqtlXplorer adds a new dimension, population genomics, to the EucGenIE and PlantGenIE environment. The resource will be of interest to researchers and molecular breeders working in Eucalyptus and other woody plant species. It provides an example of how systems genetics data can be integrated with functional genetics data to provide biological insight and formulate hypotheses. Importantly, integration within PlantGenIE enables novel comparative genomics analyses to be performed from population-scale data.

Genome-wide association and colocalization analyses identify target genes for brain arteriolosclerosis.

Brain arteriolosclerosis (B-ASC) is a small-vessel cerebrovascular pathologic change involving sclerotic thickening of arterioles in the brain. B-ASC pathology is common in aged autopsy cohorts, is associated with worse cognitive functioning, and is also associated with other neuropathologies. The genetics of dementia and other pathologic changes have been studied extensively, but to date no in-depth study has been conducted on genetic risk of autopsy-proven B-ASC. We performed the first genome-wide association study (GWAS) focusing on B-ASC using multiple cohorts. We then followed up on identified risk variants with functional analysis to investigate potential biological mechanisms. Individual GWAS and mega-analyses were conducted using data from participants in the National Alzheimer's Coordinating Center (NACC) (N = 3318) and the Religious Orders Study and Memory and Aging Project (ROSMAP) (N=1190). In total, 4,577,283 variants (MAF > 5%) were present in both datasets and passed quality control. Variants identified in GWAS were then checked for quantitative trait loci (QTL) associations in the Genotype Tissue Expression (GTEx) project. Colocalization analysis was then performed on identified QTL in brain and vascular tissues. Variants most highly associated with B-ASC were also examined for association with neuroimaging correlates of B-ASC in the Alzheimer's Disease Neuroimaging Initiative. One intergenic locus on Chromosome 6 was significantly associated with B-ASC in NACC (P = 1.4e-08) and colocalized with ELOVL4 gene expression in the brain. Other loci in both cohorts achieved less stringent thresholds of association and colocalized with multiple QTL in GTEx. One intergenic locus on Chromosome 10 near SORCS3 was suggestively associated with B-ASC in NACC and was replicated in ROSMAP. Genetic loci associated with B-ASC pathology were identified using multiple cohorts. Most loci associated with B-ASC in one cohort were not associated in the other. Multiple identified risk loci colocalized with gene expression or splicing QTL, providing evidence for potential biological mechanisms through which variants may influence disease risk. Due to the relatively small sample sizes and heterogeneity between the cohorts used, additional studies are merited to better characterize genetic risk of B-ASC.

Low Temperature Antioxidant Activity QTL Associate with Genomic Regions Involved in Physiological Cold Stress Tolerance Responses in Rice (Oryza sativa L.).

Boosting cold stress tolerance in crop plants can minimize stress-mediated yield losses. Asian rice (Oryza sativa L.), one of the most consumed cereal crops, originated from subtropical regions and is generally sensitive to low temperature environments. Within the two subspecies of rice, JAPONICA, and INDICA, the cold tolerance potential of its accessions is highly variable and depends on their genetic background. Yet, cold stress tolerance response mechanisms are complex and not well understood. This study utilized 370 accessions from the Rice Diversity Panel 1 (RDP1) to investigate and correlate four cold stress tolerance response phenotypes: membrane damage, seedling survivability, and catalase and anthocyanin antioxidative activity. Most JAPONICA accessions, and admixed accessions within JAPONICA, had lower membrane damage, higher antioxidative activity, and overall, higher seedling survivability compared to INDICA accessions. Genome-wide association study (GWAS) mapping was done using the four traits to find novel quantitative trait loci (QTL), and to validate and fine-map previously identified QTL. A total of 20 QTL associated to two or more traits were uncovered by our study. Gene Ontology (GO) term enrichment analyses satisfying four layers of filtering retrieved three potential pathways: signal transduction, maintenance of plasma membrane and cell wall integrity, and nucleic acids metabolism as general mechanisms of cold stress tolerance responses involving antioxidant activity.

Understanding the role of corneal biomechanics-associated genetic variants by bioinformatic analyses.

To analyze functions of corneal biomechanical properties (CBP)-related variants as corneal resistance factor (CRF) and corneal hysteresis (CH). Related single nucleotide polymorphisms (SNPs) and genes were identified from NHGRI-EBI GWAS catalog, GWASdb v2 and possible data in published studies. HaploReg v4.1 was used to find linkage SNPs. Functional annotations were performed by GWAVA, CADD and RegulomeDB. GTEx Portal database was used to find out expression quantitative trait locus (eQTL) association. Enrichr was used to annotate the function of GWAS gene and the associated signal pathway. STING (v11.0) database was utilized for protein interaction and network construction. The integration of 302 CH-associated and 420 CRF-associated lead SNPs has produced 531 CBP-associated lead SNPs. A total of 5,324 proxy variants identified using the HaploReg v4.1 and lead SNPs were functionally annotated. Based on the threshold (CADD ≥ 10, GWAVA ≥ 0.4 and RegulomeDB < rank 3), 23 prioritized putative regulatory SNPs were identified. Eight prioritized eQTL variants (rs75203695, rs34861673, rs846766, rs11024102, rs1377416, rs3829492, rs9934438 and rs197912) were found with strong potential of CBP regulation. It was indicated that CBP-associated genes were significantly enriched in extracellular matrix receptor interaction pathway, closely related to the phenotype of corneal dystrophy and keratoconus. COL1A1, SMAD3, BMP4 and RUNX2 occupied the core position in the co-expression network. Data integrative analysis can evaluate CBP variations and explore collagen and extracellular matrix pathways in CBP regulation, which is a promising tool to investigate biological process of corneal diseases.

A Missense Variant in the IL-6 Receptor and Protection From Peripheral Artery Disease

A Missense Variant in the IL-6 Receptor and Protection From Peripheral Artery Disease

Phenome-wide and expression quantitative trait locus associations of coronavirus disease 2019 genetic risk loci

SummaryWhile several genes and clinical traits have been associated with higher risk of severe coronavirus disease 2019 (COVID-19), how host genetic variants may interact with these parameters and contribute to severe disease is still unclear. Herein, we performed phenome-wide association study, tissue and immune-cell-specific expression quantitative trait locus (eQTL)/splicing quantitative trait locus, and colocalization analyses for genetic risk loci suggestively associated with severe COVID-19 with respiratory failure. Thirteen phenotypes/traits were associated with the severe COVID-19-associated loci at the genome-wide significance threshold, including monocyte counts, fat metabolism traits, and fibrotic idiopathic interstitial pneumonia. In addition, we identified tissue and immune subtype-specific eQTL associations affecting 48 genes, including several ones that may directly impact host immune responses, colocalized with the severe COVID-19 genome-wide association study associations, and showed altered expression in single-cell transcriptomes. Collectively, our work demonstrates that host genetic variations associated with multiple genes and traits show genetic pleiotropy with severe COVID-19 and may inform disease etiology.

Hypomethylation mediates genetic association with the major histocompatibility complex genes in Sjögren's syndrome.

Differential methylation of immune genes has been a consistent theme observed in Sjögren's syndrome (SS) in CD4+ T cells, CD19+ B cells, whole blood, and labial salivary glands (LSGs). Multiple studies have found associations supporting genetic control of DNA methylation in SS, which in the absence of reverse causation, has positive implications for the potential of epigenetic therapy. However, a formal study of the causal relationship between genetic variation, DNA methylation, and disease status is lacking. We performed a causal mediation analysis of DNA methylation as a mediator of nearby genetic association with SS using LSGs and genotype data collected from 131 female members of the Sjögren's International Collaborative Clinical Alliance registry, comprising of 64 SS cases and 67 non-cases. Bumphunter was used to first identify differentially-methylated regions (DMRs), then the causal inference test (CIT) was applied to identify DMRs mediating the association of nearby methylation quantitative trait loci (MeQTL) with SS. Bumphunter discovered 215 DMRs, with the majority located in the major histocompatibility complex (MHC) on chromosome 6p21.3. Consistent with previous findings, regions hypomethylated in SS cases were enriched for gene sets associated with immune processes. Using the CIT, we observed a total of 19 DMR-MeQTL pairs that exhibited strong evidence for a causal mediation relationship. Close to half of these DMRs reside in the MHC and their corresponding meQTLs are in the region spanning the HLA-DQA1, HLA-DQB1, and HLA-DQA2 loci. The risk of SS conferred by these corresponding MeQTLs in the MHC was further substantiated by previous genome-wide association study results, with modest evidence for independent effects. By validating the presence of causal mediation, our findings suggest both genetic and epigenetic factors contribute to disease susceptibility, and inform the development of targeted epigenetic modification as a therapeutic approach for SS.

Genetic architecture of embryo size and related traits in maize

The embryo in maize has a critical role in controlling kernel nutrition components and grain yield. We measured five embryo weight and size traits, six kernel weight and size traits, and five embryo-to-kernel ratio traits in a nested association mapping (NAM) population of 611 recombinant inbred lines (RILs) derived from four inbred lines including the high-oil, giant-embryo line BY815 as the common parent. Using three statistical methods, we identified 5–22 quantitative trait loci (QTL) for each trait, explaining 4.7%–46.7% of the phenotypic variation. The genetic architecture of maize embryo size and its related traits appeared to be dominated by multiple small-effect loci with little epistasis, and the genetic control underlying embryo size appeared to be distinct from that underlying kernel size. A trait–QTL association network included 205 nodes and 439 edges and revealed 28 key loci associated with at least three traits. Cloned maize genes including ZmUrb2, Emp12 and ZmBAM1d, maize orthologs of known rice genes that control seed size including BG1, XIAO and GS9, and 11 maize orthologs of Arabidopsis EMBRYO-DEFECTIVE (EMB) genes were identified as underlying these key loci. Further, the phenotypic and genetic relationships between embryo size and kernel size were evaluated, and genetic patterns for identified loci that control embryo size and its related traits were proposed. Our findings reveal distinct genetic architectures for embryo size, kernel size, and embryo-to-kernel ratio traits and establish a foundation for the improvement of embryo-size-mediated kernel nutrition and grain yield.

Genome survey and high-resolution genetic map provide valuable genetic resources for Fenneropenaeus chinensis

Fenneropenaeus chinensis is one of the most important aquaculture species in China. Research on its genomic and genetic structure not only helps us comprehend the genetic basis of complex economic traits, but also offers theoretical guidance in selective breeding. In the present study, a genome survey sequencing was performed to generate a rough reference genome utilized for groping preliminary genome characteristics and facilitate linkage and quantitative trait locus (QTL) mapping. Linkage mapping was conducted using a reduced-representation sequencing method 2b-RAD. In total, 36,762 SNPs were genotyped from 273 progenies in a mapping family, and a high-resolution linkage map was constructed. The consensus map contained 12,884 markers and spanned 5257.81 cM with an average marker interval of 0.41 cM, which was the first high-resolution genetic map in F. chinensis to our knowledge. QTL mapping and association analysis were carried out in 29 characters including body size, sex and disease resistance. 87 significant QTLs were detected in several traits and they were also evaluated by association analysis. Results of this study provide us valuable suggestions in genetic improvement and breeding of new varieties and also lay a basic foundation for further application of cloning of economic genes in selective breeding program and marker-assisted selection.

ERF4 affects fruit firmness through TPL4 by reducing ethylene production.

The firmness of fleshy fruit crops has a significant effect on their quality, consumer preference, shelf life and transportability. In a combined quantitative trait locus and genome-wide association studies study of apple fruit texture, we identified a mutation (C-G) in the ethylene response factor-associated amphiphilic repression (EAR) motif in the coding region of the apple ETHYLENE RESPONSE FACTOR4 (ERF4) gene. Chromatin immunoprecipitation sequencing showed that ERF4 binds to the promoter of ERF3, which is involved in regulation of ethylene biosynthesis. The EAR mutation in ERF4 results in reduced repression of ERF3 expression, which is turn promotes ethylene production and loss of fruit firmness. ERF4 acts as a transcriptional repressor whose activity is modulated by a TOPLESS co-repressor 4 (TPL4)-binding EAR repression motif. Biolayer interferometry analysis showed that the mutation in the EAR motif causes a reduction in the interaction with TPL4. Suppression of ERF4 or TPL4 promoted fruit ripening and ethylene production. Taken together, our results provide insights into how ERF4 allelic variation underlies an important fruit quality trait.

Summary-data-based Mendelian randomization prioritizes potential druggable targets for multiple sclerosis.

Multiple sclerosis is a complex autoimmune disease caused by a combination of genetic and environmental factors. Translation of Genome-Wide Association Study findings into therapeutics and effective preventive strategies has been limited to date. We used summary-data-based Mendelian randomization to synthesize findings from public expression quantitative trait locus, methylation quantitative trait locus and Multiple Sclerosis Genome-Wide Association Study datasets. By correlating the effects of methylation on multiple sclerosis, methylation on expression and expression on multiple sclerosis susceptibility, we prioritize genetic loci with evidence of influencing multiple sclerosis susceptibility. We overlay these findings onto a list of ‘druggable’ genes, i.e. genes which are currently, or could theoretically, be targeted by therapeutic compounds. We use GeNets and search tool for the retrieval of interacting genes/proteins to identify protein–protein interactions and druggable pathways enriched in our results. We extend these findings to a model of Epstein-Barr virus-infected B cells, lymphoblastoid cell lines. We conducted a systematic review of prioritized genes using the Open Targets platform to identify completed and planned trials targeting prioritized genes in multiple sclerosis and related disease areas. Expression of 45 genes in peripheral blood was strongly associated with multiple sclerosis susceptibility (False discovery rate 0.05). Of these 45 genes, 20 encode a protein which is currently targeted by an existing therapeutic compound. These genes were enriched for Gene Ontology terms pertaining to immune system function and leucocyte signalling. We refined this prioritized gene list by restricting to loci where CpG site methylation was associated with multiple sclerosis susceptibility, with gene expression and where expression was associated with multiple sclerosis susceptibility. This approach yielded a list of 15 prioritized druggable target genes for which there was evidence of a pathway linking methylation, expression and multiple sclerosis. Five of these 15 genes are targeted by existing drugs and three were replicated in a smaller expression Quantitative Trait Loci dataset (CD40, MERTK and PARP1). In lymphoblastoid cell lines, this approach prioritized 7 druggable gene targets, of which only one was prioritized by the multi-omic approach in peripheral blood (FCRL3). Systematic review of Open Targets revealed multiple early-phase trials targeting 13/20 prioritized genes in disorders related to multiple sclerosis. We use public datasets and summary-data-based Mendelian randomization to identify a list of prioritized druggable genetic targets in multiple sclerosis. We hope our findings could be translated into a platform for developing targeted preventive therapies.