Expression Quantitative Trait Loci Research Articles

Integrating single cell expression quantitative trait loci summary statistics to understand complex trait risk genes

Transcriptome-wide association study (TWAS) is a popular approach to dissect the functional consequence of disease associated non-coding variants. Most existing TWAS use bulk tissues and may not have the resolution to reveal cell-type specific target genes. Single-cell expression quantitative trait loci (sc-eQTL) datasets are emerging. The largest bulk- and sc-eQTL datasets are most conveniently available as summary statistics, but have not been broadly utilized in TWAS. Here, we present a new method EXPRESSO (EXpression PREdiction with Summary Statistics Only), to analyze sc-eQTL summary statistics, which also integrates 3D genomic data and epigenomic annotation to prioritize causal variants. EXPRESSO substantially improves existing methods. We apply EXPRESSO to analyze multi-ancestry GWAS datasets for 14 autoimmune diseases. EXPRESSO uniquely identifies 958 novel gene x trait associations, which is 26% more than the second-best method. Among them, 492 are unique to cell type level analysis and missed by TWAS using whole blood. We also develop a cell type aware drug repurposing pipeline, which leverages EXPRESSO results to identify drug compounds that can reverse disease gene expressions in relevant cell types. Our results point to multiple drugs with therapeutic potentials, including metformin for type 1 diabetes, and vitamin K for ulcerative colitis.

O03 Genetic impacts on the skin methylome

Abstract Introduction and aims DNA methylation at CpG dinucleotides is a dynamic epigenetic mechanism that has been linked to numerous skin phenotypes and diseases. Previous studies in blood show that a substantial proportion of the human DNA methylome exhibits evidence for genetic influence, and that 15–17% of narrow-sense methylation heritability is explained by genetic variants, or methylation quantitative trait loci (meQTLs). However, no genome-wide meQTL studies have previously been performed in skin. The primary purpose of this work was to identify the genetic basis of DNA methylation levels in whole skin. Methods Utilizing a twin-based heritability model, we estimated DNA methylation heritability using 414 bulk skin tissue samples from healthy older female twins from the TwinsUK cohort, profiled using the Illumina Infinium HumanMethylation450 BeadChip (Illumina, San Diego, CA, USA). We then integrated genomic data for 394 of these twins to perform genome-wide meQTL analyses. Finally, we integrated these meQTL Results with skin eQTL Results from an overlapping sample (N = 664). Results Over 8.2% of measured CpG sites showed evidence for strong heritability (h2 > 0.4), with a mean additive heritability estimate of 11.07% across all tested sites. We identified 76 442 CpG sites (18.8% of measured CpG sites) associated with a meQTL in cis (< 1 MB variant-CpG distance) and 1438 (0.35%) in trans. Both CpG sites under meQTL effects, and genetic variants underlying these meQTL effects in skin, are enriched in regions previously identified from genome-wide and epigenome-wide association studies for skin-related phenotypes, including psoriasis, melanoma and vitiligo. Additionally, some meQTLs influence gene expression via mediation of expression quantitative trait loci effects in skin, including in genes linked to melanoma and psoriasis. Conclusions Our study is the first genome-wide DNA methylation heritability and meQTL analysis of whole skin, giving insights into the regulatory landscape of epigenomic variation in skin.

MR-GGI: accurate inference of gene–gene interactions using Mendelian randomization

BackgroundResearchers have long studied the regulatory processes of genes to uncover their functions. Gene regulatory network analysis is one of the popular approaches for understanding these processes, requiring accurate identification of interactions among the genes to establish the gene regulatory network. Advances in genome-wide association studies and expression quantitative trait loci studies have led to a wealth of genomic data, facilitating more accurate inference of gene–gene interactions. However, unknown confounding factors may influence these interactions, making their interpretation complicated. Mendelian randomization (MR) has emerged as a valuable tool for causal inference in genetics, addressing confounding effects by estimating causal relationships using instrumental variables. In this paper, we propose a new statistical method, MR-GGI, for accurately inferring gene–gene interactions using Mendelian randomization.ResultsMR-GGI applies one gene as the exposure and another as the outcome, using causal cis-single-nucleotide polymorphisms as instrumental variables in the inverse-variance weighted MR model. Through simulations, we have demonstrated MR-GGI's ability to control type 1 error and maintain statistical power despite confounding effects. MR-GGI performed the best when compared to other methods using the F1 score on the DREAM5 dataset. Additionally, when applied to yeast genomic data, MR-GGI successfully identified six clusters. Through gene ontology analysis, we have confirmed that each cluster in our study performs distinct functional roles by gathering genes with specific functions.ConclusionThese findings demonstrate that MR-GGI accurately inferences gene–gene interactions despite the confounding effects in real biological environments.

Deciphering the role of a SINE-VNTR-Alu retrotransposon polymorphism as a biomarker of Parkinson’s disease progression

SINE-VNTR-Alu (SVA) retrotransposons are transposable elements which represent a source of genetic variation. We previously demonstrated that the presence/absence of a human-specific SVA, termed SVA_67, correlated with the progression of Parkinson’s disease (PD). In the present study, we demonstrate that SVA_67 acts as expression quantitative trait loci, thereby exhibiting a strong regulatory effect across the genome using whole genome and transcriptomic data from the Parkinson’s progression markers initiative cohort. We further show that SVA_67 is polymorphic for its variable number tandem repeat domain which correlates with both regulatory properties in a luciferase reporter gene assay in vitro and differential expression of multiple genes in vivo. Additionally, this variation’s utility as a biomarker is reflected in a correlation with a number of PD progression markers. These experiments highlight the plethora of transcriptomic and phenotypic changes associated with SVA_67 polymorphism which should be considered when investigating the missing heritability of neurodegenerative diseases.

High-resolution African HLA resource uncovers HLA-DRB1 expression effects underlying vaccine response.

How human genetic variation contributes to vaccine effectiveness in infants is unclear, and data are limited on these relationships in populations with African ancestries. We undertook genetic analyses of vaccine antibody responses in infants from Uganda (n = 1391), Burkina Faso (n = 353) and South Africa (n = 755), identifying associations between human leukocyte antigen (HLA) and antibody response for five of eight tested antigens spanning pertussis, diphtheria and hepatitis B vaccines. In addition, through HLA typing 1,702 individuals from 11 populations of African ancestry derived predominantly from the 1000 Genomes Project, we constructed an imputation resource, fine-mapping class II HLA-DR and DQ associations explaining up to 10% of antibody response variance in our infant cohorts. We observed differences in the genetic architecture of pertussis antibody response between the cohorts with African ancestries and an independent cohort with European ancestry, but found no in silico evidence of differences in HLA peptide binding affinity or breadth. Using immune cell expression quantitative trait loci datasets derived from African-ancestry samples from the 1000 Genomes Project, we found evidence of differential HLA-DRB1 expression correlating with inferred protection from pertussis following vaccination. This work suggests that HLA-DRB1 expression may play a role in vaccine response and should be considered alongside peptide selection to improve vaccine design.

Genetic mapping across autoimmune diseases reveals shared associations and mechanisms.

Autoimmune and inflammatory diseases are polygenic disorders of the immune system. Many genomic loci harbor risk alleles for several diseases, but the limited resolution of genetic mapping prevents determining whether the same allele is responsible, indicating a shared underlying mechanism. Here, using a collection of 129,058 cases and controls across 6 diseases, we show that ~40% of overlapping associations are due to the same allele. We improve fine-mapping resolution for shared alleles twofold by combining cases and controls across diseases, allowing us to identify more expression quantitative trait loci driven by the shared alleles. The patterns indicate widespread sharing of pathogenic mechanisms but not a single global autoimmune mechanism. Our approach can be applied to any set of traits and is particularly valuable as sample collections become depleted.

Sex affects transcriptional associations with schizophrenia across the dorsolateral prefrontal cortex, hippocampus, and caudate nucleus

Schizophrenia is a complex neuropsychiatric disorder with sexually dimorphic features, including differential symptomatology, drug responsiveness, and male incidence rate. Prior large-scale transcriptome analyses for sex differences in schizophrenia have focused on the prefrontal cortex. Analyzing BrainSeq Consortium data (caudate nucleus: n = 399, dorsolateral prefrontal cortex: n = 377, and hippocampus: n = 394), we identified 831 unique genes that exhibit sex differences across brain regions, enriched for immune-related pathways. We observed X-chromosome dosage reduction in the hippocampus of male individuals with schizophrenia. Our sex interaction model revealed 148 junctions dysregulated in a sex-specific manner in schizophrenia. Sex-specific schizophrenia analysis identified dozens of differentially expressed genes, notably enriched in immune-related pathways. Finally, our sex-interacting expression quantitative trait loci analysis revealed 704 unique genes, nine associated with schizophrenia risk. These findings emphasize the importance of sex-informed analysis of sexually dimorphic traits, inform personalized therapeutic strategies in schizophrenia, and highlight the need for increased female samples for schizophrenia analyses.

Disentangling genetic effects on transcriptional and post-transcriptional gene regulation through integrating exon and intron expression QTLs

Expression quantitative trait loci (eQTL) studies typically consider exon expression of genes and discard intronic RNA sequencing reads despite their information on RNA metabolism. Here, we quantify genetic effects on exon and intron levels of genes and their ratio in lymphoblastoid cell lines, revealing thousands of cis-QTLs of each type. While genetic effects are often shared between cis-QTL types, 7814 (47%) are not detected as top cis-QTLs at exon levels. We show that exon levels preferentially capture genetic effects on transcriptional regulation, while exon-intron-ratios better detect those on co- and post-transcriptional processes. Considering all cis-QTL types substantially increases (by 71%) the number of colocalizing variants identified by genome-wide association studies (GWAS). It further allows dissecting the potential gene regulatory processes underlying GWAS associations, suggesting comparable contributions by transcriptional (50%) and co- and post-transcriptional regulation (46%) to complex traits. Overall, integrating intronic RNA sequencing reads in eQTL studies expands our understanding of genetic effects on gene regulatory processes.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle

BackgroundCharacterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22–150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle.ResultsWe identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits.ConclusionsThis study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

Novel putative causal mutations associated with fat traits in Nellore cattle uncovered by eQTLs located in open chromatin regions

Intramuscular fat (IMF) and backfat thickness (BFT) are critical economic traits impacting meat quality. However, the genetic variants controlling these traits need to be better understood. To advance knowledge in this area, we integrated RNA-seq and single nucleotide polymorphisms (SNPs) identified in genomic and transcriptomic data to generate a linkage disequilibrium filtered panel of 553,581 variants. Expression quantitative trait loci (eQTL) analysis revealed 36,916 cis-eQTLs and 14,408 trans-eQTLs. Association analysis resulted in three eQTLs associated with BFT and 24 with IMF. Functional enrichment analysis of genes regulated by these 27 eQTLs revealed noteworthy pathways that can play a fundamental role in lipid metabolism and fat deposition, such as immune response, cytoskeleton remodeling, iron transport, and phospholipid metabolism. We next used ATAC-Seq assay to identify and overlap eQTL and open chromatin regions. Six eQTLs were in regulatory regions, four in predicted insulators and possible CCCTC-binding factor DNA binding sites, one in an active enhancer region, and the last in a low signal region. Our results provided novel insights into the transcriptional regulation of IMF and BFT, unraveling putative regulatory variants.

Abstract PO5-09-03: Multi-tissue transcriptome-wide association studies identified genes for intrinsic subtypes of breast cancer

Abstract Background: Though several genome-wide association studies (GWAS) of breast cancer (BC) have identified common variants which differ between intrinsic subtypes, the genes through which these variants act through to impact BC risk have not been fully established. Furthermore, transcriptome-wide association studies (TWAS) have thus far been primarily employed to identify genes associated with overall BC risk, while overlooking how the influence of common variation on gene expression may contribute to subtype-specific differences. Methods: In this study, we performed two complementary multi-tissue TWASs for each of the following BC intrinsic subtypes: Luminal A-like, Luminal B-like, Luminal B/HER2-negative-like, HER2-enriched-like, Triple Negative BC. These two approaches included 1) an expression-based approach that collated TWAS signals from expression quantitative trait loci (eQTLs) across multiple tissues using the Aggregated Cauchy Association Test (ACAT) and 2) a splicing-based approach that utilized two applications of ACAT to collate splicing QTLs (sQTLs) for a given gene in a tissue and then across tissues. To perform these two TWASs, we utilized e/sQTL models trained in 11 tissues from the Genotype-Tissue Expression Project including breast, ovary, uterus, vagina, EBV-transformed lymphocytes, whole blood, spleen, liver, subcutaneous adipose, visceral adipose, and cell-cultured fibroblasts. GWAS summary statistics were previously generated from 133,384 BC cases and 113,789 controls who were participants in the Breast Cancer Association Consortium (BCAC). We further performed our TWAS while conditioning e/sQTL effect sizes on nearby GWAS index SNPs. Additionally, we utilized gene-based fine-mapping of eQTLs and sQTL to identify candidate causal genes for each intrinsic subtype. Results: Overall, we identified 164 genes in 69 loci that were associated with Luminal A-like, 19 genes in 9 loci with Luminal B-like, 18 genes in 11 loci with Luminal B/HER2-negative-like, 10 genes in 7 loci with HER2-enriched-like, and 29 genes in 12 loci with TNBC. Among these genes, 17 genes had not been reported in previous TWAS of BC, and 140 genes, 1 gene, 2 genes, 2 genes, and 16 genes were uniquely associated with each of the intrinsic subtypes, respectively. Additionally, we identified one gene associated with Luminal A-like and one gene with Luminal B-like BC that were each in a locus located at least 1.4 Mb from published GWAS hits. Furthermore, we identified 106, 11, 10, 5, and 21 candidate causal genes for each of the intrinsic subtypes, respectively, that had a posterior inclusion probability of 0.9 in at least one e/sQTL. Conclusion: In summary, our multi-tissue TWAS corroborated previous GWAS loci for overall BC risk and intrinsic subtypes, while underscoring how common variation impacts BC etiology by modulating the expression and splicing of genes in multiple tissue types. Citation Format: James Li, Julian McClellan, Haoyu Zhang, Guimin Gao, Dezheng Huo. Multi-tissue transcriptome-wide association studies identified genes for intrinsic subtypes of breast cancer [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO5-09-03.

A multi-tissue, splicing-based joint transcriptome-wide association study identifies susceptibility genes for breast cancer.

Splicing-based transcriptome-wide association studies (splicing-TWASs) of breast cancer have the potential to identify susceptibility genes. However, existing splicing-TWASs test the association of individual excised introns in breast tissue only and thus have limited power to detect susceptibility genes. In this study, we performed a multi-tissue joint splicing-TWAS that integrated splicing-TWAS signals of multiple excised introns in each gene across 11 tissues that are potentially relevant to breast cancer risk. We utilized summary statistics from a meta-analysis that combined genome-wide association study (GWAS) results of 424,650 women of European ancestry. Splicing-level prediction models were trained in GTEx (v.8) data. We identified 240 genes by the multi-tissue joint splicing-TWAS at the Bonferroni-corrected significance level; in the tissue-specific splicing-TWAS that combined TWAS signals of excised introns in genes in breast tissue only, we identified nine additional significant genes. Of these 249 genes, 88 genes in 62 loci have not been reported by previous TWASs, and 17 genes in seven loci are at least 1 Mb away from published GWAS index variants. By comparing the results of our splicing-TWASs with previous gene-expression-based TWASs that used the same summary statistics and expression prediction models trained in the same reference panel, we found that 110 genes in 70 loci that are identified only by the splicing-TWASs. Our results showed that for many genes, expression quantitative trait loci (eQTL) did not show a significant impact on breast cancer risk, whereas splicing quantitative trait loci (sQTL) showed a strong impact through intron excision events.

Association of intestinal anti-inflammatory drug target genes with psychiatric Disorders: A Mendelian randomization study

IntroductionPsychiatric disorders present a substantial global public health burden with limited drug options. The gut-brain axis connects inflammatory bowel diseases and psychiatric disorders, which often have comorbidities. While some evidence hints at anti-inflammatory drugs aiding in treating psychiatric conditions, the specific effects of intestinal anti-inflammatory drugs remain unclear. ObjectivesThis study investigates the causal effect of intestinal anti-inflammatory drug targets on psychiatric disorders. We hypothesize that these drug targets may offer new insights into the treatment and prevention of such disorders. Additionally, we explore gut microbiota’s mediating role between drug target genes and psychiatric disorders. MethodsWe performed two-sample Mendelian randomization (MR) using summary data from existing expression quantitative trait loci (eQTL) and protein QTL in the brain, along with public genome-wide association studies of disease. We also explored gut microbiota’s mediating effect. The statistics encompassed six psychiatric disorders involving 9,725–500,199 individuals. Colocalization analysis enhanced the MR evidence. ResultsWe uncovered a causal link between TPMT (a target of olsalazine) expression in the amygdala and bipolar disorder (BD) risk (odds ratio [OR] = 1.08; P = 4.29 × 10−4). This association was observed even when the sigmoid colon and whole blood eQTL were considered as exposures. Colocalization analysis revealed a shared genetic variant (rs11751561) between TPMT expression and BD, with a posterior probability of 61.6 %. Interestingly, this causal effect was influenced by a decrease in the gut microbiota abundance of the genus Roseburia (effect proportion = 10.05 %). Moreover, elevated ACAT1 expression was associated with higher obsessive–compulsive disorder risk (OR = 1.62; P = 3.64 × 10−4; posterior probability = 3.1 %). ConclusionThese findings provide novel targets for the treatment of psychiatric disorders, underscore the potential of repurposing olsalazine, and emphasize the importance of TPMT and ACAT1 in future drug development.

Multi-omic Analysis of Uterine Leiomyomas in Self-Described Black and White Women: Molecular Insights into Health Disparities

Black women are at an increased risk to develop uterine leiomyomas (ULMs) and to experience worse disease prognosis compared to White women. Epidemiological and molecular factors have been identified as underlying these disparities, but there remains a paucity of deep, multi-omic analysis investigating molecular differences in ULMs from Black and White patients. To identify molecular alterations within ULM tissues correlating with patient race by multi-omic analyses of ULMs collected from cohorts of Black and White women. We performed multi-omic analysis of ULMs from Black (42) and White (47) women undergoing hysterectomy for symptomatic uterine leiomyomata. Our analysis also included application of orthogonal methods to evaluate fibroid biomechanical properties, such as second harmonic generation microscopy, uniaxial compression testing, as well as sheer wave ultrasonography analyses. We found a greater proportion of mediator complex subunit 12 (MED12) mutant ULMs from Black women (>35% increase, Mann Whitney U p = 7E-4). MED12 mutant tumors exhibited elevated abundance of extracellular matrix proteins, including several collagen isoforms, involved in regulation of the core matrisome. Histological analysis of tissue fibrosis using trichrome staining and secondary harmonic generation microscopy confirmed that MED12 mutant tumors are more fibrotic than MED12 wildtype tumors. Using Sheer Wave ultrasonography in a prospectively collected cohort, Black patients had fibroids that were firmer when compared to White patients, even when similar in size. These analyses also uncovered expression quantitative trait loci (eQTL) associated with altered allele frequencies in African and European populations that correlated with differential abundance of several proteins in ULM that are independent of MED12 mutational status, including multiple eQTL mapping to tetratricopeptide repeat protein 38. Our study shows that Black women have a higher prevalence of ULMs harboring mutations in MED12 and that this mutational status correlates with increased tissue fibrosis in comparison with wildtype ULMs. Our study provides insights into molecular alterations underlying racial disparities in ULMs and improves our understanding of the molecular etiology underlying ULM development within these populations.

Integrated analysis of single cell-RNA sequencing and Mendelian randomization identifies lactate dehydrogenase B as a target of melatonin in ischemic stroke.

Despite the success of single-cell RNA sequencing in identifying cellular heterogeneity in ischemic stroke, clarifying the mechanisms underlying these associations of differently expressed genes remains challenging. Several studies that integrate gene expression and gene expression quantitative trait loci (eQTLs) with genome wide-association study (GWAS) data to determine their causal role have been proposed. Here, we combined Mendelian randomization (MR) framework and single cell (sc) RNA sequencing to study how differently expressed genes (DEGs) mediating the effect of gene expression on ischemic stroke. The hub gene was further validated in the invitro model. We identified 2339 DEGs in 10 cell clusters. Among these DEGs, 58 genes were associated with the risk of ischemic stroke. After external validation with eQTL dataset, lactate dehydrogenase B (LDHB) is identified to be positively associated with ischemic stroke. The expression of LDHB has also been validated in sc RNA-seq with dominant expression in microglia and astrocytes, and melatonin is able to reduce the LDHB expression and activity invitro ischemic models. Our study identifies LDHB as a novel biomarker for ischemic stroke via combining the sc RNA-seq and MR analysis.

Cis- and trans-eQTL TWASs of breast and ovarian cancer identify more than 100 susceptibility genes in the BCAC and OCAC consortia.

Transcriptome-wide association studies (TWASs) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have focused on the regulatory effects of risk-associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWASs of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole-genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents an in-depth look into the role of trans eQTLs in the complex molecular mechanisms underlying these diseases.

Advancements of Breeding and Genomics in Wheat (Triticum aestivum L.): Enhancing Yield and Nutritional Value for Sustainable Agriculture

Advancements in wheat breeding and genomics presently explores the genomic interventions driving focusing on quantitative trait loci (QTL) mapping, marker-assisted selection (MAS) and genomic selection (GS). QTL mapping emerges as a pivotal method for pinpointing markers linked with desirable traits, facilitating MAS. Furthermore, genomic selection (GS) holds immense potential for crop improvement. It also delves into the current landscape of MAS and explores various prospects of GS for wheat biofortification. Looking ahead, accelerated mapping studies combined with MAS and GS schemes are poised to further enhance wheat breeding efficiency. Dense molecular maps and a large set of ESTs (Expressed Sequence Tags) have enabled genome-wide identification of gene-rich and gene-poor regions, as well as QTL, including eQTL (Expression quantitative trait loci). Additionally, markers associated with major economic traits have facilitated MAS programs in some countries and enabled map-based cloning of several genes/QTL. Resources for functional genomics, such as TILLING and RNA interference (RNAi), alongside emerging approaches like epigenetics and association mapping, are further enriching wheat genomics research. In this review, we initially present cutting-edge genome-editing technologies in crop plants, with a specific focus on wheat, addressing both functional genomics and genetic enhancement. We subsequently delineate the utilization of additional technologies, including GWAS, high-throughput genotyping and phenotyping, speed breeding, and synthetic biology, within the context of wheat breeding. We assert that integrating genome editing with other molecular breeding strategies will significantly expedite the genetic enhancement of wheat, thus contributing to sustainable global production.

Genome-wide QTL and eQTL mapping reveal genes associated with growth rate trait of the Pacific white shrimp (Litopenaeus vannamei)

BackgroundGrowth rate is a crucial economic trait for farmed animals, but the genetic regulation of this trait is largely unknown in non-model organisms such as shrimp.ResultsIn this study, we performed genome-wide phenotypic quantitative trait loci (QTL) and expression quantitative trait loci (eQTL) mapping analyses to identify genes affecting the growth rate of Pacific white shrimp (Litopenaeus vannamei), which is the most commercially-farmed crustacean worldwide. We used RNA-sequencing of 268 individuals in a mapping population, and subsequently validated our findings through gene silencing and shrimp growth experiments. We constructed a high-density genetic linkage map comprising 5533 markers spanning 44 linkage groups, with a total distance of 6205.75 cM and an average marker interval of 1.12 cM. Our analyses identified 11 QTLs significantly correlated with growth rate, and 117,525 eQTLs. By integrating QTL and eQTL data, we identified a gene (metalloreductase STEAP4) highly associated with shrimp growth rate. RNA interference (RNAi) analysis and growth experiments confirmed that STEAP4 was significantly correlated with growth rate in L. vannamei.ConclusionsOur results indicate that the comprehensive analysis of QTL and eQTL can effectively identify genes involved in complex animal traits. This is important for marker-assisted selection (MAS) of animals. Our work contributes to the development of shrimp breeding and available genetic resources.

Microglial function interacts with the environment to affect sex-specific depression risk

There is a two-fold higher incidence of depression in females compared to men with recent studies suggesting a role for microglia in conferring this sex-dependent depression risk. In this study we investigated the nature of this relation. Using GWAS enrichment, gene-set enrichment analysis and Mendelian randomization, we found minimal evidence for a direct relation between genes functionally related to microglia and sex-dependent genetic risk for depression. We then used expression quantitative trait loci and single nucleus RNA-sequencing resources to generate polygenic scores (PGS) representative of individual variation in microglial function in the adult (UK Biobank; N = 54753–72682) and fetal (ALSPAC; N = 1452) periods. The adult microglial PGS moderated the association between BMI (UK Biobank; beta = 0.001, 95 %CI 0.0009 to 0.003, P = 7.74E-6) and financial insecurity (UK Biobank; beta = 0.001, 95 %CI 0.005 to 0.015, P = 2E-4) with depressive symptoms in females. The fetal microglia PGS moderated the association between maternal prenatal depressive symptoms and offspring depressive symptoms at 24 years in females (ALSPAC; beta = 0.04, 95 %CI 0.004 to 0.07, P = 0.03). We found no evidence for an interaction between the microglial PGS and depression risk factors in males. Our results illustrate a role for microglial function in the conferral of sex-dependent depression risk following exposure to a depression risk factor.

Estimating the effect of lipid-lowering agents on novel subtypes of adult-onset diabetes.

The aims of the present study were to assess the effects of lipid-lowering drugs [HMG-CoA reductase inhibitors, proprotein convertase subtilisin/kexin type 9 inhibitors, and Niemann-Pick C1-Like 1 (NPC1L1) inhibitors] on novel subtypes of adult-onset diabetes through a Mendelian randomisation study. We first inferred causal associations between lipid-related traits [including high-density lipoprotein cholesterol, low-density lipoprotein cholesterol (LDL-C), triglycerides (TG), apolipoproteins A-I, and apolipoproteins B] and novel subtypes of adult-onset diabetes. The expression quantitative trait loci of drug target genes for three classes of lipid-lowering drugs, as well as genetic variants within or nearby drug target genes associated with LDL-C, were then utilised as proxies for the exposure of lipid-lowering drugs. Mendelian randomisation analysis was performed using summary data from genome-wide association studies of LDL-C, severe autoimmune diabetes, severe insulin-deficient diabetes (SIDD), severe insulin-resistant diabetes (SIRD), mild obesity-related diabetes (MOD), and mild age-related diabetes. There was an association between HMGCR-mediated LDL-C and the risk of SIRD [odds ratio (OR)=0.305, 95% confidence interval (CI)=0.129-0.723; p=0.007], and there was an association of PCSK9-mediated LDL-C with the risk of SIDD (OR=0.253, 95% CI=0.120-0.532; p<0.001) and MOD (OR=0.345, 95% CI=0.171-0.696; p=0.003). Moreover, NPC1L1-mediated LDL-C (OR=0.109, 95% CI=0.019-0.613; p=0.012) and the increased expression of NPC1L1 gene in blood (OR=0.727, 95% CI=0.541-0.977; p=0.034) both showed a significant association with SIRD. These results were further confirmed by sensitivity analyses. In summary, the different lipid-lowering medications have a specific effect on the increased risk of different novel subtypes of adult-onset diabetes.