Genome-wide association analysis revealed new QTL and candidate genes affecting the teat number in Dutch Large White pigs

Teat number (TNUM) is an important reproductive trait of sows, which affects the weaning survival rate of piglets. In this study, 1166 Dutch Large White pigs with TNUM phenotype were used as the research object. These pigs were genotyped by 50K SNP chip and the chip data were further imputed to the resequencing level. The estimated heritabilities of left teat number (LTN), right teat number (RTN) and total teat number (TTN) were 0.21, 0.19 and 0.3, respectively. Based on chip data, significant SNPs for RTN on SSC2, SSC5, SSC9 and SSC13 were identified using genome-wide association analysis (GWAS). Significant SNPs for TTN were identified on SSC2, SSC5 and SSC7. Based on imputed data, the GWAS identified a significant SNP (rs329158522) for LTN on SSC17, two significant SNPs (rs342855242 and rs80813115) for RTN on SSC2 and SSC9, and two significant SNPs (rs327003548 and rs326943811) for TTN on SSC5 and SSC6. Among them, four novel QTL were discovered. The Bayesian fine-mapping method was used to fine map the QTL identified in the GWAS of the imputed data, and the confidence intervals of QTL affecting LTN (SSC17: 45.22-46.20 Mb), RTN (SSC9: 122.18-122.80 Mb) and TTN (SSC5: 14.01-15.91 Mb, SSC6: 120.06-121.25 Mb) were detected. A total of 52 candidate genes were obtained. Furthermore, we identified five candidate genes, WNT10B, AQP5, FMNL3, NUAK1 and CKAP4, for the first time, which involved in breast development and other related functions by gene annotation. Overall, this study provides new molecular markers for the breeding of teat number in pigs.

New insights into the genetic loci related to egg weight and age at first egg traits in broiler breeder

Egg weight (EW) and age at first egg (AFE) are economically important traits in breeder chicken production. The genetic basis of these traits, however, is far from understood, especially for broiler breeders. In this study, genetic parameter estimation, genome-wide association analysis, meta-analysis, and selective sweep analysis were carried out to identify genetic loci associated with EW and AFE in 6,842 broiler breeders. The study found that the heritability of EW ranged from 0.42 to 0.44, while the heritability of AFE was estimated at 0.33 in the maternal line. Meta-analysis and selective sweep analysis identified two colocalized regions on GGA4 that significantly influenced EW at 32 wk (EW32W) and at 43 wk (EW43W) with both paternal and maternal lines. The genes AR, YIPF6, and STARD8 were located within the significant region (GGA4: 366.86-575.50 kb), potentially affecting EW through the regulation of follicle development, cell proliferation, and lipid transfer etc. The promising genes LCORL and NCAPG were positioned within the significant region (GGA4:75.35-75.42 Mb), potentially influencing EW through pleiotropic effects on growth and development. Additionally, 3 significant regions were associated with AFE on chromosomes GGA7, GGA19, and GGA27. All of these factors affected the AFE by influencing ovarian development. In our study, the genomic information from both paternal and maternal lines was used to identify genetic regions associated with EW and AFE. Two genomic regions and eight genes were identified as the most likely candidates affecting EW and AFE. These findings contribute to a better understanding of the genetic basis of egg production traits in broiler breeders and provide new insights into future technology development.

Identification of Novel QTL for Mercury Accumulation in Maize Using an Enlarged SNP Panel

Mercury (Hg) pollution not only poses a threat to the environment but also adversely affects the growth and development of plants, with potential repercussions for animals and humans through bioaccumulation in the food chain. Maize, a crucial source of food, industrial materials, and livestock feed, requires special attention in understanding the genetic factors influencing mercury accumulation. Developing maize varieties with low mercury accumulation is vital for both maize production and human health. In this study, a comprehensive genome-wide association study (GWAS) was conducted using an enlarged SNP panel comprising 1.25 million single nucleotide polymorphisms (SNPs) in 230 maize inbred lines across three environments. The analysis identified 111 significant SNPs within 78 quantitative trait loci (QTL), involving 169 candidate genes under the Q model. Compared to the previous study, the increased marker density and optimized statistical model led to the discovery of 74 additional QTL, demonstrating improved statistical power. Gene ontology (GO) enrichment analysis revealed that most genes participate in arsenate reduction and stress responses. Notably, GRMZM2G440968, which has been reported in previous studies, is associated with the significant SNP chr6.S_155668107 in axis tissue. It encodes a cysteine proteinase inhibitor, implying its potential role in mitigating mercury toxicity by inhibiting cysteine. Haplotype analyses provided further insights, indicating that lines carrying hap3 exhibited the lowest mercury content compared to other haplotypes. In summary, our study significantly enhances the statistical power of GWAS, identifying additional genes related to mercury accumulation and metabolism. These findings offer valuable insights into unraveling the genetic basis of mercury content in maize and contribute to the development of maize varieties with low mercury accumulation.

A genome-wide association analysis for salt tolerance during the soybean germination stage and development of KASP markers

Salt stress poses a significant challenge to crop productivity, and understanding the genetic basis of salt tolerance is paramount for breeding resilient soybean varieties. In this study, a soybean natural population was evaluated for salt tolerance during the germination stage, focusing on key germination traits, including germination rate (GR), germination energy (GE), and germination index (GI). It was seen that under salt stress, obvious inhibitions were found on these traits, with GR, GE, and GI diminishing by 32% to 54% when compared to normal conditions. These traits displayed a coefficient of variation (31.81% to 50.6%) and a substantial generalized heritability (63.87% to 86.48%). Through GWAS, a total of 1841 significant single-nucleotide polymorphisms (SNPs) were identified to be associated with these traits, distributed across chromosome 2, 5, 6, and 20. Leveraging these significant association loci, 12 candidate genes were identified to be associated with essential functions in coordinating cellular responses, regulating osmotic stress, mitigating oxidative stress, clearing reactive oxygen species (ROS), and facilitating heavy metal ion transport - all of which are pivotal for plant development and stress tolerance. To validate the candidate genes, quantitative real-time polymerase chain reaction (qRT-PCR) analysis was conducted, revealing three highly expressed genes (Glyma.02G067700, Glyma.02G068900, and Glyma.02G070000) that play pivotal roles in plant growth, development, and osmoregulation. In addition, based on these SNPs related with salt tolerance, KASP (Kompetitive Allele-Specific PCR)markers were successfully designed to genotype soybean accessions. These findings provide insight into the genetic base of soybean salt tolerance and candidate genes for enhancing soybean breeding programs in this study.

Causal relationship between inflammatory bowel disease and erythema nodosum: A two-sample bidirectional Mendelian randomization study

BACKGROUND

Previous studies have demonstrated the coexistence of erythema nodosum (EN) and inflammatory bowel disease (IBD), while the exact etiology of the co-occurrence of the two disorders remains uncertain.

METHODS

A bidirectional two-sample Mendelian randomization (MR) design was employed to determine the causal link between EN and IBD. Genetic variations associated with Crohn's disease (CD) and ulcerative colitis (UC) were derived from accessible genome-wide association studies pertaining to European ancestry. The FinnGen database was used to find the genetic variations containing EN. In the forward model, IBD was identified as the exposure, whereas in the reverse model, EN was identified as the exposure. The causal link between IBD and EN was examined using a range of different analysis techniques, the primary one being the inverse variance weighted (IVW) method, including inverse variance weighted-fixed effects (IVW-FE) and inverse-variance weighted-multiplicative random effects (IVW-MRE). To strengthen the results, assessments of sensitivity, heterogeneity, and pleiotropy were also conducted.

RESULTS

MR results showed that IBD increased the risk of EN (IVW-MRE: OR = 1.242, 95% CI = 1.068-1.443, p = 0.005). Furthermore, there was a strong correlation found between CD and a higher risk of EN (IVW-FE: OR = 1.250, 95% CI = 1.119-1.396, p = 8.036 × 10-5 ). However, UC did not appear to be linked to EN (IVW-FE: OR = 1.104, 95% CI = 0.868-1.405, p = 0.421). The reverse MR analysis findings did not imply that EN was linked to IBD. Horizontal pleiotropy did not appear to exist, and the robustness of these findings was confirmed.

CONCLUSION

The current investigation found that in European populations, IBD and its subtype CD could raise the incidence of EN.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity

STUDY QUESTION

Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins?

SUMMARY ANSWER

We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3.

WHAT IS KNOWN ALREADY

The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures.

STUDY DESIGN, SIZE, DURATION

We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls).

PARTICIPANTS/MATERIALS, SETTING, METHODS

Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation.

MAIN RESULTS AND THE ROLE OF CHANCE

This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle.

LARGE SCALE DATA

The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/).

LIMITATIONS, REASONS FOR CAUTION

Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility.

WIDER IMPLICATIONS OF THE FINDINGS

About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility.

STUDY FUNDING/COMPETING INTEREST(S)

Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Skłodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests.

TRIAL REGISTRATION NUMBER

N/A.

Unravelling novel and pleiotropic genes for cannon bone circumference and bone mineral density in Yorkshire pigs

Leg weakness is a prevalent health condition in pig farms. The augmentation of cannon bone circumference and bone mineral density can effectively improve limb strength in pigs and alleviate leg weakness. This study measured forelimb cannon bone circumference (fCBC) and rear limb cannon bone circumference (rCBC) using an inelastic tapeline and rear limb metatarsal area bone mineral density (raBMD) using a dual-energy X-ray absorptiometry bone density scanner. The samples of Yorkshire castrated boars were genotyped using a 50K single-nucleotide polymorphism (SNP) array. The SNP-chip data were imputed to the level of whole-genome sequencing data (iWGS). This study used iWGS data to perform genome-wide association studies and identified novel significant SNPs associated with fCBC on SSC6, SSC12, and SSC13, rCBC on SSC12 and SSC14, and raBMD on SSC7. Based on the high phenotypic and genetic correlations between CBC and raBMD, multi-trait meta-analysis was performed to identify pleiotropic SNPs. A significant potential pleiotropic quantitative trait locus (QTL) regulating both CBC and raBMD was identified on SSC15. Bayes fine mapping was used to establish the confidence intervals for these novel QTLs with the most refined confidence interval narrowed down to 56 kb (15.11 to 15.17 Mb on SSC12 for fCBC). Furthermore, the confidence interval for the potential pleiotropic QTL on SSC15 in the meta-analysis was narrowed down to 7.45 kb (137.55 to137.56 Mb on SSC15). Based on the biological functions of genes, the following genes were identified as novel regulatory candidates for different phenotypes: DDX42, MYSM1, FTSJ3, and MECOM for fCBC; SMURF2, and STC1 for rCBC; RGMA for raBMD. Additionally, RAMP1, which was determined to be located 23.68 kb upstream of the confidence interval of the QTL on SSC15 in the meta-analysis, was identified as a potential pleiotropic candidate gene regulating both CBC and raBMD. These findings offered valuable insights for identifying pathogenic genes and elucidating the genetic mechanisms underlying CBC and BMD.

Analysis of the Genetic Diversity in Tea Plant Germplasm in Fujian Province Based on Restriction Site-Associated DNA Sequencing

Fujian province, an important tea-producing area in China, has abundant tea cultivars. To investigate the genetic relationships of tea plant cultivars in Fujian province and the characteristics of the tea plant varieties, a total of 70 tea cultivars from Fujian and other 12 provinces in China were subjected to restriction site-associated DNA sequencing (RAD-seq). A total of 60,258,975 single nucleotide polymorphism (SNP) sites were obtained. These 70 tea plant cultivars were divided into three groups based on analyzing the phylogenetic tree, principal component, and population structure. Selection pressure analysis indicated that nucleotide diversity was high in Southern China and genetically distinct from cultivars of Fujian tea plant cultivars, according to selection pressure analysis. The selected genes have significant enrichment in pathways associated with metabolism, photosynthesis, and respiration. There were ten characteristic volatiles screened by gas chromatography-mass spectrometry (GC-MS) coupled with multivariate statistical methods, among which the differences in the contents of methyl salicylate, 3-carene, cis-3-hexen-1-ol, (E)-4-hexen-1-ol, and 3-methylbutyraldehyde can be used as reference indicators of the geographical distribution of tea plants. Furthermore, a metabolome genome-wide association study (mGWAS) revealed that 438 candidate genes were related to the aroma metabolic pathway. Further analysis showed that 31 genes of all the selected genes were screened and revealed the reasons for the genetic differences in aroma among tea plant cultivars in Fujian and Southern China. These results reveal the genetic diversity in the Fujian tea plants as well as a theoretical basis for the conservation, development, and utilization of the Fujian highly aromatic tea plant cultivars.

Genetic control of flowering time and fruit yield in citron watermelon

Flowering time and fruit yield are important traits in watermelon crop improvement. There is limited information on the inheritance and genomic loci underlying flowering time and yield performance, especially in citron watermelon. A total of 125 citron watermelon accessions were evaluated in field trials over two growing seasons for days to male and female flowers, fruit count, fruit weight, and fruit yield. The germplasm was genotyped with more than two million single-nucleotide polymorphism (SNP) markers generated via whole-genome resequencing. Trait mapping was conducted using a genome-wide association study (GWAS). Broad-sense heritability for all traits ranged from moderate to high, indicating that genetic improvement through breeding and selection is feasible. Significant marker-trait associations were uncovered for days to female flower (chromosomes Ca04, Ca05, Ca08, and Ca09), fruit count (on Ca02, Ca03, and Ca05), fruit weight (on Ca02, Ca06, Ca08, Ca10, and Ca11), and fruit yield on chromosomes Ca05, Ca07, and Ca09. The phenotypic variation explained by the significant SNPs ranged from 1.6 to 25.4, highlighting the complex genetic architecture of the evaluated traits. Candidate genes relevant to flowering time and fruit yield component traits were uncovered on chromosomes Ca02, Ca04, Ca05, Ca06, Ca09, and Ca11. These results lay a foundation for marker-assisted trait introgression of flowering time and fruit yield component traits in watermelons.

Mendelian randomization study and meta-analysis exploring the causality of age at menarche and the risk of intracerebral hemorrhage and ischemic stroke

BACKGROUND

The relationship between the age at menarche (AAM) and the risk of intracerebral hemorrhage (ICH) and ischemic stroke (IS) is still up for debate. The purpose of this study was to investigate potential causal connections between them.

METHODS

Genome-wide association analysis (GWAS) of AAM conducted by the MRC-IEU consortium was utilized for association analyses of ICH and IS by two-sample Mendelian randomization (MR) study. AAM data of the within-family GWAS consortium were used as replication phase data to verify the causal relationship between each other. Inverse variance weighting (IVW) method was the primary method used in this MR study. For additional proof, the weighted median estimation, MR-Egger regression, MR-PRESSO test, and MR-Robust Adjusted Profile Score evaluation were performed. The Cochran's Q test and the MR-PRESSO global test were used, respectively, to examine the sensitivity and pleiotropy. Random effects meta-analysis was utilized to analyze the causal data from the two consortiums to further explore the causality between AAM and ICH, IS.

RESULTS

We found that the AAM was causally linked with the risk of ICH (OR = 0.48, 95% CI: 0.28-0.80, p = 0.006). On the contrary, the causal effect from AAM to IS (OR = 0.98, 95% CI: 0.91-1.06, p = 0.64) has not been confirmed. For all subtypes of ICH, we found that nonlobar intracerebral hemorrhage (NLICH, OR = 0.41, 95% CI: 0.23-0.75, p = 0.004) but not lobar intracerebral hemorrhage (LICH, OR = 0.65, 95% CI: 0.34-1.24, p = 0.19) was associated with AAM without surprise. Similarly, we used the within-family GWAS consortium data to explore causality and found that AAM may reduce the risk of ICH (OR = 0.78, 95% CI: 0.72-0.86, p = 9.5 × 10-8 ) and NLICH (OR = 0.68, 95% CI: 0.61-0.75, p = 3.4 × 10-13 ) by IVW methods, but is not related to IS (OR = 0.97, 95% CI: 0.93-1.02, p = 0.26). These findings are further supported by the meta-analysis. Both Cochran's Q test and the MR-PRESSO global test failed to detect the presence of sensitivity.

CONCLUSION

AAM and ICH, particularly NLICH, are causally related, but not LICH, IS, or its subtypes in European population.

Inflammatory bowel disease and the risk of intracerebral hemorrhage: A Mendelian randomization study and meta-analysis

BACKGROUND

The link between inflammatory bowel disease (IBD) and intracerebral hemorrhage (ICH) is still unclear.

AIMS

We conducted a Mendelian randomization research and meta-analysis to explore the impact of IBD and its subtypes (Crohn's disease [CD], ulcerative colitis [UC]) on the risk of ICH.

METHODS

Two large genome-wide association analysis studies of International Inflammatory Bowel Disease Genetics Consortium (IIBDGC) and International Stroke Genetics Consortium as exposure (IBD, UC, and CD) and outcome (ICH) in the initial stage. IBD, CD, UC GWAS data from the FinnGen consortium were adopted for the replication phase, and ultimately, the results of the initial stage and replication phase data were combined in a meta-analysis to evaluate the causal association between IBD and its subtypes and the risk of ICH.

RESULTS

In the initial stage, we found that in the IVW (odds ratio [OR] = 0.83, 95% confidence interval [CI]: 0.71-0.96, p = .01), MR-PRESSO (OR = 0.85, 95% CI: 0.75-0.97, p = .02) and MR.RAPS (OR = 0.86, 95% CI: 0.76-0.98, p = .02) method showed that UC is associated with the risk of ICH. The causal relationship between IBD, CD, and the risk of ICH cannot be found by the IVW method. IBD and its subtypes UC, CD, and risk of ICH cannot find the presence of heterogeneity and pleiotropy. In replication stage, IBD (OR = 0.74, 95% CI: 0.59-0.94, p = .0135) related to ICH, while the IVW approach did not establish a causal link in UC and CD. The meta-analysis still indicated that UC (OR = 0.83, 95% CI: 0.72-0.93, p < .05) would lessen the risk of ICH while the causality between IBD, CD, and ICH was unable to be established.

CONCLUSION

UC was causally related to ICH, but IBD and CD are not associated with ICH. The precise pathophysiological mechanism needs to be thoroughly investigated in more detail.

Uncovering the candidate genes related to sheep body weight using multi-trait genome-wide association analysis

In sheep, body weight is an economically important trait. This study sought to map genetic loci related to weaning weight and yearling weight. To this end, a single-trait and multi-trait genome-wide association study (GWAS) was performed using a high-density 600 K single nucleotide polymorphism (SNP) chip. The results showed that 43 and 56 SNPs were significantly associated with weaning weight and yearling weight, respectively. A region associated with both weaning and yearling traits (OARX: 6.74-7.04 Mb) was identified, suggesting that the same genes could play a role in regulating both these traits. This region was found to contain three genes (TBL1X, SHROOM2 and GPR143). The most significant SNP was Affx-281066395, located at 6.94 Mb (p = 1.70 × 10-17), corresponding to the SHROOM2 gene. We also identified 93 novel SNPs elated to sheep weight using multi-trait GWAS analysis. A new genomic region (OAR10: 76.04-77.23 Mb) with 22 significant SNPs were discovered. Combining transcriptomic data from multiple tissues and genomic data in sheep, we found the HINT1, ASB11 and GPR143 genes may involve in sheep body weight. So, multi-omic anlaysis is a valuable strategy identifying candidate genes related to body weight.

Use of GWAS analysis in deciphering the inability of barley seeds to germinate after hypoxia

This article comments on:

Gómez-Álvarez EM, Tondelli A, Nghi KN, Voloboeva V, Giordano G, Valè G, Perata P, Pucciariello C. 2023. The inability of barley to germinate after submergence depends on hypoxia-induced secondary dormancy. Journal of Experimental Botany 74, 4277–4289

The basic helix-loop-helix transcription factor gene, OsbHLH38, plays a key role in controlling rice salt tolerance

The plant hormone abscisic acid (ABA) is crucial for plant seed germination and abiotic stress tolerance. However, the association between ABA sensitivity and plant abiotic stress tolerance remains largely unknown. In this study, 436 rice accessions were assessed for their sensitivity to ABA during seed germination. The considerable diversity in ABA sensitivity among rice germplasm accessions was primarily reflected by the differentiation between the Xian (indica) and Geng (japonica) subspecies and between the upland-Geng and lowland-Geng ecotypes. The upland-Geng accessions were most sensitive to ABA. Genome-wide association analyses identified four major quantitative trait loci containing 21 candidate genes associated with ABA sensitivity of which a basic helix-loop-helix transcription factor gene, OsbHLH38, was the most important for ABA sensitivity. Comprehensive functional analyses using knockout and overexpression transgenic lines revealed that OsbHLH38 expression was responsive to multiple abiotic stresses. Overexpression of OsbHLH38 increased seedling salt tolerance, while knockout of OsbHLH38 increased sensitivity to salt stress. A salt-responsive transcription factor, OsDREB2A, interacted with OsbHLH38 and was directly regulated by OsbHLH38. Moreover, OsbHLH38 affected rice abiotic stress tolerance by mediating the expression of a large set of transporter genes of phytohormones, transcription factor genes, and many downstream genes with diverse functions, including photosynthesis, redox homeostasis, and abiotic stress responsiveness. These results demonstrated that OsbHLH38 is a key regulator in plant abiotic stress tolerance.

Mining Salt Tolerance SNP Loci and Prediction of Candidate Genes in the Rice Bud Stage by Genome-Wide Association Analysis

Mining salt tolerance genes is significant for breeding high-quality salt-tolerant rice varieties in order to improve the utilization of saline-alkaline land. In this study, 173 rice accessions were measured for their germination potential (GP), germination rate (GR), seedling length (SL), root length (RL), germination potential relative to salt damage rate (GPR), germination rate relative to salt damage rate (GRR), seedling length relative to salt damage rate (SLR), relative salt damage rate at the germination stage (RSD) and comprehensive relative salt damage rate in the early seedling stage (CRS) under normal and salt stress conditions. Genome-wide association analysis was performed with 1,322,884 high-quality SNPs obtained by resequencing. Eight quantitative trait loci (QTLs) related to salt tolerance traits at the germination stage were detected in 2020 and 2021. They were related to the GPR (qGPR2) and SLR (qSLR9), which were newly discovered in this study. Three genes were predicted as salt tolerance candidate genes: LOC_Os02g40664, LOC_Os02g40810, and LOC_Os09g28310. At present, marker-assisted selection (MAS) and gene-edited breeding are becoming more widespread. Our discovery of candidate genes provides a reference for research in this field. The elite alleles identified in this study may provide a molecular basis for cultivating salt-tolerant rice varieties.

Identification of genetic biomarkers associated with pharmacokinetics and pharmacodynamics of apixaban in Chinese healthy volunteers

BACKGROUND

Apixaban is a superior direct oral anticoagulant exihibiting interindividual variability in concentration and response in the real world. The present study aimed to identify genetic biomarkers associated with pharmacokinetics (PK) and pharmacodynamics (PD) of apixaban in healthy Chinese subjects.

METHODS

This multicenter study included 181 healthy Chinese adults taking a single dose of 2.5 mg or 5 mg apixaban and assessed their PK and PD parameters. Genome-wide single nucleotide polymorphism (SNP) genotyping was performed using the Affymetrix Axiom CBC_PMRA Array. Candidate gene association analysis and genome-wide association study were conducted to identify genes with a predictive value for PK and PD parameters of apixaban.

RESULTS

Several ABCG2 variants were associated with C_max and AUC_0-t of apixaban (p < 6.12 × 10^-5) and also presented significant differences of anti-Xa_3h activity and dPT_3h according to different ABCG2 genotypes (p < 0.05). Besides, ABLIM2 variants were found to be associated with PK characteristics and F13A1 and C3 variants were associated with PD characteristics of apixaban (p < 9.46 × 10^-8).

CONCLUSION

ABCG2 variants were found to be ideal genetic biomarkers for both PK and PD characteristics of apixaban. ABLIM2, F13A1 and C3 were identified as potential candidate genes associated with inter-individual variability of apixaban. This study was registered on ClinicalTrials.gov NCT03259399.

Genetic Architecture of Parkinson's Disease

Year 2022 marks 25 years since the first mutation in familial autosomal dominant Parkinson's disease was identified. Over the years, our understanding of the role of genetic factors in the pathogenesis of familial and idiopathic forms of Parkinson's disease has expanded significantly - a number of genes for the familial form of the disease have been identified, and DNA markers for an increased risk of developing its sporadic form have been found. But, despite all the success achieved, we are far from an accurate assessment of the contribution of genetic and, even more so, epigenetic factors to the disease development. The review summarizes the information accumulated to date on the genetic architecture of Parkinson's disease and formulates issues that need to be addressed, which are primarily related to the assessment of epigenetic factors in the disease pathogenesis.

Genetic Risk of Primary Aldosteronism and Its Contribution to Hypertension: A Cross-Ancestry Meta-Analysis of Genome-Wide Association Study

BACKGROUND

Hypertension imposes substantial health and economic burden worldwide. Primary aldosteronism (PA) is one of the most common causes of secondary hypertension, causing cardiovascular events at higher risk compared with essential hypertension. However, the germline genetic contribution to the susceptibility of PA has not been well elucidated.

METHOD

We conducted a genome-wide association analysis of PA in the Japanese population and a cross-ancestry meta-analysis combined with UK Biobank and FinnGen cohorts (816 PA cases and 425 239 controls) to identify genetic variants that contribute to PA susceptibility. We also performed a comparative analysis for the risk of 42 previously established blood pressure-associated variants between PA and hypertension with the adjustment of blood pressure.

RESULTS

In the Japanese genome-wide association study, we identified 10 loci that presented suggestive evidence for the association with the PA risk (P<1.0×10^-6). In the meta-analysis, we identified 5 genome-wide significant loci (1p13, 7p15, 11p15, 12q24, and 13q12; P<5.0×10^-8), including 3 of the suggested loci in the Japanese genome-wide association study. The strongest association was observed at rs3790604 (1p13), an intronic variant of WNT2B (odds ratio, 1.50 [95% CI, 1.33-1.69]; P=5.2×10^-11). We further identified 1 nearly genome-wide significant locus (8q24, CYP11B2), which presented a significant association in the gene-based test (P=7.2×10^-7). Of interest, all of these loci were known to be associated with blood pressure in previous studies, presumably because of the prevalence of PA among individuals with hypertension. This assumption was supported by the observation that they had a significantly higher risk effect on PA than on hypertension. We also revealed that 66.7% of the previously established blood pressure-associated variants had a higher risk effect for PA than for hypertension.

CONCLUSIONS

This study demonstrates the genome-wide evidence for a genetic predisposition to PA susceptibility in the cross-ancestry cohorts and its significant contribution to the genetic background of hypertension. The strongest association with the WNT2B variants reinforces the implication of the Wnt/β-catenin pathway in the PA pathogenesis.

Genetic variations in relation to bleeding and pharmacodynamics of dabigatran in Chinese patients with nonvalvular atrial fibrillation: A nationwide multicentre prospective cohort study

INTRODUCTION

To identify the potential factors responsible for the individual variability of dabigatran, we investigated the genetic variations associated with clinical outcomes and pharmacodynamics (PD) in Chinese patients with nonvalvular atrial fibrillation (NVAF).

MATERIALS AND METHODS

Chinese patients with NVAF taking dabigatran etexilate with therapeutic doses were enrolled. The primary (bleeding events) and secondary (thromboembolic and major adverse cardiac events) outcomes for a 2-year follow-up were evaluated. Peak and trough PD parameters (anti-FIIa activity, activated partial thromboplastin time and prothrombin time) were detected. Whole-exome sequencing, genome-wide sequencing and candidate gene association analyses were performed.

RESULTS

There were 170 patients with NVAF treated with dabigatran (110 mg twice daily) who were finally included. Two single-nucleotide polymorphisms (SNPs) were significantly related with bleeding, which include UBASH3B rs2276408 (odds ratio [OR] = 8.79, 95% confidence interval [CI]: 2.99-25.83, p = 7.77 × 10-5 at sixth month visit) and FBN2 rs3805625 (OR = 8.29, 95% CI: 2.87-23.89, p = 9.08 × 10-5 at 12th month visit), as well as with increased trends at other visits (p < .05). Furthermore, minor allele carriers of 16 new SNPs increased PD levels, and those of one new SNP decreased PD values (p < 1.0 × 10-5 ). Lastly, 33 new SNPs were found to be associated with bleeding and PD among 14 candidate genes. Unfortunately, the low number of secondary outcomes precluded further association analyses.

CONCLUSIONS

Genetic variations indeed affected bleeding and PD in Chinese patients with NVAF treated with dabigatran. The functions of these suggestive genes and SNPs might further be explored and verified in more in vivo and in vitro investigations.

SNP-Target Genes Interaction Perturbing the Cancer Risk in the Post-GWAS

Cancer ranks as the second leading cause of death worldwide, and, being a genetic disease, it is highly heritable. Over the past few decades, genome-wide association studies (GWAS) have identified many risk-associated loci harboring hundreds of single nucleotide polymorphisms (SNPs). Some of these cancer-associated SNPs have been revealed as causal, and the functional characterization of the mechanisms underlying the cancer risk association has been illuminated in some instances. In this review, based on the different positions of SNPs and their modes of action, we discuss the mechanisms underlying how SNPs regulate the expression of target genes to consequently affect tumorigenesis and the development of cancer.

Genome-wide association studies for soybean epicotyl length in two environments using 3VmrMLM

Germination of soybean seed is the imminent vital process after sowing. The status of plumular axis and radicle determine whether soybean seed can emerge normally. Epicotyl, an organ between cotyledons and first functional leaves, is essential for soybean seed germination, seedling growth and early morphogenesis. Epicotyl length (EL) is a quantitative trait controlled by multiple genes/QTLs. Here, the present study analyzes the phenotypic diversity and genetic basis of EL using 951 soybean improved cultivars and landraces from Asia, America, Europe and Africa. 3VmrMLM was used to analyze the associations between EL in 2016 and 2020 and 1,639,846 SNPs for the identification of QTNs and QTN-by-environment interactions (QEIs)".A total of 180 QTNs and QEIs associated with EL were detected. Among them, 74 QTNs (ELS_Q) and 16 QEIs (ELS_QE) were identified to be associated with ELS (epicotyl length of single plant emergence), and 60 QTNs (ELT_Q) and 30 QEIs (ELT_QE) were identified to be associated with ELT (epicotyl length of three seedlings). Based on transcript abundance analysis, GO (Gene Ontology) enrichment and haplotype analysis, ten candidate genes were predicted within nine genic SNPs located in introns, upstream or downstream, which were supposed to be directly or indirectly involved in the process of seed germination and seedling development., Of 10 candidate genes, two of them (Glyma.04G122400 and Glyma.18G183600) could possibly affect epicotyl length elongation. These results indicate the genetic basis of EL and provides a valuable basis for specific functional studies of epicotyl traits.

GWAS of Reproductive Traits in Large White Pigs on Chip and Imputed Whole-Genome Sequencing Data

Total number born (TNB), number of stillborn (NSB), and gestation length (GL) are economically important traits in pig production, and disentangling the molecular mechanisms associated with traits can provide valuable insights into their genetic structure. Genotype imputation can be used as a practical tool to improve the marker density of single-nucleotide polymorphism (SNP) chips based on sequence data, thereby dramatically improving the power of genome-wide association studies (GWAS). In this study, we applied Beagle software to impute the 50 K chip data to the whole-genome sequencing (WGS) data with average imputation accuracy (R2) of 0.876. The target pigs, 2655 Large White pigs introduced from Canadian and French lines, were genotyped by a GeneSeek Porcine 50K chip. The 30 Large White reference pigs were the key ancestral individuals sequenced by whole-genome resequencing. To avoid population stratification, we identified genetic variants associated with reproductive traits by performing within-population GWAS and cross-population meta-analyses with data before and after imputation. Finally, several genes were detected and regarded as potential candidate genes for each of the traits: for the TNB trait: NOTCH2, KLF3, PLXDC2, NDUFV1, TLR10, CDC14A, EPC2, ORC4, ACVR2A, and GSC; for the NSB trait: NUB1, TGFBR3, ZDHHC14, FGF14, BAIAP2L1, EVI5, TAF1B, and BCAR3; for the GL trait: PPP2R2B, AMBP, MALRD1, HOXA11, and BICC1. In conclusion, expanding the size of the reference population and finding an optimal imputation strategy to ensure that more loci are obtained for GWAS under high imputation accuracy will contribute to the identification of causal mutations in pig breeding.

Genetic Polymorphisms Associated with Prothrombin Time and Activated Partial Thromboplastin Time in Chinese Healthy Population

(1) Background: The purpose of this study was to evaluate the effect of gene polymorphisms on prothrombin time (PT) and activated partial thromboplastin time (APTT) in a healthy Chinese population. (2) Methods: A total of 403 healthy volunteers from a series of novel oral anticoagulants (NOACs) bioequivalence trials in China were included. Coagulation tests for PT and APTT were performed in the central lab at Peking University First Hospital. Whole-exome sequencing (WES) and genome-wide association analysis were performed. (3) Results: In the correlation analysis of PT, 105 SNPs from 84 genes reached the genome-wide significance threshold (p < 1 × 10−5). Zinc Finger Protein 594 (ZNF594) rs184838268 (p = 4.50 × 10−19) was most significantly related to PT, and Actinin Alpha 1 (ACTN1) was found to interact most with other candidate genes. Significant associations with previously reported candidate genes Aurora Kinase B (AURKB), Complement C5(C5), Clock Circadian Regulator (CLOCK), and Histone Deacetylase 9(HDAC9) were detected in our dataset (p < 1 × 10−5). PiggyBac Transposable Element Derived 2(PGBD2) rs75935520 (p = 4.49 × 10−6), Bromodomain Adjacent To Zinc Finger Domain 2A(BAZ2A) rs199970765 (p = 5.69 × 10−6) and Protogenin (PRTG) rs80064850 (p = 8.69 × 10−6) were significantly correlated with APTT (p < 1 × 10−5). The heritability values of PT and APTT were 0.83 and 0.64, respectively; (4) Conclusion: The PT and APTT of healthy populations are affected by genetic polymorphisms. ZNF594 and ACTN1 variants could be novel genetic markers of PT, while PRTG polymorphisms might be associated with APTT levels. The findings could be attributed to ethnic differences, and need further investigation.

The autotetraploid potato genome provides insights into highly heterozygous species

Potato (Solanum tuberosum L.) originated in the Andes and evolved its vegetative propagation strategy through short day-dependent tuber development. Herein, we present a high-quality, chromosome-scale reference genome sequence of a tetraploid potato cultivar. The total length of this genome assembly was 2.67 Gb, with scaffold N50 and contig N50 sizes of 46.24 and 2.19 Mb, respectively. In total, 1.69 Gb repetitive sequences were obtained through de novo annotation, and long terminal repeats were the main transposable elements. A total of 126 070 protein-coding genes were annotated, of which 125 077 (99.21%) were located on chromosomes. The 48 chromosomes were classified into four haplotypes. We annotated 31 506 homologous genes, including 5913 (18.77%) genes with four homologues, 11 103 (35.24%) with three homologues, 12 177 (38.65%) with two homologues and 2313 (7.34%) with one homologue. MLH3, MSH6/7 and RFC3, which are the genes involved in the mismatch repair pathway, were found to be significantly expanded in the tetraploid potato genome relative to the diploid potato genome. Genome-wide association analysis revealed that cytochrome P450, flavonoid synthesis, chalcone enzyme, glycosyl hydrolase and glycosyl transferase genes were significantly correlated with the flesh colours of potato tuber in 150 tetraploid potatoes. This study provides valuable insights into the highly heterozygous autotetraploid potato genome and may facilitate the development of tools for potato cultivar breeding and further studies on autotetraploid crops.

Genome-Wide Identification of Candidate Genes Underlying Soluble Sugar Content in Vegetable Soybean (Glycine max L.) via Association and Expression Analysis

Soluble sugar is a major indicator of the intrinsic quality of vegetable soybean [Glycine max (L.) Merr. ]. The improvement of soluble sugar content in soybean is very important due to its healthcare functions for humans. The genetic mechanism of soluble sugar in soybean is unclear. In this study, 278 diverse soybean accessions were utilized to identify the quantitative trait nucleotides (QTNs) for total soluble sugar content in soybean seeds based on a genome-wide association study (GWAS). A total of 25,921 single-nucleotide polymorphisms (SNPs) with minor allele frequencies (MAFs) ≥ 5% and missing data ≤ 10% were selected for GWAS. Totally, thirteen QTNs associated with total soluble sugar content were identified, which were distributed on ten chromosomes. One hundred and fifteen genes near the 200-kb flanking region of these identified QTNs were considered candidate genes associated with total soluble sugar content in soybean seed. Gene-based association analysis and haplotype analysis were utilized to further identify the effect of candidate genes on total soluble sugar content. Totally, 84 SNPs from seventeen genes across four chromosomes were significantly associated with the total soluble sugar content. Among them, three SNPs from Glyma.02G292900 were identified at two locations, and other eighty-one SNPs from sixteen genes were detected at three locations. Furthermore, expression level analysis of candidate genes revealed that Glyma.02G293200 and Glyma.02G294900 were significantly positively associated with soluble sugar content and Glyma.02G294000 was significantly negatively associated with soluble sugar content. Six genes (i.e., Glyma.02G292600, Glyma.02G292700, Glyma.02G294000, Glyma.02G294300, Glyma.02G294400, and Glyma.15G264200) identified by GWAS were also detected by the analysis of differential expression genes based on soybean germplasms with higher and lower soluble sugar content. Among them, Glyma.02G294000 is the only gene that was identified by gene-based association analysis with total soluble sugar content and was considered an important candidate gene for soluble sugar content. These candidate genes and beneficial alleles would be useful for improving the soluble sugar content of soybean.

Genome-wide association study reveals the genetic architecture for calcium accumulation in grains of hexaploid wheat (Triticum aestivum L.)

BACKGROUND

Hexaploid wheat (Triticum aestivum L.) is a leading cereal crop worldwide. Understanding the mechanism of calcium (Ca) accumulation in wheat is important to reduce the risk of human micronutrient deficiencies. However, the mechanisms of Ca accumulation in wheat grain are only partly understood.

RESULTS

Here, a genome-wide association study (GWAS) was performed to dissect the genetic basis of Ca accumulation in wheat grain using an association population consisting of 207 varieties, with phenotypic data from three locations. In total, 11 non-redundant genetic loci associated with Ca concentration were identified and they explained, on average, 9.61-26.93% of the phenotypic variation. Cultivars containing more superior alleles had increased grain Ca concentrations. Notably, four non-redundant loci were mutually verified by different statistical models in at least two environments, indicating their stability across different environments. Four putative candidate genes linked to Ca accumulation were revealed from the stable genetic loci. Among them, two genes, associated with the stable genetic loci on chromosomes 4A (AX-108912427) and 3B (AX-110922471), encode the subunits of V-type Proton ATPase (TraesCS4A02G428900 and TraesCS3B02G241000), which annotated as the typical generators of a proton gradient that might be involved in Ca homeostasis in wheat grain.

CONCLUSION

To identify genetic loci associated with Ca accumulation, we conducted GWAS on Ca concentrations and detected 11 genetic loci; whereas four genetic loci were stable across different environments. A genetic loci hot spot exists at the end of chromosome 4A and associated with the putative candidate gene TraesCS4A02G428900. The candidate gene TraesCS4A02G428900 encodes V-type proton ATPase subunit e and highly expressed in wheat grains, and it possibly involved in Ca accumulation. This study increases our understanding of the genetic architecture of Ca accumulation in wheat grains, which is potentially helpful for wheat Ca biofortification pyramid breeding.

Genome-wide association study meta-analysis identifies three novel loci for circulating anti-Müllerian hormone levels in women

STUDY QUESTION

Can additional genetic variants for circulating anti-Müllerian hormone (AMH) levels be identified through a genome-wide association study (GWAS) meta-analysis including a large sample of premenopausal women?

SUMMARY ANSWER

We identified four loci associated with AMH levels at P < 5 × 10-8: the previously reported MCM8 locus and three novel signals in or near AMH, TEX41 and CDCA7.

WHAT IS KNOWN ALREADY

AMH is expressed by antral stage ovarian follicles in women, and variation in age-specific circulating AMH levels has been associated with disease outcomes. However, the physiological mechanisms underlying these AMH-disease associations are largely unknown.

STUDY DESIGN, SIZE, DURATION

We performed a GWAS meta-analysis in which we combined summary statistics of a previous AMH GWAS with GWAS data from 3705 additional women from three different cohorts.

PARTICIPANTS/MATERIALS, SETTING, METHODS

In total, we included data from 7049 premenopausal female participants of European ancestry. The median age of study participants ranged from 15.3 to 48 years across cohorts. Circulating AMH levels were measured in either serum or plasma samples using different ELISA assays. Study-specific analyses were adjusted for age at blood collection and population stratification, and summary statistics were meta-analysed using a standard error-weighted approach. Subsequently, we functionally annotated GWAS variants that reached genome-wide significance (P < 5 × 10-8). We also performed a gene-based GWAS, pathway analysis and linkage disequilibrium score regression and Mendelian randomization (MR) analyses.

MAIN RESULTS AND THE ROLE OF CHANCE

We identified four loci associated with AMH levels at P < 5 × 10-8: the previously reported MCM8 locus and three novel signals in or near AMH, TEX41 and CDCA7. The strongest signal was a missense variant in the AMH gene (rs10417628). Most prioritized genes at the other three identified loci were involved in cell cycle regulation. Genetic correlation analyses indicated a strong positive correlation among single nucleotide polymorphisms for AMH levels and for age at menopause (rg = 0.82, FDR = 0.003). Exploratory two-sample MR analyses did not support causal effects of AMH on breast cancer or polycystic ovary syndrome risk, but should be interpreted with caution as they may be underpowered and the validity of genetic instruments could not be extensively explored.

LARGE SCALE DATA

The full AMH GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/).

LIMITATIONS, REASONS FOR CAUTION

Whilst this study doubled the sample size of the most recent GWAS, the statistical power is still relatively low. As a result, we may still lack power to identify more genetic variants for AMH and to determine causal effects of AMH on, for example, breast cancer. Also, follow-up studies are needed to investigate whether the signal for the AMH gene is caused by reduced AMH detection by certain assays instead of actual lower circulating AMH levels.

WIDER IMPLICATIONS OF THE FINDINGS

Genes mapped to the MCM8, TEX41 and CDCA7 loci are involved in the cell cycle and processes such as DNA replication and apoptosis. The mechanism underlying their associations with AMH may affect the size of the ovarian follicle pool. Altogether, our results provide more insight into the biology of AMH and, accordingly, the biological processes involved in ovarian ageing.

STUDY FUNDING/COMPETING INTEREST(S)

Nurses' Health Study and Nurses' Health Study II were supported by research grants from the National Institutes of Health (CA172726, CA186107, CA50385, CA87969, CA49449, CA67262, CA178949). The UK Medical Research Council and Wellcome (217065/Z/19/Z) and the University of Bristol provide core support for ALSPAC. This publication is the work of the listed authors, who will serve as guarantors for the contents of this article. A comprehensive list of grants funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). Funding for the collection of genotype and phenotype data used here was provided by the British Heart Foundation (SP/07/008/24066), Wellcome (WT092830M and WT08806) and UK Medical Research Council (G1001357). M.C.B., A.L.G.S. and D.A.L. work in a unit that is funded by the University of Bristol and UK Medical Research Council (MC_UU_00011/6). M.C.B.'s contribution to this work was funded by a UK Medical Research Council Skills Development Fellowship (MR/P014054/1) and D.A.L. is a National Institute of Health Research Senior Investigator (NF-0616-10102). A.L.G.S. was supported by the study of Dynamic longitudinal exposome trajectories in cardiovascular and metabolic non-communicable diseases (H2020-SC1-2019-Single-Stage-RTD, project ID 874739). The Doetinchem Cohort Study was financially supported by the Ministry of Health, Welfare and Sports of the Netherlands. The funder had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript. Ansh Labs performed the AMH measurements for the Doetinchem Cohort Study free of charge. Ansh Labs was not involved in the data analysis, interpretation or reporting, nor was it financially involved in any aspect of the study. R.M.G.V. was funded by the Honours Track of MSc Epidemiology, University Medical Center Utrecht with a grant from the Netherlands Organization for Scientific Research (NWO) (022.005.021). The Study of Women's Health Across the Nation (SWAN) has grant support from the National Institutes of Health (NIH), DHHS, through the National Institute on Aging (NIA), the National Institute of Nursing Research (NINR) and the NIH Office of Research on Women's Health (ORWH) (U01NR004061; U01AG012505, U01AG012535, U01AG012531, U01AG012539, U01AG012546, U01AG012553, U01AG012554, U01AG012495). The SWAN Genomic Analyses and SWAN Legacy have grant support from the NIA (U01AG017719). The Generations Study was funded by Breast Cancer Now and the Institute of Cancer Research (ICR). The ICR acknowledges NHS funding to the NIHR Biomedical Research Centre. The content of this manuscript is solely the responsibility of the authors and does not necessarily represent official views of the funders. The Sister Study was funded by the Intramural Research Program of the National Institutes of Health (NIH), National Institute of Environmental Health Sciences (Z01-ES044005 to D.P.S.); the AMH assays were supported by the Avon Foundation (02-2012-065 to H.B. Nichols and D.P.S.). The breast cancer genome-wide association analyses were supported by the Government of Canada through Genome Canada and the Canadian Institutes of Health Research, the 'Ministère de l'Économie, de la Science et de l'Innovation du Québec' through Genome Québec and grant PSR-SIIRI-701, The National Institutes of Health (U19 CA148065, X01HG007492), Cancer Research UK (C1287/A10118, C1287/A16563, C1287/A10710) and The European Union (HEALTH-F2-2009-223175 and H2020 633784 and 634935). All studies and funders are listed in Michailidou et al. (Nature, 2017). F.J.M.B. has received fees and grant support from Merck Serono and Ferring BV. D.A.L. has received financial support from several national and international government and charitable funders as well as from Medtronic Ltd and Roche Diagnostics for research that is unrelated to this study. N.S. is scientific consultant for Ansh Laboratories. The other authors declare no competing interests.

Genetic Basis for Resistance Against Viral Nervous Necrosis: GWAS and Potential of Genomic Prediction Explored in Farmed European Sea Bass (Dicentrarchus labrax)

Viral nervous necrosis (VNN) is an infectious disease caused by the red-spotted grouper nervous necrosis virus (RGNNV) in European sea bass and is considered a serious concern for the aquaculture industry with fry and juveniles being highly susceptible. To understand the genetic basis for resistance against VNN, a survival phenotype through the challenge test against the RGNNV was recorded in populations from multiple year classes (YC2016 and YC2017). A total of 4,851 individuals from 181 families were tested, and a subset (n∼1,535) belonging to 122 families was genotyped using a ∼57K Affymetrix Axiom array. The survival against the RGNNV showed low to moderate heritability with observed scale estimates of 0.18 and 0.25 obtained using pedigree vs. genomic information, respectively. The genome-wide association analysis showed a strong signal of quantitative trait loci (QTL) at LG12 which explained ∼33% of the genetic variance. The QTL region contained multiple genes (ITPK1, PLK4, HSPA4L, REEP1, CHMP2, MRPL35, and SCUBE) with HSPA4L and/or REEP1 genes being highly relevant with a likely effect on host response in managing disease-associated symptoms. The results on the accuracy of predicting breeding values presented 20–43% advantage in accuracy using genomic over pedigree-based information which varied across model types and applied validation schemes.

Global Survey, Expressions and Association Analysis of CBLL Genes in Peanut

Cystathionine γ-synthase (CGS), methionine γ-lyase (MGL), cystathionine β-lyase (CBL) and cystathionine γ-lyase (CGL) share the Cys_Met_Meta_PP domain and play important roles in plant stress response and development. In this study, we defined the genes containing the Cys_Met_Meta_PP domain (PF01053.20) as CBL-like genes (CBLL). Twenty-nine CBLL genes were identified in the peanut genome, including 12 from cultivated peanut and 17 from wild species. These genes were distributed unevenly at the ends of different chromosomes. Evolution, gene structure, and motif analysis revealed that CBLL proteins were composed of five different evolutionary branches. Chromosome distribution pattern and synteny analysis strongly indicated that whole-genome duplication (allopolyploidization) contributed to the expansion of CBLL genes. Comparative genomics analysis showed that there were three common collinear CBLL gene pairs among peanut, Arabidopsis, grape, and soybean, but no collinear CBLL gene pairs between peanut and rice. The prediction results of cis-acting elements showed that AhCBLLs, AdCBLLs, and AiCBLLs contained different proportions of plant growth, abiotic stress, plant hormones, and light response elements. Spatial expression profiles revealed that almost all AhCBLLs had significantly higher expression in pods and seeds. All AhCBLLs could respond to heat stress, and some of them could be rapidly induced by cold, salt, submergence, heat and drought stress. Furthermore, one polymorphic site in AiCBLL7 was identified by association analysis which was closely associated with pod length (PL), pod width (PW), hundred pod weight (HPW) and hundred seed weight (HSW). The results of this study provide a foundation for further research on the function of the CBLL gene family in peanut.

Identify known and novel candidate genes associated with backfat thickness in Duroc pigs by large-scale genome-wide association analysis

Backfat thickness (BFT) is complex and economically important traits in the pig industry, since it reflects fat deposition and can be used to measure the carcass lean meat percentage in pigs. In this study, all 6,550 pigs were genotyped using the Geneseek Porcine 50K SNP Chip to identify SNPs related to BFT and to search for candidate genes through genome-wide association analysis in two Duroc populations. In total, 80 SNPs, including 39 significant and 41 suggestive SNPs, and 6 QTLs were identified significantly associated with the BFT. In addition, 9 candidate genes, including a proven major gene MC4R, 3 important candidate genes (RYR1, HMGA1, and NUDT3) which were previously described as related to BFT, and 5 novel candidate genes (SIRT2, NKAIN2, AMH, SORCS1, and SORCS3) were found based on their potential functional roles in BFT. The functions of candidate genes and gene set enrichment analysis indicate that most important pathways are related to energy homeostasis and adipogenesis. Finally, our data suggest that most of the candidate genes can be directly used for genetic improvement through molecular markers, except that the MC4R gene has an antagonistic effect on growth rate and carcass lean meat percentage in breeding. Our results will advance our understanding of the complex genetic architecture of BFT traits and laid the foundation for additional genetic studies to increase carcass lean meat percentage of pig through marker-assisted selection and/or genomic selection.

Heritability and Genomic Architecture of Episodic Exercise-Induced Collapse in Border Collies

An episodic nervous system disorder triggered by strenuous exercise, termed border collie collapse (BCC), exists in border collies and related breeds. The genetic basis of BCC is unknown but is believed to be a complex genetic disorder. Our goal was to estimate the heritability (h²_SNP) of BCC, define its underlying genetic architecture, and identify associated genomic loci using dense whole-genome single-nucleotide polymorphism (SNP) genotyping data. Genotype data were obtained for ~440,000 SNPs from 343 border collies (168 BCC cases and 175 controls). h²_SNP was calculated to be 49–61% depending on the estimated BCC prevalence. A total of 2407 SNPs across the genome accounted for nearly all the h2_SNP of BCC, with an estimated 2003 SNPs of small effect, 349 SNPs of moderate effect, and 56 SNPs of large effect. Genome-wide association analyses identified significantly associated loci on chromosomes 1, 6, 11, 20, and 28, which accounted for ~5% of the total BCC h²_SNP. We conclude that BCC is a moderately- to highly-heritable complex polygenetic disease resulting from contributions from hundreds to thousands of genetic variants with variable effect sizes. Understanding how much the BCC phenotype is determined by genetics and whether major gene mutations are likely to exist inform veterinarians and working/stock dog communities of the true nature of this condition.

Sex-Specific Causal Relations between Steroid Hormones and Obesity-A Mendelian Randomization Study

Steroid hormones act as important regulators of physiological processes including gene expression. They provide possible mechanistic explanations of observed sex-dimorphisms in obesity and coronary artery disease (CAD). Here, we aim to unravel causal relationships between steroid hormones, obesity, and CAD in a sex-specific manner. In genome-wide meta-analyses of four steroid hormone levels and one hormone ratio, we identified 17 genome-wide significant loci of which 11 were novel. Among loci, seven were female-specific, four male-specific, and one was sex-related (stronger effects in females). As one of the loci was the human leukocyte antigen (HLA) region, we analyzed HLA allele counts and found four HLA subtypes linked to 17-OH-progesterone (17-OHP), including HLA-B*14*02. Using Mendelian randomization approaches with four additional hormones as exposure, we detected causal effects of dehydroepiandrosterone sulfate (DHEA-S) and 17-OHP on body mass index (BMI) and waist-to-hip ratio (WHR). The DHEA-S effect was stronger in males. Additionally, we observed the causal effects of testosterone, estradiol, and their ratio on WHR. By mediation analysis, we found a direct sex-unspecific effect of 17-OHP on CAD while the other four hormone effects on CAD were mediated by BMI or WHR. In conclusion, we identified the sex-specific causal networks of steroid hormones, obesity-related traits, and CAD.

Identification of Novel Quantitative Trait Nucleotides and Candidate Genes for Bacterial Wilt Resistance in Tobacco (Nicotiana tabacum L.) Using Genotyping-by-Sequencing and Multi-Locus Genome-Wide Association Studies

Tobacco bacterial wilt (TBW) is a devastating soil-borne disease threatening the yield and quality of tobacco. However, its genetic foundations are not fully understood. In this study, we identified 126,602 high-quality single-nucleotide polymorphisms (SNPs) in 94 tobacco accessions using genotyping-by-sequencing (GBS) and a 94.56 KB linkage disequilibrium (LD) decay rate for candidate gene selection. The population structure analysis revealed two subpopulations with 37 and 57 tobacco accessions. Four multi-locus genome-wide association study (ML-GWAS) approaches identified 142 quantitative trait nucleotides (QTNs) in E1-E4 and the best linear unbiased prediction (BLUP), explaining 0.49-22.52% phenotypic variance. Of these, 38 novel stable QTNs were identified across at least two environments/methods, and their alleles showed significant TBW-DI differences. The number of superior alleles associated with TBW resistance for each accession ranged from 4 to 24; eight accessions had more than 18 superior alleles. Based on TBW-resistant alleles, the five best cross combinations were predicted, including MC133 × Ruyuan No. 1 and CO258 × ROX28. We identified 52 candidate genes around 38 QTNs related to TBW resistance based on homologous functional annotation and KEGG enrichment analysis, e.g., CYCD3;2, BSK1, Nitab4.5_0000641g0050, Nitab4.5_0000929g0030. To the best of our knowledge, this is the first comprehensive study to identify QTNs, superior alleles, and their candidate genes for breeding TBW-resistant tobacco varieties. The results provide further insight into the genetic architecture, marker-assisted selection, and functional genomics of TBW resistance, improving future breeding efforts to increase crop productivity.

Factors relating to bone mineral density in young and middle-aged patients with ankylosing spondylitis

BACKGROUND

Ankylosing spondylitis (AS) is a common chronic progressive rheumatic disease. The aim of this study was to explore factors influencing abnormal bone mineral density (BMD) in young and middle-aged patients with AS.

METHODS

From July 2014 to August 2018, hospitalized patients with AS and health examinees in the health examination center of our clinics, ranging in age from 20 to 50 years, were monitored. The BMD of the lumbar spine and femoral neck of AS patients and those of a healthy control group were measured using dual-energy X-ray absorption. The BMDs of AS patients were compared with respect to age, course of disease, iritis, smoking habits, sex, height, weight, body mass index (BMI), medication use, erythrocyte sedimentation rate (ESR), C-reactive protein (CRP), platelet volume, platelet count, uric acid (UA), alkaline phosphatase (AKP), and calcium ion levels. Single-nucleotide polymorphisms (SNPs) related to BMD were screened using genome-wide association analysis.

RESULTS

There was no statistical difference in the proportion of abnormal bone masses between the different body parts. The BMD of all bones in AS patients was lower than that in healthy controls (P < 0.05). Additionally, BMD was correlated with serum calcium and CRP in AS patients (P < 0.05), but not with age, platelet volume, platelet count, ESR, UA, AKP, height, weight, and BMI. The incidence of abnormal bone mass in AS patients was correlated with sex (P < 0.05), but not with medication use, iritis, or smoking. BMD of the lumbar spine in AS patients did not correlate linearly with the course of the disease, but BMD of the femoral neck correlated linearly with the course of the disease (P < 0.05). BMD was correlated with multiple SNPs in patients with AS. Lumbar BMD was correlated with rs7025373 and rs7848078. Femoral head BMD was correlated with 3:102157365, 3:102157417, rs1252202, rs1681355, rs3891857, rs7842614, and rs9870734, suggesting that genetic factors play a role in BMD in patients with AS.

CONCLUSIONS

The proportion of abnormal bone mass in AS patients was higher than that in healthy individuals of the same age. The factors related to BMD in patients with AS are gender, CRP, and blood calcium. The BMD of the femoral neck of AS patients decreases with the course of the disease, but BMD of the lumbar spine is not related to the course of the disease. BMD in AS patients is associated with multiple SNPs.

Identifying Loci Associated With Bovine Corona Virus Infection and Bovine Respiratory Disease in Dairy and Feedlot Cattle

Bovine coronavirus (BCoV) is associated with respiratory and enteric infections in both dairy and beef cattle worldwide. It is also one of a complex of pathogens associated with bovine respiratory disease (BRD), which affects millions of cattle annually. The objectives of this study were to identify loci and heritability estimates associated with BCoV infection and BRD in dairy calves and feedlot cattle. Dairy calves from California (n = 1,938) and New Mexico (n = 647) and feedlot cattle from Colorado (n = 915) and Washington (n = 934) were tested for the presence of BCoV when classified as BRD cases or controls following the McGuirk scoring system. Two comparisons associated with BCoV were investigated: (1) cattle positive for BCoV (BCoV⁺) were compared to cattle negative for BCoV (BCoV^-) and (2) cattle positive for BCoV and affected with BRD (BCoV⁺BRD⁺) were compared to cattle negative for BCoV and BRD (BCoV^-BRD^-). The Illumina BovineHD BeadChip was used for genotyping, and genome-wide association analyses (GWAA) were performed using EMMAX (efficient mixed-model association eXpedited). The GWAA for BCoV⁺ identified 51 loci (p < 1 × 10^-5; 24 feedlot, 16 dairy, 11 combined) associated with infection with BCoV. Three loci were associated with BCoV⁺ across populations. Heritability estimates for BCoV⁺ were 0.01 for dairy, 0.11 for feedlot cattle, and 0.03 for the combined population. For BCoV⁺BRD⁺, 80 loci (p < 1 × 10^-5; 26 feedlot, 25 dairy, 29 combined) were associated including 14 loci across populations. Heritability estimates for BCoV⁺BRD⁺ were 0.003 for dairy, 0.44 for feedlot cattle, and 0.07 for the combined population. Several positional candidate genes associated with BCoV and BRD in this study have been associated with other coronaviruses and respiratory infections in humans and mice. These results suggest that selection may reduce susceptibility to BCoV infection and BRD in cattle.

Genome-wide hierarchical mixed model association analysis

In genome-wide mixed model association analysis, we stratified the genomic mixed model into two hierarchies to estimate genomic breeding values (GBVs) using the genomic best linear unbiased prediction and statistically infer the association of GBVs with each SNP using the generalized least square. The hierarchical mixed model (Hi-LMM) can correct confounders effectively with polygenic effects as residuals for association tests, preventing potential false-negative errors produced with genome-wide rapid association using mixed model and regression or an efficient mixed-model association expedited (EMMAX). Meanwhile, the Hi-LMM performs the same statistical power as the exact mixed model association and the same computing efficiency as EMMAX. When the GBVs have been estimated precisely, the Hi-LMM can detect more quantitative trait nucleotides (QTNs) than existing methods. Especially under the Hi-LMM framework, joint association analysis can be made straightforward to improve the statistical power of detecting QTNs.

Functional genetics for studying the human immune system

Although small numbers of immune-mediated diseases are inherited due to rare genetic mutations, most are multifactorial diseases caused by multiple elements including genetic and environmental factors. In the case of autoimmune diseases, many disease-susceptibility genes, including several in the major histocompatibility gene complex, have been reported, and over the past 10 years, genome-wide association studies (GWAS) have been used to analyze disease-susceptibility loci in representative diseases. Furthermore, many disease susceptibility variants have been found to be related to gene expression levels. The expression of genes involved in disease pathogenesis is often cell type-specific, and this is closely related to epigenome alterations. Genomic information is present even before the onset of a disease and has a clear causal relationship to the disease (i.e., the outcome). Therefore, it is important to establish functional genetics in human immunology to understand the pathogenesis of diseases using these pieces of information. We can then apply these results to drug discovery. Here, we will review these issues, especially focusing on autoimmune diseases, and discuss current and future directions of human immune system research.

Genomic Characterization of the Fruity Aroma Gene, FaFAD1, Reveals a Gene Dosage Effect on γ-Decalactone Production in Strawberry (Fragaria × ananassa)

Strawberries produce numerous volatile compounds that contribute to the unique flavors of fruits. Among the many volatiles, γ-decalactone (γ-D) has the greatest contribution to the characteristic fruity aroma in strawberry fruit. The presence or absence of γ-D is controlled by a single locus, FaFAD1. However, this locus has not yet been systematically characterized in the octoploid strawberry genome. It has also been reported that the volatile content greatly varies among the strawberry varieties possessing FaFAD1, suggesting that another genetic factor could be responsible for the different levels of γ-D in fruit. In this study, we explored the genomic structure of FaFAD1 and determined the allele dosage of FaFAD1 that regulates variations of γ-D production in cultivated octoploid strawberry. The genome-wide association studies confirmed the major locus FaFAD1 that regulates the γ-D production in cultivated strawberry. With the hybrid capture-based next-generation sequencing analysis, a major presence-absence variation of FaFAD1 was discovered among γ-D producers and non-producers. To explore the genomic structure of FaFAD1 in the octoploid strawberry, three bacterial artificial chromosome (BAC) libraries were developed. A deletion of 8,262 bp was consistently found in the FaFAD1 region of γ-D non-producing varieties. With the newly developed InDel-based codominant marker genotyping, along with γ-D metabolite profiling data, we revealed the impact of gene dosage effect for the production of γ-D in the octoploid strawberry varieties. Altogether, this study provides systematic information of the prominent role of FaFAD1 presence and absence polymorphism in producing γ-D and proposes that both alleles of FaFAD1 are required to produce the highest content of fruity aroma in strawberry fruit.

Hierarchical mixed-model expedites genome-wide longitudinal association analysis

A hierarchical random regression model (Hi-RRM) was extended into a genome-wide association analysis for longitudinal data, which significantly reduced the dimensionality of repeated measurements. The Hi-RRM first modeled the phenotypic trajectory of each individual using a RRM and then associated phenotypic regressions with genetic markers using a multivariate mixed model (mvLMM). By spectral decomposition of genomic relationship and regression covariance matrices, the mvLMM was transformed into a multiple linear regression, which improved computing efficiency while implementing mvLMM associations in efficient mixed-model association expedited (EMMAX). Compared with the existing RRM-based association analyses, the statistical utility of Hi-RRM was demonstrated by simulation experiments. The method proposed here was also applied to find the quantitative trait nucleotides controlling the growth pattern of egg weights in poultry data.

Structural equation modeling for unraveling the multivariate genomic architecture of milk proteins in dairy cattle

The aims of this study were to investigate potential functional relationships among milk protein fractions in dairy cattle and to carry out a structural equation model (SEM) GWAS to provide a decomposition of total SNP effects into direct effects and effects mediated by traits that are upstream in a phenotypic network. To achieve these aims, we first fitted a mixed Bayesian multitrait genomic model to infer the genomic correlations among 6 milk nitrogen fractions [4 caseins (CN), namely κ-, β-, α_S1-, and α_S2-CN, and 2 whey proteins, namely β-lactoglobulin (β-LG) and α-lactalbumin (α-LA)], in a population of 989 Italian Brown Swiss cows. Animals were genotyped with the Illumina BovineSNP50 Bead Chip v.2 (Illumina Inc.). A Bayesian network approach using the max-min hill-climbing (MMHC) algorithm was implemented to model the dependencies or independence among traits. Strong and negative genomic correlations were found between β-CN and α_S1-CN (-0.706) and between β-CN and κ-CN (-0.735). The application of the MMHC algorithm revealed that κ-CN and β-CN seemed to directly or indirectly influence all other milk protein fractions. By integrating multitrait model GWAS and SEM-GWAS, we identified a total of 127 significant SNP for κ-CN, 89 SNP for β-CN, 30 SNP for α_S1-CN, and 14 SNP for α_S2-CN (mostly shared among CN and located on Bos taurus autosome 6) and 15 SNP for β-LG (mostly located on Bos taurus autosome 11), whereas no SNP passed the significance threshold for α-LA. For the significant SNP, we assessed and quantified the contribution of direct and indirect paths to total marker effect. Pathway analyses confirmed that common regulatory mechanisms (e.g., energy metabolism and hormonal and neural signals) are involved in the control of milk protein synthesis and metabolism. The information acquired might be leveraged for setting up optimal management and selection strategies aimed at improving milk quality and technological characteristics in dairy cattle.

Genetic Variation, GWAS and Accuracy of Prediction for Host Resistance to Sparicotyle chrysophrii in Farmed Gilthead Sea Bream (Sparus aurata)

Gilthead sea bream (Sparus aurata) belongs to a group of teleost which has high importance in Mediterranean aquaculture industry. However, industrial production is increasingly compromised by an elevated outbreak of diseases in sea cages, especially a disease caused by monogeneans parasite Sparicotyle chrysophrii. This parasite mainly colonizes gill tissues of host and causes considerable economical losses with mortality and reduction in growth. The aim of current study was to explore the genetics of host resistance against S. chrysophrii and investigate the potential for genomic selection to possibly accelerate genetic progress. To achieve the desired goals, a test population derived from the breeding nucleus of Andromeda Group was produced. This experimental population was established by crossing of parents mated in partial factorial crosses of ∼8 × 8 using 58 sires and 62 dams. The progeny obtained from this mating design was challenged with S. chrysophrii using a controllable cohabitation infection model. At the end of the challenge, fish were recorded for parasite count, and all the recorded fish were tissue sampled for genotyping by sequencing using 2b-RAD methodology. The initial (before challenge test) and the final body weight (after challenge test) of the fish were also recorded. The results obtained through the analysis of phenotypic records (n = 615) and the genotypic data (n = 841, 724 offspring and 117 parents) revealed that the resistance against this parasite is lowly heritable (h² = 0.147 with pedigree and 0.137 with genomic information). We observed moderately favorable genetic correlation (R_g = −0.549 to −0.807) between production traits (i.e., body weight and specific growth rate) and parasite count, which signals a possibility of indirect selection. A locus at linkage group 17 was identified that surpassed chromosome-wide Bonferroni threshold which explained 22.68% of the total genetic variance, and might be playing role in producing genetic variation. The accuracy of prediction was improved by 8% with genomic information compared to pedigree.

Genomic analyses reveal the genetic basis of early maturity and identification of loci and candidate genes in upland cotton (Gossypium hirsutum L.)

Although upland cotton (Gossypium hirsutism L.) originated in the tropics, this early maturity cotton can be planted as far north as 46°N in China due to the accumulation of numerous phenotypic and physiological adaptations during domestication. However, how the genome of early maturity cotton has been altered by strong human selection remains largely unknown. Herein, we report a cotton genome variation map generated by the resequencing of 436 cotton accessions. Whole-genome scans for sweep regions identified 357 putative selection sweeps covering 4.94% (112 Mb) of the upland cotton genome, including 5184 genes. These genes were functionally related to flowering time control, hormone catabolism, ageing and defence response adaptations to environmental changes. A genome-wide association study (GWAS) for seven early maturity traits identified 307 significant loci, 22.48% (69) of which overlapped with putative selection sweeps that occurred during the artificial selection of early maturity cotton. Several previously undescribed candidate genes associated with early maturity were identified by GWAS. This study provides insights into the genetic basis of early maturity in upland cotton as well as breeding resources for cotton improvement.

Allelic variation in rice Fertilization Independent Endosperm 1 contributes to grain width under high night temperature stress

A higher minimum (night-time) temperature is considered a greater limiting factor for reduced rice yield than a similar increase in maximum (daytime) temperature. While the physiological impact of high night temperature (HNT) has been studied, the genetic and molecular basis of HNT stress response remains unexplored. We examined the phenotypic variation for mature grain size (length and width) in a diverse set of rice accessions under HNT stress. Genome-wide association analysis identified several HNT-specific loci regulating grain size as well as loci that are common for optimal and HNT stress conditions. A novel locus contributing to grain width under HNT conditions colocalized with Fie1, a component of the FIS-PRC2 complex. Our results suggest that the allelic difference controlling grain width under HNT is a result of differential transcript-level response of Fie1 in grains developing under HNT stress. We present evidence to support the role of Fie1 in grain size regulation by testing overexpression (OE) and knockout mutants under heat stress. The OE mutants were either unaltered or had a positive impact on mature grain size under HNT, while the knockouts exhibited significant grain size reduction under these conditions.

Genome-Wide Association Mapping Identifies Novel Loci for Quantitative Resistance to Blackleg Disease in Canola

Blackleg disease, caused by the fungal pathogen Leptosphaeria maculans, continues to be a major concern for sustainable production of canola (Brassica napus L.) in many parts of the world. The deployment of effective quantitative resistance (QR) is recognized as a durable strategy in providing natural defense to pathogens. Herein, we uncover loci for resistance to blackleg in a genetically diverse panel of canola accessions by exploiting historic recombination events which occurred during domestication and selective breeding by genome-wide association analysis (GWAS). We found extensive variation in resistance to blackleg at the adult plant stage, including for upper canopy infection. Using the linkage disequilibrium and genetic relationship estimates from 12,414 high quality SNPs, GWAS identified 59 statistically significant and "suggestive" SNPs on 17 chromosomes of B. napus genome that underlie variation in resistance to blackleg, evaluated under field and shade-house conditions. Each of the SNP association accounted for up to 25.1% of additive genetic variance in resistance among diverse panel of accessions. To understand the homology of QR genomic regions with Arabidopsis thaliana genome, we searched the synteny between QR regions with 22 ancestral blocks of Brassicaceae. Comparative analyses revealed that 25 SNP associations for QR were localized in nine ancestral blocks, as a result of genomic rearrangements. We further showed that phenological traits such as flowering time, plant height, and maturity confound the genetic variation in resistance. Altogether, these findings provided new insights on the complex genetic control of the blackleg resistance and further expanded our understanding of its genetic architecture.

Dentate gyrus volume deficit in schizophrenia

BACKGROUND

Schizophrenia is associated with robust hippocampal volume deficits but subregion volume deficits, their associations with cognition, and contributing genes remain to be determined.

METHODS

Hippocampal formation (HF) subregion volumes were obtained using FreeSurfer 6.0 from individuals with schizophrenia (n = 176, mean age ± s.d. = 39.0 ± 11.5, 132 males) and healthy volunteers (n = 173, mean age ± s.d. = 37.6 ± 11.3, 123 males) with similar mean age, gender, handedness, and race distributions. Relationships between the HF subregion volume with the largest between group difference, neuropsychological performance, and single-nucleotide polymorphisms were assessed.

RESULTS

This study found a significant group by region interaction on hippocampal subregion volumes. Compared to healthy volunteers, individuals with schizophrenia had significantly smaller dentate gyrus (DG) (Cohen's d = -0.57), Cornu Ammonis (CA) 4, molecular layer of the hippocampus, hippocampal tail, and CA 1 volumes, when statistically controlling for intracranial volume; DG (d = -0.43) and CA 4 volumes remained significantly smaller when statistically controlling for mean hippocampal volume. DG volume showed the largest between group difference and significant positive associations with visual memory and speed of processing in the overall sample. Genome-wide association analysis with DG volume as the quantitative phenotype identified rs56055643 (β = 10.8, p < 5 × 10-8, 95% CI 7.0-14.5) on chromosome 3 in high linkage disequilibrium with MOBP. Gene-based analyses identified associations between SLC25A38 and RPSA and DG volume.

CONCLUSIONS

This study suggests that DG dysfunction is fundamentally involved in schizophrenia pathophysiology, that it may contribute to cognitive abnormalities in schizophrenia, and that underlying biological mechanisms may involve contributions from MOBP, SLC25A38, and RPSA.

Identification of Loci and Candidate Genes Analyses for Tocopherol Concentration of Soybean Seed

Tocopherol (Toc) occurs in soybean seeds and is extracted together with the soybean oil. Toc is utilized as an antioxidant in food and an additive in animal feed. A total of 180 representative accessions and 144 recombinant inbred lines (RILs) from the cross of 'Hefeng 25' and 'OAC Bayfield' were selected to evaluate individuals and total Toc concentrations in soybean seeds. The 180 soybean samples were sequenced by the approach of Specific Locus Amplified Fragment Sequencing (SLAF-seq). A total of 22,611 single nucleotide polymorphisms (SNPs) were developed. Nineteen quantitative trait nucleotides (QTNs) were identified associated with individual or total-Toc based on genome-wide association analysis (GWAS). Among them, three QTNs located near known QTLs, and 16 were novel. Eighteen QTLs and nine eQTLs were also detected by linkage mapping. The QTN rs9337368 on Chr.02 was colocalized according to the linkage mapping of the RILs and genome-wide association analysis and regarded as a stable locus for mining the candidate genes in association with Toc. A total of 42 candidate genes near the 200 kbp flanking region of this identified locus were found. Upon a gene-based association, 11 SNPs from five genes out of the 42 candidates were detected. Expression level analysis of five candidate genes revealed that two genes were significantly related to Toc content. The identified loci, along with the candidate genes, might be valuable for increasing the Toc concentration in soybean seeds and improving the nutritional value of soybean oil.

Mapping a Partial Andromonoecy Locus in Citrullus lanatus Using BSA-Seq and GWAS Approaches

The sexual expression of watermelon plants is the result of the distribution and occurrence of male, female, bisexual and hermaphrodite flowers on the main and secondary stems. Plants can be monoecious (producing male and female flowers), andromonoecious (producing male and hermaphrodite flowers), or partially andromonoecious (producing male, female, bisexual, and hermaphrodite flowers) within the same plant. Sex determination of individual floral buds and the distribution of the different flower types on the plant, are both controlled by ethylene. A single missense mutation in the ethylene biosynthesis gene CitACS4, is able to promote the conversion of female into hermaphrodite flowers, and therefore of monoecy (genotype MM) into partial andromonoecy (genotype Mm) or andromonoecy (genotype mm). We phenotyped and genotyped, for the M/m locus, a panel of 207 C. lanatus accessions, including five inbreds and hybrids, and found several accessions that were repeatedly phenotyped as PA (partially andromonoecious) in several locations and different years, despite being MM. A cosegregation analysis between a SNV in CitACS4 and the PA phenotype, demonstrated that the occurrence of bisexual and hermaphrodite flowers in a PA line is not dependent on CitACS4, but conferred by an unlinked recessive gene which we called pa. Two different approaches were performed to map the pa gene in the genome of C. lanatus: bulk segregant analysis sequencing (BSA-seq) and genome wide association analysis studies (GWAS). The BSA-seq study was performed using two contrasting bulks, the monoecious M-bulk and the partially andromonoecious PA-bulk, each one generated by pooling DNA from 20 F2 plants. For GWAS, 122 accessions from USDA gene bank, already re-sequenced by genotyping by sequencing (GBS), were used. The combination of the two approaches indicates that pa maps onto a genomic region expanding across 32.24-36.44 Mb in chromosome 1 of watermelon. Fine mapping narrowed down the pa locus to a 867 Kb genomic region containing 101 genes. A number of candidate genes were selected, not only for their function in ethylene biosynthesis and signalling as well as their role in flower development and sex determination, but also by the impact of the SNPs and indels differentially detected in the two sequenced bulks.

Genome-Wide Association Analyses of Equine Metabolic Syndrome Phenotypes in Welsh Ponies and Morgan Horses

Equine metabolic syndrome (EMS) is a complex trait for which few genetic studies have been published. Our study objectives were to perform within breed genome-wide association analyses (GWA) to identify associated loci in two high-risk breeds, coupled with meta-analysis to identify shared and unique loci between breeds. GWA for 12 EMS traits identified 303 and 142 associated genomic regions in 264 Welsh ponies and 286 Morgan horses, respectively. Meta-analysis demonstrated that 65 GWA regions were shared across breeds. Region boundaries were defined based on a fixed-size or the breakdown of linkage disequilibrium, and prioritized if they were: shared between breeds or across traits (high priority), identified in a single GWA cohort (medium priority), or shared across traits with no SNPs reaching genome-wide significance (low priority), resulting in 56 high, 26 medium, and seven low priority regions including 1853 candidate genes in the Welsh ponies; and 39 high, eight medium, and nine low priority regions including 1167 candidate genes in the Morgans. The prioritized regions contained protein-coding genes which were functionally enriched for pathways associated with inflammation, glucose metabolism, or lipid metabolism. These data demonstrate that EMS is a polygenic trait with breed-specific risk alleles as well as those shared across breeds.

Genome-wide association study of circulating folate one-carbon metabolites

Experimental, observational, and clinical trials support a critical role of folate one-carbon metabolism (FOCM) in colorectal cancer (CRC) development. In this report, we focus on understanding the relationship between common genetic variants and metabolites of FOCM. We conducted a genome-wide association study of FOCM biomarkers among 1,788 unaffected (without CRC) individuals of European ancestry from the Colon Cancer Family Registry. Twelve metabolites, including 5-methyltetrahydrofolate, vitamin B₂ (flavin mononucleotide and riboflavin), vitamin B₆ (4-pyridoxic acid, pyridoxal, and pyridoxamine), total homocysteine, methionine, S-adenosylmethionine, S-adenosylhomocysteine, cystathionine, and creatinine were measured from plasma using liquid chromatography-mass spectrometry (LC-MS) or LC-MS/MS. For each individual biomarker, we estimated genotype array-specific associations followed by a fixed-effect meta-analysis. We identified the variant rs35976024 (at 2p11.2 and intronic of ATOH8) associated with total homocysteine (p = 4.9 × 10^-8 ). We found a group of six highly correlated variants on chromosome 15q14 associated with cystathionine (all p < 5 × 10^-8 ), with the most significant variant rs28391580 (p = 2.8 × 10^-8 ). Two variants (rs139435405 and rs149119426) on chromosome 14q13 showed significant (p < 5 × 10^-8 ) associations with S-adenosylhomocysteine. These three biomarkers with significant associations are closely involved in homocysteine metabolism. Furthermore, when assessing the principal components (PCs) derived from seven individual biomarkers, we identified the variant rs12665366 (at 6p25.3 and intronic of EXOC2) associated with the first PC (p = 2.3 × 10^-8 ). Our data suggest that common genetic variants may play an important role in FOCM, particularly in homocysteine metabolism.