Genome-wide Association Analysis Research Articles

Genome-wide association analysis of cystatin c and creatinine kidney function in Chinese women

Genome-wide association analysis of cystatin c and creatinine kidney function in Chinese women

Genetic improvement of eating and cooking quality of rice cultivars in southern China.

The genetic improvement of rice eating and cooking quality (ECQ) is an important goal in rice breeding. It is important to understand the genetic regulation of ECQ at the genomic level for effective breeding to improve ECQ. However, the mechanisms underlying the improvement of ECQ of indica and japonica cultivars in southern China remain unclear. In this study, 290 rice cultivars (155 indica and 135 japonica cultivars) bred in southern China in the past 30 years were collected. Physicochemical indicators, namely, apparent amylose content (AAC), protein content (PC), lipid content and taste value, were measured and correlation analysis was performed. Adecrease in AAC and PC was a crucial factor for the ECQ improvement of the rice cultivars in southern China. Genome-wide association analysis and selective domestication analysis preliminarily clarified that the comprehensive utilization of major and minor genes was an important genetic basis for improvement of ECQ. An elite allele, RAmy1AA, with potential application in breeding to improve starch viscosity characteristics and ECQ, was mined. The Wxb/OsmtSSB1LT/OsDML4G/RPBFT/Du3T and Wxb/OsEro1T/Glup3G/OsNAC25G/OsBEIIbC/RAmy1AA/FLO12A gene modules, neither of which have been widely used, are proposed as the optimal allele combinations for ECQ improvement of indica and japonica cultivars in southern China. The results clarify the genetic regulation of rice ECQ improvement in southern China and provide novel genetic resources and breeding strategies for ECQ improvement in rice.

Genome-wide association of an organic naked barley diversity panel identified quantitative trait loci for disease resistance.

Foliar fungal diseases are a major limitation in organic naked barley (Hordeum vulgare L.) production. The lack of conventional fungicides in organic systems increases reliance on genetic resistance. We evaluated the severity of barley stripe rust (Puccinia striiformis f.sp. hordei Westend), leaf rust (Puccina hordei sp. hordei), spot blotch (Cochliobolus sativus, anamorph Bipolaris sorokiniana (S. Ito & Kurib.) Drechsler ex Dastur), and scald (Rhynchosporium commune Zaffarano, McDonald and Linde sp. November) on a naked barley diversity panel of 350 genotypes grown in 13 environments to identify quantitative trait loci associated with disease resistance. Genome-wide association analyses across and within environments found 10 marker trait associations for barley stripe rust, four marker trait associations for leaf rust, one marker trait association for scald, and five marker trait associations for spot blotch. Structure analysis identified six Ward groups based on genotypic diversity. Resistance to susceptible allele ratios were high for stripe rust and spot blotch, moderate for leaf rust, and low for scald. Combined phenotypic analysis values for each disease overlayed by a principal component analysis found distinct resistance and susceptibility patterns for barley stripe rust and scald. Most significant marker trait associations were previously identified in the literature, providing confirmation and potential new sources of disease resistance for genetic improvement of naked barley germplasm.

Candidate genetic loci modifying the colorectal cancer risk caused by lifestyle risk factors

Objective: 65-70% of colorectal cancer (CRC) cases are considered sporadic; they arise under the influence of environmental factors in individuals lacking a family history of CRC. Low-risk genetic variants are believed to contribute to CRC risk, in tandem with lifestyle factors. Design: 616 non-familial Swedish CRC cases with at least 1 of the following five risk factors: smoking, excessive alcohol consumption, physical inactivity, adherence to an unhealthy diet, and excess body weight were included in this study. A control group consisting of 1642 healthy individuals was used. Cases and controls were genotyped from blood samples at the Centre for Inherited Disease Research at Johns Hopkins University within the Colorectal Transdisciplinary Study (CORECT) research collaboration, using the Illumina Infinium® OncoArray-500 K BeadChip. Five separate genome-wide haplotype association analyses were performed, one for each risk factor. Logistic regression models were employed to estimate associations between haplotypes (exposure) and CRC (outcome) in cases with lifestyle risk factors versus controls. Haplotypes with an odds ratio (OR) >1 were considered candidate risk markers, denoting an area of interest in the genome. A significance threshold of p < 5×10-8 was used. Results: We found 17 haplotype regions significantly associated with CRC in cases vs controls. Several regions included genes linked to inflammation and tumor promotion. Conclusion: We concluded that having certain genetic variants was associated with an increased risk of CRC compared to healthy controls among cases with known lifestyle risk factors. The interplay of lifestyle and genetic risk factors calls for further elucidation.

QTL mapping and genome-wide association analysis reveal genetic loci and candidate gene for resistance to gray leaf spot in tropical and subtropical maize germplasm.

Using QTL mapping and GWAS, two candidate genes (Zm00001d051039 and Zm00001d051147) were consistently identified across the three different environments and BLUP values. GWAS analysis identified the candidate gene, Zm00001d044845. These genes were subsequently validated to exhibit a significant association with maize gray leaf spot (GLS) resistance. Gray leaf spot (GLS) is a major foliar disease of maize (Zea mays L.) that causes significant yield losses worldwide. Understanding the genetic mechanisms underlying gray leaf spot resistance is crucial for breeding high-yielding and disease-resistant varieties. In this study, eight tropical and subtropical germplasms were crossed with the temperate germplasm Ye107 to develop a nested association mapping (NAM) population comprising 1,653 F2:8 RILs, consisting of eight recombinant inbred line (RIL) subpopulations, using the single-seed descent method. The NAM population was evaluated for GLS resistance in three different environments, and genotyping by sequencing of the NAM population generated 593,719 high-quality single-nucleotide polymorphisms (SNPs). Linkage analysis and genome-wide association studies (GWASs) were conducted to identify candidate genes regulating GLS resistance in maize. Both analyses identified 25 QTLs and 149 SNPs that were significantly associated with GLS resistance. Candidate genes were screened 20Kb upstream and downstream of the significant SNPs, and three novel candidate genes (Zm00001d051039, Zm00001d051147, and Zm00001d044845) were identified. Zm00001d051039 and Zm00001d051147 were located on chromosome 4 and co-localized in both linkage (qGLS4-1 and qGLS4-2) and GWAS analyses. SNP-138,153,206 was located 0.499kb downstream of the candidate gene Zm00001d051039, which encodes the protein IN2-1 homolog B, a homolog of glutathione S-transferase (GST). GSTs and protein IN2-1 homolog B scavenge reactive oxygen species under various stress conditions, and GSTs are believed to protect plants from a wide range of biotic and abiotic stresses by detoxifying reactive electrophilic compounds. Zm00001d051147 encodes a probable beta-1,4-xylosyltransferase involved in the biosynthesis of xylan in the cell wall, enhancing resistance. SNP-145,813,215 was located 2.69kb downstream of the candidate gene. SNP-5,043,412 was consistently identified in three different environments and BLUP values and was located 8.788kb downstream of the candidate gene Zm00001d044845 on chromosome 9. Zm00001d044845 encodes the U-box domain-containing protein 4 (PUB4), which is involved in regulating plant immunity. qRT-PCR analysis showed that the relative expression levels of the three candidate genes were significantly upregulated in the leaves of the TML139 (resistant) parent, indicating that these three candidate genes could be associated with resistance to GLS. The findings of this study are significant for marker-assisted breeding aimed at enhancing resistance to GLS in maize and lay the foundation for further elucidation of the genetic mechanisms underlying resistance to gray leaf spot in maize and breeding of new disease-resistant varieties.

Causal links of human serum metabolites on the risk of prostate cancer: insights from genome-wide Mendelian randomization, single-cell RNA sequencing, and metabolic pathway analysis

BackgroundRecently, serum metabolites have shown potential in predicting survival outcomes and may be related to the pathogenesis of prostate cancer. Nevertheless, the precise impact concerning the genetic effect of metabolites on prostate cancer risk remains obscure. In this context, we conducted a Mendelian randomization (MR) study aiming to explore the causality between genetically determined metabolites and the risk of prostate cancer.MethodsWe conducted a two-sample MR analysis aiming to identify the underlying metabolites associated with prostate cancer. Exposure information was obtained from the largest metabolome-based genome-wide association (GWAS) data containing 7,824 Europeans. Genome-wide association analysis was utilized to detect instrumental variables (IVs) for metabolites. We applied the inverse-variance weighted (IVW) approach as the primary method, and to augment the reliability and robustness of our findings, additional analysis methods encompassing weighted median, MR-Egger, and leave-one-out analysis were utilized. MR-Egger intercept test was implemented to explore the pleiotropy. Cochran’s Q test was utilized to quantify the degree of heterogeneity. Additionally, we performed metabolic pathway analysis and single-cell RNA sequencing analysis.ResultsWe found that three serum metabolites were causally associated with prostate cancer after utilizing rigorous screening standards. Utilizing single nucleotide polymorphisms as IVs, a 1-SD increase in fructose was associated with 77% higher risk of prostate cancer (OR:1.77, 95%CI: 1.05-2.97, PIVW=0.031), a 1-SD increase in N1-methyl-3-pyridone-4-carboxamide was associated with 29% higher risk of prostate cancer (OR:1.29, 95%CI: 1.05-1.58, PIVW=0.017), and a 1-SD increase in 12-hydroxyeicosatetraenoate (12-HETE) was associated with 18% higher risk of prostate cancer (OR:1.18, 95%CI: 1.07-1.31, PIVW=0.0008). Metabolites that were causally linked to the risk of prostate cancer were mainly enriched in the valine, leucine and isoleucine biosynthesis pathway (P=0.026) and the nicotinate and nicotinamide metabolism pathway (P=0.048).ConclusionsOur MR analysis provided suggestive evidence supporting the causal relationships between three identified serum metabolites and prostate cancer, necessitating further investigation to elucidate the underlying mechanisms through which these blood metabolites and metabolic pathways may impact the initiation and progression of prostate cancer.

Genetic legacy of ancient hunter-gatherer Jomon in Japanese populations.

The tripartite ancestral structure is a recently proposed model for the genetic origin of modern Japanese, comprising indigenous Jomon hunter-gatherers and two additional continental ancestors from Northeast Asia and East Asia. To investigate the impact of the tripartite structure on genetic and phenotypic variation today, we conducted biobank-scale analyses by merging Biobank Japan (BBJ; n = 171,287) with ancient Japanese and Eurasian genomes (n = 22). We demonstrate the applicability of the tripartite model to Japanese populations throughout the archipelago, with an extremely strong correlation between Jomon ancestry and genomic variation among individuals. We also find that the genetic legacy of Jomon ancestry underlies an elevated body mass index (BMI). Genome-wide association analysis with rigorous adjustments for geographical and ancestral substructures identifies 132 variants that are informative for predicting individual Jomon ancestry. This prediction model is validated using independent Japanese cohorts (Nagahama cohort, n = 2993; the second cohort of BBJ, n = 72,695). We further confirm the phenotypic association between Jomon ancestry and BMI using East Asian individuals from UK Biobank (n = 2286). Our extensive analysis of ancient and modern genomes, involving over 250,000 participants, provides valuable insights into the genetic legacy of ancient hunter-gatherers in contemporary populations.

A genome-wide association study of the racing performance traits in Yili horses based on Blink and FarmCPU models.

Racing performance traits are the main indicators for evaluating the performance and value of sport horses. The aim of this study was to identify the key genes for racing performance traits in Yili horses by performing a genome-wide association study (GWAS). Breeding values for racing performance traits were calculated for Yili horses (n = 827) using an animal model. Genome-wide association analysis of racing performance traits in horses (n = 236) was carried out using the Blink, and FarmCPU models in GAPIT software, and genes within the significant regions were functionally annotated. The results of GWAS showed that a total of 24 significant SNP markers (P < 6.05 × 10- 9) and 22 suggestive SNP markers (P < 1.21 × 10- 7) were identified. Among them, the Blink associated 16 significant SNP loci and FarmCPU associated 12 significant SNP loci. A total of 127 candidate genes (50 significant) were annotated. Among these, CNTN6 (motor coordination), NIPA1 (neuronal development), and DCC (dopamine pathway maturation) may be the main candidate genes affecting speed traits. SHANK2 (neuronal synaptic regulation), ISCA1 (mitochondrial protein assembly), and KCNIP4 (neuronal excitability) may be the main candidate genes affecting ranking score traits. A common locus (ECA1: 22698579) was significantly associated with racing performance traits, and the function of the genes at this locus needs to be studied in depth. These findings will provide new insights into the detection and selection of genetic variants for racing performance and will help to accelerate the genetic improvement of Yili horses.

GWAS-Based Prediction of Genes Regulating the Weight of Mobilized Reserved Seeds in Sweet Corn

Seed reserve mobilization is a crucial physiological process during seed germination. Enhancing the reserve mobilization in sweet corn is vital for seed germination and seedling growth. In this study, a genome-wide association analysis (GWAS) was conducted to predict candidate genes for regulating the weight of mobilized reserved seeds (WMRS) and kernel weight (KW) in sweet corn. GWAS was performed using the BLINK model with the maize 56K SNP chip. The results indicated that there was a positive correlation between the WMRS and KW, with coefficients of variation of 68.18% and 44.63%. Association analysis identified thirteen SNPs associated with two traits, and linkage disequilibrium analysis revealed that eight of these SNPs were in strong linkage. A total of 298 candidate genes were identified within the confidence interval, of which 79 were annotated. About 20 candidate genes were identified through the comparison of homologous genes in Arabidopsis. These genes were enriched in regulating ribosome biogenesis, signal transduction, hormone synthesis, and RNA degradation processes. This study provides important insights into the genetic mechanisms governing germination traits in sweet corn, aiding further research into the localization and cloning of genes involved in the mobilization of reserve materials.

Genetic dissection of resistance to Phytophthora sojae using genome-wide association and linkage analysis in soybean [Glycine max (L.) Merr.

Two novel and one known genomic regions associated with R-gene resistance to Phytophthora sojae were identified by genome-wide association analysis and linkage analysis in soybean. Phytophthora root and stem rot (PRR) caused by Phytophthora sojae is a severe disease that causes substantial economic losses in soybean [Glycine max (L.) Merr.]. The primary approach for successful disease management of PRR is using R-gene-mediated resistance. Based on the phenotypic evaluation of 376 cultivated soybean accessions for the R-gene type resistance to P. sojae (isolate 2457), a genome-wide association analysis identified two regions on chromosomes 3 and 8. The most significant genomic region (20.7-21.3 Mbp) on chromosome 8 was a novel resistance locus where no Rps gene was previously reported. Instead, multiple copies of the UDP-glycosyltransferase superfamily protein-coding gene, associated with disease resistance, were annotated in this new locus. Another genomic region on chromosome 3 was a well-known Rps cluster. Using the Daepung × Ilpumgeomjeong RIL population, a linkage analysis confirmed these two resistance loci and identified a resistance locus on chromosome 2. A unique feature of the resistance in Ilpumgeomjeong was discovered when phenotypic distribution was projected upon eight groups of RILs carrying different combinations of resistance alleles for the three loci. Interestingly, the seven groups carrying at least one resistance allele statistically differed from the other with none, regardless of the number of resistance alleles. This suggests that the respective three different resistance genes can confer resistance to P. sojae isolate 2457. Deployment of the three regions via marker-assisted selection will facilitate effectively improving resistance to particular P. sojae isolates in soybean breeding programs.

A two-step, two-sample Mendelian randomization analysis investigating the interplay between gut microbiota, immune cells, and melanoma skin cancer.

This study aims to rigorously explore the potential causal relationships among gut microbiota (GM), immune cells, and melanoma skin cancer among participants from Europe, where this disease exhibits significant prevalence and profound societal impact. Using the genome-wide association analysis database, a double-sample Mendelian randomization (MR) analysis was drawn upon to investigate GM, immune cells, and melanoma skin cancer. The inverse variance weighted approach was applied to estimate the causal connections among these variables. A two-step MR analysis was employed to quantitatively gauge the impact of immune cells mediated GM on melanoma skin cancer. To address potential sources of bias, such as pleiotropy and heterogeneity, multiple analytical techniques were integrated. The MR analysis pinpointed 6 GM taxa related to either an augmented or declined risk of late-stage melanoma skin cancer. In the same vein, 32 immune cell phenotypes were noticed as correlates with modified risk of melanoma skin cancer. Our study also implies that the probable association between GM and melanoma could be facilitated by 5 immune cell phenotypes. The findings of our study underline certain GM taxa and immune cells as potential influencers on the onset and development of melanoma skin cancer. Importantly, our results spotlight 5 immune cell phenotypes as potential agents mediating this association.

Genome-wide meta-analysis conducted in three large biobanks expands the genetic landscape of lumbar disc herniations

Given that lumbar disc herniation (LDH) is a prevalent spinal condition that causes significant individual suffering and societal costs, the genetic basis of LDH has received relatively little research. Our aim is to increase understanding of the genetic factors influencing LDH. We perform a genome-wide association analysis (GWAS) of LDH in the FinnGen project and in Estonian and UK biobanks, followed by a genome-wide meta-analysis to combine the results. In the meta-analysis, we identify 41 loci that have not been associated with LDH in prior studies on top of the 23 known risk loci. We detect LDH-associated loci in the vicinity of genes related to inflammation, disc-related structures, and synaptic transmission. Overall, our research contributes to a deeper understanding of the genetic factors behind LDH, potentially paving the way for the development of new therapeutics, prevention methods, and treatments for symptomatic LDH in the future.

Evaluation of genetic diversity and genome-wide association studies of resistance to bacterial wilt disease in potato.

The development of novel improved varieties adapted to unstable environmental conditions is possible through the genetic diversity of breeding materials. Potato is among the most important food crops worldwide, however, there are still significant hindrances to breeding gains attributed to its autotetraploid and highly heterozygous genome. Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) is an important disease affecting potato among many economically important crops worldwide. No cultivated potato genotypes have shown a satisfactory level of resistance to bacterial wilt. Nevertheless, resistance can play a crucial role in effective integrated disease management. To understand the genetic landscape of bacterial wilt resistance in cultivated potato, we evaluated the diversity of 194 accessions from the International Potato Centre (CIP) using 9,250 single nucleotide polymorphisms (SNPs) and their associations to the response to bacterial wilt disease evaluated over two independent trials. Twenty-four accessions showed high resistance throughout both trials. Genetic diversity analysis revealed three major clusters whose subgroups were mostly represented by CIP clones derived from common parents. Genome-wide association analyses have shown six major hits: two on chromosome 8, and one on each chromosome 2, 4, 5, and 9. These results facilitate genetic dissection of bacterial wilt resistance and marker-enabled breeding in elite genotypes for potato breeding initiatives.

Genome Wide Association Analysis for Uniform Coleoptiles Emergence and Early Seedling Growth in Rice

Early seedling vigour (ESV) is a complex trait in rice. Detection of QTLs/genes controlling these traits can help us in enhancing the yield potential in rice varieties. Association mapping is a technique based on the principle of linkage disequilibrium that is used to find genes or quantitative trait loci (QTL) underlying the complex traits. In this study of haplotype breeding 281 rice genotypes were taken. ANOVA analysis showed P- value for traits and genotypes was found significant. Similarly, P-value for interaction between the traits and genotypes was also found to be highly significant (1.8663*10-208). Further, mean germination data positively correlated with mean shoot length, mean leaf number, mean culm diameter, mean shoot dry weight, and mean shoot area of 21st day among the 281 genotypes. Among the 281 number of genotypes, 111 genotypes are found to be in PCA1 and 170 genotypes are found to be in PCA2 based on the phenotypic analysis. PCA1 component constituted 29.93% and PCA2 constituted the 13.68% of total variation in the analysis. Besides, whole genome phylogenetic analysis showed three major groups of which Group 1 consists of 215 genotypes, group 2 consists of 38 genotypes and group 3 consists of 28 genotypes respectively. Especially, both the subgroups II and III comprised of the unique genotypes from the indica and aus subpopulations of rice. In this analysis, 16 significant associations (LOD Score &gt;7) for different traits were identified using the three different models (MLM, farmCPU, and blink) for GWAS studies Especially, one major QTL was identified for the mean coleoptiles’ emergence for 10 DAS on 11th chromosome (18983591) which explained 49% of the phenotypic variance. Additionally, another major QTL contributing to the shoot length variation of 29.75% was identified in the Chr02 (32954393) for shoot length trait on 28 DAS. A candidate gene namely Os02g0778400 UMP/CMP kinase A/adenylate kinase (LOC_Os02g53790) was located in the significant SNP region of the GWAS analysis. Further characterization of this gene would assist in elucidation of the mechanism regulating the early seedling length in rice under direct seeded rice.

Drosophila Toxicogenomics: genetic variation and sexual dimorphism in susceptibility to 4-Methylimidazole

Background4-methylimidazole is a ubiquitous and potentially carcinogenic environmental toxicant. Genetic factors that contribute to variation in susceptibility to its toxic effects are challenging to assess in human populations. We used the Drosophila melanogaster Genetic Reference Panel (DGRP), a living library of natural genetic variation, to identify genes with human orthologs associated with variation in susceptibility to 4-methylimidazole.ResultsWe screened 204 DGRP lines for survival following 24-hour exposure to 4-methylimidazole. We found extensive genetic variation for survival, with a broad sense heritability of 0.82; as well as genetic variation in sexual dimorphism, with a cross-sex genetic correlation of 0.59. Genome-wide association analyses identified a total of 241 candidate molecular polymorphisms in or near 273 unique genes associated with survival. These polymorphisms had either sex-specific or sex-antagonistic effects, and most had putative regulatory effects. We generated interaction networks using these candidate genes as inputs and computationally recruited genes with known physical or genetic interactions. The network genes were significantly over-represented for gene ontology terms involving all aspects of development (including nervous system development) and cellular and organismal functions as well as canonical signaling pathways, and most had human orthologs.ConclusionsThe genetic basis of variation in sensitivity to acute exposure to 4-methylimidazole in Drosophila is attributable to variation in genes and networks of genes known for their effects on multiple developmental and cellular processes, including possible neurotoxicity. Given evolutionary conservation of the underlying genes and pathways, these insights may be applicable to humans.

Genetic parameters and genome-wide association analyses for lifetime productivity in Chinese Holstein cattle

Lifetime productivity is a trait of great importance to dairy cattle populations, as it combines information from production and longevity variables. Therefore, we investigated the genetic background of lifetime productivity in high-producing dairy cattle by integrating genomics and transcriptomics datasets. A total of 3,365,612 test-day milk yield records from 134,029 Chinese Holstein cows were used to define 6 lifetime productivity traits, including lifetime milk yield covering full lifespan and 5 cumulative milk yield traits covering partial lifespan. Genetic parameters were estimated based on univariate and bivariate linear animal models and the restricted maximum likelihood method. Genome-wide association studies and weighted gene co-expression network analyses (WGCNA) were performed to identify candidate genes associated with lifetime productivity based on genomic data from 3,424 cows and peripheral blood RNA-sequencing data from 23 cows, respectively. Lifetime milk yield averaged 24,800.8 ± 14,396.6 kg (mean ± SD) across an average of 2.4 parities in Chinese Holstein population. The heritability estimates for lifetime productivity traits ranged from 0.05 (±0.01 for SE) to 0.10 (±0.02 for SE). The estimate of genetic correlation between lifetime milk yield and productive life is 0.88 (±0.3 for SE), whereas the genetic correlation with 305-d milk yield in the first lactation was 0.49 (±0.08 for SE). Absolute values for most genetic correlation estimates between lifetime productivity and type traits were lower than 0.30. Moderate genetic correlations were found between udder related traits and lifetime productivity, such as with udder depth (0.33), rear udder attachment height (0.33), and udder system (0.34). Some single nucleotide polymorphisms and gene co-expression modules significantly associated with lifetime milk yield were identified based on GWAS and WGCNA analyses, respectively. Functional enrichment analyses of the candidate genes identified revealed important pathways related to immune system, longevity, energy utilization, and metabolism, and FoxO signaling. The genes NTMT1, FNBP1, and S1PR1 were considered to be the most important candidate genes influencing lifetime productivity in Holstein cows. Overall, our findings indicate that lifetime productivity is heritable in Chinese Holstein cattle, and important candidate genes were identified by integrating genomic and transcriptomic datasets.

Genomic regions associated with embryonic loss in primiparous Holstein cows

IntroductionReproductive efficiency in dairy cattle bred by artificial insemination (AI) or via embryo transfer (ET) is important as a means of facilitating the rate of genetic improvement by utilizing superior genetics from sires and donor dams. Fertility is especially challenging in primiparous cows as they have recently calved, are still growing, and are lactating. Fertility, measured by the conception rate from the first service (CCR1), is the number of cows that become pregnant as a result of the first AI or ET compared to all cows serviced, whereas times bred (TBRD) is the number of inseminations cows bred by AI needed to achieve a successful pregnancy.MethodsThis study aimed to identify the loci associated with and gene sets enriched for the CCR1 in primiparous cows that were bred by AI (n = 882) or were ET recipients (n = 353), and to identify the loci associated with the TBRD (n= 1,851) in primiparous cows. A genome-wide association analysis and gene set enrichment analysis were performed using single-nucleotide polymorphism (SNP) data for the CCR1 and TBRD in primiparous cows. In cows bred by AI, two loci containing five positional candidate genes were associated with the CCR1 (p &amp;lt; 1 × 10−5) and one gene set with 18 leading-edge genes was enriched (NES ≥3) for the CCR1.Results and discussionIn cows that were ET recipients, five loci were associated with the CCR1, but no gene sets were enriched for the CCR1. For the second phenotype, five loci with nine positional candidate genes were associated with the TBRD and two gene sets with eight leading-edge genes were enriched for the TBRD. Identifying the loci associated with or gene sets enriched for the CCR1 and TBRD provides an opportunity to better understand embryonic loss in AI-bred and ET-recipient cows and to use this information to select for higher reproductive efficiency.

Genome-wide dissection of genes shaping inflorescence morphology in 242 Chinese south–north sorghum accessions

The inflorescences morphology (IM) of sorghum (Sorghum bicolor L. Moench) affects its resistance to pests, diseases, and grain yields. However, the specific genetic factors underlying in IM are not yet fully elucidated. Here we conducted a comprehensive genome-wide association analysis (GWAS) to identify the stable and adaptive Quantitative Trait Loci (QTL) for five IM traits (panicle length, the number of cob nodes, the number of primary branches, the largest length of the primary branch, and panicle type) in a sorghum panel, which adapted to different environments from the south to north in China. Totally, 2,015,850 high quality single nucleotide polymorphisms (SNPs) were obtained. Population structure analysis showed that two distinct genetic sub-populations were divided according to their geographic origin. Seventy-one QTLs distributed in 41 genetic regions on 9 chromosomes were identified. These regions harbored 21 high-confident candidate genes that were homologous to rice domestication genes, including 7 related to IM. Two domestication-related genes (Sobic.003G052700 and Sobic.006G247700) were located into two major QTL regions (QTL3.4721839 and QTL6.58709500) which were identified in multi-environments. Allelic variations in the two genes displayed a geographical pattern, indicating that different IM traits were selected by south and north sorghum breeders, such as south sorghums had long and loose panicles in order to adapt the hot and humid climate, while north sorghums had short and compact panicle to increase planting density and grain yield per unit area due to dry climate. This work provides new breeding strategies and resources for developing locally adapted sorghum varieties.

Hypnotic use and the risk of cardiovascular diseases in insomnia patients

Abstract Background Insomnia, commonly treated with hypnotic agents, is an independent risk factor of cardiovascular diseases. Previous research demonstrated the detrimental effect of hypnotics on the prognosis of cardiovascular patients. Purpose We aimed to examine the association between hypnotic agents and cardiovascular outcomes in general individuals with insomnia. Methods In a propensity score matched cohort of UK Biobank (UKB) participants with insomnia, Cox proportional hazard model was used to estimate the association between regular use of hypnotic agents and predetermined cardiovascular outcomes including incident coronary heart diseases (CHD), heart failure (HF), stroke, and cardiovascular death. Inverse probability of treatment weighting and competing risk models were further performed during sensitivity analysis. Drug-target Mendelian randomization (MR) analyses were employed for further evaluation of hypnotics-related genes and cardiovascular disease association. The genome-wide association data of CHD, HF, and stroke were there large-scale genome-wide association analyses. Tissue-specified expression quantitative trait loci (eQTL) data was derived from BrainSeq phase II. Results During a median follow-up of 14.3 years, the matched cohort documented a total of 1,019 CHD cases, 406 HF cases, 285 stroke cases, and 197 cardiovascular deaths. No significant association was detected between Z-meds and CHD, stroke, and cardiovascular mortality. Benzodiazepine use was significantly associated with the increased risk of CHD, HF, and cardiovascular mortality (Figure 1). The inverse probability of treatment weighting and competing risk models didn’t alter the above associations. Moreover, drug-target MR analyses corroborated the safety of Z-meds in the general population regarding cardiovascular health. However, GABRG1—a drug target for benzodiazepines—was associated with heightened risks for heart failure, and ischemic stroke (Figure 2). Conclusions Our findings suggested the heterogeneous associations between different categories of hypnotics and incident cardiovascular events in individuals with insomnia. Z-meds are safe regarding the risk of cardiovascular outcomes in the UKB population with insomnia.Figure 1Figure 2

Genome-wide association analysis reveals insights into the molecular etiology of dilated cardiomyopathy

Abstract Introduction Dilated cardiomyopathy (DCM) is a primary heart muscle disease characterized by ventricular dilatation and impaired contractility. While considered a Mendelian disorder, recent evidence highlights the important role of common genetic variation. Better characterization of the genetic architecture will provide opportunities for improved diagnosis, risk stratification and biological understanding. Methods We performed a GWAS meta-analysis of 14,256 DCM cases and 1,185,671 controls of European ancestry from 16 studies (GWASDCM). To improve discovery power, we performed multi-trait analysis of GWAS (GWASMTAG) with three genetically correlated and clinically relevant cardiac imaging traits in 36,203 individuals. To identify effector genes at loci, we used a two-step approach: 1) the nearest gene and top 3 genes prioritized by PoPS or V2G were selected as candidates; and 2) totality of evidence including 5 additional methods (coding variant, co-localization with expression, TWAS, ABC-model, and known Mendelian cardiomyopathy genes) was summarized in a single score. We performed rare variant burden analysis (MAF&amp;lt;0.1%) for predicted truncating variants (PTV) in established DCM-causing genes to highlight DCM genetic architecture, and in prioritized GWAS genes, against DCM and quantitative CMR traits in UK Biobank (UKB) and 100,000 Genomes Project (GeL), implemented by REGENIE. A polygenic score (PGS) generated using a Bayesian framework (PRS-CS) was assessed for association with the same outcomes in UKB participants, including among PTV carriers in established DCM-causing genes. Results We identified 80 independent loci, of which 65 were novel (62 loci in GWASDCM, and 54 loci in GWASMTAG) (Figure 1A). High-confidence genes were prioritized at 62 loci, including in genes with established roles in Mendelian cardiomyopathies (TTN, BAG3, FLNC, MYBPC3, FHOD3, ALPK3), and were enriched in key biological processes (Figure 1B). We characterized DCM genetic architecture across the allelic spectrum (Figure 2A). Rare variant burden analysis highlighted novel associations with DCM and quantitative traits for PTVs in NEDD4L and MAP3K7 in UKB that were replicated in GeL (Figure 2B). PGS was associated with DCM and quantitative CMR traits (OR per PGS SD 1.8, P&amp;lt;1x10-16) (Figure 2C). Individuals in the top centile had 4-fold increased risk of having DCM compared with median risk. PGS stratified penetrance of rare variants in 1,546 carriers of pathogenic variants (top quintile 7.3%, bottom quintile 1.7%, P 0.005). Conclusion We have performed the largest DCM GWAS, resulting in 80 significant loci (65 novel), and prioritized genes using a systematic approach. We identify novel DCM-causing genes (NEDD4L and MAP3K7) and generate a PGS that associates with DCM and modulates penetrance in carriers of Mendelian disease-causing rare variants.Figure 1Figure 2