BACKGROUND: Breeding polled goats is a welfare-friendly approach for horn removal in comparison to invasive methods. To gain a comprehensive understanding of the genetic basis underlying polledness in goats, we conducted whole-genome sequencing of 106 Xinong Saanen dairy goats, including 33 horned individuals, 70 polled individuals, and 3 polled intersexuality syndrome (PIS) individuals.
METHODS: The present study employed a genome-wide association study (GWAS) and linkage disequilibrium (LD) analysis to precisely map the genetic locus underlying the polled phenotype in goats.
RESULTS: The analysis conducted in our study revealed a total of 320 genome-wide significant single nucleotide polymorphisms (SNPs) associated with the horned/polled phenotype in goats. These SNPs exhibited two distinct peaks on chromosome 1, spanning from 128,817,052 to 133,005,441 bp and from 150,336,143 to 150,808,639 bp. The present study identified three genome-wide significant SNPs, namely Chr1:129789816, Chr1:129791507, and Chr1:129791577, as potential markers of PIS-affected goats. The results of our LD analysis suggested a potential association between MRPS22 and infertile intersex individuals, as well as a potential association between ERG and the polled trait in goats.
CONCLUSIONS: We have successfully identified three marker SNPs closely linked to PIS, as well as several candidate genes associated with the polled trait in goats. These results may contribute to the development of SNP chips for early prediction of PIS in goats, thereby facilitating breeding programs aimed at producing fertile herds with polled traits.

Faba bean is an important legume crop. The genetic diversity among faba bean genotypes is very important for the genetic improvement of target traits. A set of 128 fab bean genotypes that are originally from Egypt were used in this study to investigate the genetic diversity and population structure. The 128 genotypes were genotyped using the Single Primer Enrichment Technology (SPET) by which a set of 6759 SNP markers were generated after filtration. The SNP markers were distributed on all chromosomes with a range extending from 822 (Chr. 6) to 1872 (Chr.1). The SNP markers had wide ranges of polymorphic information content (PIC), gene diversity (GD), and minor allele frequency. The analysis of population structure divided the Egyptian faba bean population into five subpopulations. Considerable genetic distance was found among all genotypes, ranging from 0.1 to 0.4. The highly divergent genotype was highlighted in this study and the genetic distance among genotypes ranged from 0.1 and 0.6. Moreover, the structure of linkage disequilibrium was studied, and the analysis revealed a low level of LD in the Egyptian faba bean population. A slow LD decay at the genomic and chromosomal levels was observed. Interestingly, the distribution of haplotype blocks was presented in each chromosome and the number of haplotype block ranged from 65 (Chr. 4) to 156 (Chr. 1). Migration and genetic drift are the main reasons for the low LD in the Egyptian faba bean population. The results of this study shed light on the possibility of the genetic improvement of faba bean crop in Egypt and conducting genetic association analyses to identify candidate genes associated with target traits (e.g. protein content, grain yield, etc.) in this panel.

UNASSIGNED: The comorbidity rate of inflammatory bowel disease (IBD) and rheumatoid arthritis (RA) is high; nevertheless, the reasons behind this high rate remain unclear. Their similar genetic makeup probably contributes to this comorbidity.
UNASSIGNED: Based on data obtained from the genome-wide association study of IBD and RA, we first assessed an overall genetic association by performing the linkage disequilibrium score regression (LDSC) analysis. Further, a local correlation analysis was performed by estimating the heritability in summary statistics. Next, the causality between the two diseases was analyzed by two-sample Mendelian randomization (MR). A genetic overlap was analyzed by the conditional/conjoint false discovery rate (cond/conjFDR) method.LDSC with specific expression of gene analysis was performed to identify related tissues between the two diseases. Finally, GWAS multi-trait analysis (MTAG) was also carried out.
UNASSIGNED: IBD and RA are correlated at the genomic level, both overall and locally. The MR results suggested that IBD induced RA. We identified 20 shared loci between IBD and RA on the basis of a conjFDR of <0.01. Additionally, we identified two tissues, namely spleen and small intestine terminal ileum, which were commonly associated with both IBD and RA.
UNASSIGNED: Herein, we proved the presence of a polygenic overlap between the genetic makeup of IBD and RA and provided new insights into the genetic architecture and mechanisms underlying the high comorbidity between these two diseases.

BACKGROUND: Genome sequence variants affecting complex traits (quantitative trait loci, QTL) are enriched in functional regions of the genome, such as those marked by certain histone modifications. These variants are believed to influence gene expression. However, due to the linkage disequilibrium among nearby variants, pinpointing the precise location of QTL is challenging. We aimed to identify allele-specific binding (ASB) QTL (asbQTL) that cause variation in the level of histone modification, as measured by the height of peaks assayed by ChIP-seq (chromatin immunoprecipitation sequencing). We identified DNA sequences that predict the difference between alleles in ChIP-seq peak height in H3K4me3 and H3K27ac histone modifications in the mammary glands of cows.
RESULTS: We used a gapped k-mer support vector machine, a novel best linear unbiased prediction model, and a multiple linear regression model that combines the other two approaches to predict variant impacts on peak height. For each method, a subset of 1000 sites with the highest magnitude of predicted ASB was considered as candidate asbQTL. The accuracy of this prediction was measured by the proportion where the predicted direction matched the observed direction. Prediction accuracy ranged between 0.59 and 0.74, suggesting that these 1000 sites are enriched for asbQTL. Using independent data, we investigated functional enrichment in the candidate asbQTL set and three control groups, including non-causal ASB sites, non-ASB variants under a peak, and SNPs (single nucleotide polymorphisms) not under a peak. For H3K4me3, a higher proportion of the candidate asbQTL were confirmed as ASB when compared to the non-causal ASB sites (P < 0.01). However, these candidate asbQTL did not enrich for the other annotations, including expression QTL (eQTL), allele-specific expression QTL (aseQTL) and sites conserved across mammals (P > 0.05).
CONCLUSIONS: We identified putatively causal sites for asbQTL using the DNA sequence surrounding these sites. Our results suggest that many sites influencing histone modifications may not directly affect gene expression. However, it is important to acknowledge that distinguishing between putative causal ASB sites and other non-causal ASB sites in high linkage disequilibrium with the causal sites regarding their impact on gene expression may be challenging due to limitations in statistical power.

MICA polymorphisms have been associated with increased incidence of acute GvHD and adverse outcome in allogeneic haematopoietic stem cell transplantation (HSCT). MICB is another expressed member of MHC class I-related chain genes and its impact on HSCT outcome is yet to be fully defined. We typed a large cohort of patients and donors for MICB polymorphisms and investigated the impact of MICB matching on outcome after unrelated HSCT. 69.2% of the patients were 10/10 human leukocyte antigen (HLA) matched and 30.8% were 9/10 HLA matched. MICB typing was performed using a short amplicon-based NGS typing assay on the Illumina MiSeq platform. Differences in proteins were considered as mismatches. MICA polymorphisms were identified as possible confounder and were therefore included as parameter in the multivariate analyses. Due to the strong linkage disequilibrium with the classical HLA-genes, sub-stratification for HLA matching status was necessary, and no effect of MICB mismatches was seen in the 10/10 HLA matched group when compared to the MICB matched cases. However, in the 9/10 HLA matched group, MICB mismatched cases showed significantly worse disease free survival (DFS), GvHD and relapse free survival (GRFS) compared to the MICB matched cases (DFS: HR 1.24, p = 0.011; GRFS: HR 1.26, p = 0.002). MICA mismatches had no impact on any outcome parameter. According to our findings, effects previously attributed to MICA differences may have been confounded by MICB polymorphisms. We show that MICB differences contribute a small but relevant effect in 9/10 HLA-matched transplantations, which in turn highlights the possible usefulness of MICB typing in donor selection among similarly suitable 9/10 matched donors, especially when HLA-B mismatches have to be accepted.

Pepper agronomic traits serve as pivotal indicators for characterizing germplasm attributes and correlations. It is important to study differential genotypic variation through phenotypic differences of target traits. Whole genome resequencing was used to sequence the whole genome among different individuals of species with known reference genomes and annotations, and based on this, differential analyses of individuals or populations were carried out to identify SNPs for agronomic traits related to pepper. This study conducted a genome-wide association study encompassing 26 key agronomic traits in 182 upward-growing fruits of C. frutescens and C. annuum. The population structure (phylogenetics, population structure, population principal component analysis, genetic relationship) and linkage disequilibrium analysis were realized to ensure the accuracy and reliability of GWAS results, and the optimal statistical model was determined. A total of 929 SNPs significantly associated with 26 agronomic traits, were identified, alongside the detection of 519 candidate genes within 100 kb region adjacent to these SNPs. Additionally, through gene annotation and expression pattern scrutiny, genes such as GAUT1, COP10, and DDB1 correlated with fruit traits in Capsicum frutescens and Capsicum annuum were validated via qRT-PCR. In the CH20 (Capsicum annuum) and YB-4 (Capsicum frutescens) cultivars, GAUT1 and COP10 were cloned with cDNA lengths of 1065 bp and 561 bp, respectively, exhibiting only a small number of single nucleotide variations and nucleotide deletions. This validation provides a robust reference for molecular marker-assisted breeding of pepper agronomic traits, offering both genetic resources and theoretical foundations for future endeavors in molecular marker-assisted breeding for pepper.

The growth hormone (GH) gene plays a vital role in regulating animal metabolism and body size, making it a potential candidate for influencing livestock performance. This study aimed to investigate the polymorphisms within the GH gene and their associations with 10 biometric traits in the Sumbawa cattle population of Indonesia. Biometric trait data and blood samples were collected from 112 Sumbawa cattle individuals, and their GH gene sequences were analyzed using two sets of primers for amplification. Seven single nucleotide polymorphisms (SNPs) were identified in the GH gene: g.442C>T, g.446G>C, g.558C>T, g.649C>A, g.1492C>A, g.1510C>A, and g.1578G>A. All SNPs were located in the intronic region except for SNP g.558C>T, which was found in the coding sequence (CDS) region. The SNP g.558C>T is classified as a synonymous variant. Haplotype analysis revealed a strong linkage disequilibrium between SNPs g.558C>T and g.649C>A. Distributions of genotypes and alleles of all SNPs were in agreement with the Hardy-Weinberg equilibrium (p > 0.05, χ2 < 15.56), except for SNPs g.446G>C and g.1492C>A. The association study showed that the SNP g.442C>T significantly (p < 0.05) affected HL, BL, SH, and PH traits in Sumbawa cattle. Additionally, the g.446G>C and g.558C>T were also found to be associated with PH and CC traits, respectively. The polymorphisms detected in the GH gene could have implications for selection programs to enhance desired biometric traits in Sumbawa cattle. Improving livestock productivity can be done by understanding genetic diversity and its relationship with phenotypic characteristics.

UNASSIGNED: Gender incongruence (GI) is characterized by a marked incongruence between an individual\'s experienced/expressed gender and the assigned sex at birth. It includes strong displeasure about his or her sexual anatomy and secondary sex characteristics. In some people, this condition produces a strong distress with anxiety and depression named gender dysphoria (GD). This condition appears to be associated with genetic, epigenetics, hormonal as well as social factors. Given that L-glutamate is the major excitatory neurotransmitter in the central nervous system, also associated with male sexual behavior as well as depression, we aimed to determine whether metabotropic glutamate receptors are involved in GD.
UNASSIGNED: We analyzed 74 single nucleotide polymorphisms located at the metabotropic glutamate receptors (mGluR1, mGluR3, mGluR4, mGluR5, mGluR7 and mGluR8) in 94 transgender versus 94 cisgender people. The allele and genotype frequencies were analyzed by c2 test contrasting male and female cisgender and transgender populations. The strength of the associations was measured by binary logistic regression, estimating the odds ratio (OR) for each genotype. Measurement of linkage disequilibrium, and subsequent measurement of haplotype frequencies were also performed considering three levels of significance: P ≤ 0.05, P ≤ 0.005 and P ≤ 0.0005. Furthermore, false positives were controlled with the Bonferroni correction (P ≤ 0.05/74 = 0.00067).
UNASSIGNED: After analysis of allele and genotypic frequencies, we found twenty-five polymorphisms with significant differences at level P ≤ 0.05, five at P ≤ 0.005 and two at P ≤ 0.0005. Furthermore, the only two polymorphisms (rs9838094 and rs1818033) that passed the Bonferroni correction were both related to the metabotropic glutamate receptor 7 (mGluR7) and showed significant differences for multiple patterns of inheritance. Moreover, the haplotype T/G [OR=0.34 (0.19-0.62); P<0.0004] had a lower representation in the transgender population than in the cisgender population, with no evidence of sex cross-interaction.
UNASSIGNED: We provide genetic evidence that the mGluR7, and therefore glutamatergic neurotransmission, may be involved in GI and GD.

UNASSIGNED: Epidemiologic studies have suggested co-morbidity between hypothyroidism and psychiatric disorders. However, the shared genetic etiology and causal relationship between them remain currently unclear.
UNASSIGNED: We assessed the genetic correlations between hypothyroidism and psychiatric disorders [anxiety disorders (ANX), schizophrenia (SCZ), major depressive disorder (MDD), and bipolar disorder (BIP)] using summary association statistics from genome-wide association studies (GWAS). Two disease-associated pleiotropic risk loci and genes were identified, and pathway enrichment, tissue enrichment, and other analyses were performed to determine their specific functions. Furthermore, we explored the causal relationship between them through Mendelian randomization (MR) analysis.
UNASSIGNED: We found significant genetic correlations between hypothyroidism with ANX, SCZ, and MDD, both in the Linkage disequilibrium score regression (LDSC) approach and the high-definition likelihood (HDL) approach. Meanwhile, the strongest correlation was observed between hypothyroidism and MDD (LDSC: rg=0.264, P=7.35×10-12; HDL: rg=0.304, P=4.14×10-17). We also determined a significant genetic correlation between MDD with free thyroxine (FT4) and thyroid-stimulating hormone (TSH) levels. A total of 30 pleiotropic risk loci were identified between hypothyroidism and psychiatric disorders, of which the 15q14 locus was identified in both ANX and SCZ (P values are 6.59×10-11 and 2.10×10-12, respectively) and the 6p22.1 locus was identified in both MDD and SCZ (P values are 1.05×10-8 and 5.75×10-14, respectively). Sixteen pleiotropic risk loci were identified between MDD and indicators of thyroid function, of which, four loci associated with MDD (1p32.3, 6p22.1, 10q21.1, 11q13.4) were identified in both FT4 normal level and Hypothyroidism. Further, 79 pleiotropic genes were identified using Magma gene analysis (P<0.05/18776 = 2.66×10-6). Tissue-specific enrichment analysis revealed that these genes were highly enriched into six brain-related tissues. The pathway analysis mainly involved nucleosome assembly and lipoprotein particles. Finally, our two-sample MR analysis showed a significant causal effect of MDD on the increased risk of hypothyroidism, and BIP may reduce TSH normal levels.
UNASSIGNED: Our findings not only provided evidence of a shared genetic etiology between hypothyroidism and psychiatric disorders, but also provided insights into the causal relationships and biological mechanisms that underlie their relationship. These findings contribute to a better understanding of the pleiotropy between hypothyroidism and psychiatric disorders, while having important implications for intervention and treatment goals for these disorders.

BACKGROUND: The genotype-by-environment interaction (GxE) in beef cattle can be investigated using reaction norm models to assess environmental sensitivity and, combined with genome-wide association studies (GWAS), to map genomic regions related to animal adaptation. Including genetic markers from whole-genome sequencing in reaction norm (RN) models allows us to identify high-resolution candidate genes across environmental gradients through GWAS. Hence, we performed a GWAS via the RN approach using whole-genome sequencing data, focusing on mapping candidate genes associated with the expression of reproductive and growth traits in Nellore cattle. For this purpose, we used phenotypic data for age at first calving (AFC), scrotal circumference (SC), post-weaning weight gain (PWG), and yearling weight (YW). A total of 20,000 males and 7,159 females genotyped with 770k were imputed to the whole sequence (29 M). After quality control and linkage disequilibrium (LD) pruning, there remained ∼ 2.41 M SNPs for SC, PWG, and YW and ∼ 5.06 M SNPs for AFC.
RESULTS: Significant SNPs were identified on Bos taurus autosomes (BTA) 10, 11, 14, 18, 19, 20, 21, 24, 25 and 27 for AFC and on BTA 4, 5 and 8 for SC. For growth traits, significant SNP markers were identified on BTA 3, 5 and 20 for YW and PWG. A total of 56 positional candidate genes were identified for AFC, 9 for SC, 3 for PWG, and 24 for YW. The significant SNPs detected for the reaction norm coefficients in Nellore cattle were found to be associated with growth, adaptative, and reproductive traits. These candidate genes are involved in biological mechanisms related to lipid metabolism, immune response, mitogen-activated protein kinase (MAPK) signaling pathway, and energy and phosphate metabolism.
CONCLUSIONS: GWAS results highlighted differences in the physiological processes linked to lipid metabolism, immune response, MAPK signaling pathway, and energy and phosphate metabolism, providing insights into how different environmental conditions interact with specific genes affecting animal adaptation, productivity, and reproductive performance. The shared genomic regions between the intercept and slope are directly implicated in the regulation of growth and reproductive traits in Nellore cattle raised under different environmental conditions.