Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium

Polygenic Indices (a.k.a. Polygenic Scores) in Social Science: A Guide for Interpretation and Evaluation

Polygenic indices (PGI)-the new recommended label for polygenic scores (PGS) in social science-are genetic summary scales often used to represent an individual's liability for a disease, trait, or behavior based on the additive effects of measured genetic variants. Enthusiasm for linking genetic data with social outcomes and the inclusion of premade PGIs in social science datasets have facilitated increased uptake of PGIs in social science research-a trend that will likely continue. Yet, most social scientists lack the expertise to interpret and evaluate PGIs in social science research. Here, we provide a primer on PGIs for social scientists focusing on key concepts, unique statistical genetic considerations, and best practices in calculation, estimation, reporting, and interpretation. We summarize our recommended best practices as a checklist to aid social scientists in evaluating and interpreting studies with PGIs. We conclude by discussing the similarities between PGIs and standard social science scales and unique interpretative considerations.

Malignancy following solid organ transplantation: Current techniques for determination of donor versus recipient origin

Among the post-transplantation complications that patients may encounter, the transmission of a donor-derived malignant neoplasm is uncommon but potentially life threatening. The determination of donor versus recipient origin is essential particularly in the setting of multiple transplant recipients from the donor. Advances in molecular biology now allow accurate discrimination utilizing routine tissue samples in a timely and cost-effective manner. The techniques are routinely performed in hospital molecular biology laboratories and are also available in commercial labs. The current methodologies are discussed and future possibilities are presented for clinicians caring for solid organ recipients.

Unraveling the genomic landscape of Campylorhynchus wrens along western Ecuador's precipitation gradient: Insights into hybridization, isolation by distance, and isolation by the environment

Environmental gradients have the potential to influence genetic differentiation among populations ultimately leading to allopatric speciation. However, environmental gradients can also facilitate hybridization between closely related taxa. We investigated a putative hybrid zone in western Ecuador, involving two polytypic wren species (Aves: Troglodytidae), Campylorhynchus zonatus and C. fasciatus. Our study addressed two primary questions: (1) Is there evidence of population structure and genetic admixture between these taxa in western Ecuador? and (2) What are the relative contributions of isolation by distance and isolation by the environment to the observed genetic differentiation along the environmental gradient in this region? We analyzed 4409 single-nucleotide polymorphisms (SNPs) from 112 blood samples sequenced using ddRadSeq and a de novo assembly. The optimum number of genetic clusters ranged from 2 to 4, aligning with geographic origins, known phylogenetics, and physical or ecological constraints. We observed notable transitions in admixture proportions along the environmental gradient in western Ecuador between C. z. brevirostris and the northern and southern genetic clusters of C. f. pallescens. Genetic differentiation between the two C. f. pallescens populations could be attributed to an unreported potential physical barrier in central western Ecuador, where the proximity of the Andes to the coastline restricts lowland habitats, limiting dispersal and gene flow, especially among dry-habitat specialists. The observed admixture in C. f. pallescens suggests that this subspecies may be a hybrid between C. z. brevirostris and C. fasciatus, with varying degrees of admixture in western Ecuador and northwestern Peru. We found evidence of isolation by distance, while isolation by the environment was less pronounced but still significant for annual mean precipitation and precipitation seasonality. This study enhances our understanding of avian population genomics in tropical regions.

Population-based nanopore sequencing of the HIV-1 pangenome to identify drug resistance mutations

HIV-1 drug resistance genotypic tests have primarily been performed by Sanger sequencing of gene segments encoding different drug target proteins. Since the number of targets has increased with the addition of a new class of antiretroviral drugs, a simple high-throughput system for assessing nucleotide sequences throughout the HIV-1 genome is required. Here, we developed a new solution using nanopore sequencing of viral pangenomes amplified by PCR. Benchmark tests using HIV-1 molecular clones demonstrated an accuracy of up to 99.9%. In addition, validation tests of our protocol in 106 clinical samples demonstrated high concordance of drug resistance and tropism genotypes (92.5% and 98.1%, respectively) between the nanopore sequencing-based results and archived clinical determinations made based on Sanger sequencing data. These results suggest that our new approach will be a powerful solution for the comprehensive survey of HIV-1 drug resistance mutations in clinical settings.

An approach for searching genes in signaling pathways or gene-gene interaction networks related to Hypertension in the Mexican population

INTRODUCTION

The selection of single nucleotide polymorphisms (SNPs) to evaluate the genetic susceptibility in complex traits is often conducted in isolation, without considering the entire set of genes. Incorporating signaling pathways or gene-gene interaction search may provide a more comprehensive approach to selecting SNP candidates for further study.

OBJECTIVE

To propose a systematic procedure for identifying SNPs candidates with complex traits such as hypertension and blood pressure.

METHODS

Sequential stages to SNPs selection: 1) literature review to identify SNPs, following the PRISMA methodology, 2) identification and selection of signaling pathways and selection of gene-gene interaction networks using the STRING software, and 3) application of specific criteria for SNPs candidates, including: a) SNPs with minor allele frequency > 5% in the target population, b) SNPs located within genes involved in three or more signaling pathways, and c) SNPs that are not in linkage disequilibrium, with a D'or r2 value < 0.8.

RESULTS

Stage 1) A total of 44 publications were selected, providing information on 230 genes evaluated with blood pressure. Stage 2) Using the STRING software, we selected 7 signaling pathways with a false discovery rate < 0.0001 and strength ≥ 0.8; and we identified 16 genes belonging to gene-gene interaction networks, six of them share ≥ 3 signaling pathways. Stage 3) Finally, 7 SNPs were selected for genotyping in the Health Workers Cohort Study. We observed a positive association between SNPs with hypertension incidence in males (rs1130214, rs3807989) and females (rs5051, rs2493123).

CONCLUSION

Our methodological proposal may be a reliable way for selecting SNPs candidates to study complex traits.

Identification of eQTLs using different sets of single nucleotide polymorphisms associated with carcass and body composition traits in pigs

BACKGROUND

Mapping expression quantitative trait loci (eQTLs) in skeletal muscle tissue in pigs is crucial for understanding the relationship between genetic variation and phenotypic expression of carcass traits in meat animals. Therefore, the primary objective of this study was to evaluate the impact of different sets of single nucleotide polymorphisms (SNP), including scenarios removing SNPs pruned for linkage disequilibrium (LD) and SNPs derived from SNP chip arrays and RNA-seq data from liver, brain, and skeletal muscle tissues, on the identification of eQTLs in the Longissimus lumborum tissue, associated with carcass and body composition traits in Large White pigs. The SNPs identified from muscle mRNA were combined with SNPs identified in the brain and liver tissue transcriptomes, as well as SNPs from the GGP Porcine 50 K SNP chip array. Cis- and trans-eQTLs were identified based on the skeletal muscle gene expression level, followed by functional genomic analyses and statistical associations with carcass and body composition traits in Large White pigs.

RESULTS

The number of cis- and trans-eQTLs identified across different sets of SNPs (scenarios) ranged from 261 to 2,539 and from 29 to 13,721, respectively. Furthermore, 6,180 genes were modulated by eQTLs in at least one of the scenarios evaluated. The eQTLs identified were not significantly associated with carcass and body composition traits but were significantly enriched for many traits in the "Meat and Carcass" type QTL. The scenarios with the highest number of cis- (n = 304) and trans- (n = 5,993) modulated genes were the unpruned and LD-pruned SNP set scenarios identified from the muscle transcriptome. These genes include 84 transcription factor coding genes.

CONCLUSIONS

After LD pruning, the set of SNPs identified based on the transcriptome of the skeletal muscle tissue of pigs resulted in the highest number of genes modulated by eQTLs. Most eQTLs are of the trans type and are associated with genes influencing complex traits in pigs, such as transcription factors and enhancers. Furthermore, the incorporation of SNPs from other genomic regions to the set of SNPs identified in the porcine skeletal muscle transcriptome contributed to the identification of eQTLs that had not been identified based on the porcine skeletal muscle transcriptome alone.

The mediation effect of glucocorticoid receptor gene (NR3C1) methylation between childhood maltreatment and depressive symptoms in Chinese adolescents: A 2-year longitudinal study

This three-wave longitudinal study examined whether methylation alterations in promoter exon 1F of a stress-related gene-NR3C1 (NR3C1-1F)-explained the longitudinal associations between childhood maltreatment and adolescent depressive symptoms. A total of 370 Han Chinese adolescents (Mage  = 16.31 ± 1.28 years; 51.4% girls) recruited from Shandong, China were tracked from 2018 to 2020. The results showed that the severity of childhood maltreatment, especially that of emotional abuse and physical neglect, conferred risk for adolescent depressive symptoms via reducing NR3C1-1F methylation levels. These mediation effects of NR3C1-1F methylation did not vary between adolescent sex or NR3C1 BclI and Tth111I polymorphisms. The findings highlight how childhood maltreatment contributes to psychopathology development at a biological level.

Influence of Admixture on Phenotypes

Individuals of European descent have historically been the focus of genetic studies and possess relatively homogenous genomes. As a result, analytical methods have been developed and optimized with such genomes in mind. African-descent and Latino individuals generally possess genomes of greater architectural complexity due to mosaic genomic ancestry, which can extensively and intricately impact phenotypic expression. As such, genetic analyses of admixed individuals require that genetic admixture be quantified to accurately model the impact of genetic variation on phenotypic expression. In this overview, we explore how fundamental genetic concepts such as linkage disequilibrium and differential allele frequency interact with genetic admixture to uniquely influence phenotypes in admixed individuals. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC.

Whole chloroplast genome-specific non-synonymous SNPs reveal the presence of substantial diversity in the pigeonpea mini-core collection

To unravel the plastid genome diversity among the cultivated groups of the pigeonpea germplasm, we characterized the SNP occurrence and distribution of 142 pigeonpea mini-core collections based on their reference-based assembly of the chloroplast genome. A total of 8921 SNPs were found, which were again filtered and finally 3871 non-synonymous SNPs were detected and used for diversity estimates. These 3871 SNPs were classified into 12 groups and were present in only 44 of the 125 genes, demonstrating the presence of a precise mechanism for maintaining the whole chloroplast genome throughout evolution. The Acetyl-CoA carboxylase D gene possesses the maximum number of SNPs (12.29%), but the Adenosine Tri-Phosphate synthatase cluster genes (atpA, atpB, atpE, atpF, atpH, and atpI) altogether bear 43.34% of the SNPs making them most diverse. Various diversity estimates, such as the number of effective alleles (1.013), Watterson's estimate (0.19), Tajima's D ( - 3.15), Shannon's information index (0.036), suggest the presence of less diversity in the cultivated gene pool of chloroplast genomes. The genetic relatedness estimates based on pairwise correlations were also in congruence with these diversity descriptors and indicate the prevalence of rare alleles in the accessions. Interestingly, no stratification was observed either through STRUCTURE, PCoA, or phylogenetic analysis, indicating the common origin of the chloroplast in all the accessions used, irrespective of their geographical distribution. Further 6194 Cleaved Amplified Polymorphic Sequences (CAPS) markers for 531 SNPs were developed and validated in a selected set of germplasm. Based on these results, we inferred that all of the cultivated gene pools of pigeonpea have a common origin for the chloroplast genome and they possess less diversity in protein-coding regions, indicating a stable and evolved plastid genome. At the same time, all diversity analysis indicates the occurrence of rare alleles, suggesting the suitability of the mini-core collection in future pigeonpea improvement programs. In addition, the development of chloroplast genome-based CAPS markers would have utility in pigeonpea breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03785-8.

Microsatellites used in forensics are in regions enriched for trait-associated variants

The 20 short tandem repeat (STR) loci of the combined DNA index system (CODIS) are the basis of the vast majority of forensic genetics in the United States. One argument for permissive rules about the collection of CODIS genotypes is that the CODIS loci are thought to contain little information about ancestry or traits. However, in the past 20 years, a growing field has identified hundreds of thousands of genotype-trait associations. Here, we conduct a survey of the landscape of such associations surrounding the CODIS loci as compared with non-CODIS STRs. Although this study cannot establish or quantify associations between CODIS genotypes and phenotypes, we find that the regions around the CODIS loci are enriched for both known pathogenic variants (> 90th percentile) and for trait-associated SNPs identified in genome-wide association studies (GWAS) (≥ 95th percentile in 10kb and 100kb flanking regions), compared with other random sets of autosomal tetranucleotide-repeat STRs.

Haplotype blocks for genomic prediction: a comparative evaluation in multiple crop datasets

In modern plant breeding, genomic selection is becoming the gold standard for selection of superior genotypes. The basis for genomic prediction models is a set of phenotyped lines along with their genotypic profile. With high marker density and linkage disequilibrium (LD) between markers, genotype data in breeding populations tends to exhibit considerable redundancy. Therefore, interest is growing in the use of haplotype blocks to overcome redundancy by summarizing co-inherited features. Moreover, haplotype blocks can help to capture local epistasis caused by interacting loci. Here, we compared genomic prediction methods that either used single SNPs or haplotype blocks with regards to their prediction accuracy for important traits in crop datasets. We used four published datasets from canola, maize, wheat and soybean. Different approaches to construct haplotype blocks were compared, including blocks based on LD, physical distance, number of adjacent markers and the algorithms implemented in the software "Haploview" and "HaploBlocker". The tested prediction methods included Genomic Best Linear Unbiased Prediction (GBLUP), Extended GBLUP to account for additive by additive epistasis (EGBLUP), Bayesian LASSO and Reproducing Kernel Hilbert Space (RKHS) regression. We found improved prediction accuracy in some traits when using haplotype blocks compared to SNP-based predictions, however the magnitude of improvement was very trait- and model-specific. Especially in settings with low marker density, haplotype blocks can improve genomic prediction accuracy. In most cases, physically large haplotype blocks yielded a strong decrease in prediction accuracy. Especially when prediction accuracy varies greatly across different prediction models, prediction based on haplotype blocks can improve prediction accuracy of underperforming models. However, there is no "best" method to build haplotype blocks, since prediction accuracy varied considerably across methods and traits. Hence, criteria used to define haplotype blocks should not be viewed as fixed biological parameters, but rather as hyperparameters that need to be adjusted for every dataset.

Polymorphisms in the Nonhomologous End-joining DNA Repair Pathway are Associated with HPV Integration in Cervical Dysplasia

Previous evidence indicates that human papillomavirus (HPV) integration status may be associated with cervical cancer development and progression. However, host genetic variation within genes that may play important roles in the viral integration process is understudied. The aim of this study was to examine the association between HPV16 and HPV18 viral integration status and SNPs in nonhomologous-end-joining (NHEJ) DNA repair pathway genes on cervical dysplasia. Women enrolled in two large trials of optical technologies for cervical cancer detection and positive for HPV16 or HPV18 were selected for HPV integration analysis and genotyping. Associations between SNPs and cytology (normal, low-grade, or high-grade lesions) were evaluated. Among women with cervical dysplasia, polytomous logistic regression models were used to evaluate the effect of each SNP on viral integration status. Of the 710 women evaluated [149 high-grade squamous intraepithelial lesion (HSIL), 251; low-grade squamous intraepithelial lesion (LSIL, 310 normal)], 395 (55.6%) were positive for HPV16 and 192 (27%) were positive for HPV18. Tag-SNPs in 13 DNA repair genes, including RAD50, WRN, and XRCC4, were significantly associated with cervical dysplasia. HPV16 integration status was differential across cervical cytology, but overall, most participants had a mix of both episomal and integrated HPV16. Four tag-SNPs in the XRCC4 gene were found to be significantly associated with HPV16 integration status. Our findings indicate that host genetic variation in NHEJ DNA repair pathway genes, specifically XRCC4, are significantly associated with HPV integration, and that these genes may play an important role in determining cervical cancer development and progression.

PREVENTION RELEVANCE

HPV integration in premalignant lesions and is thought to be an important driver of carcinogenesis. However, it is unclear what factors promote integration. The use of targeted genotyping among women presenting with cervical dysplasia has the potential to be an effective tool in assessing the likelihood of progression to cancer.

RASA3 is a candidate gene in sickle cell disease-associated pulmonary hypertension and pulmonary arterial hypertension

Pulmonary hypertension (PH) is associated with significant morbidity and mortality. RASA3 is a GTPase activating protein integral to angiogenesis and endothelial barrier function. In this study, we explore the association of RASA3 genetic variation with PH risk in patients with sickle cell disease (SCD)-associated PH and pulmonary arterial hypertension (PAH). Cis-expression quantitative trait loci (eQTL) were queried for RASA3 using whole genome genotype arrays and gene expression profiles derived from peripheral blood mononuclear cells (PBMC) of three SCD cohorts. Genome-wide single nucleotide polymorphisms (SNPs) near or in the RASA3 gene that may associate with lung RASA3 expression were identified, reduced to 9 tagging SNPs for RASA3 and associated with markers of PH. Associations between the top RASA3 SNP and PAH severity were corroborated using data from the PAH Biobank and analyzed based on European or African ancestry (EA, AA). We found that PBMC RASA3 expression was lower in patients with SCD-associated PH as defined by echocardiography and right heart catheterization and was associated with higher mortality. One eQTL for RASA3 (rs9525228) was identified, with the risk allele correlating with PH risk, higher tricuspid regurgitant jet velocity and higher pulmonary vascular resistance in patients with SCD-associated PH. rs9525228 associated with markers of precapillary PH and decreased survival in individuals of EA but not AA. In conclusion, RASA3 is a novel candidate gene in SCD-associated PH and PAH, with RASA3 expression appearing to be protective. Further studies are ongoing to delineate the role of RASA3 in PH.

Predicting late-stage age-related macular degeneration by integrating marginally weak SNPs in GWA studies

Introduction: Age-related macular degeneration (AMD) is a progressive neurodegenerative disease and the leading cause of blindness in developed countries. Current genome-wide association studies (GWAS) for late-stage age-related macular degeneration are mainly single-marker-based approaches, which investigate one Single-Nucleotide Polymorphism (SNP) at a time and postpone the integration of inter-marker Linkage-disequilibrium (LD) information in the downstream fine mappings. Recent studies showed that directly incorporating inter-marker connection/correlation into variants detection can help discover novel marginally weak single-nucleotide polymorphisms, which are often missed in conventional genome-wide association studies, and can also help improve disease prediction accuracy. Methods: Single-marker analysis is performed first to detect marginally strong single-nucleotide polymorphisms. Then the whole-genome linkage-disequilibrium spectrum is explored and used to search for high-linkage-disequilibrium connected single-nucleotide polymorphism clusters for each strong single-nucleotide polymorphism detected. Marginally weak single-nucleotide polymorphisms are selected via a joint linear discriminant model with the detected single-nucleotide polymorphism clusters. Prediction is made based on the selected strong and weak single-nucleotide polymorphisms. Results: Several previously identified late-stage age-related macular degeneration susceptibility genes, for example, BTBD16, C3, CFH, CFHR3, HTARA1, are confirmed. Novel genes DENND1B, PLK5, ARHGAP45, and BAG6 are discovered as marginally weak signals. Overall prediction accuracy of 76.8% and 73.2% was achieved with and without the inclusion of the identified marginally weak signals, respectively. Conclusion: Marginally weak single-nucleotide polymorphisms, detected from integrating inter-marker linkage-disequilibrium information, may have strong predictive effects on age-related macular degeneration. Detecting and integrating such marginally weak signals can help with a better understanding of the underlying disease-development mechanisms for age-related macular degeneration and more accurate prognostics.

Donor and recipient polymorphisms of MAPK signaling pathway genes influence post-transplant liver function in Chinese liver transplant patients taking tacrolimus

Tacrolimus (TAC) is an immunosuppressive drug that is widely used for patients who underwent liver transplantation. In addition to inhibiting the action of calcineurin, TAC also exerts its immunosuppressive effects by interfering with mitogen activated protein kinase (MAPK) pathway. In this study, we investigated the impact of both recipient and donor genetic polymorphisms of MAPK kinase kinase (MAP3K) genes on clinical events in Han Chinese liver transplantation recipients taking TAC. Fifty-seven tag SNPs of 11 genes (MEKK1, MEKK2, MEKK4, MLK1, MLK3, ASK1, TAO1, TAO2, Tpl2, TAK1 and ZAK1) in the MAPK pathway were detected by MALDI-TOF MS assay in 175 TAC-treated liver transplant recipients. The associations of SNPs with incidence of acute rejection, TAC-induced acute nephrotoxicity, and post-transplantation liver and kidney function were explored using Kaplan-Meier survival analysis, Cox-proportional hazard model and linear mixed model, respectively. For the sites significantly associated with clinical events, the dual-luciferase reporter gene system was used to perform preliminary function verification. The results showed that (1) Donor-recipient combinational (D-R) MEKK1 rs62355944 and D-R MLK1 rs8006424 genotypes were significant influence factors of post-transplantation γ-glutamyl transpeptidase (GGT) level (P < 0.0001); (2) D-R MLK1 rs8006424 genotypes were found to significantly affect the alkaline phosphatase (ALP) level after transplantation (P < 0.0001). The results of the dual luciferase reporter gene system demonstrated that the luciferase activity of the pGL3-rs8006424A was significantly higher than that of pGL3-rs8006424G (3.47 ± 0.10 vs 2.97 ± 0.08, P = 0.002). Therefore, MEKK1 rs62355944 and MLK1 rs8006424 might serve as biomarkers to predict post-transplant liver function in liver transplant patients.

Association of CTNND2 gene polymorphism with schizophrenia: Two-sample case-control study in Chinese Han population

OBJECTIVES

Genetic factors play an important role in the etiology of schizophrenia (SZ). Catenin Delta 2 (CTNND2) is one of the genes regulating neuronal development in the brain. It is unclear whether CTNND2 is involved in SZ. With the hypothesis that CTNND2 may be a risk gene for SZ, we performed a case-control association analysis to investigate if CTNND2 gene single nucleotide polymorphisms (SNPs) are implicated in SZ in a Han Chinese northern population.

MATERIALS AND METHODS

We recruited subjects from 2010 to 2022 from the Han population of northern Henan and divided them into two case-control samples, including a discovery sample (SZ = 528 and control = 528) and replication sample (SZ = 2458 and control = 6914). Twenty-one SNPs were genotyped on the Illumina BeadStation 500G platform using GoldenGate technology and analyzed by PLINK. Positive and Negative Syndrome Scale (PANSS) was used to assess clinical symptoms.

RESULTS

Rs16901943, rs7733427, and rs2168878 SNPs were associated with SZ (Chi² = 7.484, 11.576, and 5.391, respectively, df = 1; p = 0.006, 0.00067, and 0.02, respectively) in two samples. Rs10058868 was associated with SZ in male patients in the discovery sample (Chi² = 6.264, df = 1, p = .044). Only rs7733427 survived Bonferroni correction. Linkage disequilibrium block three haplotypes were associated with SZ in the discovery and total sample. PANSS analysis of the four SNPs implicated rs10058868 and rs2168878 with symptoms of depression and excitement, respectively, in the SZ patients.

CONCLUSION

Four SNPs were identified as being correlated with SZ. The CTNND2 gene may be involved in susceptibility to SZ.

Application potential of chicken DNA chip in domestic pigeon species - Preliminary results

Introducing the SNP technology to pigeon breeding will enhance the competitiveness of a sector that produces one of the healthiest and best quality meats. The present study aimed to test the applicability of the Illumina Chicken_50K_CobbCons array on 24 domestic pigeon individuals from the Mirthys hybrids and Racing pigeon breeds. A total of 53,313 SNPs were genotyped. Principal component analysis shows a significant overlap between the two groups. The chip performed poorly in this data set, with a call rate per sample of 0.474 (49%). The low call rate was likely due to an increase in the evolutionary distance. A total of 356 SNPs were retained after a relatively strict quality control. We have demonstrated that it is technically feasible to use a chicken microarray chip on pigeon samples. Presumably, with a larger sample size and by assigning phenotypic data, efficiency would be improved, allowing more thorough analyses, such as genome-wide association studies.

Assessing effectiveness of many-objective evolutionary algorithms for selection of tag SNPs

BACKGROUND

Studies on genome-wide associations help to determine the cause of many genetic diseases. Genome-wide associations typically focus on associations between single-nucleotide polymorphisms (SNPs). Genotyping every SNP in a chromosomal region for identifying genetic variation is computationally very expensive. A representative subset of SNPs, called tag SNPs, can be used to identify genetic variation. Small tag SNPs save the computation time of genotyping platform, however, there could be missing data or genotyping errors in small tag SNPs. This study aims to solve Tag SNPs selection problem using many-objective evolutionary algorithms.

METHODS

Tag SNPs selection can be viewed as an optimization problem with some trade-offs between objectives, e.g. minimizing the number of tag SNPs and maximizing tolerance for missing data. In this study, the tag SNPs selection problem is formulated as a many-objective problem. Nondominated Sorting based Genetic Algorithm (NSGA-III), and Multi-Objective Evolutionary Algorithm based on Decomposition (MOEA/D), which are Many-Objective evolutionary algorithms, have been applied and investigated for optimal tag SNPs selection. This study also investigates different initialization methods like greedy and random initialization. optimization.

RESULTS

The evaluation measures used for comparing results for different algorithms are Hypervolume, Range, SumMin, MinSum, Tolerance rate, and Average Hamming distance. Overall MOEA/D algorithm gives superior results as compared to other algorithms in most cases. NSGA-III outperforms NSGA-II and other compared algorithms on maximum tolerance rate, and SPEA2 outperforms all algorithms on average hamming distance.

CONCLUSION

Experimental results show that the performance of our proposed many-objective algorithms is much superior as compared to the results of existing methods. The outcomes show the advantages of greedy initialization over random initialization using NSGA-III, SPEA2, and MOEA/D to solve the tag SNPs selection as many-objective optimization problem.

A comprehensive evaluation of polygenic score and genotype imputation performances of human SNP arrays in diverse populations

Regardless of the overwhelming use of next-generation sequencing technologies, microarray-based genotyping combined with the imputation of untyped variants remains a cost-effective means to interrogate genetic variations across the human genome. This technology is widely used in genome-wide association studies (GWAS) at bio-bank scales, and more recently, in polygenic score (PGS) analysis to predict and stratify disease risk. Over the last decade, human genotyping arrays have undergone a tremendous growth in both number and content making a comprehensive evaluation of their performances became more important. Here, we performed a comprehensive performance assessment for 23 available human genotyping arrays in 6 ancestry groups using diverse public and in-house datasets. The analyses focus on performance estimation of derived imputation (in terms of accuracy and coverage) and PGS (in terms of concordance to PGS estimated from whole-genome sequencing data) in three different traits and diseases. We found that the arrays with a higher number of SNPs are not necessarily the ones with higher imputation performance, but the arrays that are well-optimized for the targeted population could provide very good imputation performance. In addition, PGS estimated by imputed SNP array data is highly correlated to PGS estimated by whole-genome sequencing data in most cases. When optimal arrays are used, the correlations of PGS between two types of data are higher than 0.97, but interestingly, arrays with high density can result in lower PGS performance. Our results suggest the importance of properly selecting a suitable genotyping array for PGS applications. Finally, we developed a web tool that provides interactive analyses of tag SNP contents and imputation performance based on population and genomic regions of interest. This study would act as a practical guide for researchers to design their genotyping arrays-based studies. The tool is available at: https://genome.vinbigdata.org/tools/saa/ .

Breast cancer-associated SNP rs72755295 is a cis-regulatory variation for human EXO1

Breast cancer is the most common malignant tumor in women. A previous genome-wide association study reports that rs72755295, a SNP locating at intron of EXO1 (exonuclease 1), is associated with breast cancer. Due to the complete linkage disequilibrium between rs72755295 and rs4149909, a nonsynonymous mutation for EXO1, rs4149909 is supposed to be the causal SNP. Since EXO1 is overexpressed in breast carcinoma samples, we hypothesized that the genetic variations in this locus might confer breast cancer risk by regulating EXO1 expression. To substantiate this, a functional genomics study was performed. The dual luciferase assay indicated that G of rs72755295 presents significantly higher relative enhancer activity than A, thus verifying that this SNP can influence gene expression in breast cell. Through chromosome conformation capture it was disclosed that the enhancer containing rs72755295 can interact with the EXO1 promoter. RNA-seq analysis indicated that EXO1 expression is dependent on the rs72755295 genotype. By chromatin immunoprecipitation, the transcription factor PAX6 (paired box 6) was recognized to bind the region spanning rs72755295. In electrophoretic mobility shift assay, G of rs72755295 displays obviously higher binding affinity with nuclear protein than A. Our results indicated that rs72755295 is a cis-regulatory variation for EXO1 and might confer breast cancer risk besides rs4149909.

Identification of rs2736099 as a novel cis-regulatory variation for TERT and implications for tumorigenesis and cell proliferation

PURPOSE

Lung cancer is a malignant tumor with obvious genetic predisposition. Association studies have proposed that rs2853677, a SNP localizing at intron region of TERT (telomerase reverse transcriptase), is significantly associated with TERT expression, telomere length and eventually lung cancer risk. However, functional genomics work indicates that rs2853677 is not with the ability to alter gene expression. All these facts make us hypothesize that some other genetic variation(s) are in linkage disequilibrium (LD) with rs2853677 and influence TERT expression.

METHODS

LD pattern in rs2853677 nearby region was analyzed based on 1000 genomes data for three representative populations in the world and functional genomics research was performed for this locus.

RESULTS

Only one SNP, rs2736099, is in strong LD with rs2853677 in East Asian. Dual-luciferase reporter assay verifies that rs2736099 can regulate gene expression and should be the causal SNP for this disease. Through chromosome conformation capture assay, it is disclosed that the enhancer surrounding rs2736099 can interact with TERT promoter. Through chromatin immunoprecipitation, the transcription factor SP1 (Sp1 transcription factor) is recognized for the chromatin segment spanning rs2736099.

CONCLUSIONS

Our results provide the missing piece between genetic variation at this locus and lung cancer risk, which is also applied to tumorigenesis in other tissues and cell proliferation.

LmTag: functional-enrichment and imputation-aware tag SNP selection for population-specific genotyping arrays

Despite the rapid development of sequencing technology, single-nucleotide polymorphism (SNP) arrays are still the most cost-effective genotyping solutions for large-scale genomic research and applications. Recent years have witnessed the rapid development of numerous genotyping platforms of different sizes and designs, but population-specific platforms are still lacking, especially for those in developing countries. SNP arrays designed for these countries should be cost-effective (small size), yet incorporate key information needed to associate genotypes with traits. A key design principle for most current platforms is to improve genome-wide imputation so that more SNPs not included in the array (imputed SNPs) can be predicted. However, current tag SNP selection methods mostly focus on imputation accuracy and coverage, but not the functional content of the array. It is those functional SNPs that are most likely associated with traits. Here, we propose LmTag, a novel method for tag SNP selection that not only improves imputation performance but also prioritizes highly functional SNP markers. We apply LmTag on a wide range of populations using both public and in-house whole-genome sequencing databases. Our results show that LmTag improved both functional marker prioritization and genome-wide imputation accuracy compared to existing methods. This novel approach could contribute to the next generation genotyping arrays that provide excellent imputation capability as well as facilitate array-based functional genetic studies. Such arrays are particularly suitable for under-represented populations in developing countries or non-model species, where little genomics data are available while investment in genome sequencing or high-density SNP arrays is limited. $\textrm{LmTag}$ is available at: https://github.com/datngu/LmTag.

Informative SNP Selection Based on a Fuzzy Clustering and Improved Binary Particle Swarm Optimization Algorithm

Single-nucleotide polymorphism (SNP) involves the replacement of a single nucleotide in a deoxyribonucleic acid (DNA) sequence and is often linked to the development of specific diseases. Although current genotyping methods can tag SNP loci within biological samples to provide accurate genetic information for a disease associated, they have limited prediction accuracy. Furthermore, they are complex to perform and may result in the prediction of an excessive number of tag SNP loci, which may not always be associated with the disease. Therefore in this manuscript, we aimed to evaluate the impact of a newly optimized fuzzy clustering and binary particle swarm optimization algorithm (FCBPSO) on the accuracy and running time of informative SNP selection. Fuzzy clustering and FCBPSO were first applied to identify the equivalence relation and the candidate tag SNP set to reduce the redundancy between loci. The FCBPSO algorithm was then optimized and used to obtain the final tag SNP set. The prediction performance and running time of the newly developed model were compared with other traditional methods, including NMC, SPSO, and MCMR. The prediction accuracy of the FCBPSO algorithm was always higher than that of the other algorithms especially as the number of tag SNPs increased. However, when the number of tag SNPs was low, the prediction accuracy of FCBPSO was slightly lower than that of MCMR (add prediction accuracy values for each algorithm). However, the running time of the FCBPSO algorithm was always lower than that of MCMR. FCBPSO not only reduced the size and dimension of the optimization problem but also simplified the training of the prediction model. This improved the prediction accuracy of the model and reduced the running time when compared with other traditional methods.

Genetic LGALS1 Variants Are Associated with Heterogeneity in Galectin-1 Serum Levels in Patients with Early Arthritis

Galectin 1 (Gal1) exerts immunomodulatory effects leading to therapeutic effects in autoimmune animal models. Patients with rheumatoid arthritis have been reported to show higher Gal1 serum levels than the healthy population. Our study aimed to find genetic variants on the Gal1 gene (LGALS1) modulating its expression and/or clinical features in patients with early arthritis (EA). LGALS1 was sequenced in 53 EA patients to characterize all genetic variants. Then, we genotyped rs9622682, rs929039, and rs4820293, which covered the main genetic variation in LGALS1, in 532 EA patients. Gal1 and IL-6 serum levels were measured by ELISA and Gal1 also by western blot (WB) in lymphocytes from patients with specific genotypes. Once disease activity improved with treatment, patients with at least one copy of the minor allele in rs9622682 and rs929039 or those with GG genotype in rs4820293 showed significantly higher Gal1 serum levels (p < 0.05). These genotypic combinations were also associated with higher Gal1 expression in lymphocytes by WB and lower IL-6 serum levels in EA patients. In summary, our study suggests that genetic variants studied in LGALS1 can explain heterogeneity in Gal1 serum levels showing that patients with higher Gal1 levels have lower serum IL-6 levels.

Depressive Symptoms Predict Clinical Recurrence of Inflammatory Bowel Disease

BACKGROUND

Inflammatory bowel disease (IBD) patients are at high risk for depression, and depression has been shown to affect disease course. We examined interrelations between depression, genetic risk factors for depression, and IBD flares.

METHOD

In 1973 patients (1137 Crohn's disease, 836 ulcerative colitis) of the Swiss IBD Cohort Study (SIBDCS), depressive status (hospital anxiety and depression subscale for depression, HADS-D ≥11) was assessed on a yearly basis. We investigated the impact of depression on IBD-relevant clinical outcomes in Cox proportional hazards models. We used active disease (CDAI ≥150 or MTWAI ≥10) and 2 published composite flare definitions-FNCE (physician-reported flare, nonresponse to therapy, new complication, or extraintestinal manifestation) and AFFSST (active disease, physician-reported flare, fistula, stenosis, and new systemic therapy)-as clinical end points. Additionally, 62 preselected single nucleotide polymorphisms (SNPs) were screened for cross-sectional associations with depression, and if present, their predictive value for future depression and clinical deterioration was assessed.

RESULTS

Depression was a strong risk factor for disease-related end points, including active disease (adjusted hazard ratio [aHR], 3.55; P < 0.001), AFFSST (aHR, 1.62; P < 0.001), and FNCE (aHR, 1.35; P = 0.019). The SNP rs2522833 was significantly associated with depression at enrollment (q = 0.059). The TC allele of rs588765 was negatively associated with the presence of depression at enrollment (q = 0.050) and after enrollment (aHR, 0.67; P = 0.035) and with fewer active disease states (aHR, 0.72; P = 0.045) during follow-up.

CONCLUSION

In IBD, depressive symptoms and inflammatory activity are intimately related. Depressive symptoms were a strong predictor of clinical deterioration, and genetic markers may play a role in this relationship.

Effects of cyclophosphamide related genetic variants on clinical outcomes of adult hematopoietic cell transplant patients

PURPOSE

Genetic variants may influence the pharmacokinetics and pharmacodynamics (PKPD) of cyclophosphamide (CY). CY plays a critical role in conditioning chemotherapy for hematopoietic cell transplantation (HCT), but its use is limited by toxicity. We explored the effect of genetic variants, potentially affecting PKPD of CY, and outcomes after HCT.

METHODS

This observational pharmacogenomic study included 85 adults with hematologic malignancies who received reduced intensity conditioning with CY, fludarabine, and total body irradiation. We collected recipient DNA prior to HCT and evaluated 97 candidate variants in 66 genes and 3 metabolism phenotypes potentially involved in PKPD pathways of CY. In multivariable analysis we investigated the association between the genotypes and four clinical outcomes: Day 180 non-relapse mortality (NRM) and day 180 overall survival (OS), acute graft-versus-host-disease (aGVHD) grades 2-4, and engraftment. p values were not adjusted for multiple testing.

RESULTS

The median recipient age was 63 years (range 21-75). Acute myeloid leukemia was the most common diagnosis (34%; n = 29). In multivariable analysis adjusted for exposure to phosphoramide mustard, the final active metabolite of CY, we identified 6 variants in 6 genes associated with at least one of the clinical outcomes. An ABCC4 variant (rs9561778) was associated with poor Day 180 NRM (p < 0.01), MUTYH variant (rs3219484) with higher Day 180 NRM and aGVHD (both p < 0.01), and SYNE1 variant (rs4331993) with better Day 180 OS and engraftment (both p ≤ 0.01).

CONCLUSION

The present study suggests that genetic variants influencing the PKPD of CY may help identify patients at risk for inferior outcomes after HCT using CY-based reduced-intensity conditioning.

Novel Design of Imputation-Enabled SNP Arrays for Breeding and Research Applications Supporting Multi-Species Hybridization

Array-based single nucleotide polymorphism (SNP) genotyping platforms have low genotype error and missing data rates compared to genotyping-by-sequencing technologies. However, design decisions used to create array-based SNP genotyping assays for both research and breeding applications are critical to their success. We describe a novel approach applicable to any animal or plant species for the design of cost-effective imputation-enabled SNP genotyping arrays with broad utility and demonstrate its application through the development of the Illumina Infinium Wheat Barley 40K SNP array Version 1.0. We show that the approach delivers high quality and high resolution data for wheat and barley, including when samples are jointly hybridised. The new array aims to maximally capture haplotypic diversity in globally diverse wheat and barley germplasm while minimizing ascertainment bias. Comprising mostly biallelic markers that were designed to be species-specific and single-copy, the array permits highly accurate imputation in diverse germplasm to improve the statistical power of genome-wide association studies (GWAS) and genomic selection. The SNP content captures tetraploid wheat (A- and B-genome) and Aegilops tauschii Coss. (D-genome) diversity and delineates synthetic and tetraploid wheat from other wheat, as well as tetraploid species and subgroups. The content includes SNP tagging key trait loci in wheat and barley, as well as direct connections to other genotyping platforms and legacy datasets. The utility of the array is enhanced through the web-based tool, Pretzel (https://plantinformatics.io/) which enables the content of the array to be visualized and interrogated interactively in the context of numerous genetic and genomic resources to be connected more seamlessly to research and breeding. The array is available for use by the international wheat and barley community.

Novel association between a transient receptor potential cation channel subfamily M member 5 expression quantitative trait locus rs35197079 and decreased susceptibility of gestational diabetes mellitus in a Chinese population

AIMS/INTRODUCTION

Emerging evidence suggests that expression quantitative trait loci (eQTLs) are more likely to associate with complex diseases. Transient receptor potential cation channel subfamily M member 5 (TRPM5) is a ubiquitously expressed voltage-gated cation channel that acts indispensably to trigger insulin secretion in pancreatic β-cells. The present study evaluated the association between TRPM5 eQTL single-nucleotide polymorphisms and the risk of gestational diabetes mellitus (GDM) in a Chinese population.

MATERIALS AND METHODS

A total of 380 unrelated Chinese pregnant women including 241 GDM patients and 139 controls were included in this study. The eQTL single-nucleotide polymorphisms of TRPM5 were obtained from the GTEx eQTL Browser, and were subsequently genotyped using the Agena MassARRAY iPLEX platform.

RESULTS

Logistic regression analysis and linear regression analysis showed that rs35197079 and rs74848824 were significantly associated with reduced GDM risk and lower fasting plasma glucose levels after adjusting confounder factors in dominant genetic models. Stratification analysis based on pre-pregnancy body mass index validated a strong association between rs35197079 and GDM susceptibility in underweight and normal weight individuals. Luciferase and electrophoretic mobility shift assays carried out in rat pancreatic β-cells showed that rs35197079 was functional.

CONCLUSIONS

The TRPM5 eQTL single-nucleotide polymorphism rs35197079 was associated with decreased GDM susceptibility in a Chinese population, especially in underweight and normal weight pregnant women, and it was functional in modulating gene transcription.

Genetic control of arbuscular mycorrhizal colonization by Rhizophagus intraradices in Helianthus annuus (L.)

Plant symbiosis with arbuscular mycorrhizal (AM) fungi provides many benefits, including increased nutrient uptake, drought tolerance, and belowground pathogen resistance. To develop a better understanding of the genetic architecture of mycorrhizal symbiosis, we conducted a genome-wide association study (GWAS) of this plant-fungal interaction in cultivated sunflower. A diversity panel of cultivated sunflower (Helianthus annuus L.) was phenotyped for root colonization under inoculation with the AM fungus Rhizophagus intraradices. Using a mixed linear model approach with a high-density genetic map, we identified genomic regions that are likely associated with R. intraradices colonization in sunflower. Additionally, we used a set of twelve diverse lines to assess the effect that inoculation with R. intraradices has on dried shoot biomass and macronutrient uptake. Colonization among lines in the mapping panel ranged from 0-70% and was not correlated with mycorrhizal growth response, shoot phosphorus response, or shoot potassium response among the Core 12 lines. Association mapping yielded three single-nucleotide polymorphisms (SNPs) that were significantly associated with R. intraradices colonization. This is the first study to use GWAS to identify genomic regions associated with AM colonization in an Asterid eudicot species. Three genes of interest identified from the regions containing these SNPs are likely related to plant defense.

The Bnapus50K array: a quick and versatile genotyping tool for Brassica napus genomic breeding and research

Rapeseed is a globally cultivated commercial crop, primarily grown for its oil. High-density single nucleotide polymorphism (SNP) arrays are widely used as a standard genotyping tool for rapeseed research, including for gene mapping, genome-wide association studies, germplasm resource analysis, and cluster analysis. Although considerable rapeseed genome sequencing data have been released, DNA arrays are still an attractive choice for providing additional genetic data in an era of high-throughput whole-genome sequencing. Here, we integrated re-sequencing DNA array data (32,216, 304 SNPs) from 505 inbred rapeseed lines, allowing us to develop a sensitive and efficient genotyping DNA array, Bnapus50K, with a more consistent genetic and physical distribution of probes. A total of 42,090 high-quality probes were filtered and synthesized, with an average distance between adjacent SNPs of 8 kb. To improve the practical application potential of this array in rapeseed breeding, we also added 1,618 functional probes related to important agronomic traits such as oil content, disease resistance, male sterility, and flowering time. The additional probes also included those specifically for detecting genetically modified material. These probes show a good detection efficiency and are therefore useful for gene mapping, along with crop variety improvement and identification. The novel Bnapus50K DNA array developed in this study could prove to be a quick and versatile genotyping tool for B. napus genomic breeding and research.

Genome-wide association study identifies new loci associated with noise-induced tinnitus in Chinese populations

Background

Tinnitus is an auditory phantom sensation in the absence of an acoustic stimulus, which affects nearly 15% of the population. Excessive noise exposure is one of the main causes of tinnitus. To now, the knowledge of the genetic determinants of susceptibility to tinnitus remains limited.

Results

We performed a two-stage genome-wide association study (GWAS) and identified that two single nucleotide polymorphisms (SNPs), rs2846071 located in the intergenic region at 11q13.5 (odds ratio [OR] = 2.14, 95% confidence interval [CI] = 1.96–3.40, combined P = 4.89 × 10^− 6) and rs4149577 located in the intron of TNFRSF1A gene at 12p13.31 (OR = 2.05, 95% CI = 1.89–2.51, combined P = 6.88 × 10^− 6), are significantly associated with the susceptibility to noise-induced tinnitus. Furthermore, the expression quantitative trait loci (eQTL) analyses revealed that rs2846071 is significantly correlated with the expression of WNT11 gene, and rs4149577 with the expression of TNFRSF1A gene in multiple brain tissues (all P < 0.05). The newly identified candidate gene WNT11 is involved in Wnt pathway, and TNFRSF1A in the tumor necrosis factor pathway, respectively. Pathway enrichment analyses also showed that these two pathways are closely relevant to tinnitus.

Conclusions

Our findings highlight two novel loci at 11q13.5 and 12p13.31 conferring susceptibility to noise-induced tinnitus. and suggest that the WNT11 and TNFRSF1A genes might be the candidate causal targets of 11q13.5 and 12p13.31 loci, respectively.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00987-y.

Is BRD7 associated with spermatogenesis impairment and male infertility in humans? A case-control study in a Han Chinese population

Background

Bromodomain-containing protein 7 (BRD7), a member of the bromodomain-containing protein family, plays important roles in chromatin modification and transcriptional regulation. A recent model of Brd7-knockout mice presented azoospermia and male infertility, implying the potential role of BRD7 in spermatogenic failure in humans. This case-control study aimed to explore the association of the BRD7 gene with spermatogenic efficiency and the risk of spermatogenic defects in humans.

Results

A total of six heterozygous variants were detected in the coding and splicing regions of the BRD7 gene in patients with azoospermia. For each of four rare variants predicted to potentially damage BRD7 function, we further identified these four variants in oligozoospermia and normozoospermia as well. However, no difference in the allele and genotype frequencies of rare variants were observed between cases with spermatogenic failure and controls with normozoospermia; the sperm products of variant carriers were similar to those of noncarriers. Moreover, similar distribution of the alleles, genotypes and haplotypes of seven tag single nucleotide polymorphisms (tagSNPs) was observed between the cases with azoospermia and oligozoospermia and controls with normozoospermia; associations of tagSNP-distinguished BRD7 alleles with sperm products were not identified.

Conclusions

The lack of an association of BRD7-linked rare and common variants with spermatogenic failure implied a limited contribution of the BRD7 gene to spermatogenic efficiency and susceptibility to male infertility in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12610-021-00139-3.

Genome-Wide Association Studies for Fusarium Head Blight Resistance and Its Trade-Off With Grain Yield in Soft Red Winter Wheat

Identification of quantitative trait loci for Fusarium head blight (FHB) resistance from different sources and pyramiding them into cultivars could provide effective protection against FHB. The objective of this study was to characterize a soft red winter wheat (SRWW) breeding population that has been subjected to intense germplasm introduction and alien introgression for FHB resistance in the past. The population was evaluated under misted FHB nurseries inoculated with Fusarium graminearum-infested corn spawn for two years. Phenotypic data included disease incidence (INC), disease severity (SEV), Fusarium damaged kernels (FDK), FHB index (FHBdx), and deoxynivalenol concentration (DON). Genome-wide association studies using 13,784 SNP markers identified 25 genomic regions at -logP ≥ 4.0 that were associated with five FHB-related traits. Of these 25, the marker trait associations that explained more than 5% phenotypic variation were localized on chromosomes 1A, 2B, 3B, 5A, 7A, 7B, and 7D, and from diverse sources including adapted SRWW lines such as Truman and Bess, and unadapted common wheat lines such as Ning7840 and Fundulea 201R. Furthermore, individuals with favorable alleles at the four loci Fhb1, Qfhb.nc-2B.1 (Q_2B.1), Q_7D.1, and Q_7D.2 showed better FDK and DON scores (but not INC, SEV, and FHBdx) compared with other allelic combinations. Our data also showed while pyramiding multiple loci provides protection against FHB disease, it has a significant trade-off with grain yield.

A current guide to candidate gene association studies

Important factors contribute to a gained momentum in candidate gene association studies (CGASs), including the generalized use of next-generation sequencing (NGS), growing opportunities for hospital-based research, and the availability of open-source databases and bioinformatics tools. This article summarizes the general principles and analytical methods as a guide to CGASs in today's favorable context.

Genome Divergence and Dynamics in the Thin-Tailed Desert Sheep From Sudan

With climate change bound to affect food and feed production, emphasis will shift to resilient and adapted indigenous livestock to sustain animal production. However, indigenous livestock comprise several varieties, strains and ecotypes whose genomes are poorly characterized. Here, we investigated genomic variation in an African thin-tailed Desert Sheep sampled in Sudan, using 600K genotype data generated from 92 individuals representing five ecotypes. We included data from 18 fat-tailed and 45 thin-tailed sheep from China, to investigate shared ancestry and perform comparative genomic analysis. We observed a clear genomic differentiation between the African thin-tailed Desert Sheep and the Chinese thin-tailed and fat-tailed sheep, suggesting a broad genetic structure between the fat-tailed and thin-tailed sheep in general, and that at least two autosomal gene pools comprise the genome profile of the thin-tailed sheep. Further analysis detected two distinct genetic clusters in both the African thin-tailed Desert Sheep and the Chinese thin-tailed sheep, suggesting a fine-scale and complex genome architecture in thin-tailed sheep. Selection signature analysis suggested differences in adaptation, production, reproduction and morphology likely underly the fine-scale genetic structure in the African thin-tailed Desert Sheep. This may need to be considered in designing breeding programs and genome-wide association studies.

Weak Association Between the Glutamate Decarboxylase 1 Gene (GAD1) and Schizophrenia in Han Chinese Population

Objectives

Schizophrenia (SZ) is a complex psychiatric disorder with high heritability, and genetic components are thought to be pivotal risk factors for this illness. The glutamate decarboxylase 1 gene (GAD1) was hypothesized to be a candidate risk locus for SZ given its crucial role in the GABAergic neurotransmission system, and previous studies have examined the associations of single nucleotide polymorphisms (SNPs) spanning the GAD1 gene with SZ. However, inconsistent results were obtained. We hence examined the associations between GAD1 SNPs and SZ in two independent case-control samples of Han Chinese ancestry.

Materials and Methods

Two Han Chinese SZ case-control samples, referred as the discovery sample and the replication sample, respectively, were recruited for the current study. The discovery sample comprised of 528 paranoid SZ cases (with age of first onset ≥ 18) and 528 healthy controls; the independent replication sample contained 1,256 early onset SZ cases (with age of first onset < 18) and 2,661 healthy controls. Logistic regression analysis was performed to examine the associations between GAD1 SNPs and SZ.

Results

Ten SNPs covering GAD1 gene were analyzed in the discovery sample, and two SNPs showed nominal associations with SZ (rs2241165, P = 0.0181, OR = 1.261; rs2241164, P = 0.0225, OR = 1.219). SNP rs2241164 was also nominally significant in the independent replication sample (P = 0.0462, OR = 1.110), and the significance became stronger in a subsequent meta-analysis combining both discovery and replication samples (P = 0.00398, OR = 1.138). Nevertheless, such association could not survive multiple corrections, although the effect size of rs2241164 was comparable with other SZ risk loci identified in genome-wide association studies (GWAS) in Han Chinese population. We also examined the associations between GAD1 SNPs and SZ in published datasets of SZ GWAS in East Asians and Europeans, and no significant associations were observed.

Conclusion

We observed weak associations between GAD1 SNPs and risk of SZ in Han Chinese populations. Further analyses in larger Han Chinese samples with more detailed phenotyping are necessary to elucidate the genetic correlation between GAD1 SNPs and SZ.

Japonica Array NEO with increased genome-wide coverage and abundant disease risk SNPs

Ethnic-specific SNP arrays are becoming more important to increase the power of genome-wide association studies in diverse population. In the Tohoku Medical Megabank Project, we have been developing a series of Japonica Arrays (JPA) for genotyping participants based on reference panels constructed from whole-genome sequence data of the Japanese population. Here, we designed a novel version of the SNP array for the Japanese population, called Japonica Array NEO (JPA NEO), comprising a total of 666,883 markers. Among them, 654,246 tag SNPs of autosomes and X chromosome were selected from an expanded reference panel of 3,552 Japanese, 3.5KJPNv2, using pairwise r² of linkage disequilibrium measures. Additionally, 28,298 markers were included for the evaluation of previously identified disease risk markers from the literature and databases, and those present in the Japanese population were extracted using the reference panel. Through genotyping 286 Japanese samples, we found that the imputation quality r² and INFO score in the minor allele frequency bin >2.5–5% were >0.9 and >0.8, respectively, and >12 million markers were imputed with an INFO score >0.8. From these results, JPA NEO is a promising tool for genotyping the Japanese population with genome-wide coverage, contributing to the development of genetic risk scores.

Genetic Overlap Between Attention Deficit/Hyperactivity Disorder and Autism Spectrum Disorder in SHANK2 Gene

Background: Recent findings indicated a high comorbidity between attention-deficit/hyperactivity disorder (ADHD) and autism spectrum disorder (ASD), as well as shared genetic influences on them. The latter might contribute at least partly to the former clinical scenario. This study aimed at investigating whether SHANK genes were potential pleiotropic genes to the two said disorders, underlying their genetic overlap.

Methods: This study recruited 298 boys with ADHD (including 256 family trios of 1 ADHD boy and his 2 biological parents), 134 boys with ASD, 109 boys with both ADHD and ASD, and 232 typically developing boys as community controls. They were aged between 6 and 11 years old.

Results: There was no significant difference in allele frequency of a number of single nucleotide polymorphisms (SNPs) in SHANK2/SHANK3 between the three clinical groups (ADHD, ASD, and ADHD + ASD) and between the two control groups (community controls and pseudo-controls), respectively. The three clinical groups and the two control groups were thus, respectively, combined. A comparison between the two aggregated samples identified significant evidence of disease association for three SHANK2 SNPs with both ADHD and ASD, even after multiple testing correction: rs11236616 (OR = 0.762, permuted p = 0.0376), rs7106631 (OR = 0.720, permuted p = 0.0034), and rs9888288 (OR = 0.770, permuted p = 0.0407). Comparisons among individual groups pointed to a similar trend of findings.

Conclusion:SHANK2 could be considered a potential pleiotropic gene underlying the genetic overlap between ADHD and ASD. This might contribute partly to their high comorbidity in the afflicted children.

Single nucleotide polymorphisms (SNPs) in prostate cancer: its implications in diagnostics and therapeutics

Prostate cancer is one of the most frequently diagnosed malignancies in developed countries and approximately 248,530 new cases of prostate cancer are likely to be diagnosed in the United States in 2021. During the late 1990s and 2000s, the prostate cancer-related death rate has decreased by 4% per year on average because of advancements in prostate-specific antigen (PSA) testing. However, the non-specificity of PSA to distinguish between benign and malignant forms of cancer is a major concern in the management of prostate cancer. Despite other risk factors in the pathogenesis of prostate cancer, recent advancement in molecular genetics suggests that genetic heredity plays a crucial role in prostate carcinogenesis. Approximately, 60% of heritability and more than 100 well-recognized single-nucleotide-polymorphisms (SNPs) have been found to be associated with prostate cancer and constitute a major risk factor in the development of prostate cancer. Recent findings revealed that a low to moderate effect on the progression of prostate cancer of individual SNPs was observed compared to a strong progressive effect when SNPs were in combination. Here, in this review, we made an attempt to critically analyze the role of SNPs and associated genes in the development of prostate cancer and their implications in diagnostics and therapeutics. A better understanding of the role of SNPs in prostate cancer susceptibility may improve risk prediction, enhance fine-mapping, and furnish new insights into the underlying pathophysiology of prostate cancer.

NF-κB1 intronic region polymorphisms as risk factor for head and neck cancer in HPV-infected population from Pakistan

BACKGROUND

Head and neck cancer (HNC) developed due to the number of risk factors, including infection of Human Papillomavirus (HPV). The genetic predisposition also plays an important role in deregulating the NF-κB pathway, and certain polymorphisms are reported to affect the pathway genes.

OBJECTIVES

The present study was conducted for the detection of HPV and polymorphisms in the NF-κB1 gene of HNC patients in the Pakistani population.

METHODS

Genomic DNA from HNC tumors samples were extracted using the Exgene SV DNA extraction Kit. Allele-specific PCR and direct sequencing were done for analysis of NF-κB1 SNPs, 94ins/del (rs28362491), rs1598858, and rs4648068.

RESULTS

The genotypes AG (36.2%/ 12%) of rs1598858, and AG (28.3%/ 12%) and GG (28.3%/ 22%) of rs4648068 were associated with significantly (p≤0.05) increased risk of head and neck cancer in studied population. Furthermore, among the HNC cases, genotypes AGrs1598858 (p≤0.014) and GGrs4648068 (p≤0.001) had increased risk of HPV related cancers. Tobacco use (OR-3.158442; [1.140, 8.754]), lymph nodes involvement (OR 4.05128; [1.854, 8.852]), and poorly differentiated tumors (OR 1.997155; [0.940, 4.245]) were positively associated with HPV induced cancers.

CONCLUSION

It was the first comprehensive study from Pakistan, to evaluate the polymorphic variants of NF-κB1. Genotypes AGrs4648068, GGrs4648068, and AGrs1598858 of NF-κB1 gene are associated with increased risk of head and neck cancers in the Pakistani population. It can be concluded that HPV infection, lymph nodes and tobacco use can act synergetic to each other and add up in modulating HNC when present together with intronic SNPs of NF-κB1 gene.

Grouping of genomic markers in populations with family structure

BACKGROUND

Linkage and linkage disequilibrium (LD) between genome regions cause dependencies among genomic markers. Due to family stratification in populations with non-random mating in livestock or crop, the standard measures of population LD such as [Formula: see text] may be biased. Grouping of markers according to their interdependence needs to account for the actual population structure in order to allow proper inference in genome-based evaluations.

RESULTS

Given a matrix reflecting the strength of association between markers, groups are built successively using a greedy algorithm; largest groups are built at first. As an option, a representative marker is selected for each group. We provide an implementation of the grouping approach as a new function to the R package hscovar. This package enables the calculation of the theoretical covariance between biallelic markers for half- or full-sib families and the derivation of representative markers. In case studies, we have shown that the number of groups comprising dependent markers was smaller and representative SNPs were spread more uniformly over the investigated chromosome region when the family stratification was respected compared to a population-LD approach. In a simulation study, we observed that sensitivity and specificity of a genome-based association study improved if selection of representative markers took family structure into account.

CONCLUSIONS

Chromosome segments which frequently recombine in the underlying population can be identified from the matrix of pairwise dependence between markers. Representative markers can be exploited, for instance, for dimension reduction prior to a genome-based association study or the grouping structure itself can be employed in a grouped penalization approach.

Genome-wide association analysis unveils novel QTLs for seminal root system architecture traits in Ethiopian durum wheat

BACKGROUND

Genetic improvement of root system architecture is essential to improve water and nutrient use efficiency of crops or to boost their productivity under stress or non-optimal soil conditions. One hundred ninety-two Ethiopian durum wheat accessions comprising 167 historical landraces and 25 modern cultivars were assembled for GWAS analysis to identify QTLs for root system architecture (RSA) traits and genotyped with a high-density 90 K wheat SNP array by Illumina.

RESULTS

Using a non-roll, paper-based root phenotyping platform, a total of 2880 seedlings and 14,947 seminal roots were measured at the three-leaf stage to collect data for total root length (TRL), total root number (TRN), root growth angle (RGA), average root length (ARL), bulk root dry weight (RDW), individual root dry weight (IRW), bulk shoot dry weight (SDW), presence of six seminal roots per seedling (RT6) and root shoot ratio (RSR). Analysis of variance revealed highly significant differences between accessions for all RSA traits. Four major (- log₁₀P ≥ 4) and 34 nominal (- log₁₀P ≥ 3) QTLs were identified and grouped in 16 RSA QTL clusters across chromosomes. A higher number of significant RSA QTL were identified on chromosome 4B particularly for root vigor traits (root length, number and/or weight).

CONCLUSIONS

After projecting the identified QTLs on to a high-density tetraploid consensus map along with previously reported RSA QTL in both durum and bread wheat, fourteen nominal QTLs were found to be novel and could potentially be used to tailor RSA in elite lines. The major RGA QTLs on chromosome 6AL detected in the current study and reported in previous studies is a good candidate for cloning the causative underlining sequence and identifying the beneficial haplotypes able to positively affect yield under water- or nutrient-limited conditions.

A pharmacogenetics study of platinum-based chemotherapy in lung cancer: ABCG2 polymorphism and its genetic interaction with SLC31A1 are associated with response and survival

Objective: The expression and function of platinum transporters affect drug tissue concentration and therapeutic effects. We had previously characterized functional variant of platinum intake transporter SLC31A1 gene. We aimed to investigate the association of platinum efflux transporter gene ABCG2 polymorphism and combined ABCG2 and SLC31A1 polymorphisms with clinical outcomes of NSCLC patients receiving platinum-based chemotherapy. Methods: We genotyped thirteen tagging and functional SNPs of ABCG2 in 1004 patients, and assessed their association with response, toxicity and survival using unconditional logistic regression and Cox proportional hazards regression analyses respectively. Results: Nonsynonymous rs2231142 (odds ratio [OR] 2.07; 95 % confidence interval [CI] 1.26-3.63), rs1871744 (OR 0.60; 95 % CI 0.42-0.87) and their haplotype and diplotype were associated with objective response. Rs4148157 was associated with shorter overall survival (Log-rank P = 0.002; hazard ratio [HR] 1.22; 95 % CI 1.05-1.42). Furthermore, the combined SLC31A1 rs2233914 and ABCG2 rs1871744 genotype was significantly associated with poor response (OR 0.31; 95 % CI 0.17-0.56; P _interaction = 0.003). And the combined genotypes of the functional rs10759637 of SLC31A1 and the nonsynonymous rs2231142 (Log-rank P = 5.20×10^-5; HR 1.47; 95 % CI 1.19-1.81; P _interaction = 0.007) or linked rs4148157 of ABCG2 were significantly associated with poor survival. Conclusion: This study reveals divergent association of ABCG2 polymorphism with response and survival of NSCLC patients receiving platinum-based chemotherapy, demonstrates the combined effects of functional variants of ABCG2 and SLC31A1 on clinical outcomes, and highlights pharmacogenetic relevance of platinum transporter genes interaction.

Engulfment and cell motility 1 (ELMO1) and apolipoprotein A1 (APOA1) as candidate genes for sickle cell nephropathy

Sickle cell disease (SCD) and apolipoprotein L1 (APOL1) G1/G2 variants increase chronic kidney disease (CKD) risk in African Americans by poorly understood mechanisms. We applied bioinformatics to identify new candidate genes associated with SCD-related CKD. An interaction network demonstrated APOA1 connecting haemoglobin subunit β (HBB) and APOL1 with 36 other candidate genes. Gene expression revealed upregulation of engulfment and cell motility 1 (ELMO1) and downregulation of APOA1 in the kidney cortex of SCD versus non-SCD mice. Analysis of candidate genes identified ELMO1 rs10951509 to be associated with albuminuria and APOA1 rs11216132 with haemoglobinuria in patients with SCD. A bioinformatic approach highlights ELMO1 and APOA1 as potentially associated with SCD nephropathy.

Associations between Plasma Choline Metabolites and Genetic Polymorphisms in One-Carbon Metabolism in Postmenopausal Women: The Women's Health Initiative Observational Study

BACKGROUND

Choline plays an integral role in one-carbon metabolism in the body, but it is unclear whether genetic polymorphisms are associated with variations in plasma choline and its metabolites.

OBJECTIVES

This study aimed to evaluate the association of genetic variants in choline and one-carbon metabolism with plasma choline and its metabolites.

METHODS

We analyzed data from 1423 postmenopausal women in a case-control study nested within the Women's Health Initiative Observational Study. Plasma concentrations of choline, betaine, dimethylglycine (DMG), and trimethylamine N-oxide were determined in 12-h fasting blood samples collected at baseline (1993-1998). Candidate and tagging single-nucleotide polymorphisms (SNPs) were genotyped in betaine-homocysteine S-methyltransferase (BHMT), BHMT2, 5,10-methylenetetrahydrofolate reductase (MTHFR), methylenetetrahydrofolate dehydrogenase (NADP+ dependent 1) (MTHFD1), 5-methyltetrahydrofolate-homocysteine methyltransferase (MTR), and 5-methyltetrahydrofolate-homocysteine methyltransferase reductase (MTRR). Linear regression was used to derive percentage difference in plasma concentrations per variant allele, adjusting for confounders, including B-vitamin biomarkers. Potential effect modification by plasma vitamin B-12, vitamin B-6, and folate concentrations and folic-acid fortification periods was examined.

RESULTS

The candidate SNP BHMT R239Q (rs3733890) was associated with lower concentrations of plasma betaine and DMG concentrations (-4.00% and -6.75% per variant allele, respectively; both nominal P < 0.05). Another candidate SNP, BHMT2 rs626105 A>G, was associated with higher plasma DMG concentration (13.0%; P < 0.0001). Several tagSNPs in these 2 genes were associated with plasma concentrations after correction for multiple comparisons. Vitamin B-12 status was a significant effect modifier of the association between the genetic variant BHMT2 rs626105 A>G and plasma DMG concentration.

CONCLUSIONS

Genetic variations in metabolic enzymes were associated with plasma concentrations of choline and its metabolites. Our findings contribute to the knowledge on the variation in blood nutrient concentrations in postmenopausal women.

Investigation of mechanisms underlying chaotic genetic patchiness in the intertidal marbled crab Pachygrapsus marmoratus (Brachyura: Grapsidae) across the Ligurian Sea

Background

Studies on marine community dynamics and population structures are limited by the lack of exhaustive knowledge on the larval dispersal component of connectivity. Genetic data represents a powerful tool in understanding such processes in the marine realm. When dealing with dispersion and connectivity in marine ecosystems, many evidences show patterns of genetic structure that cannot be explained by any clear geographic trend and may show temporal instability. This scenario is usually referred to as chaotic genetic patchiness, whose driving mechanisms are recognized to be selection, temporal shifts in local population dynamics, sweepstakes reproductive success and collective dispersal.

In this study we focused on the marbled crab Pachygrapsus marmoratus that inhabits the rocky shores of the Mediterranean Sea, Black Sea and East Atlantic Ocean, and disperses through planktonic larvae for about 1 month. P. marmoratus exhibits unexpectedly low connectivity levels at local scale, although well-defined phylogeographic patterns across the species’ distribution range were described. This has been explained as an effect of subtle geographic barriers or due to sweepstake reproductive success. In order to verify a chaotic genetic patchiness scenario, and to explore mechanisms underlying it, we planned our investigation within the Ligurian Sea, an isolated basin of the western Mediterranean Sea, and we genotyped 321 individuals at 11 microsatellite loci.

Results

We recorded genetic heterogeneity among our Ligurian Sea samples with the occurrence of genetic clusters not matching the original populations and a slight inter-population divergence, with the geographically most distant populations being the genetically most similar ones. Moreover, individuals from each site were assigned to all the genetic clusters. We also recorded evidences of self-recruitment and a higher than expected within-site kinship.

Conclusions

Overall, our results suggest that the chaotic genetic patchiness we found in P. marmoratus Ligurian Sea populations is the result of a combination of differences in reproductive success, en masse larval dispersion and local larval retention. This study defines P. marmoratus as an example of marine spawner whose genetic pool is not homogenous at population level, but rather split in a chaotic mosaic of slightly differentiated genetic patches derived from complex and dynamic ecological processes.

Temporal Variation in Genetic Composition of Migratory Helicoverpa Zea in Peripheral Populations

Migrant populations of Helicoverpa zea (Boddie) captured during 2002, 2005, 2016, and 2018 from Landisville and Rock Springs in Pennsylvania, USA were genotyped using 85 single nucleotide polymorphism (SNP) markers. Samples (n = 702) genotyped were divided into 16 putative populations based on collection time and site. Fixation indices (F-statistics), analysis of molecular variance, and discriminant analysis of principal components were used to examine within and among population genetic variation. The observed and expected heterozygosity in putative populations ranged from 0.317-0.418 and 0.320-0.359, respectively. Broad range of FST (0.0-0.2742) and FIS (0.0-0.2330) values indicated different genotype frequencies between and within the populations, respectively. High genetic diversity within and low genetic differentiation between populations was found in 2002 and 2005. Interestingly, high genetic differentiation between populations from two collection sites observed in 2018 populations was not evident in within-site comparisons of putative populations collected on different dates during the season. The shift of H. zea population genetic makeup in 2018 may be influenced by multiple biotic and abiotic factors including tropical storms. Continued assessment of these peripheral populations of H. zea will be needed to assess the impacts of genetic changes on pest control and resistance management tactics.