Quantitative Trait Loci Research Articles

Population-wide DNA methylation polymorphisms at single-nucleotide resolution in 207 cotton accessions reveal epigenomic contributions to complex traits.

DNA methylation plays multiple regulatory roles in crop development. However, the relationships of methylation polymorphisms with genetic polymorphisms, gene expression, and phenotypic variation in natural crop populations remain largely unknown. Here, we surveyed high-quality methylomes, transcriptomes, and genomes obtained from the 20-days-post-anthesis (DPA) cotton fibers of 207 accessions and extended the classical framework of population genetics to epigenetics. Over 287 million single methylation polymorphisms (SMPs) were identified, 100 times more than the number of single nucleotide polymorphisms (SNPs). These SMPs were significantly enriched in intragenic regions while depleted in transposable elements. Association analysis further identified a total of 5,426,782 cis-methylation quantitative trait loci (cis-meQTLs), 5078 cis-expression quantitative trait methylation (cis-eQTMs), and 9157 expression quantitative trait loci (eQTLs). Notably, 36.39% of cis-eQTM genes were not associated with genetic variation, indicating that a large number of SMPs associated with gene expression variation are independent of SNPs. In addition, out of the 1715 epigenetic loci associated with yield and fiber quality traits, only 36 (2.10%) were shared with genome-wide association study (GWAS) loci. The construction of multi-omics regulatory networks revealed 43 cis-eQTM genes potentially involved in fiber development, which cannot be identified by GWAS alone. Among these genes, the role of one encoding CBL-interacting protein kinase 10 in fiber length regulation was successfully validated through gene editing. Taken together, our findings prove that DNA methylation data can serve as an additional resource for breeding purposes and can offer opportunities to enhance and expedite the crop improvement process.

EQTL in diseased colon tissue identifies novel target genes associated with IBD.

Genome-wide association studies (GWAS) have identified over 300 loci associated with the inflammatory bowel diseases (IBD), but putative causal genes for most are unknown. We conducted the largest disease-focused expression quantitative trait loci (eQTL) analysis using colon tissue from 252 IBD patients to determine genetic effects on gene expression and potential contribution to IBD. Combined with two non-IBD colon eQTL studies, we identified 194 potential target genes for 108 GWAS loci. eQTL in IBD tissue were enriched for IBD GWAS loci colocalizations, provided novel evidence for IBD-associated genes such as ABO and TNFRSF14 , and identified additional target genes compared to non-IBD tissue eQTL. IBD-associated eQTL unique to diseased tissue had distinct regulatory and functional characteristics with increased effect sizes. Together, these highlight the importance of eQTL studies in diseased tissue for understanding functional consequences of genetic variants, and elucidating molecular mechanisms and regulation of key genes involved in IBD.

A regulatory variant rs9379874 in T1D risk region 6p22.2 affects BTN3A1 expression regulating T cell function.

Genome-wide association studies (GWAS) have identified that 6p22.2 region is associated with type 1 diabetes (T1D) risk in the Chinese Han population. This study aims to reveal associations between this risk region and T1D subgroups and related clinical features, and further identify causal variant(s) and target gene(s) in this region. 2608 T1D and 4814 healthy controls were recruited from East, Central, and South China. Baseline data and genotyping for rs4320356 were collected. The most likely causal variant and gene were identified by bioinformatics analysis, dual-luciferase reporter assays, expression quantitative trait loci (eQTL), and functional annotation of the non-coding region within the 6p22.2 region. The leading variant rs4320356 in the 6p22.2 region was associated with T1D risk in the Chinese and Europeans. However, this variant was not significantly associated with islet function or autoimmunity. In silico analysis suggested rs9379874 was the most potential causal variant for T1D risk among thymus, spleen, and T cells, overlapping with the enhancer-related histone mark in multiple T cell subsets. Dual luciferase reporter assay and eQTL showed that the T allele of rs9379874 increased BTN3A1 expression by binding to FOXA1. Public single-cell RNA sequencing analysis indicated that BTN3A1 was related to T-cell activation, ATP metabolism, and cytokine metabolism pathways, which might contribute to T1D development. This study indicates that a functional variant rs9379874 regulates BTN3A1 expression, expanding the genomic landscape of T1D risk and offering a potential target for developing novel therapies.

Enhancing rice aroma through innovative approaches

Aromatic rice is used extensively in many different cuisines around the world for its wonderful aroma and cooking qualities. Aromatic rice varieties such as Basmati and non-Basmati fragrant rice have gained popularity in both domestic and foreign markets, despite their origins being predominantly in Southeast Asia and the Indian subcontinent. The primary gene responsible for rice aroma is the fgr/Badh2/Os2-AP, situated on chromosome 8 and encodes betaine aldehyde dehydrogenase 2 (Badh2). Key aroma compounds are attributed to over 500 volatiles. The primary aromatic molecule in rice, 2-acetyl-1-pyrroline (2-AP), accumulates as a result of mutations in this gene and gives rice its distinctive scent. Aroma is not decided by single compound rather it is decided by volatile profile and also by environmental factors. The identification of Quantitative Trait Loci (QTLs) linked to fragrance features on different chromosomes has improved our comprehension of the genetic processes behind rice scent. Advances in genetic engineering, particularly CRISPR/Cas9 and TALEN have facilitated the manipulation of the Badh2 gene, enhancing aroma profiles in rice. Additionally, gene silencing and introgression techniques have also proven in increasing 2-AP content. The review explores the biochemical properties and advancement of aromatic rice, emphasizing its complex inheritance patterns and potential for breeding improvement.

Genetic variations in NER pathway gene polymorphisms and Wilms tumor risk: A six-center case-control study in East China.

The nucleotide excision repair (NER) system is one of the main ways to protect organisms from DNA damage caused by endogenous and exogenous carcinogens. NER deficiency increases genome mutations, chromosomal aberrations, and cancer viability. However, the genetic association between Wilms tumor and NER pathway gene polymorphisms needs to be further validated. We assessed the associations between 19 NER gene polymorphisms and Wilms tumor susceptibility in 416 cases and 936 controls from East China via the TaqMan method. We found that xeroderma pigmentosum group D (XPD) rs238406 and rs13181 significantly decreased the risk of Wilms tumor [adjusted odds ratio (OR) = 0.59, 95% confidence interval (CI) = 0.46-0.75, p <.0001; adjusted OR = 0.63, 95% CI = 0.44-0.89, p = .009, respectively]. Furthermore, the rs751402 and rs2296147 polymorphisms in the xeroderma pigmentosum group G (XPG) gene were significantly correlated with an increased risk for Wilms tumor (adjusted OR = 1.47, 95% CI = 1.03-2.09, p = .034; adjusted OR = 2.14, 95% CI = 1.29-3.56, p = .003, respectively). Expression quantitative trait loci (eQTL) analysis revealed that these four polymorphisms may affect the expression of genes that are adjacent to XPD and XPG. Our study provides evidence that XPD and XPG gene polymorphisms are associated with Wilms tumor risk. Nonetheless, these findings should be confirmed in a larger sample size.

Heat shock responsive genes in Brassicaceae: genome-wide identification, phylogeny, and evolutionary associations within and between genera.

Heat stress affects the growth and development of Brassicaceae crops. Plant breeders aim to mitigate the effects of heat stress by selecting for heat stress tolerance, but the genes responsible for heat stress in Brassicaceae remain largely unknown. During heat stress, heat shock proteins (HSPs) function as molecular chaperones to aid in protein folding, and heat shock transcription factors (HSFs) serve as transcriptional regulators of HSP expression. We identified 5002 heat shock related genes, including HSPs and HSFs, across 32 genomes in Brassicaceae. Among these, 3347 genes were duplicated, with segmented duplication primarily contributing to their expansion. We identified 466 physical gene clusters, including 240 homogenous clusters and 226 heterogeneous clusters, shedding light on the organization of heat shock related genes. Notably, 37 genes were co-located with published thermotolerance quantitative trait loci, which supports their functional role in conferring heat stress tolerance. This study provides a comprehensive resource for the identification of functional Brassicaceae heat shock related genes, elucidates their clustering and duplication patterns and establishes the genomic foundation for future heat tolerance research. We hypothesise that genetic variants in HSP and HSF genes in certain species have potential for improving heat stress tolerance in Brassicaceae crops.

Genome-wide association studies in a diverse strawberry collection unveil loci controlling agronomic and fruit quality traits.

Strawberries (Fragaria sp.) are cherished for their organoleptic properties and nutritional value. However, breeding new cultivars involves the simultaneous selection of many agronomic and fruit quality traits, including fruit firmness and extended postharvest life. The strawberry germplasm collection here studied exhibited extensive phenotypic variation in 26 agronomic and fruit quality traits across three consecutive seasons. Phenotypic correlations and principal component analysis revealed relationships among traits and accessions, emphasizing the impact of plant breeding on fruit weight and firmness to the detriment of sugar or vitamin C content. Genetic diversity analysis on 124 accessions using 44,408 markers denoted a population structure divided into six subpopulations still retaining considerable diversity. Genome-wide association studies for the 26 traits unveiled 121 significant marker-trait associations distributed across 95 quantitative trait loci (QTLs). Multiple associations were detected for fruit firmness, a key breeding target, including a prominent locus on chromosome 6A. The candidate gene FaPG1, controlling fruit softening and postharvest shelf life, was identified within this QTL region. Differential expression of FaPG1 confirmed its role as the primary contributor to natural variation in fruit firmness. A kompetitive allele-specific PCR assay based on the single nucleotide polymorphism (SNP) AX-184242253, associated with the 6A QTL, predicts a substantial increase in fruit firmness, validating its utility for marker-assisted selection. In essence, this comprehensive study provides insights into the phenotypic and genetic landscape of the strawberry collection and lays a robust foundation for propelling the development of superior strawberry cultivars through precision breeding.

The Tetratricopeptide repeat protein TaTPR-B1 regulates spike compactness in bread wheat.

Spike compactness (SC) is strongly associated with wheat (Triticum aestivum L.) grain yield. In this study, we conducted a quantitative trait locus (QTL) analysis using a doubled haploid (DH) population derived from a cross between two common wheat varieties with contrasting spike morphology, revealing 16 stable QTLs associated with SC. The effect of a major QTL, QSc.cau-6B.1, was validated in 231 F7 recombinant inbred lines (RILs) derived from the same cross as the DH population. Using two residual heterozygous lines (RHLs), we delimited QSc.cau-6B.1 to an approximately 0.5-Mbp physical interval containing four high-confidence genes. The tetratricopeptide repeat-TraesCS6B03G1214400 (TaTPR-B1) was the priority candidate gene according to sequence and expression variations between near-isogenic lines. Accordingly, TaTPR-B1 knockout in the common wheat variety 'CB037' significantly increased SC compared to the wild type (WT). Conversely, TaTPR-B1 overexpression in the common wheat variety 'Fielder' significantly decreased SC compared to the WT. Moreover, we developed a PCR-based marker targeting the 32-bp insertion/deletion (InDel) between the two TaTPR-B1 alleles, which could be practical and valuable in modern wheat breeding programs for diagnostic purposes. Collectively, these findings provide insight into the genetic basis of SC in common wheat and present a valuable target with a breeder-friendly diagnostic marker for gene pyramid breeding.

Salt tolerance candidate genes identified by QTL mapping, RNA-seq, and functional analysis in tilapia

Salt tolerance in fish is crucial for aquaculture as it enhances survival and productivity in varying salinity conditions, thus expanding the range of viable farming environments and improving economic sustainability. Through QTL mapping and GWAS in a hybrid F2 family of Mozambique and Nile tilapia, two large-effect QTL on chromosomes 11 and 18 were identified respectively. These two QTL explained a total of 39.9 % of the phenotypic variance. The identification of a QTL on LG11 suggests the presence of a previously unrecognized genetic factor contributing to salt tolerance. Comparative transcriptomic analysis of gill and kidney tissues between susceptible and tolerant tilapias highlights the importance of osmotic balance in regulation of salt tolerance in tilapia. Integration of QTL and RNA-seq data identified two candidate genes: acyl-coenzyme A thioesterase 5 (acot5) and sodium- and chloride-dependent taurine transporter (slc6a6) likely playing critical roles in such process. Functional analysis showed that over-expressing acot5 or slc6a6 in grouper kidney cells increased viability under salt stress by 4.46 % and 17.53 %, respectively. Subcellular localization revealed nuclear presence of ACOT5 and stress-induced nuclear translocation of SLC6A6. These findings highlight acot5 and slc6a6 as candidates for genetic manipulation and selection to enhance salt tolerance in tilapia, guiding genetic improvement efforts and promoting sustainable practices.

Implication of digestive functions and microbiota in the establishment of muscle glycogen differences between divergent lines for ultimate pH

Both the quality of chicken meat and the quality of chicks are influenced by the level of breast muscle glycogen reserves. In order to study the role of digestive metabolism in establishing this muscular phenotype, we compared two divergent chicken lines for the ultimate pH (pHu) of the breast meat, a proxy for glycogen reserves. Males aged 4 weeks had twice the breast muscle glycogen content in the pHu- line (low pHu) than in the pHu + line (high pHu). The increase in glycogen reserves (pHu-) was associated with a higher relative weight of the proventriculus and gizzard, as well as better apparent ileal digestibility of nitrogen and calcium. The diversity of the cecal microbiota was comparable, but three bacterial genera (Lachnospira, Lachnospiraceae UCG-010, Caproiciproducens) varied between the lines. The differences observed could lead to down-regulation of carbon fixation in prokaryotes and of the citrate cycle in the pHu + line. RNA-seq analysis of the jejunum, the major site of nutrient absorption, revealed 149 genes differentially expressed (DE) between the lines, including several genes linked to immunity, hormonal response and circadian rhythms that are less expressed in pHu + animals. Others involved in cell migration and proliferation, and more generally tissue morphogenesis, also differed between the lines. Among the DE genes, several co-localized with Quantitative Trait Loci (QTL) controlling pHu and selection signatures identified in the divergent lines, such as the gene coding for ghrelin, a hormone regulating appetite.

Discovery of key volatile compound and regulatory gene governing the ‘taro-like’ aroma in wax gourd (Benincasa hispida)

The distinctive ‘taro-like’ aroma in wax gourd has become a crucial flavor characteristic, significantly enhancing its market value. Despite its importance, the volatile compounds responsible for this aroma and their genetic regulatory mechanisms have not been sufficiently studied. Utilizing headspace-gas chromatography-quadrupole time of flight-mass spectrometry (HS-GC-QTOF-MS), this research examined the volatile compounds in wax gourd germplasms with and without the ‘taro-like’ aroma, identifying 2-acetyl-1-pyrroline (2-AP) as the primary volatile responsible for this aroma. To uncover the genetic basis of this trait, segregating populations and recombinant inbred lines (RILs) were developed from MY-1 x GX-71 germplasms, which exhibit significant differences in aroma expression. Genetic analysis revealed that a single recessive gene governs the ‘taro-like’ aroma trait in wax gourd. Through bulked segregant analysis sequencing (BSA-seq) and quantitative trait loci (QTL) mapping, the candidate region for this aroma gene was localized to a 285.29 kb interval on Chromosome 2. Further examination of candidate gene expression and sequence variations highlighted the BhBADH gene's critical role in the formation of the ‘taro-like’ aroma. Additionally, RNA-seq analysis revealed significant enrichment of differentially expressed genes (DEGs) associated with metabolic pathways involved in environmental information processing, which contribute to 2-AP synthesis in wax gourd with a ‘taro-like’ aroma. Our research provides a theoretical foundation for the molecular breeding to enhance the aroma quality of wax gourd.

Genetic dissection of Meloidogyne incognita resistance genes based on VIGS functional analysis in Cucumis metuliferus

The southern root-knot nematode, Meloidogyne incognita, is a highly serious plant parasitic nematode species that causes significant economic losses in various crops, including cucumber (Cucumis sativus L.). Currently, there are no commercial cultivars available with resistance to M. incognita in cucumber. However, the African horned melon (Cucumis metuliferus Naud.), a semi-wild relative of cucumber, has shown high resistance to M. incognita. In this study, we constructed an ultrahigh-density genetic linkage bin-map using low-coverage sequences from an F2 population generated through the cross between C. metuliferus inbred lines CM3 and CM27. Finally, we identified a QTL (quantitative trait locus, QTL3.1) with a LOD (logarithm of the odds) score of 3.84, explaining 8.4% of the resistance variation. Subsequently, by combining the results of qPCR (quantitative PCR) and VIGS (virus-induced gene silencing), we identified two genes, EVM0025394 and EVM0006042, that are potentially involved in the resistance to M. incognita in CM3. The identification of QTLs and candidate genes in this study serve as a basis for further functional analysis and lay the groundwork for harnessing this resistance trait.

Quantitative Trait Locus Mapping Combined with RNA Sequencing Identified Candidate Genes for Resistance to Powdery Mildew in Bitter Gourd (Momordica charantia L.)

Improving the powdery mildew resistance of bitter gourd is highly important for achieving high yield and high quality. To better understand the genetic basis of powdery mildew resistance in bitter gourd, this study analyzed 300 lines of recombinant inbred lines (RILs) formed by hybridizing the powdery mildew-resistant material MC18 and the powdery mildew-susceptible material MC402. A high-density genetic map of 1222.04 cM was constructed via incorporating 1,996,505 SNPs generated by resequencing data from 180 lines, and quantitative trait locus (QTL) positioning was performed using phenotypic data at different inoculation stages. A total of seven QTLs related to powdery mildew resistance were identified on four chromosomes, among which qPm-3-1 was detected multiple times and at multiple stages after inoculation. By selecting 18 KASP markers that were evenly distributed throughout the region, 250 lines and parents were genotyped, and the interval was narrowed to 207.22 kb, which explained 13.91% of the phenotypic variation. Through RNA-seq analysis of the parents, 11,868 differentially expressed genes (DEGs) were screened. By combining genetic analysis, gene coexpression, and sequence comparison analysis of extreme materials, two candidate genes controlling powdery mildew resistance in bitter gourd were identified (evm.TU.chr3.2934 (C3H) and evm.TU.chr3.2946 (F-box-LRR)). These results represent a step forward in understanding the genetic regulatory network of powdery mildew resistance in bitter gourd and lay a molecular foundation for the genetic improvement in powdery mildew resistance.

Parental assigned chromosomes for cultivated cacao provides insights into genetic architecture underlying resistance to vascular streak dieback.

Diseases of Theobroma cacao L. (Malvaceae) disrupt cocoa bean supply and economically impact growers. Vascular streak dieback (VSD), caused by Ceratobasidium theobromae, is a new encounter disease of cacao currently contained to southeast Asia and Melanesia. Resistance to VSD has been tested with large progeny trials in Sulawesi, Indonesia, and in Papua New Guinea with the identification of informative quantitative trait loci (QTLs). Using a VSD susceptible progeny tree (clone 26), derived from a resistant and susceptible parental cross, we assembled the genome to chromosome-level and discriminated alleles inherited from either resistant or susceptible parents. The parentally phased genomes were annotated for all predicted genes and then specifically for resistance genes of the nucleotide-binding site leucine-rich repeat class (NLR). On investigation, we determined the presence of NLR clusters and other potential disease response gene candidates in proximity to informative QTLs. We identified structural variants within NLRs inherited from parentals. We present the first diploid, fully scaffolded, and parentally phased genome resource for T. cacao L. and provide insights into the genetics underlying resistance and susceptibility to VSD.

Proteo-genomic analyses in relatively lean Chinese adults identify proteins and pathways that affect general and central adiposity levels

Adiposity is an established risk factor for multiple diseases, but the causal relationships of different adiposity types with circulating protein biomarkers have not been systematically investigated. We examine the causal associations of general and central adiposity with 2923 plasma proteins among 3977 Chinese adults (mean BMI = 23.9 kg/m²). Genetically-predicted body mass index (BMI), body fat percentage (BF%), waist circumference (WC), and waist-to-hip ratio (WHR) are significantly (FDR < 0.05) associated with 399, 239, 436, and 283 proteins, respectively, with 80 proteins associated with all four and 275 with only one adiposity trait. WHR is associated with the most proteins (n = 90) after adjusting for other adiposity traits. These associations are largely replicated in Europeans (mean BMI = 27.4 kg/m²). Two-sample Mendelian randomisation (MR) analyses in East Asians using cis-protein quantitative trait locus (cis-pQTLs) identified in GWAS find 30/2 proteins significantly affect levels of BMI/WC, respectively, with 10 showing evidence of colocalisation, and seven (inter-alpha-trypsin inhibitor heavy chain H3, complement factor B, EGF-containing fibulin-like extracellular matrix protein 1, thioredoxin domain-containing protein 15, alpha-2-antiplasmin, fibronectin, mimecan) are replicated in separate MR using different cis-pQTLs identified in Europeans. These findings identified potential novel mechanisms and targets, to our knowledge, for improved treatment and prevention of obesity and associated diseases.

HBI: a hierarchical Bayesian interaction model to estimate cell-type-specific methylation quantitative trait loci incorporating priors from cell-sorted bisulfite sequencing data

Methylation quantitative trait loci (meQTLs) quantify the effects of genetic variants on DNA methylation levels. However, most published studies utilize bulk methylation datasets composed of different cell types and limit our understanding of cell-type-specific methylation regulation. We propose a hierarchical Bayesian interaction (HBI) model to infer cell-type-specific meQTLs, which integrates a large-scale bulk methylation data and a small-scale cell-type-specific methylation data. Through simulations, we show that HBI enhances the estimation of cell-type-specific meQTLs. In real data analyses, we demonstrate that HBI can further improve the functional annotation of genetic variants and identify biologically relevant cell types for complex traits.

Breeding for flooding tolerance in rice: Advancements and future perspectives

Rice, the resilient grain, is cultivated in different types of ecosystems, from seacoast to hilly areas. Unfortunately, due to climate change, it frequently suffers from submergence stress during its growth period. Anoxic stress at the germination phase, flash flooding during the vegetative phase and water stagnation in low-lying areas are the major types of flooding in rice. When floodwaters rise, rice adapts. Some varieties stretch their stems toward the surface, gasping for air. Others remain dormant, conserving energy, like the FR13A landrace, SUB1 equips rice with underwater endurance. It orchestrates a genetic symphony, fine-tuning metabolic pathways and signalling survival. Some of the promising quantitative trait locus (QTLs) identified are qAG-9-2, which is responsible for anaerobic germination tolerance; qSUB1 for vegetative stage submergence tolerance, and qTIL12 for deep water adaptation. Identifying other novel QTLs and donors helps to breed varieties tolerant to different types of submergence stress. Along with Swarna SUB1, many mega-submergence-tolerant varieties have been developed and released for cultivation in Asia. As we cultivate these versatile survivors, we sow hope for a food-secure future.

A high-throughput protocol for testing heat-stress tolerance in pollen

AbstractViable pollen is crucial for fertilization, but pollen is generally highly susceptible to heat stress. A quick, reliable method for testing the heat-stress tolerance of pollen is needed to improve the heat-stress tolerance in plants, but current methods require considerable space and labor. In addition, many such methods only test tolerance to a single constant temperature, making it time-consuming to screen heat tolerance over a wide temperature range and to examine the dynamics of pollen viability at different temperatures. To address this issue, we aimed to: (1) develop an easy, reliable method for measuring pollen viability at different temperatures; and (2) identify the best temperature range for screening pollen with high heat-stress tolerance. We harvested mature pollen from wheat (Triticum aestivum) plants and transferred it to a 96-well plate filled with liquid medium containing sucrose. We placed the plate in a PCR machine operating under a gradient PCR program to simultaneously test a range of temperatures. After incubating the pollen for 4 h, at temperatures ranging from 21.9 to 47 °C, we examined the pollen grains under a light microscope and employed a specific image analysis pipeline to assess the effects of temperature on pollen morphology, germination, and tube growth. This method facilitated the high-throughput screening of many pollen samples, enabling rapid, reliable, and precise analysis of pollen viability in response to temperature. Our approach should be applicable to other plant species and could be used to identify quantitative trait loci or genes influencing heat stress tolerance in pollen for breeding programs.

Multiomics integration unravels genotype‐microbiome interactions shaping the conjunctival transcriptome

AbstractGene expression is a molecular phenotype shaped by the interplay between genotype and environment. The microbiome represents a critical environmental exposure for the host. However, the genotype‐microbiome interactions (GMIs) shaping the host transcriptome remain largely unexplored. Here, we integrated conjunctival multiomics data from 120 pairs of twins to investigate GMIs. We identified 272,972 expression quantitative trait loci associated with 5946 genes and 241,073 genotype‐controlled correlations between gene expression and microbial abundance. We developed a modeling approach, MicroGenix, that screens for GMIs from host genome, transcriptome, and microbiome data and identifies GMIs associated with the disease through gene‐based association tests. We applied MicroGenix to the conjunctival data set and found that incorporating GMIs into gene expression prediction models significantly increased prediction accuracy. The genes with increased accuracy were overrepresented by those encoding cell adhesion molecules. We further used MicroGenix to predict the transcriptome from conjunctival metagenomes and identified GMIs associated with ocular surface disease. Our work provides a resource for studying host‐microbiome interactions at the conjunctiva and a computational approach to investigate GMIs using multiomics data.

A network-based systems genetics framework identifies pathobiology and drug repurposing in Parkinson's disease.

Parkinson's disease (PD) is the second most prevalent neurodegenerative disorder. However, current treatments are directed at symptoms and lack ability to slow or prevent disease progression. Large-scale genome-wide association studies (GWAS) have identified numerous genomic loci associated with PD, which may guide the development of disease-modifying treatments. We presented a systems genetics approach to identify potential risk genes and repurposable drugs for PD. First, we leveraged non-coding GWAS loci effects on multiple human brain-specific quantitative trait loci (xQTLs) under the protein-protein interactome (PPI) network. We then prioritized a set of PD likely risk genes (pdRGs) by integrating five types of molecular xQTLs: expression (eQTLs), protein (pQTLs), splicing (sQTLs), methylation (meQTLs), and histone acetylation (haQTLs). We also integrated network proximity-based drug repurposing and patient electronic health record (EHR) data observations to propose potential drug candidates for PD treatments. We identified 175 pdRGs from QTL-regulated GWAS findings, such as SNCA , CTSB , LRRK2, DGKQ , CD38 and CD44 . Multi-omics data validation revealed that the identified pdRGs are likely to be druggable targets, differentially expressed in multiple cell types and impact both the parkin ubiquitin-proteasome and alpha-synuclein (a-syn) pathways. Based on the network proximity-based drug repurposing followed by EHR data validation, we identified usage of simvastatin as being significantly associated with reduced incidence of PD (fall outcome: hazard ratio (HR) = 0.91, 95% confidence interval (CI): 0.87-0.94; for dementia outcome: HR = 0.88, 95% CI: 0.86-0.89), after adjusting for 267 covariates. Our network-based systems genetics framework identifies potential risk genes and repurposable drugs for PD and other neurodegenerative diseases if broadly applied.