Genome-wide association study reveals 14 new SNPs and confirms two structural variants highly associated with the horned/polled phenotype in goats

Abstract Background There is a long-term interest in investigating the genetic basis of the horned/polled phenotype in domestic goats. Here, we report a genome-wide association study (GWAS) to detect the genetic loci affecting the polled phenotype in goats. Results We obtained a total of 13,980,209 biallelic SNPs, using the genotyping-by-sequencing data from 45 Jintang Black (JT) goats, which included 32 female and nine male goats, and four individuals with the polled intersex syndrome (PIS). Using a mixed-model based GWAS, we identified two association signals, which were located at 150,334,857–150,817,260 bp (P = 5.15 × 10− 119) and 128,286,704–131,306,537 bp (P = 2.74 × 10− 15) on chromosome 1. The genotype distributions of the 14 most significantly associated SNPs were completely correlated with horn status in goats, based on the whole-genome sequencing (WGS) data from JT and two other Chinese horned breeds. However, variant annotation suggested that none of the detected SNPs within the associated regions were plausible causal mutations. Via additional read-depth analyses and visual inspections of WGS data, we found a 10.1-kb deletion (CHI1:g. 129424781_129434939del) and a 480-kb duplication (CHI1:150,334,286–150,818,098 bp) encompassing two genes KCNJ15 and ERG in the associated regions of polled and PIS-affected goats. Notably, the 10.1-kb deletion also served as the insertion site for the 480-kb duplication, as validated by PCR and Sanger sequencing. Our WGS genotyping showed that all horned goats were homozygous for the reference alleles without either the structural variants (SVs), whereas the PIS-affected goats were homozygous for both the SVs. We also demonstrated that horned, polled, and PIS-affected individuals among 333 goats from JT and three other Chinese horned breeds can be accurately classified via PCR amplification and agarose gel electrophoresis of two fragments in both SVs. Conclusion Our results revealed that two genomic regions on chromosome 1 are major loci affecting the polled phenotypes in goats. We provided a diagnostic PCR to accurately classify horned, polled, and PIS-affected goats, which will enable a reliable genetic test for the early-in-life prediction of horn status in goats.

Download Full-text

Integration of genome-wide association study and expression quantitative trait locus mapping for identification of endometriosis-associated genes

Scientific Reports ◽

10.1038/s41598-020-79515-4 ◽

2021 ◽

Vol 11 (1) ◽

Author(s):

Ya-Ching Chou ◽

Ming-Jer Chen ◽

Pi-Hua Chen ◽

Ching-Wen Chang ◽

Mu-Hsien Yu ◽

...

Keyword(s):

Association Study ◽

Quantitative Trait ◽

Genome Wide Association Study ◽

Chromosome 9 ◽

Genome Wide Association ◽

Chromosome 1 ◽

Eqtl Analysis ◽

Taiwanese Population ◽

Genome Wide ◽

A Genome

AbstractTo determine whether genetic predisposition to endometriosis varies depending on ethnicity and in association with expression quantitative trait loci (eQTL) in a Taiwanese population. We conducted a genome-wide association study (GWAS) and replicated it in 259 individuals with laparoscopy-confirmed stage III or IV endometriosis (cases) and 171 women without endometriosis (controls). Their genomic DNA was extracted from blood and evaluated by the GWAS of Taiwan Biobank Array. Novel genetic variants that predispose individuals to endometriosis were identified using GWAS and replication, including rs10739199 (P = 6.75 × 10−5) and rs2025392 (P = 8.01 × 10−5) at chromosome 9, rs1998998 (P = 6.5 × 10−6) at chromosome 14, and rs6576560 (P = 9.7 × 10−6) at chromosome 15. After imputation, strong signals were exhibited by rs10822312 (P = 1.80 × 10−7) at chromosome 10, rs58991632 (P = 1.92 × 10−6) and rs2273422 (P = 2.42 × 10−6) at chromosome 20, and rs12566078 (P = 2.5 × 10−6) at chromosome 1. We used the Genotype-Tissue Expression (GTEx) database to observe eQTL. Among these SNPs, the cis-eQTL rs13126673 of inturned planar cell polarity protein (INTU) showed significant association with INTU expression (P = 5.1 × 10–33). Moreover, the eQTL analysis was performed on endometriotic tissues from women with endometriosis. The expression of INTU in 78 endometriotic tissue of women with endometriosis is associated with rs13126673 genotype (P = 0.034). To our knowledge, this is the first GWAS to link endometriosis and eQTL in a Taiwanese population.

Download Full-text

Genome-wide association study of phenotypes measuring progression from first cocaine or opioid use to dependence reveals novel risk genes

Exploration of Medicine ◽

10.37349/emed.2020.00032 ◽

2021 ◽

Author(s):

Richard Sherva ◽

Congcong Zhu ◽

Leah Wetherill ◽

Howard J. Edenberg ◽

Emma Johnson ◽

...

Keyword(s):

Association Study ◽

Genome Wide Association Study ◽

Meta Analysis ◽

Chromosome 21 ◽

Genome Wide Association ◽

Chromosome 1 ◽

Opioid Use ◽

Discovery Sample ◽

Genome Wide ◽

A Genome

Aim: Substance use disorders (SUD) result in substantial morbidity and mortality worldwide. Opioids, and to a lesser extent cocaine, contribute to a large percentage of this health burden. Despite their high heritability, few genetic risk loci have been identified for either opioid or cocaine dependence (OD or CD, respectively). A genome-wide association study of OD and CD related phenotypes reflecting the time between first self-reported use of these substances and a first DSM-IV dependence diagnosis was conducted. Methods: Cox proportional hazards regression in a discovery sample of 6,188 African-Americans (AAs) and 6,835 European-Americans (EAs) participants in a genetic study of multiple substance dependence phenotypes were used to test for association between genetic variants and these outcomes. The top findings were tested for replication in two independent cohorts. Results: In the discovery sample, three independent regions containing variants associated with time to dependence at P < 5 x 10-8 were identified, one (rs61835088 = 1.03 x 10-8) for cocaine in the combined EA-AA meta-analysis in the gene FAM78B on chromosome 1, and two for opioids in the AA portion of the sample in intergenic regions of chromosomes 4 (rs4860439, P = 1.37 x 10-8) and 9 (rs7032521, P = 3.30 x 10-8). After meta-analysis with data from the replication cohorts, the signal at rs61835088 improved (HR = 0.87, P = 3.71 x 10-9 and an intergenic SNP on chromosome 21 (rs2825295, HR = 1.14, P = 2.57 x 10-8) that missed the significance threshold in the AA discovery sample became genome-wide significant (GWS) for CD. Conclusions: Although the two GWS variants are not in genes with obvious links to SUD biology and have modest effect sizes, they are statistically robust and show evidence for association in independent samples. These results may point to novel pathways contributing to disease progression and highlight the utility of related phenotypes to better understand the genetics of SUDs.

Download Full-text

Identification of genetic loci and candidate genes related to soybean flowering through genome wide association study

BMC Genomics ◽

10.1186/s12864-019-6324-7 ◽

2019 ◽

Vol 20 (1) ◽

Author(s):

Minmin Li ◽

Ying Liu ◽

Yahan Tao ◽

Chongjing Xu ◽

Xin Li ◽

...

Keyword(s):

Association Study ◽

Flowering Time ◽

Candidate Genes ◽

Genome Wide Association Study ◽

Growth Period ◽

Snp Markers ◽

Genome Wide Association ◽

Chromosome 1 ◽

Genome Wide ◽

A Genome

Abstract Background As a photoperiod-sensitive and self-pollinated species, the growth periods traits play important roles in the adaptability and yield of soybean. To examine the genetic architecture of soybean growth periods, we performed a genome-wide association study (GWAS) using a panel of 278 soybean accessions and 34,710 single nucleotide polymorphisms (SNPs) with minor allele frequencies (MAF) higher than 0.04 detected by the specific-locus amplified fragment sequencing (SLAF-seq) with a 6.14-fold average sequencing depth. GWAS was conducted by a compressed mixed linear model (CMLM) involving in both relative kinship and population structure. Results GWAS revealed that 37 significant SNP peaks associated with soybean flowering time or other growth periods related traits including full bloom, beginning pod, full pod, beginning seed, and full seed in two or more environments at -log10(P) > 3.75 or -log10(P) > 4.44 were distributed on 14 chromosomes, including chromosome 1, 2, 3, 5, 6, 9, 11, 12, 13, 14, 15, 17, 18, 19. Fourteen SNPs were novel loci and 23 SNPs were located within known QTLs or 75 kb near the known SNPs. Five candidate genes (Glyma.05G101800, Glyma.11G140100, Glyma.11G142900, Glyma.19G099700, Glyma.19G100900) in a 90 kb genomic region of each side of four significant SNPs (Gm5_27111367, Gm11_10629613, Gm11_10950924, Gm19_34768458) based on the average LD decay were homologs of Arabidopsis flowering time genes of AT5G48385.1, AT3G46510.1, AT5G59780.3, AT1G28050.1, and AT3G26790.1. These genes encoding FRI (FRIGIDA), PUB13 (plant U-box 13), MYB59, CONSTANS, and FUS3 proteins respectively might play important roles in controlling soybean growth periods. Conclusions This study identified putative SNP markers associated with soybean growth period traits, which could be used for the marker-assisted selection of soybean growth period traits. Furthermore, the possible candidate genes involved in the control of soybean flowering time were predicted.

Download Full-text

Genome-wide association study of phenotypes measuring progression from first cocaine or opioid use to dependence reveals novel risk genes

Exploration of Medicine ◽

10.37349/emed.2021.00032 ◽

2021 ◽

Author(s):

Richard Sherva ◽

Congcong Zhu ◽

Leah Wetherill ◽

Howard J. Edenberg ◽

Emma Johnson ◽

...

Keyword(s):

Association Study ◽

Genome Wide Association Study ◽

Meta Analysis ◽

Chromosome 21 ◽

Genome Wide Association ◽

Chromosome 1 ◽

Opioid Use ◽

Discovery Sample ◽

Genome Wide ◽

A Genome

Aim: Substance use disorders (SUD) result in substantial morbidity and mortality worldwide. Opioids, and to a lesser extent cocaine, contribute to a large percentage of this health burden. Despite their high heritability, few genetic risk loci have been identified for either opioid or cocaine dependence (OD or CD, respectively). A genome-wide association study of OD and CD related phenotypes reflecting the time between first self-reported use of these substances and a first DSM-IV dependence diagnosis was conducted. Methods: Cox proportional hazards regression in a discovery sample of 6,188 African-Americans (AAs) and 6,835 European-Americans (EAs) participants in a genetic study of multiple substance dependence phenotypes were used to test for association between genetic variants and these outcomes. The top findings were tested for replication in two independent cohorts. Results: In the discovery sample, three independent regions containing variants associated with time to dependence at P < 5 × 10−8 were identified, one (rs61835088 = 1.03 × 10−8) for cocaine in the combined EA-AA meta-analysis in the gene FAM78B on chromosome 1, and two for opioids in the AA portion of the sample in intergenic regions of chromosomes 4 (rs4860439, P = 1.37 × 10−8) and 9 (rs7032521, P = 3.30 × 10−8). After meta-analysis with data from the replication cohorts, the signal at rs61835088 improved (HR = 0.87, P = 3.71 × 10−9 and an intergenic SNP on chromosome 21 (rs2825295, HR = 1.14, P = 2.57 × 10−8) that missed the significance threshold in the AA discovery sample became genome-wide significant (GWS) for CD. Conclusions: Although the two GWS variants are not in genes with obvious links to SUD biology and have modest effect sizes, they are statistically robust and show evidence for association in independent samples. These results may point to novel pathways contributing to disease progression and highlight the utility of related phenotypes to better understand the genetics of SUDs.

Download Full-text