Association of Cadherin-Related Family Member 1 with Traumatic Brain Injury

The cadherin family plays a pivotal role in orchestrating synapse formation in the central nervous system. Cadherin-related family member 1 (CDHR1) is a photoreceptor-specific calmodulin belonging to the expansive cadherin superfamily. However, its role in traumatic brain injury (TBI) remains largely unknown. CDHR1 expression across various brain tissue sites was analyzed using the GSE104687 dataset. Employing a summary-data-based Mendelian Randomization (SMR) approach, integrated analyses were performed by amalgamating genome-wide association study abstracts from TBI with public data on expressed quantitative trait loci and DNA methylation QTL from both blood and diverse brain tissues. CDHR1 expression and localization in different brain tissues were meticulously delineated using western blotting, immunohistochemistry, and enzyme-linked immunosorbent assay. CDHR1 expression was consistently elevated in the TBI group compared to that in the sham group across multiple tissues. The inflammatory response emerged as a crucial biological mechanism, and pro-inflammatory and anti-inflammatory factors were not expressed in either group. Integrated SMR analyses encompassing both blood and brain tissues substantiated the heightened CDHR1 expression profiles, with methylation modifications emerging as potential contributing factors for increased TBI risk. This was corroborated by western blotting and immunohistochemistry, confirming augmented CDHR1 expression following TBI. This multi-omics-based genetic association study highlights the elevated TBI risk associated with CDHR1 expression coupled with putative methylation modifications. These findings provide compelling evidence for future targeted investigations and offer promising avenues for developing interventional therapies for TBI.

Mendelian randomization of circulating proteome identifies IFN-γ as a druggable target in aplastic anemia

BACKGROUND

Aplastic anemia (AA) is a kind of bone marrow failure (BMF) characterized by pancytopenia with hypoplasia/aplasia of bone marrow. Immunosuppressive therapy and bone marrow transplantation are effective methods to treat severe aplastic anemia. However, the efficacy is limited by complications and the availability of suitable donors. This study aimed to determine whether any circulating druggable protein levels may have causal effects on AA and provide potential novel drug targets for AA.

METHODS

Genetic variants strongly associated with circulating druggable protein levels to perform Mendelian randomization (MR) analyses were used. The effect of these druggable protein levels on AA risk was measured using the summary statistics from a large-scale proteomic genome-wide association study (GWAS) and FinnGen database ( https://www.finngen.fi/en/access_results ). Multivariable MR analyses were performed to statistically adjust for potential confounders, including platelet counts, reticulocyte counts, neutrophil counts, and proportions of hematopoietic stem cells.

RESULTS

The data showed that higher level of circulating IFN-γ levels was causally associated with AA susceptibility. The causal effects of circulating IFN-γ levels on the AA were broadly consistent, when adjusted for platelet counts, reticulocyte counts, neutrophil counts and proportions of hematopoietic stem cells.

CONCLUSIONS

High levels of circulating IFN-γ levels might increase the risk of AA and might provide a potential novel target for AA.

Study on the causal relationship between educational attainment and delirium: A two-sample Mendelian randomization study

This study aimed to investigate whether there is a causal relationship between educational attainment and delirium at the genetic level using the Mendelian randomization method, and provide new evidence for studies in this field. We found a causal relationship between educational attainment and delirium at the genetic level after excluding confounders using Mendelian randomization. The inverse variance weighting method of random effects was the main analysis method. The weighted median and Mendelian Randomization-Egger methods, as well as simple, and weighted modes were used as supplementary analysis methods. Additionally, horizontal pleiotropy tests were conducted, including the Mendelian Randomization-Egger intercept test and Mendelian Randomization Pleiotropy RESidual Sum and Outlier. Cochran's Q statistic was used to assess the size of heterogeneity. We retrieved all second single nucleotide polymorphism features and performed multivariate Mendelian randomization to adjust for the effect of potential confounders on our results. The inverse variance weighting suggested a negative correlation between genetically predicted educational attainment and delirium (0.67[0.49-0.92], p = 0.013); Mendelian Randomization Pleiotropy RESidual Sum and Outlier (0.67[0.49-0.92], p = 0.013) and multivariate Mendelian randomization (0.52[0.33-0.82], p = 0.005) results were generally consistent with the inverse variance weighting method. The Mendelian Randomization-Egger, simple, and weighted mode results were consistent with the inverse variance weighting results. Our results were not affected by pleiotropy or heterogeneity (p > 0.05, for both pleiotropy and heterogeneity). In addition, the "leave-one-out" analysis showed that the results of our Mendelian randomization analysis were not influenced by individual single nucleotide polymorphisms. Studies have found a causal relationship between educational attainment and delirium at the genetic level; higher educational attainment may be a protective factor against delirium. Clinically, more attention should be paid to patients at a high risk of delirium with low educational attainment.

Limbic-predominant age-related TDP-43 encephalopathy (LATE-NC): Co-pathologies and genetic risk factors provide clues about pathogenesis

Limbic-predominant age-related TDP-43 encephalopathy neuropathologic change (LATE-NC) is detectable at autopsy in more than one-third of people beyond age 85 years and is robustly associated with dementia independent of other pathologies. Although LATE-NC has a large impact on public health, there remain uncertainties about the underlying biologic mechanisms. Here, we review the literature from human studies that may shed light on pathogenetic mechanisms. It is increasingly clear that certain combinations of pathologic changes tend to coexist in aging brains. Although "pure" LATE-NC is not rare, LATE-NC often coexists in the same brains with Alzheimer disease neuropathologic change, brain arteriolosclerosis, hippocampal sclerosis of aging, and/or age-related tau astrogliopathy (ARTAG). The patterns of pathologic comorbidities provide circumstantial evidence of mechanistic interactions ("synergies") between the pathologies, and also suggest common upstream influences. As to primary mediators of vulnerability to neuropathologic changes, genetics may play key roles. Genes associated with LATE-NC include TMEM106B, GRN, APOE, SORL1, ABCC9, and others. Although the anatomic distribution of TDP-43 pathology defines the condition, important cofactors for LATE-NC may include Tau pathology, endolysosomal pathways, and blood-brain barrier dysfunction. A review of the human phenomenology offers insights into disease-driving mechanisms, and may provide clues for diagnostic and therapeutic targets.

Genome-wide association study and development of molecular markers for yield and quality traits in peanut (Arachis hypogaea L.)

BACKGROUND

This study aims to decipher the genetic basis governing yield components and quality attributes of peanuts, a critical aspect for advancing molecular breeding techniques. Integrating genotype re-sequencing and phenotypic evaluations of seven yield components and two grain quality traits across four distinct environments allowed for the execution of a genome-wide association study (GWAS).

RESULTS

The nine phenotypic traits were all continuous and followed a normal distribution. The broad heritability ranged from 88.09 to 98.08%, and the genotype-environment interaction effects were all significant. There was a highly significant negative correlation between protein content (PC) and oil content (OC). The 10× genome re-sequencing of 199 peanut accessions yielded a total of 631,988 high-quality single nucleotide polymorphisms (SNPs), with 374 significant SNP loci identified in association with the nine traits of interest. Notably, 66 of these pertinent SNPs were detected in multiple environments, and 48 of them were linked to multiple traits of interest. Five loci situated on chromosome 16 (Chr16) exhibited pleiotropic effects on yield traits, accounting for 17.64-32.61% of the observed phenotypic variation. Two loci on Chr08 were found to be strongly associated with protein and oil contents, accounting for 12.86% and 14.06% of their respective phenotypic variations, respectively. Linkage disequilibrium (LD) block analysis of these seven loci unraveled five nonsynonymous variants, leading to the identification of one yield-related candidate gene and two quality-related candidate genes. The correlation between phenotypic variation and SNP loci in these candidate genes was validated by Kompetitive allele-specific PCR (KASP) marker analysis.

CONCLUSIONS

Overall, molecular markers were developed for genetic loci associated with yield and quality traits through a GWAS investigation of 199 peanut accessions across four distinct environments. These molecular tools can aid in the development of desirable peanut germplasm with an equilibrium of yield and quality through marker-assisted breeding.

Genome-wide analyses reveal shared genetic architecture and novel risk loci between opioid use disorder and general cognitive ability

BACKGROUND

Opioid use disorder (OUD), a serious health burden worldwide, is associated with lower cognitive function. Recent studies have demonstrated a negative genetic correlation between OUD and general cognitive ability (COG), indicating a shared genetic basis. However, the specific genetic variants involved, and the underlying molecular mechanisms remain poorly understood. Here, we aimed to quantify and identify the genetic basis underlying OUD and COG.

METHODS

We quantified the extent of genetic overlap between OUD and COG using a bivariate causal mixture model (MiXeR) and identified specific genetic loci applying conditional/conjunctional FDR. Finally, we investigated biological function and expression of implicated genes using available resources.

RESULTS

We estimated that ~94% of OUD variants (4.8k out of 5.1k variants) also influence COG. We identified three novel OUD risk loci and one locus shared between OUD and COG. Loci identified implicated biological substrates in the basal ganglia.

CONCLUSION

We provide new insights into the complex genetic risk architecture of OUD and its genetic relationship with COG.

Comparative pangenomic analysis of Campylobacter fetus isolated from Spanish bulls and other mammalian species

Campylobacter fetus comprises two closely related mammal-associated subspecies: Campylobacter fetus subsp. fetus (Cff) and Campylobacter fetus subsp. venerealis (Cfv). The latter causes bovine genital campylobacteriosis, a sexually-transmitted disease endemic in Spain that results in significant economic losses in the cattle industry. Here, 33 C. fetus Spanish isolates were whole-genome sequenced and compared with 62 publicly available C. fetus genomes from other countries. Genome-based taxonomic identification revealed high concordance with in silico PCR, confirming Spanish isolates as Cff (n = 4), Cfv (n = 9) and Cfv biovar intermedius (Cfvi, n = 20). MLST analysis assigned the Spanish isolates to 6 STs, including three novel: ST-76 and ST-77 for Cfv and ST-78 for Cff. Core genome SNP phylogenetic analysis of the 95 genomes identified multiple clusters, revealing associations at subspecies and biovar level between genomes with the same ST and separating the Cfvi genomes from Spain and other countries. A genome-wide association study identified pqqL as a Cfv-specific gene and a potential candidate for more accurate identification methods. Functionality analysis revealed variations in the accessory genome of C. fetus subspecies and biovars that deserve further studies. These results provide valuable information about the regional variants of C. fetus present in Spain and the genetic diversity and predicted functionality of the different subspecies.

Genome-wide association study and polygenic risk scores of retinal thickness across the cognitive continuum: data from the NORFACE cohort

BACKGROUND

Several studies have reported a relationship between retinal thickness and dementia. Therefore, optical coherence tomography (OCT) has been proposed as an early diagnosis method for Alzheimer's disease (AD). In this study, we performed a genome-wide association study (GWAS) aimed at identifying genes associated with retinal nerve fiber layer (RNFL) and ganglion cell inner plexiform layer (GCIPL) thickness assessed by OCT and exploring the relationships between the spectrum of cognitive decline (including AD and non-AD cases) and retinal thickness.

METHODS

RNFL and GCIPL thickness at the macula were determined using two different OCT devices (Triton and Maestro). These determinations were tested for association with common single nucleotide polymorphism (SNPs) using adjusted linear regression models and combined using meta-analysis methods. Polygenic risk scores (PRSs) for retinal thickness and AD were generated.

RESULTS

Several genetic loci affecting retinal thickness were identified across the genome in accordance with previous reports. The genetic overlap between retinal thickness and dementia, however, was weak and limited to the GCIPL layer; only those observable with all-type dementia cases were considered.

CONCLUSIONS

Our study does not support the existence of a genetic link between dementia and retinal thickness.

The role of leptin in ischemic stroke and its subtypes: A Mendelian randomization study

BACKGROUND AND AIMS

Observational studies have suggested a relationship between leptin and risk of stroke. However, evidence for the association remains inconsistent, and whether the association reflects a causal relationship remains to be established. To clarify this relationship, we adopted a two-sample Mendelian randomization (MR) analysis to investigate whether leptin plays a causal role in the risk of stroke and its subtypes.

METHODS AND RESULTS

Five independent single-nucleotide polymorphisms (SNPs) associated with the leptin level from genome-wide association studies (GWASs) of European individuals were selected. We performed an MR analysis using the inverse-variance-weighted (IVW) as primary method to examine the causal effects of leptin on ischemic stroke (IS). Moreover, MR-Egger intercept and Cochran's Q statistic were also performed to detect the pleiotropy or heterogeneity of our MR results. Genetically predicted circulating leptin level was not associated with ischemic stroke [odds ratio (OR): 1.48, 95% confidence interval (CI): 0.78-2.8, P = 0.22], large artery stroke (OR: 1.44, 95% CI: 0.39-5.25, P = 0.57), cardioembolic stroke (OR:1.33, 95% CI: 0.55-3.22, P = 0.52), and small vessel stroke (OR: 1.48, 95% CI: 0.39-5.63, P = 0.56) using the IVW method. Likewise, there is no convincing evidence for the associations between leptin levels and cardiovascular diseases (CVD) risk factors.

CONCLUSIONS

This study did not provide evidence that leptin levels are associated with increased risk of stroke and its subtypes.

Genetic causal association between frozen shoulder and carpal tunnel syndrome: a two-sample mendelian randomization

OBJECTIVE

Observational studies have suggested an association between frozen shoulder (FS) and carpal tunnel syndrome (CTS). However, due to challenges in establishing a temporal sequence, the causal relationship between these two conditions remains elusive. This study, based on aggregated data from large-scale population-wide genome-wide association studies (GWAS), investigates the genetic causality between FS and CTS.

METHODS

Initially, a series of quality control measures were employed to select single nucleotide polymorphisms (SNPs) closely associated with the exposure factors. Two-sample Mendelian randomization (MR) was utilized to examine the genetic causality between FS and CTS, employing methods including Inverse-Variance Weighted (IVW), MR-Egger, Weighted Median, Simple Mode, and Weighted Mode approaches. Subsequently, sensitivity analyses were conducted to assess the robustness of the MR analysis results.

RESULTS

IVW analysis results indicate a positive causal relationship between CTS and FS (p < 0.05, OR > 1), while a negative causal relationship between the two conditions was not observed. Heterogeneity tests suggest minimal heterogeneity in our IVW analysis results (p > 0.05). Multivariable MR testing also indicates no pleiotropy in our IVW analysis (p > 0.05), and stepwise exclusion tests demonstrate the reliability and stability of the MR analysis results. Gene Ontology (GO) pathway analysis reveals enrichment of genes regulated by the associated SNPs in the TGFβ-related pathways.

CONCLUSION

This study provides evidence of the genetic causal association between frozen shoulder and carpal tunnel syndrome and provides new insights into the genetics of fibrotic disorders.

FORGEdb: a tool for identifying candidate functional variants and uncovering target genes and mechanisms for complex diseases

The majority of disease-associated variants identified through genome-wide association studies are located outside of protein-coding regions. Prioritizing candidate regulatory variants and gene targets to identify potential biological mechanisms for further functional experiments can be challenging. To address this challenge, we developed FORGEdb ( https://forgedb.cancer.gov/ ; https://forge2.altiusinstitute.org/files/forgedb.html ; and https://doi.org/10.5281/zenodo.10067458 ), a standalone and web-based tool that integrates multiple datasets, delivering information on associated regulatory elements, transcription factor binding sites, and target genes for over 37 million variants. FORGEdb scores provide researchers with a quantitative assessment of the relative importance of each variant for targeted functional experiments.

Schizophrenia in the genetic era: a review from development history, clinical features and genomic research approaches to insights of susceptibility genes

Schizophrenia is a devastating neuropsychiatric disorder affecting 1% of the world population and ranks as one of the disorders providing the most severe burden for society. Schizophrenia etiology remains obscure involving multi-risk factors, such as genetic, environmental, nutritional, and developmental factors. Complex interactions of genetic and environmental factors have been implicated in the etiology of schizophrenia. This review provides an overview of the historical origins, pathophysiological mechanisms, diagnosis, clinical symptoms and corresponding treatment of schizophrenia. In addition, as schizophrenia is a polygenic, genetic disorder caused by the combined action of multiple micro-effective genes, we further detail several approaches, such as candidate gene association study (CGAS) and genome-wide association study (GWAS), which are commonly used in schizophrenia genomics studies. A number of GWASs about schizophrenia have been performed with the hope to identify novel, consistent and influential risk genetic factors. Finally, some schizophrenia susceptibility genes have been identified and reported in recent years and their biological functions are also listed. This review may serve as a summary of past research on schizophrenia genomics and susceptibility genes (NRG1, DISC1, RELN, BDNF, MSI2), which may point the way to future schizophrenia genetics research. In addition, depending on the above discovery of susceptibility genes and their exact function, the development and application of antipsychotic drugs will be promoted in the future.

Exploring the causal association of rheumatoid arthritis with atrial fibrillation: a Mendelian randomization study

BACKGROUND

It has been proved that rheumatoid arthritis (RA) patients have high incidence of atrial fibrillation (AF). Nevertheless, whether they have causal relevance is uncertain. This study aimed to explore and verify the authenticity of causal relationship between RA and AF using Mendelian randomization (MR).

METHODS

The genome-wide association study (GWAS) summary data from Biobank Japan Project (BBJ) (RA, 4199 cases and 208,254 controls) were regarded as exposure data and the GWAS data from European Bio-informatics Institute database (EBI) (AF, 15,979 cases and 102,776 controls) as outcome data. The causal effect was appraised by the inverse variance weighted (IVW) method, MR-Egger regression, and weighted median estimator. MR-robust adjusted profile score (MR-RAPS) method was delivered to examine the robustness of causal relationship and MR Pleiotropy Residual Sum and Outlier (MR-PRESSO) method to control horizontal (directional) pleiotropy.

RESULTS

The results indicated that RA increased the risk of AF (IVW, the odds ratio (OR) = 1.060; 95% confidence interval (CI), 1.028 to 1.092; p = 1.411 × 10-4; weighted median, OR = 1.046, 95% CI, 1.002 to 1.093, p = 0.047). The MR analysis also showed this causal effect through all four IVW methods with various statistical algorithms. Both MR-RAPS and MR-PRESSO supported the causality of RA and AF. Also, the MR-PRESSO result indicated the absence of apparent pleiotropy.

CONCLUSION

There is a causal association between RA and AF. RA patients are genetically more vulnerable to AF. This study may contribute to further exploring early clinical prevention and fundamental mechanism of AF in patients with RA. Key Points • We provided some genetic evidence for the causal link between rheumatoid arthritis (RA) and atrial fibrillation (AF) with multiple Mendelian randomization (MR) methods. • RA patients were genetically more vulnerable to AF. • This study partly shed light on latent fundamental mechanisms underlying RA-induced AF and inspired future studies on RA-AF relationship.

Genome-wide QTL mapping across three tissues highlights several Alzheimer's and Parkinson's disease loci potentially acting via DNA methylation

DNA methylation (DNAm) is an epigenetic mark with essential roles in disease development and predisposition. Here, we created genome-wide maps of methylation quantitative trait loci (meQTL) in three peripheral tissues and used Mendelian randomization (MR) analyses to assess the potential causal relationships between DNAm and risk for two common neurodegenerative disorders, i.e. Alzheimer's disease (AD) and Parkinson's disease (PD). Genome-wide single nucleotide polymorphism (SNP; ~5.5M sites) and DNAm (~850K CpG sites) data were generated from whole blood (n=1,058), buccal (n=1,527) and saliva (n=837) specimens. We identified between 11 and 15 million genome-wide significant (p<10-14) SNP-CpG associations in each tissue. Combining these meQTL GWAS results with recent AD/PD GWAS summary statistics by MR strongly suggests that the previously described associations between PSMC3, PICALM, and TSPAN14 and AD may be founded on differential DNAm in or near these genes. In addition, there is strong, albeit less unequivocal, support for causal links between DNAm at PRDM7 in AD as well as at KANSL1/MAPT in AD and PD. Our study adds valuable insights on AD/PD pathogenesis by combining two high-resolution "omics" domains, and the meQTL data shared along with this publication will allow like-minded analyses in other diseases.

Genome-wide association study reveals HSF2, GJA1 and TRIM36 as susceptibility genes for preeclampsia: a community-based population study in Tianjin, China

BACKGROUND

Preeclampsia (PE) mainly occurs in pregnant women and is hereditary. Several genome-wide association studies (GWAS) on Caucasian samples have reported some gene loci that are associated with preeclampsia. However, these studies have not reached consistent conclusions. No previous GWAS has examined preeclampsia in the Chinese Han population.

METHOD

This study aimed to identify common genetic variations associated with preeclampsia in the Chinese Han population through two-stage case‒control studies. The discovery cohort included 92 patients with severe preeclampsia and 187 healthy controls. The validation cohort included 52 patients with preeclampsia and 104 controls. A genome-wide association study was performed to identify putative preeclampsia genes in the discovery cohort, with validation in the validation cohort.

RESULTS

In the discovery cohort, GWAS demonstrated that 19 single-nucleotide polymorphisms (SNPs) were associated with preeclampsia (P < 10-5). The pathway analysis revealed that these 19 SNP representative genes were mainly enriched in the adenylyl cyclase-inhibiting G-protein coupled receptor signaling pathway. After validation in the validation cohort, rs13176432 and rs13210237 remained closely related to preeclampsia (P<0.05). In the combined data set, the frequency of the G allele in rs13176432 was significantly higher in cases with preeclampsia than in controls (P = 5 × 10-6). The frequency of the A allele in rs13210237 was higher in the preeclampsia group (P = 8 × 10-6). The rs13210237 representative genes include HSF2 and GJA1, while the rs13176432 representative gene is TRIM36. There were no differences in genotype distribution between the early-onset and late-onset preeclampsia groups (P > 0.05). Furthermore, rs13210237 and rs13176432 were related to preeclampsia in the adjusted regression model (P < 0.000).

CONCLUSION

In this study of two independent cohorts, we found that rs13210237 and rs13176432 might be novel preeclampsia-susceptible genetic factors in the Han population in China. However, there was no association between the onset of preeclampsia and these genotypes.

Zebrafish as a Model for Osteoporosis: Functional Validations of Genome-Wide Association Studies

PURPOSE OF REVIEW

GWAS, as a largely correlational analysis, requires in vitro or in vivo validation. Zebrafish (Danio rerio) have many advantages for studying the genetics of human diseases. Since gene editing in zebrafish has been highly valuable for studying embryonic skeletal developmental processes that are prenatally or perinatally lethal in mammalian models, we are reviewing pros and cons of this model.

RECENT FINDINGS

The true power for the use of zebrafish is the ease by which the genome can be edited, especially using the CRISPR/Cas9 system. Gene editing, followed by phenotyping, for complex traits such as BMD, is beneficial, but the major physiological differences between the fish and mammals must be considered. Like mammals, zebrafish do have main bone cells; thus, both in vivo stem cell analyses and in vivo imaging are doable. Yet, the "long" bones of fish are peculiar, and their bone cavities do not contain bone marrow. Partial duplication of the zebrafish genome should be taken into account. Overall, small fish toolkit can provide unmatched opportunities for genetic modifications and morphological investigation as a follow-up to human-first discovery.

Genome-wide identification of m6A-associated single nucleotide polymorphisms in complex diseases of nervous system

N6-methyladenosine (m6A) is one of the most abundant chemical modifications on RNA and can affect the occurrence and development of diseases. Some studies have shown that the expressions of some m6A-related genes are significantly regulated by single nucleotide variants (SNV). However, the function of m6A-associated single nucleotide polymorphisms (m6A-SNP) remains unclear in multiple sclerosis (MS), Alzheimer's disease (AD) and Parkinson's disease (PD). Here, we identified the disease-associated m6A-SNPs by integrating genome-wide association study (GWAS) and m6A-SNPs from the RMVar database, and confirmed the relationship between these identified m6A-SNPs and their target genes in eQTL analysis and gene differential expression analysis. Finally, 26 genes corresponding to 20 m6A-SNPs with eQTL signals were identified and differentially expressed (P < 0.05) in MS, 15 genes corresponding to 12 m6A-SNPs (P < 1e-04) were differentially expressed in AD, and 27 PD-associated m6A-SNPs that regulated the expression of 31 genes were identified. There were 5 HLA genes with eQTL signals (HLA-DQB1, HLA-DRB1, HLA-DQA1, HLA-DQA2 and HLA-DQB1-AS1) to be detected in the three diseases. In summary, our study provided new insights into understanding the potential roles of these m6A-SNPs in disease pathogenesis as well as therapeutic target.

Genomic region and origin for selected traits during differentiation of small-fruit cultivars in Japanese apricot (Prunus mume)

The Japanese apricot (Prunus mume) is a popular fruit tree in Japan. However, the genetic factors associated with fruit trait variations are poorly understood. In this study, we investigated nine fruit-associated traits, including harvesting time, fruit diameter, fruit shape, fruit weight, stone (endocarp) weight, ratio of stone weight to fruit weight, and rate of fruit gumming, using 110 Japanese apricot accessions over four years. A genome-wide association study (GWAS) was performed for these traits and strong signals were detected on chromosome 6 for harvesting time and fruit diameters. These peaks were shown to undergo strong artificial selection during the differentiation of small-fruit cultivars. The genomic region defined by the GWAS and XP-nSL analyses harbored several candidate genes associated with plant hormone regulation. Furthermore, the alleles of small-fruit cultivars in this region were shown to have genetic proximity to some Chinese cultivars of P. mume. These results indicate that the small-fruit trait originated in China; after being introduced into Japan, it was preferred and selected by the Japanese people, resulting in the differentiation of small-fruit cultivars.

Genome-wide association studies to identify quantitative trait loci and positional candidate genes affecting meat quality-related traits in pigs

Meat quality comprises a set of key traits such as pH, meat color, water-holding capacity, tenderness and marbling. These traits are complex because they are affected by multiple genetic and environmental factors. The aim of this study was to investigate the molecular genetic basis underlying nine meat quality-related traits in a Yorkshire pig population using a genome-wide association study (GWAS) and subsequent biological pathway analysis. In total, 45,926 single nucleotide polymorphism (SNP) markers from 543 pigs were selected for the GWAS after quality control. Data were analyzed using a genome-wide efficient mixed model association (GEMMA) method. This linear mixed model-based approach identified two quantitative trait loci (QTLs) for meat color (b*) on chromosome 2 (SSC2) and one QTL for shear force on chromosome 8 (SSC8). These QTLs acted additively on the two phenotypes and explained 3.92%-4.57% of the phenotypic variance of the traits of interest. The genes encoding HAUS8 on SSC2 and an lncRNA on SSC8 were identified as positional candidate genes for these QTLs. The results of the biological pathway analysis revealed that positional candidate genes for meat color (b*) were enriched in pathways related to muscle development, muscle growth, intramuscular adipocyte differentiation, and lipid accumulation in muscle, whereas positional candidate genes for shear force were overrepresented in pathways related to cell growth, cell differentiation, and fatty acids synthesis. Further verification of these identified SNPs and genes in other independent populations could provide valuable information for understanding the variations in pork quality-related traits.

CaHMA1 promotes Cd accumulation in pepper fruit

The main planting areas for pepper (Capsicum sp.) are high in cadmium (Cd), which is the most prevalent heavy metal pollutant worldwide. Breeding pepper cultivars with low Cd levels can promote sustainable agricultural production and ensure the safety of pepper products. To identify breeding targets for reducing Cd accumulation in pepper fruits, we performed a genome-wide association study on 186 accessions. Polymorphisms were associated with fruit Cd content in a genomic region containing a homolog of Arabidopsis (Arabidopsis thaliana) Heavy metal-transporting ATPase 1 (HMA1) encoding a P-type ATPase. In two cultivars with contrasting Cd accumulation, transcriptome analysis revealed differentially expressed genes enriched for carbohydrate metabolism and photosynthesis in fruits with high Cd accumulation, and a Cd2+/Zn2+-exporting ATPase gene (HMA). Heterologous expression of CaHMA1 in yeast increases Cd sensitivity. Overexpression of CaHMA1 conferred a severe increase in Cd content in Arabidopsis plants, whereas reduced CaHMA1 expression in pepper fruits decreased Cd content. We propose that CaHMA1 expression may be an important component of the high Cd accumulation in pepper plants.

Multivariate GWAS of Alzheimer's disease CSF biomarker profiles implies GRIN2D in synaptic functioning

BACKGROUND

Genome-wide association studies (GWAS) of Alzheimer's disease (AD) have identified several risk loci, but many remain unknown. Cerebrospinal fluid (CSF) biomarkers may aid in gene discovery and we previously demonstrated that six CSF biomarkers (β-amyloid, total/phosphorylated tau, NfL, YKL-40, and neurogranin) cluster into five principal components (PC), each representing statistically independent biological processes. Here, we aimed to (1) identify common genetic variants associated with these CSF profiles, (2) assess the role of associated variants in AD pathophysiology, and (3) explore potential sex differences.

METHODS

We performed GWAS for each of the five biomarker PCs in two multi-center studies (EMIF-AD and ADNI). In total, 973 participants (n = 205 controls, n = 546 mild cognitive impairment, n = 222 AD) were analyzed for 7,433,949 common SNPs and 19,511 protein-coding genes. Structural equation models tested whether biomarker PCs mediate genetic risk effects on AD, and stratified and interaction models probed for sex-specific effects.

RESULTS

Five loci showed genome-wide significant association with CSF profiles, two were novel (rs145791381 [inflammation] and GRIN2D [synaptic functioning]) and three were previously described (APOE, TMEM106B, and CHI3L1). Follow-up analyses of the two novel signals in independent datasets only supported the GRIN2D locus, which contains several functionally interesting candidate genes. Mediation tests indicated that variants in APOE are associated with AD status via processes related to amyloid and tau pathology, while markers in TMEM106B and CHI3L1 are associated with AD only via neuronal injury/inflammation. Additionally, seven loci showed sex-specific associations with AD biomarkers.

CONCLUSIONS

These results suggest that pathway and sex-specific analyses can improve our understanding of AD genetics and may contribute to precision medicine.

Causality between heart failure and epigenetic age: a bidirectional Mendelian randomization study

AIMS

Heart failure (HF) is a prevalent age-related cardiovascular disease with poor prognosis in the elderly population. This study aimed to establish the causal relationship between ageing and HF by conducting a bidirectional Mendelian randomization (MR) analysis on epigenetic age (a marker of ageing) and HF.

METHODS AND RESULTS

Genome-wide association study data for epigenetic age (GrimAge, HorvathAge, HannumAge, and PhenoAge) and HF were collected and assessed for significant genetic variables. A bidirectional MR analysis was carried out using the random-effects inverse-variance weighted (IVW) method as the primary approach, while other methods (MR-Egger, weighted median, simple mode, and weighted mode) and multiple sensitivity analyses (heterogeneity analysis, leave-one-out sensitivity analysis, and horizontal pleiotropy analysis) were employed to evaluate the impact of epigenetic age on HF and vice versa. Bidirectional MR analysis of two samples revealed that the epigenetic PhenoAge clock increased the risk of HF [IVW odds ratio (OR) 1.015, 95% confidence interval (CI) 1.002-1.028, P = 0.028 and weighted median OR 1.020, 95% CI 1.001-1.038, P = 0.039]. Other results were not statistically significant.

CONCLUSIONS

The bidirectional MR analysis demonstrated a causal link between genetically predicted epigenetic age and HF in individuals of European descent. Further research into epigenetic age in other populations and additional genetic information related to HF is warranted.

Genome-wide association study of cooking-caused grain expansion in rice (Oryza sativa L.)

Cooking-caused rice grain expansion (CCRGE) is a critical trait for evaluating the cooking quality of rice. Previous quantitative trait locus (QTL) mapping studies on CCRGE have been limited to bi-parental populations, which restrict the exploration of natural variation and mapping resolution. To comprehensively and precisely dissect the genetic basis of CCRGE, we performed a genome-wide association study (GWAS) on three related indices: grain breadth expansion index (GBEI), grain length expansion index (GLEI), and grain length-breadth ratio expansion index (GREI), using 345 rice accessions grown in two years (environments) and 193,582 SNP markers. By analyzing each environment separately using seven different methods (3VmrMLM, mrMLM, FASTmrMLM, FASTmrEMMA, pLARmEB, pKWmEB, ISIS EM-BLASSO), we identified a total of 32, 19 and 27 reliable quantitative trait nucleotides (QTNs) associated with GBEI, GLEI and GREI, respectively. Furthermore, by jointly analyzing the two environments using 3VmrMLM, we discovered 19, 22 and 25 QTNs, as well as 9, 5 and 7 QTN-by-environment interaction (QEIs) associated with GBEI, GLEI and GREI, respectively. Notably, 12, 9 and 15 QTNs for GBEI, GLEI and GREI were found within the intervals of previously reported QTLs. In the vicinity of these QTNs or QEIs, based on analyses of mutation type, gene ontology classification, haplotype, and expression pattern, we identified five candidate genes that are related to starch synthesis and endosperm development. The five candidate genes, namely, LOC_Os04g53310 (OsSSIIIb, near QTN qGREI-4.5s), LOC_Os05g02070 (OsMT2b, near QTN qGLEI-5.1s), LOC_Os06g04200 (wx, near QEI qGBEI-6.1i and QTNs qGREI-6.1s and qGLEI-6.1t), LOC_Os06g12450 (OsSSIIa, near QTN qGLEI-6.2t), and LOC_Os08g09230 (OsSSIIIa, near QTN qGBEI-8.1t), are predicted to be involved in the process of rice grain starch synthesis and to influence grain expansion after cooking. Our findings provide valuable insights and will facilitate genetic research and improvement of CCRGE.

Genome-wide association study reveals candidate genes controlling root system architecture under low phosphorus supply at seedling stage in Brassica napus

Optimal root system architecture (RSA) is essential for vigorous growth and yield in crops. Plants have evolved adaptive mechanisms in response to low phosphorus (LP) stress, and one of those is changes in RSA. Here, more than five million single-nucleotide polymorphisms (SNPs) obtained from whole-genome re-sequencing data (WGR) of an association panel of 370 oilseed rape (Brassica napus L.) were used to conduct a genome-wide association study (GWAS) of RSA traits of the panel at LP in "pouch and wick" system. Fifty-two SNPs were forcefully associated with lateral root length (LRL), total root length (TRL), lateral root density (LRD), lateral root number (LRN), mean lateral root length (MLRL), and root dry weight (RDW) at LP. There were significant correlations between phenotypic variation and the number of favorable alleles of the associated loci on chromosomes A06 (chrA06_20030601), C03 (chrC03_3535483), and C07 (chrC07_42348561), respectively. Three candidate genes (BnaA06g29270D, BnaC03g07130D, and BnaC07g43230D) were detected by combining transcriptome, candidate gene association analysis, and haplotype analysis. Cultivar carrying "CCGC" at BnaA06g29270DHap1, "CAAT" at BnaC03g07130DHap1, and "ATC" at BnaC07g43230DHap1 had greater LRL, LRN, and RDW than lines carrying other haplotypes at LP supply. The RSA of a cultivar harboring the three favorable haplotypes was further confirmed by solution culture experiments. These findings define exquisite insights into genetic architectures underlying B. napus RSA at LP and provide valuable gene resources for root breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01411-2.

Genetic diversity and genome-wide association study of 13 agronomic traits in 977 Beta vulgaris L. germplasms

BACKGROUND

Sugar beet (Beta vulgaris L.) is an economically essential sugar crop worldwide. Its agronomic traits are highly diverse and phenotypically plastic, influencing taproot yield and quality. The National Beet Medium-term Gene Bank in China maintains more than 1700 beet germplasms with diverse countries of origin. However, it lacks detailed genetic background associated with morphological variability and diversity.

RESULTS

Here, a comprehensive genome-wide association study (GWAS) of 13 agronomic traits was conducted in a panel of 977 sugar beet accessions. Almost all phenotypic traits exhibited wide genetic diversity and high coefficient of variation (CV). A total of 170,750 high-quality single-nucleotide polymorphisms (SNPs) were obtained using the genotyping-by-sequencing (GBS). Neighbour-joining phylogenetic analysis, principal component analysis, population structure and kinship showed no obvious relationships among these genotypes based on subgroups or regional sources. GWAS was carried out using a mixed linear model, and 159 significant associations were detected for these traits. Within the 25 kb linkage disequilibrium decay of the associated markers, NRT1/PTR FAMILY 6.3 (BVRB_5g097760); nudix hydrolase 15 (BVRB_8g182070) and TRANSPORT INHIBITOR RESPONSE 1 (BVRB_8g181550); transcription factor MYB77 (BVRB_2g023500); and ethylene-responsive transcription factor ERF014 (BVRB_1g000090) were predicted to be strongly associated with the taproot traits of root groove depth (RGD); root shape (RS); crown size (CS); and flesh colour (FC), respectively. For the aboveground traits, UDP-glycosyltransferase 79B6 (BVRB_9g223780) and NAC domain-containing protein 7 (BVRB_5g097990); F-box protein At1g10780 (BVRB_6g140760); phosphate transporter PHO1 (BVRB_3g048660); F-box protein CPR1 (BVRB_8g181140); and transcription factor MYB77 (BVRB_2g023500) and alcohol acyltransferase 9 (BVRB_2g023460) might be associated with the hypocotyl colour (HC); plant type (PT); petiole length (PL); cotyledon size (C); and fascicled leaf type (FLT) of sugar beet, respectively. AP-2 complex subunit mu (BVRB_5g106130), trihelix transcription factor ASIL2 (BVRB_2g041790) and late embryogenesis abundant protein 18 (BVRB_5g106150) might be involved in pollen quantity (PQ) variation. The candidate genes extensively participated in hormone response, nitrogen and phosphorus transportation, secondary metabolism, fertilization and embryo maturation.

CONCLUSIONS

The genetic basis of agronomical traits is complicated in heterozygous diploid sugar beet. The putative valuable genes found in this study will help further elucidate the molecular mechanism of each phenotypic trait for beet breeding.

Exploring the novel SNPs in neuroticism and birth weight based on GWAS datasets

OBJECTIVES

Epidemiological studies have confirmed that low birth weight (BW) is related to neuroticism and they may have a common genetic mechanism based on phenotypic correlation research. We conducted our study on a European population with 159,208 neuroticism and 289,142 birth weight samples. In this study, we aimed to identify new neuroticism single nucleotide polymorphisms (SNPs) and pleiotropic SNPs associated with neuroticism and BW and to provide more theoretical basis for the pathogenesis of the disease.

METHODS

We estimated the pleiotropic enrichment between neuroticism and BW in two independent Genome-wide association studies (GWAS) when the statistical thresholds were Conditional False Discovery Rate (cFDR) < 0.01 and Conjunctional Conditional False Discovery Rate (ccFDR) < 0.05. We performed gene annotation and gene functional analysis on the selected significant SNPs to determine the biological role of gene function and pathogenesis. Two-sample Mendelian Randomization (TSMR) analysis was performed to explore the causal relationship between the neuroticism and BW.

RESULTS

The conditional quantile-quantile plots (Q-Q plot) indicated that neuroticism and BW have strong genetic pleiotropy enrichment trends. With the threshold of cFDR < 0.001, we identified 126 SNPs related to neuroticism and 172 SNPs related to BW. With the threshold of ccFDR < 0.05, we identified 62 SNPs related to both neuroticism and BW. Among these SNPs, rs8039305 and rs35755513 have eQTL (expressed quantitative trait loci) and meQTL (methylation quantitative trait loci) effects simultaneously. Through GO enrichment analysis we also found that the two pathways of positive regulation of "mesenchymal cell proliferation" and "DNA-binding transcription factor activity" were significantly enriched in neuroticism and BW. Mendelian randomization analysis results indicate that there is no obvious causal relationship between neuroticism and birth weight.

CONCLUSION

We found 126 SNPs related to neuroticism, 172 SNPs related to BW and 62 SNPs associated with both neuroticism and BW, which provided a theoretical basis for their genetic mechanism and novel potential targets for treatment/intervention development.

Association between intake of the n-3 polyunsaturated fatty acid docosahexaenoic acid (n-3 PUFA DHA) and reduced risk of ovarian cancer: A systematic Mendelian Randomization study

BACKGROUND & AIMS

Whether the intake of docosahexaenoic acid (DHA), an n-3 polyunsaturated fatty acid, is beneficial for ovarian cancer (OC) remains controversial and we hope to disentangle this puzzle using genetic data from large-scale populations in European and Asian.

METHODS

We employed, for the first time, a systematic Mendelian randomization (MR) design to comprehensively evaluate the causal effect of plasma DHA levels, an objective biomarker of DHA intake, on OC risk in European and then verified the extrapolation of the results in the Asian. Data in the analysis included genetic association data obtained from large-scale genome-wide association studies with 13,499 individuals for plasma DHA measurements and 66,450 individuals for OC in the European population, and 1361 individuals for plasma DHA measurements and 61,457 individuals for OC in the Asian population. The causal relationship between DHA and OC was estimated using the inverse-variance weighted approach, together with extensive validation and sensitivity analyses to verify the main results.

RESULTS

In the European population, MR evidence suggested a causal relationship between higher plasma DHA levels and lower OC risk (OR, 0.89 for OC per one-SD increment in DHA; 95% CI, 0.83 to 0.96; P = 0.003). Subgroup analysis by histological type of OC indicated that this observed association was stronger among endometrioid ovarian cancer (EOC) (OR, 0.82; 95% CI, 0.69 to 0.96; P = 0.014). A similar causal association of borderline significance was reached in the Asian replication set. The above results were consistently supported by a series of validation and sensitivity analyses.

CONCLUSION

Our study provided robust genetic evidence for a protective association between plasma DHA levels and lower risk of OC, especially EOC, in the European population. These findings may inform prevention strategies and interventions directed towards DHA intake and OC.

Dissecting the genetic heterogeneity of gastric cancer

BACKGROUND

Gastric cancer (GC) is clinically heterogenous according to location (cardia/non-cardia) and histopathology (diffuse/intestinal). We aimed to characterize the genetic risk architecture of GC according to its subtypes. Another aim was to examine whether cardia GC and oesophageal adenocarcinoma (OAC) and its precursor lesion Barrett's oesophagus (BO), which are all located at the gastro-oesophageal junction (GOJ), share polygenic risk architecture.

METHODS

We did a meta-analysis of ten European genome-wide association studies (GWAS) of GC and its subtypes. All patients had a histopathologically confirmed diagnosis of gastric adenocarcinoma. For the identification of risk genes among GWAS loci we did a transcriptome-wide association study (TWAS) and expression quantitative trait locus (eQTL) study from gastric corpus and antrum mucosa. To test whether cardia GC and OAC/BO share genetic aetiology we also used a European GWAS sample with OAC/BO.

FINDINGS

Our GWAS consisting of 5816 patients and 10,999 controls highlights the genetic heterogeneity of GC according to its subtypes. We newly identified two and replicated five GC risk loci, all of them with subtype-specific association. The gastric transcriptome data consisting of 361 corpus and 342 antrum mucosa samples revealed that an upregulated expression of MUC1, ANKRD50, PTGER4, and PSCA are plausible GC-pathomechanisms at four GWAS loci. At another risk locus, we found that the blood-group 0 exerts protective effects for non-cardia and diffuse GC, while blood-group A increases risk for both GC subtypes. Furthermore, our GWAS on cardia GC and OAC/BO (10,279 patients, 16,527 controls) showed that both cancer entities share genetic aetiology at the polygenic level and identified two new risk loci on the single-marker level.

INTERPRETATION

Our findings show that the pathophysiology of GC is genetically heterogenous according to location and histopathology. Moreover, our findings point to common molecular mechanisms underlying cardia GC and OAC/BO.

FUNDING

German Research Foundation (DFG).

Genetic association study identifies genetic variants for non-alcoholic fatty liver without comorbidities in the Korean population

BACKGROUND

Non-alcoholic fatty liver (NAFL) refers to a disease in which fat builds up in the liver, similar to what occurs for those who drink a lot of alcohol, even in cases of not drinking alcohol at all or only in a small amount. Along with non-alcoholic steatohepatitis (NASH), NAFL is a type of non-alcoholic fatty liver disease (NAFLD). Currently, the prevalence of NAFLD is increasing worldwide. A wide range of comorbidities that can increase the risk of NAFLD includes obesity, type 2 diabetes, dyslipidemia, and metabolic syndrome.

OBJECTIVE

This study aimed to discover genetic variants for NAFL in the Korean population.

METHODS

Differing from previous studies, we conducted a genome-wide association study for NAFL in the selected subjects without comorbidities to rule out bias due to the inclusion of confounding effects of comorbidities. We grouped 424 NAFL cases and 5,402 controls from the Korean Genome and Epidemiology Study (KoGES) subjects without comorbidities such as dyslipidemia, type 2 diabetes, and metabolic syndrome. All subjects including cases and controls did not consume alcohol at all, or consumed less than 20 g/day for men and less than 10 g/day for women.

RESULTS

The logistic association analysis adjusting for sex, age, BMI, and waist circumference identified one novel genome-wide significant variant (rs7996045, P = 2.3 × 10^-8) for NAFL. This variant was located in the intron of CLDN10 and was not detected using previous conventional approaches in which confounding effects resulting from comorbidities were not considered in the study design stage. In addition, we detected several genetic variants showing suggestive association for NAFL (P < 10^-5).

CONCLUSION

The unique strategy in our association analysis of excluding major confounding factors provides, for the first time, an insight into the genuine genetic basis influencing NAFL.

Genome-wide association study of leaf-related traits in tea plant in Guizhou based on genotyping-by-sequencing

BACKGROUND

Studying the genetic characteristics of tea plant (Camellia spp.) leaf traits is essential for improving yield and quality through breeding and selection. Guizhou Plateau, an important part of the original center of tea plants, has rich genetic resources. However, few studies have explored the associations between tea plant leaf traits and single nucleotide polymorphism (SNP) markers in Guizhou.

RESULTS

In this study, we used the genotyping-by-sequencing (GBS) method to identify 100,829 SNP markers from 338 accessions of tea germplasm in Guizhou Plateau, a region with rich genetic resources. We assessed population structure based on high-quality SNPs, constructed phylogenetic relationships, and performed genome-wide association studies (GWASs). Four inferred pure groups (G-I, G-II, G-III, and G-IV) and one inferred admixture group (G-V), were identified by a population structure analysis, and verified by principal component analyses and phylogenetic analyses. Through GWAS, we identified six candidate genes associated with four leaf traits, including mature leaf size, texture, color and shape. Specifically, two candidate genes, located on chromosomes 1 and 9, were significantly associated with mature leaf size, while two genes, located on chromosomes 8 and 11, were significantly associated with mature leaf texture. Additionally, two candidate genes, located on chromosomes 1 and 2 were identified as being associated with mature leaf color and mature leaf shape, respectively. We verified the expression level of two candidate genes was verified using reverse transcription quantitative polymerase chain reaction (RT-qPCR) and designed a derived cleaved amplified polymorphism (dCAPS) marker that co-segregated with mature leaf size, which could be used for marker-assisted selection (MAS) breeding in Camellia sinensis.

CONCLUSIONS

In the present study, by using GWAS approaches with the 338 tea accessions population in Guizhou, we revealed a list of SNPs markers and candidate genes that were significantly associated with four leaf traits. This work provides theoretical and practical basis for the genetic breeding of related traits in tea plant leaves.

Genetic variation in BnGRP1 contributes to low phosphorus tolerance in Brassica napus

Lack of phosphorus (P) is a major environmental factor affecting rapeseed (Brassica napus. L) root growth and development. For breeding purposes, it is crucial to identify the molecular mechanisms of root system architecture (RSA) traits underlying low P tolerance in rapeseed. The natural variations in the glycine-rich protein gene, BnGRP1, were analyzed in the natural population of 400 rapeseed cultivars under low P stress through genome-wide association study (GWAS) and transcriptome analyses. Based on 11 SNP mutations in BnGRP1 sequence, ten types of haplotypes (Hap) were formed. Compared with the other types, the cultivar of BnGRP1Hap1 type in the panel demonstrated the longest root length and heaviest root weight. BnGRP1Hap1 overexpression in rapeseed depicted the ability to enhance its resistance in response to low P tolerance. CRISPR/Cas9-derived BnGRP1Hap4 knockout mutations in rapeseed can lead to sensitivity to low P stress. Furthermore, BnGRP1Hap1 influences the expression of phosphate transporter 1 (PHT1) genes associated with P absorption. Overall, the findings of this study highlight new mechanisms of GRP1 genes in enhancing low P tolerance in rapeseed.

Genome-Wide Association Studies (GWAS)

Most of the breeding targets are quantitative traits. In exploring the quantitative trait locus (QTL) system of a trait, linkage mapping was established using sparse polymerase chain reaction (PCR) markers. With the genome-wide sequencing technology advanced, genome-wide association study (GWAS) was developed for natural (germplasm) populations using dense genomic markers, which facilitates the identification of the complete QTL system with their multiple alleles on genomic locations. GWAS makes use of the linkage disequilibrium (LD) due to historical saturate recombination and high-density genomic markers to detect QTLs through statistical test for the association between molecular markers and phenotypes. However, due to inbreeding and mixture of source populations, the germplasm population often has complex and unknown structure, which leads to false positives/negatives in GWAS. Various GWAS methods have been proposed to reduce false positives/negatives, including those of the general linear model and the mixed linear model, which focused mainly on finding a handful of major QTLs under single-locus model for major gene cloning and could not detect directly the multiple alleles using bi-allelic single-nucleotide polymorphism (SNP) marker. As a relatively thorough detection of QTLs with their multiple alleles is required for germplasm population, the restricted two-stage multi-locus multi-allele GWAS (RTM-GWAS) procedure was proposed for identifying the QTL system with varying multiple alleles. From the RTM-GWAS results, a QTL-allele matrix is constructed as a compact form of the population genetic constitution, which can be further used for crop genetic and breeding studies, including major gene mining, population evolution, and breeding by genetic design.

Genome-wide association study of quality traits and starch pasting properties of maize kernels

BACKGROUND

Starch are the main nutritional components of maize (Zea mays L.), and starch pasting properties are widely used as essential indicators for quality estimation. Based on the previous studies, various genes related to pasting properties have been identified in maize. However, the loci underlying variations in starch pasting properties in maize inbred lines remain to be identified.

RESULTS

To investigate the genetic architecture of these traits, the starch pasting properties were examined based on 292 maize inbred lines, which were genotyped with the MaizeSNP50 BeadChip composed of 55,126 evenly spaced, random SNPs. A genome-wide association study (GWAS) implemented in the software package FarmCPU was employed to identify genomic loci for the starch pasting properties. 48 SNPs were found to be associated with pasting properties. Moreover, 37 candidate genes were correlated with pasting properties. Among the candidate genes, GRMZM2G143646 and GRMZM2G166407 were associated with breakdown and final viscosity significantly, and both genes encode PPR (Pentatricopeptide repeat) protein. We used GWAS to explore candidate genes of maize starch pasting properties in this study. The identified candidate genes will be useful for further understanding of the genetic architecture of starch pasting properties in maize.

CONCLUSION

This study showed a complex regulation network about maize quality trait and starch pasting properties. It may provide some useful markers for marker assisted selection and a basis for cloning the genes behind these SNPs.

Transcriptome-wide association study identifies novel genes associated with bone mineral density and lean body mass in children

OBJECTIVE

To identify novel candidate genes whose expression is associated with bone mineral density (BMD) and body lean mass (LM) in children.

METHODS

A tissue-specific transcriptome-wide association study (TWAS) was conducted utilizing a large-scale genome-wide association study (GWAS) dataset associated with BMD and LM and involving 10,414 participants. The measurement of BMD and LM phenotypes was made based on total-body dual-energy X-ray absorptiometry (TB-DXA) scans. TWAS was conducted by using FUSION software. Reference panels for muscle skeleton (MS), peripheral blood (NBL) and whole blood (YBL) were used for TWAS analysis. Functional enrichment and protein-protein interaction (PPI) analyses of the genes identified by TWAS were performed by using the online tool Metascape ( http://metascape.org ).

RESULTS

For BMD, we identified 174 genes with P < 0.05, such as IKZF1 (P = 1.46 × 10^-9) and CHKB (P = 8.31 × 10^-7). For LM, we identified 208 genes with P < 0.05, such as COPS5 (P = 3.03 × 10^-12) and MRPS33 (P = 5.45 × 10^-10). Gene ontology (GO) enrichment analysis of the BMD-associated genes revealed 200 GO terms, such as protein catabolic process (Log P = -5.09) and steroid hormone-mediated signaling pathway (Log P = -3.13). GO enrichment analysis of the LM-associated genes detected 287 GO terms, such as the apoptotic signaling pathway (Log P = -8.08) and lipid storage (Log P = -3.55).

CONCLUSION

This study identified several candidate genes for BMD and LM in children, providing novel clues to the genetic mechanisms underlying the development of childhood BMD and LM.

Single-Nuclei RNA Sequencing of 5 Regions of the Human Prenatal Brain Implicates Developing Neuron Populations in Genetic Risk for Schizophrenia

BACKGROUND

While a variety of evidence supports a prenatal component in schizophrenia, there are few data regarding the cell populations involved. We sought to identify cells of the human prenatal brain mediating genetic risk for schizophrenia by integrating cell-specific gene expression measures generated through single-nuclei RNA sequencing with recent large-scale genome-wide association study (GWAS) and exome sequencing data for the condition.

METHODS

Single-nuclei RNA sequencing was performed on 5 brain regions (frontal cortex, ganglionic eminence, hippocampus, thalamus, and cerebellum) from 3 fetuses from the second trimester of gestation. Enrichment of schizophrenia common variant genetic liability and rare damaging coding variation was assessed in relation to gene expression specificity within each identified cell population.

RESULTS

Common risk variants were prominently enriched within genes with high expression specificity for developing neuron populations within the frontal cortex, ganglionic eminence, and hippocampus. Enrichments were largely independent of genes expressed in neuronal populations of the adult brain that have been implicated in schizophrenia through the same methods. Genes containing an excess of rare damaging variants in schizophrenia had higher expression specificity for developing glutamatergic neurons of the frontal cortex and hippocampus that were also enriched for common variant liability.

CONCLUSIONS

We found evidence for a distinct contribution of prenatal neuronal development to genetic risk for schizophrenia, involving specific populations of developing neurons within the second-trimester fetal brain. Our study significantly advances the understanding of the neurodevelopmental origins of schizophrenia and provides a resource with which to investigate the prenatal antecedents of other psychiatric and neurologic disorders.

Cellular Genome-wide Association Study Identifies Common Genetic Variation Influencing Lithium-Induced Neural Progenitor Proliferation

BACKGROUND

Bipolar disorder is a highly heritable neuropsychiatric condition affecting more than 1% of the human population. Lithium salts are commonly prescribed as a mood stabilizer for individuals with bipolar disorder. Lithium is clinically effective in approximately half of treated individuals, and their genetic backgrounds are known to influence treatment outcomes. While the mechanism of lithium's therapeutic action is unclear, it stimulates adult neural progenitor cell proliferation, similar to some antidepressant drugs.

METHODS

To identify common genetic variants that modulate lithium-induced proliferation, we conducted an EdU incorporation assay in a library of 80 genotyped human neural progenitor cells treated with lithium. These data were used to perform a genome-wide association study to identify common genetic variants that influence lithium-induced neural progenitor cell proliferation. We manipulated the expression of a putatively causal gene using CRISPRi/a (clustered regularly interspaced short palindromic repeats interference/activation) constructs to experimentally verify lithium-induced proliferation effects.

RESULTS

We identified a locus on chr3p21.1 associated with lithium-induced proliferation. This locus is also associated with bipolar disorder risk, schizophrenia risk, and interindividual differences in intelligence. We identified a single gene, GNL3, whose expression temporally increased in an allele-specific fashion following lithium treatment. Experimentally increasing the expression of GNL3 led to increased proliferation under baseline conditions, while experimentally decreasing GNL3 expression suppressed lithium-induced proliferation.

CONCLUSIONS

Our experiments reveal that common genetic variation modulates lithium-induced neural progenitor proliferation and that GNL3 expression is necessary for the full proliferation-stimulating effects of lithium. These results suggest that performing genome-wide associations in genetically diverse human cell lines is a useful approach to discover context-specific pharmacogenomic effects.

Exome-wide variation in a diverse barley panel reveals genetic associations with ten agronomic traits in Eastern landraces

Barley (Hordeum vulgare ssp. vulgare) was one of the first crops to be domesticated and is adapted to a wide range of environments. Worldwide barley germplasm collections possess valuable allelic variations that could further improve barley productivity. Although barley genomics has offered a global picture of allelic variation among varieties and its association with various agronomic traits, polymorphisms from East Asian varieties remain scarce. In this study, we analyzed exome polymorphisms in a panel of 274 barley varieties collected worldwide, including 137 varieties from East Asian countries and Ethiopia. We revealed the underlying population structure and conducted genome-wide association studies for ten agronomic traits. Moreover, we examined genome-wide associations for traits related to grain size such as awn length and glume length. Our results demonstrate the value of diverse barley germplasm panels containing Eastern varieties, highlighting their distinct genomic signatures relative to Western subpopulations.

The genetic structure of pain in depression patients: A genome-wide association study and proteome-wide association study

BACKGROUND

Comparing with the general population, the pain in depression patients has more complex biological mechanism. We aim to explore the etiological mechanism of pain in depression patients from the perspective of genetics.

METHODS

Utilizing the UK Biobank samples with self-reported depression status or PHQ score ≥10, we conducted genome-wide association studies (GWAS) of seven pain traits (N = 1,133-58,349). Additionally, we used FUSION pipeline to perform proteome-wide association study (PWAS) and transcriptome-wide association study (TWAS) by integrating GWAS summary data with two different proteome reference weights (ROS/MAP and Banner) and Rnaseq gene expression reference weights, respectively.

RESULTS

GWAS identified 3 significant genes associated with different pain traits in depression patients, including TRIOBP (P_GWAS = 4.48 × 10^-8) for stomach or abdominal pain, SLC9A9(P_GWAS = 2.77 × 10^-8) for multisite chronic pain (MCP) and ADGRF1 (P_GWAS = 1.51 × 10^-8) for neck or shoulder pain. In addition, PWAS and TWAS analysis also identified multiple candidate genes associated with different pain traits in depression patients, such as TPRG1L (P_PWAS-Banner = 3.38 × 10^-2) and SIRPA (P_PWAS-Banner = 3.65 × 10^-2) for MCP, etc. Notably, when comparing the results of PWAS and TWAS analysis, we found overlapping candidate genes in these pain traits, such as GSTM3 (P_PWAS-_Banner = 3.38 × 10^-2, P_TWAS = 6.92 × 10^-3) in the stomach or abdominal pain phenotype, ATG7 (P_PWAS-_Rosmap = 3.15 × 10^-2, P_TWAS = 2.98 × 10^-2) in the MCP, etc. CONCLUSIONS: We identified multiple novel candidate genes for pain traits in depression patients from different perspectives of genetics, which provided novel clues for understanding the genetic mechanisms underlying the pain in depression patients.

Genome-wide association study for resistance in bread wheat (Triticum aestivum L.) to stripe rust (Puccinia striiformis f. sp. tritici) races in Argentina

BACKGROUND

Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most devastating diseases of the wheat crop. It causes significant reductions in both grain yield and grain quality. In recent years, new and more virulent races have overcome many of the known resistance genes in Argentinian germplasm. In order to identify loci conferring resistance to the local races of Pst for effective utilization in future breeding programs, a genome-wide association study (GWAS) was performed using a collection of 245 bread wheat lines genotyped with 90 K SNPs.

RESULTS

To search for adult plant resistance (APR) the panel was evaluated for disease severity (DS) and area under disease progress curve (AUDPC) in field trials during two years under natural infection conditions. To look for seedling or all-stage resistance (ASR) the panel was evaluated to determine infection type (IT) under greenhouse conditions against two prevalent races in Argentina. The phenotypic data showed that the panel possessed enough genetic variability for searching for sources of resistance to Pst. Significant correlations between years were observed for Pst response in the field and high heritability values were found for DS (H² = 0.89) and AUDPC (H² = 0.93). Based on GWAS, eight markers associated with Pst resistance (FDR < 0.01) were identified, of these, five were associated with ASR (on chromosomes 1B, 2A, 3A and 5B) and three with APR (on chromosomes 3B and 7A). These markers explained between 2% and 32.62% of the phenotypic variation. Five of the markers corresponded with previously reported Yr genes/QTL, while the other three (QYr.Bce.1B.sd.1, QYr.Bce.3A.sd and QYr.Bce.3B.APR.2) might be novel resistance loci.

CONCLUSION

Our results revealed high genetic variation for resistance to Argentinian stripe rust races in the germplasm used here. It constitutes a very promising step towards the improvement of Pst resistance of bread wheat in Argentina. Also, the identification of new resistance loci would represent a substantial advance for diversifying the current set of resistance genes and to advance in the improvement of the durable resistance to the disease.

rs2013278 in the multiple immunological-trait susceptibility locus CD28 regulates the production of non-functional splicing isoforms

Background

Ligation of CD28 with ligands such as CD80 or CD86 provides a critical second signal alongside antigen presentation by class II major histocompatibility complex expressed on antigen-presenting cells through the T cell antigen receptor for naïve T cell activation. A number of studies suggested that CD28 plays an important role in the pathogenesis of various human diseases. Recent genome-wide association studies (GWASs) identified CD28 as a susceptibility locus for lymphocyte and eosinophil counts, multiple sclerosis, ulcerative colitis, celiac disease, rheumatoid arthritis, asthma, and primary biliary cholangitis. However, the primary functional variant and molecular mechanisms of disease susceptibility in this locus remain to be elucidated. This study aimed to identify the primary functional variant from thousands of genetic variants in the CD28 locus and elucidate its functional effect on the CD28 molecule.

Results

Among the genetic variants exhibiting stronger linkage disequilibrium (LD) with all GWAS-lead variants in the CD28 locus, rs2013278, located in the Rbfox binding motif related to splicing regulation, was identified as a primary functional variant related to multiple immunological traits. Relative endogenous expression levels of CD28 splicing isoforms (CD28i and CD28Δex2) compared with full-length CD28 in allele knock-in cell lines generated using CRISPR/Cas9 were directly regulated by rs2013278 (P < 0.05). Although full-length CD28 protein expressed on Jurkat T cells showed higher binding affinity for CD80/CD86, both CD28i and CD28Δex2 encoded loss-of-function isoforms.

Conclusion

The present study demonstrated for the first time that CD28 has a shared disease-related primary functional variant (i.e., rs2013278) that regulates the CD28 alternative splicing that generates loss-of-function isoforms. They reduce disease risk by inducing anergy of effector T cells that over-react to autoantigens and allergens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-022-00419-7.

Genome-wide association study of aromatase inhibitor discontinuation due to musculoskeletal symptoms

OBJECTIVE

Aromatase inhibitors (AIs) are commonly used to treat hormone receptor positive (HR +) breast cancer. AI-induced musculoskeletal syndrome (AIMSS) is a common toxicity that causes AI treatment discontinuation. The objective of this genome-wide association study (GWAS) was to identify genetic variants associated with discontinuation of AI therapy due to AIMSS and attempt to replicate previously reported associations.

METHODS

In the Exemestane and Letrozole Pharmacogenetics (ELPh) study, postmenopausal patients with HR + non-metastatic breast cancer were randomized to letrozole or exemestane. Genome-wide genotyping of germline DNA was conducted followed by imputation. Each imputed variant was tested for association with time-to-treatment discontinuation due to AIMSS using a Cox proportional hazards model assuming additive genetic effects and adjusting for age, baseline pain score, prior taxane treatment, and AI arm. Secondary analyses were conducted within each AI arm and analyses of candidate variants previously reported to be associated with AIMSS risk.

RESULTS

Four hundred ELPh participants were included in the combined analysis. Two variants surpassed the genome-wide significance level in the primary analysis (p value < 5 × 10-8), an intronic variant (rs79048288) within CCDC148 (HR = 4.42, 95% CI: 2.67-7.33) and an intergenic variant (rs912571) upstream of PPP1R14C (HR = 0.30, 95% CI: 0.20-0.47). In the secondary analysis, rs74418677, which is known to be associated with expression of SUPT20H, was significantly associated with discontinuation of letrozole therapy due to AIMSS (HR = 5.91, 95% CI: 3.16-11.06). We were able to replicate associations for candidate variants previously reported to be associated with AIMSS in this cohort, but were not able to replicate associations for any other variants previously reported in other patient cohorts.

CONCLUSIONS

Our GWAS findings identify several candidate variants that may be associated with AIMSS risk from AI generally or letrozole specifically. Validation of these associations in independent cohorts is needed before translating these findings into clinical practice to improve treatment outcomes in patients with HR + breast cancer.

Post-traumatic stress disorder in the Canadian Longitudinal Study on Aging: A genome-wide association study

OBJECTIVE

Canada exhibits one of highest lifetime prevalence for post-traumatic stress disorder (PTSD), but the etiology of this debilitating mental health condition still remains largely unknown. This study aims to examine the genetics of PTSD in the Canadian Longitudinal Study on Aging (CLSA) to identify potential genetic factors involved in the development of PTSD.

METHOD

The CLSA sample was screened for primary (PTSD status) and secondary outcomes (avoidance, detachment, guardedness, and nightmares) based on the Primary Care PTSD Screen Scale (PC-PTSD). After GWAS quality control and whole-genome imputation, single-marker, gene-based, and polygenic risk score (PRS) analyses were performed.

RESULTS

Based on available genotype and phenotype data, N = 16,535 individuals were selected for the analyses. While genome-wide analyses did not show significant findings for our primary and secondary outcomes, PRS analyses showed variable levels of association between PC-PTSD items with trauma, major depressive disorder, schizophrenia, bipolar disorder, educational attainment, and insomnia (p &lt; 5e-4).

CONCLUSION

This is the first GWAS of PTSD status and individual PC-PTSD items in a population sample of older adults from Canada. This study was also able to replicate findings from previous studies. Genetic investigations into individual symptom components of PTSD may help untangle the complex genetic architecture of PTSD.

Patterns of genetic variation and QTLs controlling grain traits in a collection of global wheat germplasm revealed by high-quality SNP markers

BACKGROUND

Establish a molecular breeding program involved assembling a diverse germplasm collection and generating accurate genotypes to characterize their genetic potential and associate them with agronomic traits. In this study, we acquired over eight hundred wheat accessions from international gene banks and assessed their genetic relatedness using high-quality SNP genotypes. Understanding the scope of genomic variation in this collection allows the breeders to utilize the genetic resources efficiently while improving wheat yield and quality.

RESULTS

A wheat diversity panel comprising 39 durum wheat, 60 spelt wheat, and 765 bread wheat accessions was genotyped on iSelect 90 K wheat SNP arrays. A total of 57,398 SNP markers were mapped to IWGSC RefSeq v2.1 assembly, over 30,000 polymorphic SNPs in the A, B, D genomes were used to analyze population structure and diversity, the results revealed the separation of the three species and the differentiation of CIMMYT improved breeding lines and landraces or widely grown cultivars. In addition, several chromosomal regions under selection were detected. A subset of 280 bread wheat accessions was evaluated for grain traits, including grain length, width, surface area, and color. Genome-wide association studies (GWAS) revealed that several chromosomal regions were significantly linked to known quantitative trait loci (QTL) controlling grain-related traits. One of the SNP peaks at the end of chromosome 7A was in strong linkage disequilibrium (LD) with WAPO-A1, a gene that governs yield components.

CONCLUSIONS

Here, the most updated and accurate physical positions of SNPs on 90 K genotyping array are provided for the first time. The diverse germplasm collection and associated genotypes are available for the wheat researchers to use in their molecular breeding program. We expect these resources to broaden the genetic basis of original breeding and pre-breeding materials and ultimately identify molecular markers associated with important agronomic traits which are evaluated in diverse environmental conditions.

SSR-based genome-wide association study in turkish durum wheat germplasms revealed novel QTL of accumulated platinum

BACKGROUND

Durum wheat has a genetic capacity to accumulate toxic metals that can exceed the safety limit of the international standards, which may seriously affect human health. Identifying germplasms with low, nontoxic accumulated metal contents is important to select and develop new varieties. Thus, the objective of this study is to identify the levels of accumulated platinum in durum wheat and detect novel QTL.

METHODS AND RESULTS

Platinum contents were determined using 130 durum genotypes. Results generally showed low values of accumulated Pt and significantly less than the maximum grain's Pt content determined by international standards. Pt contents among genotypes varied from ≤ 0.001 to 0.72 µg/kg with an average of 0.02. Landraces showed the lowest average accumulated Pt. GWAS was then performed with 780 SSR markers. Five QTL were detected and explained 14.4-23.1% of the total phenotypic variation. Chromosomes 3 A, 3B, and 5B appear to be hotspots and may play a crucial role in accumulated Pt and were harbored in 1, 3, and 1 QTL, respectively.

CONCLUSIONS

This assessment of accumulated Pt within a unique panel included accessions mostly from Turkish regions, and GWAS used is the first study regarding accumulated Pt indices to reveal novel QTL. It will allow breeders to accelerate their selection of proper genotypes according to desired alleles and offer an opportunity to apply MAS to minimize Pt toxicity in durum wheat. Results indicated that the significance of genome (B) regions are likely related to the inheritance control of Pt content and may play a pivotal role regarding durum wheat's Pt contents. Nonetheless, these novel QTL should be validated in independent populations in numerous environments.

Functional genomic analysis delineates regulatory mechanisms of GWAS-identified bipolar disorder risk variants

Background

Genome-wide association studies (GWASs) have identified multiple risk loci for bipolar disorder (BD). However, pinpointing functional (or causal) variants in the reported risk loci and elucidating their regulatory mechanisms remain challenging.

Methods

We first integrated chromatin immunoprecipitation sequencing (ChIP-Seq) data from human brain tissues (or neuronal cell lines) and position weight matrix (PWM) data to identify functional single-nucleotide polymorphisms (SNPs). Then, we verified the regulatory effects of these transcription factor (TF) binding–disrupting SNPs (hereafter referred to as “functional SNPs”) through a series of experiments, including reporter gene assays, allele-specific expression (ASE) analysis, TF knockdown, CRISPR/Cas9-mediated genome editing, and expression quantitative trait loci (eQTL) analysis. Finally, we overexpressed PACS1 (whose expression was most significantly associated with the identified functional SNPs rs10896081 and rs3862386) in mouse primary cortical neurons to investigate if PACS1 affects dendritic spine density.

Results

We identified 16 functional SNPs (in 9 risk loci); these functional SNPs disrupted the binding of 7 TFs, for example, CTCF and REST binding was frequently disrupted. We then identified the potential target genes whose expression in the human brain was regulated by these functional SNPs through eQTL analysis. Of note, we showed dysregulation of some target genes of the identified TF binding–disrupting SNPs in BD patients compared with controls, and overexpression of PACS1 reduced the density of dendritic spines, revealing the possible biological mechanisms of these functional SNPs in BD.

Conclusions

Our study identifies functional SNPs in some reported risk loci and sheds light on the regulatory mechanisms of BD risk variants. Further functional characterization and mechanistic studies of these functional SNPs and candidate genes will help to elucidate BD pathogenesis and develop new therapeutic approaches and drugs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01057-3.

Comprehensive evaluation of Chinese peanut mini-mini core collection and QTL mapping for aflatoxin resistance

BACKGROUND

Aflatoxin contamination caused by Aspergillus fungi has been a serious factor affecting food safety of peanut (Arachis hypogaea L.) because aflatoxins are highly harmful for human and animal health. As three mechanisms of resistance to aflatoxin in peanut including shell infection resistance, seed infection resistance and aflatoxin production resistance exist among naturally evolved germplasm stocks, it is highly crucial to pyramid these three resistances for promoting peanut industry development and protecting consumers' health. However, less research effort has been made yet to investigate the differentiation and genetic relationship among the three resistances in diversified peanut germplasm collections.

RESULTS

In this study, the Chinese peanut mini-mini core collection selected from a large basic collection was systematically evaluated for the three resistances against A. flavus for the first time. The research revealed a wide variation among the diversified peanut accessions for all the three resistances. Totally, 14 resistant accessions were identified, including three with shell infection resistance, seven with seed infection resistance and five with aflatoxin production resistance. A special accession, Zh.h1312, was identified with both seed infection and aflatoxin production resistance. Among the five botanic types of A. hypogaea, the var. vulgaris (Spanish type) belonging to subspecies fastigiata is the only one which possessed all the three resistances. There was no close correlation between shell infection resistance and other two resistances, while there was a significant positive correlation between seed infection and toxin production resistance. All the three resistances had a significant negative correlation with pod or seed size. A total of 16 SNPs/InDels associated with the three resistances were identified through genome-wide association study (GWAS). Through comparative analysis, Zh.h1312 with seed infection resistance and aflatoxin production resistance was also revealed to possess all the resistance alleles of associated loci for seed infection index and aflatoxin content.

CONCLUSIONS

This study provided the first comprehensive understanding of differentiation of aflatoxin resistance in diversified peanut germplasm collection, and would further contribute to the genetic enhancement for resistance to aflatoxin contamination.

Experimental demonstration and pan-structurome prediction of climate-associated riboSNitches in Arabidopsis

Background

Genome-wide association studies (GWAS) aim to correlate phenotypic changes with genotypic variation. Upon transcription, single nucleotide variants (SNVs) may alter mRNA structure, with potential impacts on transcript stability, macromolecular interactions, and translation. However, plant genomes have not been assessed for the presence of these structure-altering polymorphisms or “riboSNitches.”

Results

We experimentally demonstrate the presence of riboSNitches in transcripts of two Arabidopsis genes, ZINC RIBBON 3 (ZR3) and COTTON GOLGI-RELATED 3 (CGR3), which are associated with continentality and temperature variation in the natural environment. These riboSNitches are also associated with differences in the abundance of their respective transcripts, implying a role in regulating the gene's expression in adaptation to local climate conditions. We then computationally predict riboSNitches transcriptome-wide in mRNAs of 879 naturally inbred Arabidopsis accessions. We characterize correlations between SNPs/riboSNitches in these accessions and 434 climate descriptors of their local environments, suggesting a role of these variants in local adaptation. We integrate this information in CLIMtools V2.0 and provide a new web resource, T-CLIM, that reveals associations between transcript abundance variation and local environmental variation.

Conclusion

We functionally validate two plant riboSNitches and, for the first time, demonstrate riboSNitch conditionality dependent on temperature, coining the term “conditional riboSNitch.” We provide the first pan-genome-wide prediction of riboSNitches in plants. We expand our previous CLIMtools web resource with riboSNitch information and with 1868 additional Arabidopsis genomes and 269 additional climate conditions, which will greatly facilitate in silico studies of natural genetic variation, its phenotypic consequences, and its role in local adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02656-4.

Screening of stable resistant accessions and identification of resistance loci to Barley yellow mosaic virus disease

Background

The disease caused by Barley yellow mosaic virus (BaYMV) infection is a serious threat to autumn-sown barley (Hordeum vulgare L.) production in Europe, East Asia and Iran. Due to the rapid diversification of BaYMV strains, it is urgent to discover novel germplasm and genes to assist breeding new varieties with resistance to different BaYMV strains, thus minimizing the effect of BaYMV disease on barley cropping.

Methods

A natural population consisting of 181 barley accessions with different levels of resistance to BaYMV disease was selected for field resistance identification in two separate locations (Yangzhou and Yancheng, Jiangsu Province, China). Additive main effects and multiplicative interaction (AMMI) analysis was used to identify accessions with stable resistance. Genome-wide association study (GWAS) of BaYMV disease resistance was broadly performed by combining both single nucleotide polymorphisms (SNPs) and specific molecular markers associated with the reported BaYMV disease resistance genes. Furthermore, the viral protein genome linked (VPg) sequences of the virus were amplified and analyzed to assess the differences between the BaYMV strains sourced from the different experimental sites.

Results

Seven barley accessions with lower standardized Area Under the Disease Progress Steps (sAUDPS) index in every environment were identified and shown to have stable resistance to BaYMV disease in each assessed location. Apart from the reported BaYMV disease resistance genes rym4 and rym5, one novel resistance locus explaining 24.21% of the phenotypic variation was identified at the Yangzhou testing site, while two other novel resistance loci that contributed 19.23% and 19.79% of the phenotypic variation were identified at the Yancheng testing site, respectively. Further analysis regarding the difference in the VPg sequence of the predominant strain of BaYMV collected from these two testing sites may explain the difference of resistance loci differentially identified under geographically distinct regions. Our research provides novel genetic resources and resistance loci for breeding barley varieties for BaMYV disease resistance.

Innovative computational approaches shed light on genetic mechanisms underlying cognitive impairment among children born extremely preterm

BACKGROUND

Although survival rates for infants born extremely preterm (gestation < 28 weeks) have improved significantly in recent decades, neurodevelopmental impairment remains a major concern. Children born extremely preterm remain at high risk for cognitive impairment from early childhood to adulthood. However, there is limited evidence on genetic factors associated with cognitive impairment in this population.

METHODS

First, we used a latent profile analysis (LPA) approach to characterize neurocognitive function at age 10 for children born extremely preterm. Children were classified into two groups: (1) no or low cognitive impairment, and (2) moderate-to-severe cognitive impairment. Second, we performed TOPMed-based genotype imputation on samples with genotype array data (n = 528). Third, we then conducted a genome-wide association study (GWAS) for LPA-inferred cognitive impairment. Finally, computational analysis was conducted to explore potential mechanisms underlying the variant x LPA association.

RESULTS

We identified two loci reaching genome-wide significance (p value < 5e-8): TEA domain transcription factor 4 (TEAD4 at rs11829294, p value = 2.40e-8) and syntaxin 18 (STX18 at rs79453226, p value = 1.91e-8). Integrative analysis with brain expression quantitative trait loci (eQTL), chromatin conformation, and epigenomic annotations suggests tetraspanin 9 (TSPAN9) and protein arginine methyltransferase 8 (PRMT8) as potential functional genes underlying the GWAS signal at the TEAD4 locus.

CONCLUSIONS

We conducted a novel computational analysis by utilizing an LPA-inferred phenotype with genetics data for the first time. This study suggests that rs11829294 and its LD buddies have potential regulatory roles on genes that could impact neurocognitive impairment for extreme preterm born children.