Genome-wide Association Analysis Research Articles

Genome-wide association analyses of breast cancer in women of African ancestry identify new susceptibility loci and improve risk prediction.

We performed genome-wide association studies of breast cancer including 18,034 cases and 22,104 controls of African ancestry. Genetic variants at 12 loci were associated with breast cancer risk (P < 5 × 10-8), including associations of a low-frequency missense variant rs61751053 in ARHGEF38 with overall breast cancer (odds ratio (OR) = 1.48) and a common variant rs76664032 at chromosome 2q14.2 with triple-negative breast cancer (TNBC) (OR = 1.30). Approximately 15.4% of cases with TNBC carried six risk alleles in three genome-wide association study-identified TNBC risk variants, with an OR of 4.21 (95% confidence interval = 2.66-7.03) compared with those carrying fewer than two risk alleles. A polygenic risk score (PRS) showed an area under the receiver operating characteristic curve of 0.60 for the prediction of breast cancer risk, which outperformed PRS derived using data from females of European ancestry. Our study markedly increases the population diversity in genetic studies for breast cancer and demonstrates the utility of PRS for risk prediction in females of African ancestry.

Sweet corn association panel and genome-wide association analysis reveal loci for chilling-tolerant germination

Sweet corn is highly susceptible to the deleterious effects of low temperatures during the initial stages of growth and development. Employing a 56K chip, high-throughput single-nucleotide polymorphism (SNP) sequencing was conducted on 100 sweet corn inbred lines. Subsequently, six germination indicators—germination rate, germination index, germination time, relative germination rate, relative germination index, and relative germination time—were utilized for genome-wide association analysis. Candidate genes were identified via comparative analysis of homologous genes in Arabidopsis and rice, and their functions were validated using quantitative real-time polymerase chain reaction (qRT-PCR). The results revealed 35,430 high-quality SNPs, 16 of which were significantly correlated. Within 50 kb upstream and downstream of the identified SNPs, 46 associated genes were identified, of which six were confirmed as candidate genes. Their expression patterns indicated that Zm11ΒHSDL5 and Zm2OGO likely play negative and positive regulatory roles, respectively, in the low-temperature germination of sweet corn. Thus, we determined that these two genes are responsible for regulating the low-temperature germination of sweet corn. This study contributes valuable theoretical support for improving sweet corn breeding and may aid in the creation of specific germplasm resources geared toward enhancing low-temperature tolerance in sweet corn.

Large-scale transcriptomic and genomic analyses reveal a novel functional gene SERPINB6 for chicken carcass traits

BackgroundCarcass traits are crucial indicators of meat production efficiency. However, the molecular regulatory mechanisms associated with these traits remain unclear.ResultsIn this study, we conducted comprehensive transcriptomic and genomic analyses on 399 Tiannong partridge chickens to identify key genes and variants associated with carcass traits and to elucidate the underlying regulatory mechanisms. Based on association analyses with the elastic net (EN) model, we identified 12 candidate genes (AMY1A, AP3B2, CEBPG, EEF2, EIF4EBP1, FGFR1, FOXD3, GOLM1, LOC107052698, PABPC1, SERPINB6 and TBC1D16) for 4 carcass-related traits, namely live weight, dressed weight, eviscerated weight, and breast muscle weight. SERPINB6 was identified as the only overlapping gene by 3 analyses, EN model analysis, weighted gene co-expression network analysis and differential expression analysis. Cell-level experiments confirmed that SERPINB6 promotes the proliferation of chicken DF1 cells and primary myoblasts. Further expression genome-wide association study and association analysis indicated that rs317934171 is the critical site that enhances SERPINB6 expression. Furthermore, a dual-luciferase reporter assay proved that gga-miR-1615 targets the 3′UTR of SERPINB6.ConclusionsCollectively, our findings reveal that SERPINB6 serves as a novel gene for chicken carcass traits by promoting fibroblast and myoblast proliferation. Additionally, the downstream variant rs317934171 regulates SERPINB6 expression. These results identify a new target gene and molecular marker for the molecular mechanisms of chicken carcass traits.

Predictive genome-wide association and economic analyses of pooled DNA in a commercial feedlot setting

Understanding the genetic profile of cattle allows producers to make predictions of performance and disease. However, it is often fiscally and logistically unfeasible for commercial cattle producers to individually genotype animals. As such, our objec­tive was to utilize prospectively-pooled DNA from commercial feedlot cattle for a genome-wide association study (GWAS) of ec­onomically-important carcass traits and measure the economic efficiency of this technique.

Genetic Diversity and Genome-Wide Association Analysis of the Hulled/Naked Trait in a Barley Collection from Shanghai Agricultural Gene Bank

Genetic Diversity and Genome-Wide Association Analysis of the Hulled/Naked Trait in a Barley Collection from Shanghai Agricultural Gene Bank

Genomic epidemiology of Staphylococcus aureus from the Iberian Peninsula highlights the expansion of livestock associated-CC398 towards wildlife

Staphylococcus aureus is a versatile pathobiont, exhibiting a broad host range, including humans, other mammals, and avian species. Host specificity determinants, virulence, and antimicrobial resistance genes are often shared by strains circulating at the animal-human interface. While transmission dynamics studies have shown strain exchange between humans and livestock, knowledge of the source, genetic diversification, and transmission drivers of S. aureus in wildlife lag behind.In this work, we explore a wide array of S. aureus genomes from different sources in the Iberian Peninsula to understand population structure, gene content and niche adaptation at the human-livestock-wildlife nexus. Through Bayesian inference, we address the hypothesis that S. aureus strains in wildlife originate from humanized landscapes, either from contact with humans or through interactions with livestock. Phylogenetic reconstruction applied to whole genome sequence data was completed with a dataset of 450 isolates featuring multiple clones from the 1990–2022 period and a subset of CC398 strains representing the 2008–2022 period. Phylodynamic signatures of S. aureus from the Iberian Peninsula suggest widespread circulation of most clones among humans before jumping to other hosts. The number of transitions of CC398 strains within each host category (human, livestock, wildlife) was high (88.26 %), while the posterior probability of transitions from livestock to wildlife was remarkably high (0.99). Microbial genome-wide association analysis did not evidence genome rearrangements nor biomarkers suggesting S. aureus niche adaptation to wildlife, thus supporting recent spill overs. Altogether, our findings indicate that S. aureus isolates collected in the past years from wildlife most likely represent multiple introduction events from livestock. The clonal origin of CC398 and its potential to disseminate and evolve through different animal host species are highlighted, calling for management practices at the livestock-wildlife axis to improve biosecurity and thus restrict S. aureus transmission and niche expansion along gradients of human influence.

Exploring the relationship between admixture and genetic susceptibility to attention deficit hyperactivity disorder in two Latin American cohorts.

Contemporary research on the genomics of Attention Deficit Hyperactivity Disorder (ADHD) often underrepresents admixed populations of diverse genomic ancestries, such as Latin Americans. This study explores the relationship between admixture and genetic associations for ADHD in Colombian and Mexican cohorts. Some 546 participants in two groups, ADHD and Control, were genotyped with Infinium PsychArray®. Global ancestry levels were estimated using overall admixture proportions and principal component analysis, while local ancestry was determined using a method to estimate ancestral components along the genome. Genome-wide association analysis (GWAS) was conducted to identify significant associations. Differences between Colombia and Mexico were evaluated using appropriate statistical tests. 354 Single-nucleotide polymorphisms (SNPs) and Single-nucleotide variants (SNVs) related to some genes and intergenic regions exhibited suggestive significance (p-value < 5*10e-5) in the GWAS. None of the variants revealed genome-wide significance (p-value < 5*10e-8). The study identified a significant relationship between risk SNPs and the European component of admixture, notably observed in the LOC105379109 gene. Despite differences in risk association loci, such as FOXP2, our findings suggest a possible homogeneity in genetic variation's impact on ADHD between Colombian and Mexican populations. Current reference datasets for ADHD predominantly consist of samples with high European ancestry, underscoring the need for further research to enhance the representation of reference populations and improve the identification of ADHD risk traits in Latin Americans.

Genome wide association analysis reveals the genetic sites of heat resistance in largemouth bass (Micropterus salmoides)

Largemouth bass (Micropterus salmoides) is an important commercially cultured freshwater fish. High temperature is a common challenge encountered in the summer cultivation of largemouth bass. In this study, we conducted a genome-wide association study using general linear models and linear mixed models to identify significant single nucleotide polymorphisms (SNPs) and candidate genes associated with loss of equilibrium (LOE) and weight gain (WG) during high temperature challenge. We collected fin tissue from 353 individuals of largemouth bass exposed to high water temperature (37.5 °C) and performed whole-genome resequencing. Through association analysis, we successfully identified 10 SNPs and 38 candidate genes that were significantly associated with LOE, and four SNPs and 19 candidate genes that were significantly associated with WG. Using quantitative PCR and Western blot, we tested the expression of eight genes related to LOE in heat-sensitive and heat-tolerant individuals, and identified significant differences in the expression of the following five genes: Tmub2, Ngfr, Hgf, Hdac6, and Sypl. For WG, we assessed the expression of three candidate genes, among which Mrc2l displayed significant differences in expression between the high WG group and low WG group. We discuss the functions and potential biological processes of candidate genes for LOE and WG. Our results provide basic genetic information for breeding largemouth bass with high temperature tolerance.

Variation in WIDTH OF LEAF AND GRAIN contributes to grain and leaf size by controlling LARGE2 stability in rice.

Grain and flag leaf size are two important agronomic traits that influence grain yield in rice (Oryza sativa). Many QTLs and genes that regulate these traits individually have been identified, however, few QTLs and genes that simultaneously control these two traits have been identified. In this study, we conducted a genome-wide association analysis in rice and detected a major locus, WIDTH OF LEAF AND GRAIN (WLG), that associated with both grain and flag leaf width. WLG encodes a RING-domain E3 ubiquitin ligase. WLGhap.B, which possesses five SNP variations compared to WLGhap.A, encodes a protein with enhanced ubiquitination activity that confers increased rice leaf width and grain size, whereas mutation of WLG leads to narrower leaves and smaller grains. Both WLGhap.A and WLGhap.B interact with LARGE2, a HETC-type E3 ligase, however, WLGhap.B exhibits stronger interaction with LARGE2, thus higher ubiquitination activity towards LARGE2 compared with WLGhap.A. Lysine1021 is crucial for the ubiquitination of LARGE2 by WLG. Loss-of-function of LARGE2 in wlg-1 phenocopies large2-c in grain and leaf width, suggesting that WLG acts upstream of LARGE2. These findings reveal the genetic and molecular mechanism by which the WLG-LARGE2 module mediates grain and leaf size in rice, and suggest the potential of WLGhap.B in improving rice yield.

Genome-Wide Association Study Identifies IFIH1 and HLA-DQB1*05:02 Loci Associated With Anti-NMDAR Encephalitis.

Anti-N-methyl-d-aspartate receptor (NMDAR) encephalitis is a rare autoimmune neurologic disorder, the genetic etiology of which remains poorly understood. Our study aims to investigate the genetic basis of this disease in the Chinese Han population. We performed a genome-wide association study and fine-mapping study within the major histocompatibility complex (MHC) region of 413 Chinese patients with anti-NMDAR encephalitis recruited from 6 large tertiary hospitals and 7,127 healthy controls. Our genome-wide association analysis identified a strong association at the IFIH1 locus on chromosome 2q24.2 (rs3747517, p = 1.06 × 10-8, OR = 1.55, 95% CI, 1.34-1.80), outside of the human leukocyte antigen (HLA) region. Furthermore, through a fine-mapping study of the MHC region, we discovered associations for 3 specific HLA class I and II alleles. Notably, HLA-DQB1*05:02 (p = 1.43 × 10-12; OR, 2.10; 95% CI 1.70-2.59) demonstrates the strongest association among classical HLA alleles, closely followed by HLA-A*11:01 (p = 4.36 × 10-7; OR, 1.52; 95% CI 1.29-1.79) and HLA-A*02:07 (p = 1.28 × 10-8; OR, 1.87; 95% CI 1.50-2.31). In addition, we uncovered 2 main HLA amino acid variation associated with anti-NMDAR encephalitis including HLA-DQβ1-126H (p = 1.43 × 10-12; OR, 2.10; 95% CI 1.70-2.59), exhibiting a predisposing effect, and HLA-B-97R (p = 3.40 × 10-8; OR, 0.63; 95% CI 0.53-0.74), conferring a protective effect. Computational docking analysis suggested a close relationship between the NR1 subunit of NMDAR and DQB1*05:02. Our findings indicate that genetic variation in IFIH1, involved in the type I interferon signaling pathway and innate immunity, along with variations in the HLA class I and class II genes, has substantial implications for the susceptibility to anti-NMDAR encephalitis in the Chinese Han population.

Integrative multi-omics analyses to identify the genetic and functional mechanisms underlying ovarian cancer risk regions

To identify credible causal risk variants (CCVs) associated with different histotypes of epithelial ovarian cancer (EOC), we performed genome-wide association analysis for 470,825 genotyped and 10,163,797 imputed SNPs in 25,981 EOC cases and 105,724 controls of European origin. We identified five histotype-specific EOC risk regions (p value <5×10-8) and confirmed previously reported associations for 27 risk regions. Conditional analyses identified an additional 11 signals independent of the primary signal at six risk regions (pvalue <10-5). Fine mapping identified 4,008 CCVs in these regions, of which 1,452 CCVs were located in ovarian cancer-related chromatin marks with significant enrichment in active enhancers, active promoters, and active regions for CCVs from each EOC histotype. Transcriptome-wide association and colocalization analyses across histotypes using tissue-specific and cross-tissue datasets identified 86 candidate susceptibility genes in known EOC risk regions and 32 genes in 23 additional genomic regions that may represent novel EOC risk loci (false discovery rate <0.05). Finally, by integrating genome-wide HiChIP interactome analysis with transcriptome-wide association study (TWAS), variant effect predictor, transcription factor ChIP-seq, and motifbreakR data, we identified candidate gene-CCV interactions at each locus. This included risk loci where TWAS identified one or more candidate susceptibility genes (e.g., HOXD-AS2, HOXD8, and HOXD3 at 2q31) and other loci where no candidate gene was identified (e.g., MYC and PVT1 at 8q24) by TWAS. In summary, this study describes a functional framework and provides a greater understanding of the biological significance of risk alleles and candidate gene targets at EOC susceptibility loci identified by a genome-wide association study.

Causal relationship between intervertebral disc degeneration and osteoporosis: a bidirectional two-sample Mendelian randomization study.

The relationship between intervertebral disc degeneration (IVDD) and osteoporosis (OP), diagnosed primarily using bone mineral density (BMD), remains unclear so far. The present study, therefore, aimed to investigate the potential relationship between osteoporosis and intervertebral disc degeneration using Mendelian randomization and genome-wide association analyses. Specifically, the impact of bone mineral density on the development of intervertebral disc degeneration was evaluated. The genome-wide association studies (GWAS) summary data of OP/BMDs and IVDD were collected from the FinnGen consortium, the GEFOS consortium, and MRC-IEU. The relationship between IVDD and OP was then explored using TSMR. The inverse-variance weighted (IVW) method was adopted as the primary effect estimate, and the reliability and stability of the results were validated using various methods, including MR-Egger, weighted median, simple mode, weighted mode, and MR-PRESSO. No significant causal relationship was observed between OP and IVDD (IVW, P > 0.05) or between femoral neck BMD (FA-BMD) and IVDD when OP and FA-BMD were used as exposures. However, increased levels of total body BMD (TB-BMD) and lumbar spine BMD (LS-BMD) were revealed as significant risk factors for IVDD (TB-BMD: IVW, OR = 1.201, 95% CI: 1.123-1.284, P = 8.72 × 10-8; LS-BMD: IVW, OR = 1.179, 95% CI: 1.083-1.284, P = 1.43 × 10-4). Interestingly, both heel BMD (eBMD) and femur neck BMD (FN-BMD) exhibited potential causal relationships (eBMD: IVW, OR = 1.068, 95% CI: 1.008-1.131, P = 0.0248; FN-BMD, IVW, OR = 1.161, 95% CI: 1.041-1.295, P = 0.0074) with the risk of IVDD. The reverse MR analysis revealed no statistically causal impact of IVDD on OP and the level of BMD (P > 0.05). OP and the level of FA-BMD were revealed to have no causal relationship with IVDD. The increased levels of TB-BMD and LS-BMD could promote the occurrence of IVDD. Both eBMD and FN-BMD have potential causal relationships with the risk of IVDD. No significant relationship exists between IVDD and the risk of OP. Further research is warranted to comprehensively comprehend the molecular mechanisms underlying the impact of OP and BMD on IVDD and vice versa.

Genome-wide association studies identifies genetic loci related to fatty acid and branched-chain amino acid metabolism and histone modifications under varying nitrogen treatments in safflower (Carthamus tinctorius).

Effective identification and usage of genetic variation are prerequisites for developing nutrient-efficient cultivars. A collection of 94 safflower (Carthamus tinctorius ) genotypes (G) was investigated for important morphological and photosynthetic traits at four nitrogen (N) treatments. We found significant variation for all the studied traits except chlorophyll b (chl b ) among safflower genotypes, nitrogen treatments and G×N interaction. The examined traits showed a 2.82-50.00% increase in response to N application. Biological yield (BY) reflected a significantly positive correlation with fresh shoot weight (FSW), root length (RL), fresh root weight (FRW) and number of leaves (NOL), while a significantly positive correlation was also observed among carotenoids (C), chlorophyll a (chl a ), chl b and total chlorophyll content (CT) under all treatments. Superior genotypes with respect to plant height (PH), FSW, NOL, RL, FRW and BY were clustered into Group 3, while genotypes with better mean performance regarding chl a , chl b C and CT were clustered into Group 2 as observed in principal component analysis. The identified eight best-performing genotypes could be useful to develop improved nitrogen efficient cultivars. Genome-wide association analysis resulted in 32 marker-trait associations (MTAs) under four treatments. Markers namely DArT-45481731 , DArT-17812864 , DArT-15670279 and DArT-45482737 were found consistent. Protein-protein interaction networks of loci associated with MTAs were related to fatty acid and branched-chain amino acid metabolism and histone modifications.

The role of dietary preferences in osteoarthritis: a Mendelian randomization study using genome-wide association analysis data from the UK Biobank.

To understand the impact of individual preferences for specific dietary items on OA, and to help inform the development of effective and targeted OA prevention and management strategies, we performed a Mendelian randomization analysis between dietary preferences and osteoarthritis. This study utilized genetic data from the UK Biobank to investigate the association between OA and 21 different common dietary items. Instrumental variables representing European populations were carefully selected based on their genetic significance and linkage disequilibrium. In cases where a dietary item had few relevant genetic markers, a more lenient selection threshold was applied. To prevent bias, the analysis excluded single nucleotide polymorphisms (SNPs) associated with factors such as body mass index (BMI) and cholesterol. Using inverse-variance weighting (IVW) and Mendelian randomization, significant associations were detected between certain dietary items and OA. Using Mendelian randomization to examine the relationship between 21 different dietary items and OA, significant associations were found for coffee, peas, watercress, and cheese, where the first two had a promoting effect and the last two an inhibiting effect on OA. Due to heterogeneity in the test results for cheese, a random IVW representation was used. The results of sensitivity analysis showed no significant heterogeneity or horizontal pleiotropy in the selected SNPS, demonstrating the reliability of Mendelian randomization analysis. This study identified coffee, peas, watercress, and cheese as food items that may have significant dietary effects on osteoarthritis. This information may be useful to consider in the development of OA management strategies.

An autosomal recessive variant in PYGM causes myophosphorylase deficiency in Red Angus composite cattle

BackgroundBetween 2020 and 2022, eight calves in a Nebraska herd (composite Simmental, Red Angus, Gelbvieh) displayed exercise intolerance during forced activity. In some cases, the calves collapsed and did not recover. Available sire pedigrees contained a paternal ancestor within 2–4 generations in all affected calves. Pedigrees of the calves’ dams were unavailable, however, the cows were ranch-raised and retained from prior breeding seasons, where bulls used for breeding occasionally had a common ancestor. Therefore, it was hypothesized that a de novo autosomal recessive variant was causative of exercise intolerance in these calves.ResultsA genome-wide association analysis utilizing SNP data from 6 affected calves and 715 herd mates, followed by whole-genome sequencing of 2 affected calves led to the identification of a variant in the gene PYGM (BTA29:g.42989581G > A). The variant, confirmed to be present in the skeletal muscle transcriptome, was predicted to produce a premature stop codon (p.Arg650*). The protein product of PYGM, myophosphorylase, breaks down glycogen in skeletal muscle. Glycogen concentrations were fluorometrically assayed as glucose residues demonstrating significantly elevated glycogen concentrations in affected calves compared to cattle carrying the variant and to wild-type controls. The absence of the PYGM protein product in skeletal muscle was confirmed by immunohistochemistry and label-free quantitative proteomics analysis; muscle degeneration was confirmed in biopsy and necropsy samples. Elevated skeletal muscle glycogen persisted after harvest, resulting in a high pH and dark-cutting beef, which is negatively perceived by consumers and results in an economic loss to the industry. Carriers of the variant did not exhibit differences in meat quality or any measures of animal well-being.ConclusionsMyophosphorylase deficiency poses welfare concerns for affected animals and negatively impacts the final product. The association of the recessive genotype with dark-cutting beef further demonstrates the importance of genetics to not only animal health but to the quality of their product. Although cattle heterozygous for the variant may not immediately affect the beef industry, identifying carriers will enable selection and breeding strategies to prevent the production of affected calves.

Genetic dissection of yield and quality-related traits in a Colombian Andigenum potato collection, revealed by genome-wide association and genomic prediction analyses

Genetic dissection of yield and quality-related traits in a Colombian Andigenum potato collection, revealed by genome-wide association and genomic prediction analyses

Integrated genome-wide association and transcriptomic analysis to identify receptor kinase genes to stripe rust resistance in wheat germplasm from southwestern China

Stripe rust of wheat, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most important diseases of wheat worldwide. Identification of new and elite Pst-resistance loci or genes has the potential to enhance overall resistance to this pathogen. Here, we conducted an integrated genome-wide association study (GWAS) and transcriptomic analysis to screen for loci associated with resistance to stripe rust in 335 accessions from Yunnan, including 311 landraces and 24 cultivars. Based on the environmental phenotype, we identified 113 protein kinases significantly associated with Pst resistance using mixed linear model (MLM) and generalized linear model (GLM) models. Transcriptomic analysis revealed that 52 of 113 protein kinases identified by GWAS were up and down regulated in response to Pst infection. Among these genes, a total of 15 receptor kinase genes were identified associated with Pst resistance. 11 candidate genes were newly discovered in Yunnan wheat germplasm. Our results revealed that resistance alleles to stripe rust were accumulated in Yunnan wheat germplasm, implying direct or indirect selection for improving stripe rust resistance in elite wheat breeding programs.

Genome-wide meta-analysis of short-tandem repeats for Parkinson's disease risk using genotype imputation

Abstract Idiopathic Parkinson’s disease is determined by a combination of genetic and environmental factors. Recently, the first genome-wide association study on short-tandem repeats in Parkinson’s disease reported on eight suggestive short-tandem repeat-based risk loci (α=5.3×10-6), of which four were novel, i.e. they had not been implicated in Parkinson’s disease risk by genome-wide association analyses of single-nucleotide polymorphisms before. Here, we tested these eight candidate short-tandem repeats in a large, independent Parkinson’s disease case-control dataset (n=4,757). Furthermore, we combined the results from both studies by meta-analysis resulting in the largest Parkinson’s disease genome-wide association study of short-tandem repeats to date (n=43,844). Lastly, we investigated whether leading short-tandem repeat risk variants exert functional effects on gene expression regulation based on methylation quantitative trait locus data in human post-mortem brain (n=142). None of the eight previously reported short-tandem repeats were significantly associated with Parkinson’s disease in our independent dataset after multiple testing correction (α=6.25x10-3). However, we observed modest support for short tandem repeats near CCAR2 and NCOR1 in the updated meta-analyses of all available data. While the genome-wide meta-analysis did not reveal additional study-wide significant (α=6.3×10-7) short-tandem repeat signals, we identified seven novel suggestive Parkinson’s disease short-tandem repeat risk loci (α=5.3×10-6). Of these, especially a short-tandem repeat near MEIOSIN showed consistent evidence for association across datasets. CCAR2, NCOR1 and one novel suggestive locus identified here (LINC01012) emerged from colocalization analyses showing evidence for a shared causal short-tandem repeat variant affecting both Parkinson’s disease risk and cis DNA methylation in brain. Larger studies, ideally using short-tandem repeats called from whole-sequencing data, are needed to more fully investigate their role in Parkinson’s disease.

Genome-wide association analyses reveal genotype-by-environment interactions of growth and organ weights in gilthead seabream (Sparus aurata)

Gilthead seabream is one of the most important fish species for Mediterranean aquaculture. More than 80% of the genetically improved fingerlings of this species originate from a single country, Greece. A previous study revealed moderate genotype-by-environment (G x E) interactions for growth traits and organ weights between commercial production sites in Greece and Spain. The main objectives of this study were to identify QTLs and genes associated with growth traits and organ weights in gilthead seabream and to identify the biological processes and pathways controlling these traits under different temperature environments. For growth traits, we found fourteen genome-wide suggestive SNP effects in Spain. A strong peak between 1.97 Mb to 4.69 Mb on chromosome (chr) 22 was shared among three growth traits (harvest weight, fillet weight and thermal growth coefficient). Fourteen genome-wide suggestive SNP effects were also identified in Greece. However, none overlap with those found in Spain. Two SNPs on chr5 (chr5:5729959, chr5:2336903) and two SNPs on chr6 (chr6:1957732, chr6:3659811) were associated with all three growth traits. For organ weights, a total of 15 SNPs were associated above the suggestive threshold in Spain. Four SNPs (chr22:1966940, chr22:6094383, chr22:67838, chr15:16798259) are shared between viscera weight and liver weight, while one SNP (chr22:67838) is shared among viscera, liver and heart weight. Interestingly, three of these four SNPs (chr22:1966940, chr22:67838, and chr15:16798259) also overlap with significant SNPs associated with growth traits. None of the SNPs associated with viscera weight, liver weight and heart weight are shared with Greece. GO and KEGG functional enrichment analyses showed that light absorption and ECM-receptor interaction are significantly enriched with growth traits QTLs in Spain but not in Greece. For organ weights, response to stimulus is the most prominent process in Spain, while cell adhesion, and regulation of phosphorylation are more prominent in Greece. Genomic architectures for growth traits and organ weights differ between environments. These findings not only explain part of the G x E interactions, but also give insight into how genetics determines differences in how environmental conditions affect growth and organ weights.

Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder.

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.