Evaluation of heritability partitioning approaches in livestock populations.

BACKGROUND: Heritability partitioning approaches estimate the contribution of different functional classes, such as coding or regulatory variants, to the genetic variance. This information allows a better understanding of the genetic architecture of complex traits, including complex diseases, but can also help improve the accuracy of genomic selection in livestock species. However, methods have mainly been tested on human genomic data, whereas livestock populations have specific characteristics, such as high levels of relatedness, small effective population size or long-range levels of linkage disequilibrium. RESULTS: Here, we used data from 14,762 cows, imputed at the whole-genome sequence level for 11,537,240 variants, to simulate traits in a typical livestock population and evaluate the accuracy of two state-of-the-art heritability partitioning methods, GREML and a Bayesian mixture model. In simulations where a single functional class had increased contribution to heritability, we observed that the estimators were unbiased but had low precision. When causal variants were enriched in variants with low (< 0.05) or high (> 0.20) minor allele frequency or low (below 1st quartile) or high (above 3rd quartile) linkage disequilibrium scores, it was necessary to partition the genetic variance into multiple classes defined on the basis of allele frequencies or LD scores to obtain unbiased results. When multiple functional classes had variable contributions to heritability, estimators showed higher levels of variation and confounding between certain categories was observed. In addition, estimators from small categories were particularly imprecise. However, the estimates and their ranking were still informative about the contribution of the classes. We also demonstrated that using methods that estimate the contribution of a single category at a time, a commonly used approach, results in an overestimation. Finally, we applied the methods to phenotypes for muscular development and height and estimated that, on average, variants in open chromatin regions had a higher contribution to the genetic variance (> 45%), while variants in coding regions had the strongest individual effects (> 25-fold enrichment on average). Conversely, variants in intergenic or intronic regions showed lower levels of enrichment (0.2 and 0.6-fold on average, respectively). CONCLUSIONS: Heritability partitioning approaches should be used cautiously in livestock populations, in particular for small categories. Two-component approaches that fit only one functional category at a time lead to biased estimators and should not be used.

Effective population size in field pea.

BACKGROUND: Effective population size (Ne) is a pivotal parameter in population genetics as it can provide information on the rate of inbreeding and the contemporary status of genetic diversity in breeding populations. The population with smaller Ne can lead to faster inbreeding, with little potential for genetic gain making selections ineffective. The importance of Ne has become increasingly recognized in plant breeding, which can help breeders monitor and enhance the genetic variability or redesign their selection protocols. Here, we present the first Ne estimates based on linkage disequilibrium (LD) in the pea genome. RESULTS: We calculated and compared Ne using SNP markers from North Dakota State University (NDSU) modern breeding lines and United States Department of Agriculture (USDA) diversity panel. The extent of LD was highly variable not only between populations but also among different regions and chromosomes of the genome. Overall, NDSU had a higher and longer-range LD than the USDA that could extend up to 500 Kb, with a genome-wide average r2 of 0.57 (vs 0.34), likely due to its lower recombination rates and the selection background. The estimated Ne for the USDA was nearly three-fold higher (Ne = 174) than NDSU (Ne = 64), which can be confounded by a high degree of population structure due to the selfing nature of pea. CONCLUSIONS: Our results provided insights into the genetic diversity of the germplasm studied, which can guide plant breeders to actively monitor Ne in successive cycles of breeding to sustain viability of the breeding efforts in the long term.

Association between OX40L polymorphism and type 2 diabetes mellitus in Iranians.

INTRODUCTION: Diabetes mellitus (DM) is one of the leading causes of morbidity and mortality worldwide. It is a multifactorial disease that genetic and environmental factors contribute to its development. The aim of the study was to investigate the association of OX40L promoter gene polymorphisms with type 2 diabetes mellitus (T2DM) in Iranians. MATERIALS AND METHODS: Three hundred and sixty-eight subjects including 184 healthy subjects and 184 T2DM patients were enrolled in our study. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) was applied to detect genotype and allele frequencies of rs3850641, rs1234313 and rs10912580. In addition, SNPStats web tool was applied to estimate haplotype frequency and linkage disequilibrium (LD). RESULTS: The distribution of tested polymorphisms was statistically different between the T2DM patients and healthy subjects (P < 0.01). rs1234313 AG (OR = 0.375, 95% CI = 0.193-0.727, P = 0.004) and rs10912580 AG (OR = 0.351, 95% CI = 0.162-0.758, P = 0.008) genotypes were associated with the decreased risk of T2DM in Iranians. Moreover, our prediction revealed that AAG (OR = 0.46, 95% CI= (0.28-0.76), P = 0.0028) and GAG (OR = 0.24, 95% CI= (0.13-0.45), P < 0.0001) haplotypes were related to the reduced risk of the disease. However, the tested polymorphisms had no effect on biochemical parameters and body mass index (BMI) in the patient group (P > 0.05). CONCLUSION: Our findings revealed that OX40L promoter gene polymorphisms are associated with T2DM. Moreover, genotype and allelic variations were related to the decreased risk of T2DM in Iranians. Further studies are recommended to show whether these polymorphic variations could affect OX40/OX40L interaction or OX40L phenotype.

CoPheScan: phenome-wide association studies accounting for linkage disequilibrium.

Phenome-wide association studies (PheWAS) facilitate the discovery of associations between a single genetic variant with multiple phenotypes. For variants which impact a specific protein, this can help identify additional therapeutic indications or on-target side effects of intervening on that protein. However, PheWAS is restricted by an inability to distinguish confounding due to linkage disequilibrium (LD) from true pleiotropy. Here we describe CoPheScan (Coloc adapted Phenome-wide Scan), a Bayesian approach that enables an intuitive and systematic exploration of causal associations while simultaneously addressing LD confounding. We demonstrate its performance through simulation, showing considerably better control of false positive rates than a conventional approach not accounting for LD. We used CoPheScan to perform PheWAS of protein-truncating variants and fine-mapped variants from disease and pQTL studies, in 2275 disease phenotypes from the UK Biobank. Our results identify the complexity of known pleiotropic genes such as APOE, and suggest a new causal role for TGM3 in skin cancer.

Genotyping-by-sequencing targets genic regions and improves resolution of genome-wide association studies in autotetraploid potato.

KEY MESSAGE: De novo genotyping in potato using methylation-sensitive GBS discovers SNPs largely confined to genic or gene-associated regions and displays enhanced effectiveness in estimating LD decay rates, population structure and detecting GWAS associations over 'fixed' SNP genotyping platform. Study also reports the genetic architectures including robust sequence-tagged marker-trait associations for sixteen important potato traits potentially carrying higher transferability across a wider range of germplasm. This study deploys recent advancements in polyploid analytical approaches to perform complex trait analyses in cultivated tetraploid potato. The study employs a 'fixed' SNP Infinium array platform and a 'flexible and open' genome complexity reduction-based sequencing method (GBS, genotyping-by-sequencing) to perform genome-wide association studies (GWAS) for several key potato traits including the assessment of population structure and linkage disequilibrium (LD) in the studied population. GBS SNPs discovered here were largely confined (~ 90%) to genic or gene-associated regions of the genome demonstrating the utility of using a methylation-sensitive restriction enzyme (PstI) for library construction. As compared to Infinium array SNPs, GBS SNPs displayed enhanced effectiveness in estimating LD decay rates and discriminating population subgroups. GWAS using a combined set of 30,363 SNPs identified 189 unique QTL marker-trait associations (QTL-MTAs) covering all studied traits. The majority of the QTL-MTAs were from GBS SNPs potentially illustrating the effectiveness of marker-dense de novo genotyping platforms in overcoming ascertainment bias and providing a more accurate correction for different levels of relatedness in GWAS models. GWAS also detected QTL 'hotspots' for several traits at previously known as well as newly identified genomic locations. Due to the current study exploiting genome-wide genotyping and de novo SNP discovery simultaneously on a large tetraploid panel representing a greater diversity of the cultivated potato gene pool, the reported sequence-tagged MTAs are likely to have higher transferability across a wider range of potato germplasm and increased utility for expediting genomics-assisted breeding for the several complex traits studied.

Genome-wide association study provided insights into the polled phenotype and polled intersex syndrome (PIS) in goats.

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. CONCLUSION: 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.

Genetic causality and metabolite pathway identifying the relationship of blood metabolites and psoriasis.

BACKGROUND: Psoriasis is a chronic inflammatory disease that causes significant disability. However, little is known about the underlying metabolic mechanisms of psoriasis. Our study aims to investigate the causality of 975 blood metabolites with the risk of psoriasis. MATERIALS AND METHODS: We mainly applied genetic analysis to explore the possible associations between 975 blood metabolites and psoriasis. The inverse variance weighted (IVW) method was used as the primary analysis to assess the possible association of blood metabolites with psoriasis. Moreover, generalized summary-data-based Mendelian randomization (GSMR) was used as a supplementary analysis. In addition, linkage disequilibrium score regression (LDSC) was used to investigate their genetic correction further. Metabolic pathway analysis of the most suggested metabolites was also performed using MetaboAnalyst 5.0. RESULTS: In our primary analysis, 17 metabolites, including unsaturated fatty acids, phospholipids, and triglycerides traits, were selected as potential factors in psoriasis, with odd ratios (OR) ranging from 0.986 to 1.01. The GSMR method confirmed the above results (ß = 0.001, p < 0.05). LDSC analysis mainly suggested the genetic correlation of psoriasis with genetic correlations (rg) from 0.088 to 0.155. Based on the selected metabolites, metabolic pathway analysis suggested seven metabolic pathways including ketone body that may be prominent pathways for metabolites in psoriasis. CONCLUSION: Our study supports the causal role of unsaturated fatty acid properties and lipid traits with psoriasis. These properties may be regulated by the ketone body metabolic pathway.

Drug resistance markers in Plasmodium vivax isolates from a Kanchanaburi province, Thailand between January to May 2023.

BACKGROUND: Plasmodium vivax has become the predominant species in the border regions of Thailand. The emergence and spread of antimalarial drug resistance in P. vivax is one of the significant challenges for malaria control. Continuous surveillance of drug resistance is therefore necessary for monitoring the development of drug resistance in the region. This study aims to investigate the prevalence of the mutation in the P. vivax multidrug resistant 1 (Pvmdr1), dihydrofolate reductase (Pvdhfr), and dihydropteroate synthetase (Pvdhps) genes conferred resistance to chloroquine (CQ), pyrimethamine (P) and sulfadoxine (S), respectively. METHOD: 100 P. vivax isolates were obtained between January to May 2023 from a Kanchanaburi province, western Thailand. Nucleotide sequences of Pvmdr1, Pvdhfr, and Pvdhps genes were amplified and sequenced. The frequency of single nucleotide polymorphisms (SNPs)-haplotypes of drug-resistant alleles was assessed. The linkage disequilibrium (LD) tests were also analyzed. RESULTS: In Pvmdr1, T958M, Y976F, and F1076L, mutations were detected in 100%, 21%, and 23% of the isolates, respectively. In Pvdhfr, the quadruple mutant allele (I57R58M61T117) prevailed in 84% of the samples, followed by (L57R58M61T117) in 11%. For Pvdhps, the double mutant allele (G383G553) was detected (48%), followed by the triple mutant allele (G383M512G553) (47%) of the isolates. The most prevalent combination of Pvdhfr (I57R58M61T117) and Pvdhps (G383G553) alleles was sextuple mutated haplotypes (48%). For LD analysis, the association in the SNPs pairs was found between the intragenic and intergenic regions of the Pvdhfr and Pvdhps genes. CONCLUSION: The study has recently updated the high prevalence of three gene mutations associated with CQ and SP resistance. Genetic monitoring is therefore important to intensify in the regions to further assess the spread of drug resistant. Our data also provide evidence on the distribution of drug resistance for the early warning system, thereby threatening P. vivax malaria treatment policy decisions at the national level.

GWAS elucidated grain yield genetics in Indian spring wheat under diverse water conditions.

KEY MESSAGE: Underpinned natural variations and key genes associated with yield under different water regimes, and identified genomic signatures of genetic gain in the Indian wheat breeding program. A novel KASP marker for TKW under water stress was developed and validated. A comprehensive genome-wide association study was conducted on 300 spring wheat genotypes to elucidate the natural variations associated with grain yield and its eleven contributing traits under fully irrigated, restricted water, and simulated no water conditions. Utilizing the 35K Wheat Breeders' Array, we identified 1155 quantitative trait nucleotides (QTNs), with 207 QTNs exhibiting stability across diverse conditions. These QTNs were further delimited into 539 genomic regions using a genome-wide LD value of 3.0 Mbp, revealing pleiotropic control across traits and conditions. Sub-genome A was significantly associated with traits under irrigated conditions, while sub-genome B showed more QTNs under water stressed conditions. Favourable alleles with significantly associated QTNs were delineated, with a notable pyramiding effect for enhancing trait performance. Additionally, allele of only 921 QTNs significantly affected the population mean. Allele profiling highlighted C-306 as a most potential source of drought tolerance. Moreover, 762 genes overlapping significant QTNs were identified, narrowing down to 27 putative candidate genes overlapping 29 novel and functional SNPs expressing (≥ 0.5 tpm) relevance across various growth conditions. A new KASP assay was developed, targeting a gene TraesCS2A03G1123700 regulating thousand kernel weight under severe drought condition. Genomic selection models (GBLUP, BayesB, MxE, and R-Norm) demonstrated an average prediction accuracy of 0.06-0.58 across environments, indicating potential for trait selection. Retrospective analysis of the Indian wheat breeding program supported a genetic gain in GY at the rate of ca. 0.56% per breeding cycle, since 1960, supporting the identification of genomic signatures driving trait selection and genetic gain. These findings offer insight into improving the rate of genetic gain in wheat breeding programs globally.

Cold tolerance SNPs and candidate gene mining in the soybean germination stage based on genome-wide association analysis.

KEY MESSAGE: Three QTLs associated with low-temperature tolerance were identified by genome-wide association analysis, and 15 candidate genes were identified by haplotype analysis and gene expression analyses. Low temperature is a critical factor affecting the geographical distribution, growth, development, and yield of soybeans, with cold stress during seed germination leading to substantial productivity loss. In this study, an association panel comprising 260 soybean accessions was evaluated for four germination traits and four cold tolerance index traits, revealing extensive variation in cold tolerance. Genome-wide association study (GWAS) identified 10 quantitative trait nucleotides (QTNs) associated with cold tolerance, utilizing 30,799 single nucleotide polymorphisms (SNPs) and four GWAS models. Linkage disequilibrium (LD) analysis positioned these QTNs within three cold-tolerance quantitative trait loci (QTL) and, with QTL19-1, was positioned by three multi-locus models, underscoring its importance as a key QTL. Integrative haplotype analysis, supplemented by transcriptome analysis, uncovered 15 candidate genes. The haplotypes within the genes Glyma.18G044200, Glyma.18G044300, Glyma.18G044900, Glyma.18G045100, Glyma.19G222500, and Glyma.19G222600 exhibited significant phenotypic variations, with differential expression in materials with varying cold tolerance. The QTNs and candidate genes identified in this study offer substantial potential for marker-assisted selection and gene editing in breeding cold-tolerant soybeans, providing valuable insights into the genetic mechanisms underlying cold tolerance during soybean germination.

Common and novel haplotype structures between different types of cancer.

BACKGROUND: Background: Genome-wide association studies (GWAS) have identified hundreds of genetic variants associated with cancer risk. GWAS data are important for cancer prevention and understanding the underlying mechanisms of cancer. AIMS: This study aimed to investigate the genetic association between different types of cancer using GWAS data and a bioinformatics approach. METHODS AND RESULTS: The significant GWAS variants associated with more than one cancer type were identified. Common linkage disequilibrium (LD) variants between different types of cancer were identified by 1000 genomes phase 3 LD data. Haplotype blocks were identified by analyzing 1000 Genomes phase 3 genotyping data in the GWAS populations. Subsequent analyses included functional SNP analyses and TCGA gene expression. The results associated with significant GWAS variants (P<5E-8) showed the following haplotype associations in European population: GT rs4808075-rs8170 haplotype on BABAM1 with breast and ovarian cancers, GC rs16857609-rs11693806 haplotype on DIRC3 with breast and thyroid cancers, GCG rs380286-rs401681-rs31487 haplotype on CLPTM1L with skin and lung cancers, GGG rs4430796-rs11651052-rs11263763 haplotype on HNF1B with prostate and endometrial cancers, and GT rs10505477-rs6983267 haplotype on CASC8 associated with colorectal and prostate cancers. All these genes had significantly different expressions in tumor tissues (P<1E-3). In addition, the rs11693806 variant is located in the hsa-miR-873-5p binding site and has an enhancing effect on the hsa-miR-873-5p:DIRC3 interaction. CONCLUSION: These novel haplotype structures and miRNA:lncRNA interactions are important for understanding the common genetic link between cancers. These results can potentially be used in genetic panels.

LDER-GE estimates phenotypic variance component of gene-environment interactions in human complex traits accurately with GE interaction summary statistics and full LD information.

Gene-environment (GE) interactions are essential in understanding human complex traits. Identifying these interactions is necessary for deciphering the biological basis of such traits. In this study, we review state-of-art methods for estimating the proportion of phenotypic variance explained by genome-wide GE interactions and introduce a novel statistical method Linkage-Disequilibrium Eigenvalue Regression for Gene-Environment interactions (LDER-GE). LDER-GE improves the accuracy of estimating the phenotypic variance component explained by genome-wide GE interactions using large-scale biobank association summary statistics. LDER-GE leverages the complete Linkage Disequilibrium (LD) matrix, as opposed to only the diagonal squared LD matrix utilized by LDSC (Linkage Disequilibrium Score)-based methods. Our extensive simulation studies demonstrate that LDER-GE performs better than LDSC-based approaches by enhancing statistical efficiency by ~23%. This improvement is equivalent to a sample size increase of around 51%. Additionally, LDER-GE effectively controls type-I error rate and produces unbiased results. We conducted an analysis using UK Biobank data, comprising 307 259 unrelated European-Ancestry subjects and 966 766 variants, across 217 environmental covariate-phenotype (E-Y) pairs. LDER-GE identified 34 significant E-Y pairs while LDSC-based method only identified 23 significant E-Y pairs with 22 overlapped with LDER-GE. Furthermore, we employed LDER-GE to estimate the aggregated variance component attributed to multiple GE interactions, leading to an increase in the explained phenotypic variance with GE interactions compared to considering main genetic effects only. Our results suggest the importance of impacts of GE interactions on human complex traits.

Micronutrient-Associated Single Nucleotide Polymorphism and Mental Health: A Mendelian Randomization Study.

PURPOSE: Previous studies have demonstrated the link between micronutrients and mental health. However, it remains uncertain whether this connection is causal. We aim to investigate the potential causal effects of micronutrients on mental health based on linkage disequilibrium score (LDSC) regression and Mendelian randomization (MR) analysis. METHODS: Utilizing publicly available genome-wide association study (GWAS) summary datasets, we performed LDSC and MR analysis to identify candidate micronutrients with potential causal effects on mental health. Single nucleotide polymorphisms (SNPs) significantly linked with candidate micronutrients with a genome-wide significance level (p < 5 × 10-8) were selected as instrumental variables (IVs). To estimate the causal effect of candidate micronutrients on mental health, we employed inverse variance weighted (IVW) regression. Additionally, two sensitivity analyses, MR-Egger and weighted median, were performed to validate our results. RESULTS: We found evidence supporting significant causal associations between micronutrients and mental health. LDSC detected several candidate micronutrients, including serum iron (genetic correlation = -0.134, p = 0.032) and vitamin C (genetic correlation = -0.335, p < 0.001) for attention-deficit/hyperactivity disorder (ADHD), iron-binding capacity (genetic correlation = 0.210, p = 0.037) for Alzheimer's disease (AD), and vitamin B12 (genetic correlation = -0.178, p = 0.044) for major depressive disorder (MDD). Further MR analysis suggested a potential causal relationship between vitamin B12 and MDD (b = -0.139, p = 0.009). There was no significant heterogeneity or pleiotropy, indicating the validity of the findings. CONCLUSION: In this study, we identified underlying causal relationships between micronutrients and mental health. Notably, more research is necessary to clarify the underlying biological mechanisms by which micronutrients affect mental health.

A robust cis-Mendelian randomization method with application to drug target discovery.

Mendelian randomization (MR) uses genetic variants as instrumental variables (IVs) to investigate causal relationships between traits. Unlike conventional MR, cis-MR focuses on a single genomic region using only cis-SNPs. For example, using cis-pQTLs for a protein as exposure for a disease opens a cost-effective path for drug target discovery. However, few methods effectively handle pleiotropy and linkage disequilibrium (LD) of cis-SNPs. Here, we propose cisMR-cML, a method based on constrained maximum likelihood, robust to IV assumption violations with strong theoretical support. We further clarify the severe but largely neglected consequences of the current practice of modeling marginal, instead of conditional genetic effects, and only using exposure-associated SNPs in cis-MR analysis. Numerical studies demonstrated our method's superiority over other existing methods. In a drug-target analysis for coronary artery disease (CAD), including a proteome-wide application, we identified three potential drug targets, PCSK9, COLEC11 and FGFR1 for CAD.

Worldwide distribution of genetic factors related to severity of COVID-19 infection.

BACKGROUND: Genome-wide association studies of COVID-19 severity have been carried out mostly on European or East Asian populations with small representation of other world regions. Here we explore the worldwide distribution and linkage disequilibrium (LD) patterns of genetic variants previously associated with COVID-19 severity. METHODS: We followed up the results of a large Spanish genome-wide meta-analysis on 26 populations from the 1000 Genomes Project by calculating allele frequencies and LD scores of the nine most significant SNPs. We also used the entire set of summary statistics to compute polygenic risk scores (PRSs) and carried out comparisons at the population and continental level. RESULTS: We observed the strongest differences among continental regions for the five top SNPs in chromosome 3. European, American, and South Asian populations showed similar LD patterns. Average PRSs in South Asian and American populations were consistently higher than those observed in Europeans. While PRS distributions were similar among South Asians, the American populations showed striking differences among them. CONCLUSIONS: Considering the caveats of PRS transferability across ethnicities, our analysis showed that American populations present the highest genetic risk score, hence potentially higher propensity, for COVID-19 severity. Independent validation is warranted with additional summary statistics and phenotype data.

A multi-ethnic reference panel to impute HLA classical and non-classical class I alleles in admixed samples: Testing imputation accuracy in an admixed sample from Brazil.

The MHC class I region contains crucial genes for the innate and adaptive immune response, playing a key role in susceptibility to many autoimmune and infectious diseases. Genome-wide association studies have identified numerous disease-associated SNPs within this region. However, these associations do not fully capture the immune-biological relevance of specific HLA alleles. HLA imputation techniques may leverage available SNP arrays by predicting allele genotypes based on the linkage disequilibrium between SNPs and specific HLA alleles. Successful imputation requires diverse and large reference panels, especially for admixed populations. This study employed a bioinformatics approach to call SNPs and HLA alleles in multi-ethnic samples from the 1000 genomes (1KG) dataset and admixed individuals from Brazil (SABE), utilising 30X whole-genome sequencing data. Using HIBAG, we created three reference panels: 1KG (n = 2504), SABE (n = 1171), and the full model (n = 3675) encompassing all samples. In extensive cross-validation of these reference panels, the multi-ethnic 1KG reference exhibited overall superior performance than the reference with only Brazilian samples. However, the best results were achieved with the full model. Additionally, we expanded the scope of imputation by developing reference panels for non-classical, MICA, MICB and HLA-H genes, previously unavailable for multi-ethnic populations. Validation in an independent Brazilian dataset showcased the superiority of our reference panels over the Michigan Imputation Server, particularly in predicting HLA-B alleles among Brazilians. Our investigations underscored the need to enhance or adapt reference panels to encompass the target population's genetic diversity, emphasising the significance of multiethnic references for accurate imputation across different populations.

Shared genetic links between hypothyroidism and psychiatric disorders: evidence from a comprehensive genetic analysis.

Background: Epidemiologic studies have suggested co-morbidity between hypothyroidism and psychiatric disorders. However, the shared genetic etiology and causal relationship between them remain currently unclear. Methods: 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. Results: 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. Conclusions: 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.

Cost-effective genomic prediction of critical economic traits in sturgeons through low-coverage sequencing.

Low-coverage whole-genome sequencing (LCS) offers a cost-effective alternative for sturgeon breeding, especially given the lack of SNP chips and the high costs associated with whole-genome sequencing. In this study, the efficiency of LCS for genotype imputation and genomic prediction was assessed in 643 sequenced Russian sturgeons (∼13.68×). The results showed that using BaseVar+STITCH at a sequencing depth of 2× with a sample size larger than 300 resulted in the highest genotyping accuracy. In addition, when the sequencing depth reached 0.5× and SNP density was reduced to 50 K through linkage disequilibrium pruning, the prediction accuracy was comparable to that of whole sequencing depth. Furthermore, an incremental feature selection method has the potential to improve prediction accuracy. This study suggests that the combination of LCS and imputation can be a cost-effective strategy, contributing to the genetic improvement of economic traits and promoting genetic gains in aquaculture species.

Molecular genetic diversity and linkage disequilibrium structure of the Egyptian faba bean using Single Primer Enrichment Technology (SPET).

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.

Temporal genomics in Hawaiian crickets reveals compensatory intragenomic coadaptation during adaptive evolution.

Theory predicts that compensatory genetic changes reduce negative indirect effects of selected variants during adaptive evolution, but evidence is scarce. Here, we test this in a wild population of Hawaiian crickets using temporal genomics and a high-quality chromosome-level cricket genome. In this population, a mutation, flatwing, silences males and rapidly spread due to an acoustically-orienting parasitoid. Our sampling spanned a social transition during which flatwing fixed and the population went silent. We find long-range linkage disequilibrium around the putative flatwing locus was maintained over time, and hitchhiking genes had functions related to negative flatwing-associated effects. We develop a combinatorial enrichment approach using transcriptome data to test for compensatory, intragenomic coevolution. Temporal changes in genomic selection were distributed genome-wide and functionally associated with the population's transition to silence, particularly behavioural responses to silent environments. Our results demonstrate how 'adaptation begets adaptation'; changes to the sociogenetic environment accompanying rapid trait evolution can generate selection provoking further, compensatory adaptation.