Genetic variants in microRNA and microRNA biogenesis pathway genes and breast cancer risk among women of African ancestry

MicroRNAs (miRNA) regulate breast biology by binding to specific RNA sequences, leading to RNA degradation and inhibition of translation of their target genes. While germline genetic variations may disrupt some of these interactions between miRNAs and their targets, studies assessing the relationship between genetic variations in the miRNA network and breast cancer risk are still limited, particularly among women of African ancestry. We systematically put together a list of 822 and 10,468 genetic variants among primary miRNA sequences and 38 genes in the miRNA biogenesis pathway, respectively; and examined their association with breast cancer risk in the ROOT consortium which includes women of African ancestry. Findings were replicated in an independent consortium. Logistic regression was used to estimate the odds ratio (OR) and 95 % confidence intervals (CI). For overall breast cancer risk, three single-nucleotide polymorphisms (SNPs) in miRNA biogenesis genes DROSHA rs78393591 (OR = 0.69, 95 % CI: 0.55–0.88, P = 0.003), ESR1 rs523736 (OR = 0.88, 95 % CI: 0.82–0.95, P = 3.99 × 10^?4), and ZCCHC11 rs114101502 (OR = 1.33, 95 % CI: 1.11–1.59, P = 0.002), and one SNP in primary miRNA sequence (rs116159732 in miR-6826, OR = 0.74, 95 % CI: 0.63–0.89, P = 0.001) were found to have significant associations in both discovery and validation phases. In a subgroup analysis, two SNPs were associated with risk of estrogen receptor (ER)-negative breast cancer, and three SNPs were associated with risk of ER-positive breast cancer. Several variants in miRNA and miRNA biogenesis pathway genes were associated with breast cancer risk. Risk associations varied by ER status, suggesting potential new mechanisms in etiology.     

Sex influences eQTL effects of SLE and Sjögren’s syndrome-associated genetic polymorphisms

Background

Systemic lupus erythematosus (SLE) and primary Sjögren’s syndrome (pSS) are autoimmune disorders characterized by autoantibodies, dysregulated B cells, and notably high female-to-male incidence ratios. Genome-wide association studies have identified several susceptibility SNPs for both diseases. Many SNPs in the genome are expression quantitative trait loci (eQTLs), with context-dependent effects. Assuming that sex is a biological context, we investigated whether SLE/pSS SNPs act as eQTLs in B cells and used a disease-targeted approach to understand if they display sex-specific effects.

Methods

We used genome-wide genotype and gene expression data from primary B cells from 125 males and 162 females. The MatrixEQTL R package was used to identify eQTLs within a genomic window of 2 Mb centered on each of 22 established SLE and/or pSS susceptibility SNPs. To find sex-specific eQTLs, we used a linear model with a SNP * sex interaction term.

Results

We found ten SNPs affecting the expression of 16 different genes (FDR <?0.05). rs7574865-INPP1, rs7574865-MYO1B, rs4938573-CD3D, rs11755393-SNRPC, and rs4963128-PHRF1 were novel observations for the immune compartment and B cells. By analyzing the SNP * sex interaction terms, we identified six genes with differentially regulated expression in females compared to males, depending on the genotype of SLE/pSS-associated SNPs: SLC39A8 (BANK1 locus), CD74 (TNIP1 locus), PXK, CTSB (BLK/FAM167A locus), ARCN1 (CXCR5 locus), and DHX9 (NCF2 locus).

Conclusions

We identified several unknown sex-specific eQTL effects of SLE/pSS-associated genetic polymorphisms and provide novel insight into how gene-sex interactions may contribute to the sex bias in systemic autoimmune diseases.

     

Genome-wide association study of body weight in Wenshang Barred chicken based on the SLAF-seq technology

Chicken body weight (BW) is an economically important trait, and many studies have been conducted on genetic selection for BW. However, previous studies have detected functional chromosome mutations or regions using gene chips. The present study used the specific-locus amplified fragment sequencing (SLAF-seq) technology to perform a genome-wide association study (GWAS) on purebred Wengshang Barred chicken. A total of 1,286,715 single-nucleotide polymorphisms (SNPs) were detected, and 175,211 SNPs were selected as candidate SNPs for genome-wide association analysis using TASSEL general linear models. Six SNP markers reached genome-wide significance. Of these, rs732048524, rs735522839, rs738991545, and rs15837818 were significantly associated with body weight at 28 days (BW28), while rs314086457 and rs315694878 were significantly associated with BW120. These SNPs are close to seven genes (PRSS23, ME3, FAM181B, NABP1, SDPR, TSSK6L2, and RBBP8). Moreover, 24 BW-associated SNPs reached “suggestive” genome-wide significance. Of these, 6, 13, 1, and 4 SNPs were associated with BW28, BW56, BW80, and BW120, respectively. These results would enrich the studies on BW and promote the use of Chinese chicken, especially the Wenshang Barred chicken.     

An integrated approach of bioinformatic prediction and in vitro analysis identified that miR-34a targets <Emphasis Type="Italic">MET</Emphasis> and <Emphasis Type="Italic">AXL</Emphasis> in triple-negative breast cancer

Background

Breast cancer is the most prevalent cancer among women, and AXL and MET are the key genes in the PI3K/AKT/mTOR pathway as critical elements in proliferation and invasion of cancer cells. MicroRNAs (miRNAs) are small non-coding RNAs regulating the expression of genes.

Methods

Bioinformatic approaches were used to find a miRNA that simultaneously targets both AXL and MET 3′-UTRs. The expression of target miRNA was evaluated in triple-negative (MDA-MB-231) and HER2-overexpressing (SK-BR-3) breast cancer cell lines as well as normal breast cells, MCF-10A, using quantitative real-time PCR. Then, the miRNA was overexpressed in normal and cancer cell lines using a lentiviral vector system. Afterwards, effects of overexpressed miRNA on the expression of AXL and MET genes were evaluated using quantitative real-time PCR.

Results

By applying bioinformatic software and programs, miRNAs that target the 3′-UTR of both AXL and MET mRNAs were determined, and according to the scores, miR-34a was selected for further analyses. The expression level of miR-34a in MDA-MB-231 and SK-BR-3 was lower than that of MCF-10A. Furthermore, AXL and MET expression in SK-BR-3 and MDA-MB-231 was lower and higher, respectively, than that of MCF-10A. After miR-34a overexpression, MET and AXL were downregulated in MDA-MB-231. In addition, MET was downregulated in SK-BR-3, while AXL was upregulated in this cell line.

Conclusions

These findings may indicate that miR-34a is an oncogenic miRNA, downregulated in the distinct breast cancer subtypes. It also targets MET and AXL 3′-UTRs in triple-negative breast cancer. Therefore, it can be considered as a therapeutic target in this type of breast cancer.

     

Association analysis of ILVBL gene polymorphisms with aspirin-exacerbated respiratory disease in asthma

Background

We previously reported that the ILVBL gene on chromosome 19p13.1 was associated with the risk for aspirin-exacerbated respiratory disease (AERD) and the percent decline of forced expired volume in one second (FEV1) after an oral aspirin challenge test. In this study, we confirmed the association between polymorphisms and haplotypes of the ILVBL gene and the risk for AERD and its phenotype.

Methods

We recruited 141 AERD and 995 aspirin-tolerant asthmatic (ATA) subjects. All study subjects underwent an oral aspirin challenge (OAC). Nine single nucleotide polymorphisms (SNPs) with minor allele frequencies above 0.05, which were present in the region from 2 kb upstream to 0.5 kb downstream of ILVBL in Asian populations, were selected and genotyped.

Results

In an allelic association analysis, seven of nine SNPs were significantly associated with the risk for AERD after correction for multiple comparisons. In a codominant model, the five SNPs making up block2 (rs2240299, rs7507755, rs1468198, rs2074261, and rs13301) showed significant associations with the risk for AERD (corrected P?=?0.001–0.004, OR?=?0.59–0.64). Rs1468198 was also significantly associated with the percent decline in FEV1 in OAC tests after correction for multiple comparisons in the codominant model (corrected P?=?0.033), but the other four SNPs in hapblock2 were not.

Conclusion

To the best of our knowledge, this is the first report of an association between SNPs on ILVBL and AERD. SNPs on ILVBL could be promising genetic markers of this condition.

     

Role of Gene-Gene Interactions in the Chromosomal Instability in Workers at Coal Thermal Power Plants

The 23 polymorphic variants in genes encoding the enzymes of xenobiotics biotransformation (CYP1A1 (rs4646903), CYP1A2 (rs762551), GSTP1 (rs1138272, rs1695), GSTM1 (del), and GSTT1 (del)), DNA repair (XRCC1 (rs25489, rs25487), APEX1 (rs1130409), hOGG1 (rs1052133), ADPRT (rs1136410), XPD (rs13181), XPG (rs17655), XPC (rs2228001), ATM (rs1801516), NBS (rs1805794), XRCC2 (rs3218536), and XRCC3 (rs861539)), antioxidant system (MnSOD (rs4880) and GPx1 (rs1050450)), cell cycle control and apoptosis (TP53 (rs1042522)), DNA methylation (MTHFR (rs1801133) and MTR (rs1805087)), and chromosomal aberrations in lymphocytes in the workers at thermal power plants were analyzed. We found that allelic variants in the CYP1A1 (rs4646903), hOGG1 (rs1052133), XRCC1 (rs25487), and APEX1 (rs1130409) genes were associated with an increased level of chromosomal aberrations in workers. Informative models of gene-gene interactions including CYP1A1 (rs4646903, T>C), CYP1A2 (rs762551, C>A), GSTT1 (del); XRCC1 (rs25487, G>A), MTHFR (rs1801133, C>T), GSTT1 (del); XRCC1 (rs25487, G>A), APEX1 (rs1130409, T>G), TP53 (rs1042522, G>C) determining the formation of the increased frequency of chromosomal aberrations in the workers at coal thermal power plants were discovered.     

<Emphasis Type="Italic">EPSTI1</Emphasis> polymorphisms are associated with systemic lupus erythematosus

Epithelial stromal interaction 1 (EPSTI1) is an interferon (IFN) response gene, which was originally identified as a stromal fibroblast-induced gene in breast cancer. Our previous study using a customized SNP chip found evidence of an association between EPSTI1 and susceptibility to the chronic inflammatory disease, systematic lupus erythematosus (SLE). This study aimed to validate whether polymorphisms in EPSTI1 are associated with susceptibility to SLE. We analyzed genotype and allele frequencies of SNPs at EPSTI1 using genomic DNA from 119 patients with SLE and 512 healthy controls. We found that the genotype frequencies of rs1044856 and rs1359184 in patients with SLE were significantly different from those found in the control group (P?=?0.03 and P?=?0.01, respectively). In addition, we found that genotype and allele frequencies of rs1359184 in female patients with SLE were significantly different from those found in female controls (P?=?0.02 and P?=?0.04, respectively). We identified two major haplotypes in EPSTI1 that were significantly different between patients with SLE and healthy controls (P?=?0.01 and P?=?0.05, respectively). Furthermore, we found that rs1359184 and rs1044856 in EPSTI1 were associated with antinuclear antibody (ANA) and erythrocyte sedimentation rate (ESR) levels in patients with SLE (P?=?0.0035 and P?=?0.021, respectively). Our findings indicate that polymorphisms in EPSTI1 are associated with susceptibility to SLE and that haplotypes at EPSTI1 may be useful genetic markers for SLE.     

Association of polymorphisms in intron 2 of FGFR2 and breast cancer risk in Chinese women

Recent genome-wide association studies (GWAS) demonstrated that genetic variation in intron 2 of fibroblast growth factor receptor 2 (FGFR2) was a novel risk for breast cancer. We investigated whether two SNPs rs1219648 and rs2981582 in intron 2 of FGFR2 were associated with the risk of breast cancer in Chinese women. A total of 340 female breast cancer patients and 400 normal age-matched controls were recruited. Two SNPs were genotyped using matrix-assisted laser desorption/ionization mass spectrometry. The two SNPs rs1219648 and rs2981582 showed no association with the risk of breast cancer. A subgroup analysis by menopausal status demonstrated that the distribution of rs2981582 T alleles, including CT and TT genotypes, was significantly higher in premenopausal patients compared with postmenopausal patients. The TT genotype in rs2981582 was more strongly associated with ER-positive than with ER-negative tumors by ER status analysis. Analysis by haplotypes showed that no haplotypes associated with breast cancer. The results showed no association between two SNPs, rs1219648 and rs2981582 and breast cancer risk, although in a stratified analysis rs2981582 strongly associated with premenopausal and ER-positive breast cancer patients in Chinese women.     

Genetic Variations in MicroRNA Processing Genes Are Associated with Susceptibility in Depression

MicroRNAs (miRNAs) are new prominent gene expression regulators that have critical roles in neural development by regulating synaptic functions, and miRNA biogenesis may play an important role in psychiatric disorders. Despite emerging evidences demonstrating that single-nucleotide polymorphisms in the miRNA processing genes were associated with cancer and cardiovascular disorders, evidences about association between variants of the genes and depression are lacking. This study aims to find the association between miRNA processing gene variants and depression. We genotyped three polymorphisms from three miRNA processing genes in a case-control study including 314 patients and 252 matched healthy controls. The high-resolution melting method was used to genotype the three loci. Frequencies of genotypes and alleles showed significant difference between patients with depression and healthy controls in DGCR8 rs3757 and AGO1 rs636832. An allele frequency was significantly higher in rs3757 and lower in rs636832, respectively. Variant allele of DGCR8 rs3757 was associated with increased risk of suicidal tendency and improvement response to antidepressant treatment, whereas the variant of AGO1 rs636832 showed decreased risk of suicidal tendency, suicidal behavior, and recurrence. Besides allele frequency showed significant difference when compared patients with remission to controls, no significant differences were found in GEMIN4 rs7813 between patients and healthy controls. DGCR8 rs3757 and AGO1 rs636832 were found to have significant association with depression, and GEMIN4 rs7813 did not affect susceptibility to depression. These observations suggested that miRNA processing polymorphisms may affect depression risk and treatment.     

Association study of genetic markers of schizophrenia and its cognitive endophenotypes

A replicative analysis of associations of 15 SNPs located in the regions of 11 genes (TCF4, VRK2, NOTCH4, ZNF804A, AGBL1, RELN, ZFP64P1, KCNB2, CSMD1, CPVL, NRIP1) and three intergenic regions (SLCO6A1/LINCOO491, LOC105376248/LOC105376249, SPA17/NRGN) with schizophrenia was conducted in the Russian population of the Siberian region. These SNPs were previously identified in genome-wide association studies (GWAS) of schizophrenia and cognitive abnormalities. The present study confirmed associations of KCNB2 rs2247572, CSMD1 rs2616984, and intergenic rs12807809 located in SPA17/NRGN with schizophrenia. It was established that the frequency of the CSMD1 rs2616984 G/G genotype was higher in patients compared to the control group (OR = 1.73; CI: 1.14–2.62; р = 0.0337). The frequencies of the KCNB2 rs2247572 TT genotype (OR = 0.41; CI: 0.20–0.87; р = 0.0485) and intergenic rs12807809 CT genotype located in SPA17/NRGN (OR = 0.70; CI: 0.53–0.94; р = 0.0464) were significantly decreased in patients compared to the control group.     

Bayesian logistic regression in detection of gene–steroid interaction for cancer at <Emphasis Type="Italic">PDLIM5</Emphasis> locus

The PDZ and LIM domain 5 (PDLIM5) gene may play a role in cancer, bipolar disorder, major depression, alcohol dependence and schizophrenia; however, little is known about the interaction effect of steroid and PDLIM5 gene on cancer. This study examined 47 single-nucleotide polymorphisms (SNPs) within the PDLIM5 gene in the Marshfield sample with 716 cancer patients (any diagnosed cancer, excluding minor skin cancer) and 2848 noncancer controls. Multiple logistic regression model in PLINK software was used to examine the association of each SNP with cancer. Bayesian logistic regression in PROC GENMOD in SAS statistical software, ver. 9.4 was used to detect gene–steroid interactions influencing cancer. Single marker analysis using PLINK identified 12 SNPs associated with cancer (P < 0.05); especially, SNP rs6532496 revealed the strongest association with cancer (P = 6.84 × 10^?3); while the next best signal was rs951613 (P = 7.46 × 10^?3). Classic logistic regression in PROC GENMOD showed that both rs6532496 and rs951613 revealed strong gene–steroid interaction effects (OR = 2.18, 95% CI = 1.31?3.63 with P = 2.9 × 10^?3 for rs6532496 and OR = 2.07, 95% CI = 1.24 ?3.45 with P = 5.43 × 10^?3 for rs951613, respectively). Results from Bayesian logistic regression showed stronger interaction effects (OR = 2.26, 95% CI = 1.2 ?3.38 for rs6532496 and OR = 2.14, 95% CI = 1.14 ?3.2 for rs951613, respectively). All the 12 SNPs associated with cancer revealed significant gene–steroid interaction effects (P < 0.05); whereas 13 SNPs showed gene–steroid interaction effects without main effect on cancer. SNP rs4634230 revealed the strongest gene–steroid interaction effect (OR = 2.49, 95% CI = 1.5 ?4.13 with P = 4.0 × 10^?4 based on the classic logistic regression and OR = 2.59, 95% CI = 1.4 ?3.97 from Bayesian logistic regression; respectively). This study provides evidence of common genetic variants within the PDLIM5 gene and interactions between PLDIM5 gene polymorphisms and steroid use influencing cancer.     

Association of genes of different functional classes with type 1 diabetes

The genetic structure of susceptibility to type 1 diabetes (T1D) in the population of Tomsk was studied. We had a group of T1D patients (N = 285) and a population sample (N = 300) and we studied 58 SNPs localized in the 47 genes which products are involved in various metabolic pathways and processes as fibrogenesis, endothelial dysfunction, and inflammation. Genotyping was performed by mass spectrometry using the Sequenom MassARRAY system (United States). We compared the group of T1D patients and the population sample and found an association with the predisposition to disease for seven markers: rs3765124 of the ADAMDEC1 gene, genotype AA (p = 0.004), allele A (p = 0.033); rs1007856 of the ITGB5 gene, genotype TT (p = 0.015), allele T (p = 0.036); rs20579 of the LIG1 gene, genotype CC (p = 0.004), allele C (p = 0.002); rs12980602 of the IFNL2 gene, allele C (p = 0.029); rs4986819 of the PARP4 gene, allele C (p = 0.044); rs1143674 of the ITGA4 gene genotype GG (p = 0.002); rs679620 of the MMP3 gene, genotype AA (p = 0.008). Thus, the products of genes associated with T1D belong to different molecular classes: metalloproteases (ADAMDEC1, MMP3), cytokines (IL28A), cell surface receptors (ITGA4), adhesion molecules (ITGB5), DNA ligases (LIG1), and ribosyltransferase enzymes (PARP4). The ADAMDEC1, ITGA4, and ITGB5 genes belong to two biological processes: cell communication and signal transduction. The LIG1 and PARP4 genes regulate the metabolism of nucleic acids, MMP3 is involved in the regulation of protein metabolism, and the IFNL2 is involved in the immune response.     

Polymorphisms in the 3′-UTR of <Emphasis Type="Italic">SCD5</Emphasis> gene are associated with hepatocellular carcinoma in Korean population

The purpose of the study was to assess the relationship between polymorphisms of the SCD5 and MMP1 gene and hepatocellular carcinoma (HCC). The gene polymorphisms with a minor allele frequency (MAF)?>?0.05 were selected eight SNPs (rs6840, rs1065403, rs3821974, and rs3733230 in 3?-UTR; rs4693472, rs3733227, rs1848067, and rs6535374 in intron region) of SCD5 gene and two SNPs (rs1799750 and rs1144393 in promoter region) of MMP1 gene. The genotype of SCD5 and MMP1 gene SNPs were determined by direct sequencing and pyrosequencing, respectively. One hundred sixty-two patients with HCC and two hundred twenty-five control subjects were recruited in Korean male population. In terms of genotype frequencies, SCD5 genotype TC, GA, AG, and CG of rs6840, rs1065403, rs3821974, and rs3733230, respectively were higher in control group than HCC. In addition, these genotype decreased the risk (rs6840; OR 0.55, 95% CI 0.31–0.99; rs1065403; OR 0.46, 95% CI 0.26–0.83; rs3821974; OR 0.56, 95% CI 0.31–0.99; rs3733230; OR 0.62, 95% CI 0.34–1.12) of HCC, which may work as a prevention of HCC. At least one minor allele carrier of SCD5 gene polymorphisms were related to decreased risk of HCC for AFP cut-point levels >?200 or >?400 ng/ml, respectively. Our results indicate that polymorphisms in the 3?-UTR of the SCD5 gene may associated with HCC in the Korean male population.