A phenome-wide association study identifies effects of copy-number variation of VNTRs and multicopy genes on multiple human traits

The human genome contains tens of thousands of large tandem repeats and hundreds of genes that show common and highly variable copy-number changes. Due to their large size and repetitive nature, these variable number tandem repeats (VNTRs) and multicopy genes are generally recalcitrant to standard genotyping approaches and, as a result, this class of variation is poorly characterized. However, several recent studies have demonstrated that copy-number variation of VNTRs can modify local gene expression, epigenetics, and human traits, indicating that many have a functional role. Here, using read depth from whole-genome sequencing to profile copy number, we report results of a phenome-wide association study (PheWAS) of VNTRs and multicopy genes in a discovery cohort of ∼35,000 samples, identifying 32 traits associated with copy number of 38 VNTRs and multicopy genes at 1% FDR. We replicated many of these signals in an independent cohort and observed that VNTRs showing trait associations were significantly enriched for expression QTLs with nearby genes, providing strong support for our results. Fine-mapping studies indicated that in the majority (∼90%) of cases, the VNTRs and multicopy genes we identified represent the causal variants underlying the observed associations. Furthermore, several lie in regions where prior SNV-based GWASs have failed to identify any significant associations with these traits. Our study indicates that copy number of VNTRs and multicopy genes contributes to diverse human traits and suggests that complex structural variants potentially explain some of the so-called "missing heritability" of SNV-based GWASs.

Analysis of Interleukin-1 Receptor Antagonist Variable Number Tandem Repeat Variant in Recurrent Aphthous Stomatitis

BACKGROUND/AIMS

Recurrent aphthous stomatitis (RAS) is one of the common oral inflammatory diseases. As immunological and genetic factors have been held responsible for the pathogenesis of RAS, the objective of this study was to determine whether the interleukin-1 receptor antagonist (IL-1Ra) gene variable number tandem repeat (VNTR) variant is a risk factor for the development of RAS in Turkish patients and to define its contribution to the increased risk.

METHODS

The IL-1Ra VNTR variant was evaluated in 169 RAS patients and 171 healthy controls by the polymerase chain reaction (PCR) method.

RESULTS

No statistically significant difference was found in the genotype distributions and allele frequencies of IL-1Ra VNTR variant between RAS patients and healthy controls.

CONCLUSION

Lack of association between IL-1Ra VNTR variant and RAS could indicate that IL-1Ra has no significant role in the pathophysiology of RAS. However, it still appears to be very worthwhile to continue to search for cytokine gene variants in order to predict the development of such disease.

A PERIOD3 variable number tandem repeat polymorphism modulates melatonin treatment response in delayed sleep-wake phase disorder

We examined whether a polymorphism of the PERIOD3 gene (PER3; rs57875989) modulated the sleep-promoting effects of melatonin in Delayed Sleep-Wake Phase Disorder (DSWPD). One hundred and four individuals (53 males; 29.4 ±10.0 years) with DSWPD and a delayed dim light melatonin onset (DLMO) collected buccal swabs for genotyping (PER3^4/4 n = 43; PER3 5 allele [heterozygous and homozygous] n = 60). Participants were randomised to placebo or 0.5 mg melatonin taken 1 hour before desired bedtime (or ~1.45 hours before DLMO), with sleep attempted at desired bedtime (4 weeks; 5-7 nights/week). We assessed sleep (diary and actigraphy), Pittsburgh Sleep Quality Index (PSQI), Insomnia Severity Index (ISI), Patient-Reported Outcomes Measurement Information System (PROMIS: Sleep Disturbance, Sleep-Related Impairment), Sheehan Disability Scale (SDS) and Patient- and Clinician-Global Improvement (PGI-C, CGI-C). Melatonin treatment response on actigraphic sleep onset time did not differ between genotypes. For PER3^4/4 carriers, self-reported sleep onset time was advanced by a larger amount and sleep onset latency (SOL) was shorter in melatonin-treated patients compared to those receiving placebo (P = .008), while actigraphic sleep efficiency in the first third of the sleep episode (SE T1) did not differ. For PER3 5 carriers, actigraphic SOL and SE T1 showed a larger improvement with melatonin (P < .001). Melatonin improved ISI (P = .005), PROMIS sleep disturbance (P < .001) and sleep-related impairment (P = .017), SDS (P = .019), PGI-C (P = .028) and CGI-C (P = .016) in PER3^4/4 individuals only. Melatonin did not advance circadian phase. Overall, PER3^4/4 DSWPD patients have a greater response to melatonin treatment. PER3 genotyping may therefore improve DSWPD patient outcomes.

Enhancer Function of MicroRNA-3681 Derived from Long Terminal Repeats Represses the Activity of Variable Number Tandem Repeats in the 3' UTR of SHISA7

microRNAs (miRNAs) are non-coding RNA molecules involved in the regulation of gene expression. miRNAs inhibit gene expression by binding to the 3' untranslated region (UTR) of their target gene. miRNAs can originate from transposable elements (TEs), which comprise approximately half of the eukaryotic genome and one type of TE, called the long terminal repeat (LTR) is found in class of retrotransposons. Amongst the miRNAs derived from LTR, hsa-miR-3681 was chosen and analyzed using bioinformatics tools and experimental analysis. Studies on hsa-miR-3681 have been scarce and this study provides the relative expression analysis of hsa-miR-3681-5p from humans, chimpanzees, crab-eating monkeys, and mice. Luciferase assay for hsa-miR-3681-5p and its target gene SHISA7 supports our hypothesis that the number of miRNA binding sites affects target gene expression. Especially, the variable number tandem repeat (VNTR) and hsa-miR-3681-5p share the binding sites in the 3' UTR of SHISA7, which leads the enhancer function of hsa-miR-3681-5p to inhibit the activity of VNTR. In conclusion, hsa-miR-3681-5p acts as a super-enhancer and the enhancer function of hsa-miR-3681-5p acts as a repressor of VNTR activity in the 3' UTR of SHISA7.

Investigation of the relationship between ischemic stroke and endothelial nitric oxide synthase gene polymorphisms [G894T, intron 4 VNTR and T786C]

Background/aim

We aimed to investigate the associations between endothelial nitric oxide synthase (eNOS) gene polymorphisms [G894T (rs1799983)], intron 4 (27-bpTR) variable number tandem repeat (VNTR) and T786C (rs2070744), and ischemic stroke in the Anatolian population.

Materials and methods

This case-control study included 112 patients with “stroke of undetermined etiology” and 160 controls. Real-time polymerase chain reaction (RT-PCR) analysis was used to analyze these polymorphisms. Between-group frequencies of alleles and genotypes were compared using binary logistic regression analysis.

Results

No significant difference was observed between the two groups in terms of the genotype and allele distributions of the eNOS G894T (rs1799983) polymorphism (P > 0.05). The a alleles and the 4b/a and 4a/a genotypes of the intron 4 (27-bpTR) VNTR polymorphism had significantly higher frequencies in the patient group than in the control group (OR: 2.715, P < 0.001; OR: 3.396, P < 0.001; OR: 10.631, P = 0.016, respectively). On the contrary, the TC genotype and C alleles of the T786C (rs2070744) polymorphism had a significantly lower frequency in the patient group than in the control group (OR: 0.244, P < 0.001, OR: 0.605, P = 0.006, respectively).

Conclusion

Our findings indicate that the eNOS G894T and T786C [rs2070744] polymorphisms are not associated with the risk of ischemic stroke, whereas the intron 4 [27-bpTR] VNTR may be a risk factor in the Anatolian population.

Lack of association between interleukin-1 receptor antagonist gene 86-bp VNTR polymorphism and ischemic stroke: A meta-analysis

OBJECTIVE

The results of published studies which examined the association between variable number tandem repeat (VNTR) polymorphism of interleukin-1 receptor antagonist (IL-1RN) and ischemic stroke (IS) are conflicting. Thus, we performed a meta-analysis to examine the potential association between IL-1RN VNTR polymorphism and IS risk.

METHODS

A systematic literature search of PubMed, Embase, Medline, Web of Science, Cochrane Library, Chinese Biomedical Literature Database, Chinese National Knowledge Infrastructure, CQVIP, and WANFANG Database identified 10 studies with 2331 cases and 3335 controls. The pooled odds ratio (OR) with 95% confidence interval (95% CI) was calculated to evaluate the strength of the association. Subgroup analysis and meta-regression analysis were used to investigate the potential sources of heterogeneity. Begg funnel plots were used to explore the publication bias.

RESULTS

In this study, no enough proof was found to prove the association between IL-1RN 86-bp VNTR polymorphism and IS risk with random-effects model in the homozygous model (1/1 vs 2/2, OR = 0.97, 95% CI = 0.50-1.87, Pheterogeneity = .00), the heterozygous model (1/2 vs 2/2, OR = 0.64, 95% CI = 0.41-1.01, Pheterogeneity = .10), the dominant model (1/1 + 1/2 vs 2/2, OR = 0.85, 95% CI = 0.51-1.42, Pheterogeneity = .02), the recessive model (1/1 vs 1/2 + 2/2, OR = 0.69, 95% CI = 0.46-1.03, Pheterogeneity = .00), and allelic model (1 vs 2, OR = 1.24, 95% CI = 0.89-1.74, Pheterogeneity = .00). A marginally significant negative association was observed between IL-1RN 86-bp VNTR polymorphism and IS risk in the heterozygous model in the fixed-effects model (1/2 vs 2/2, OR = 0.71, 95% CI = 0.53-0.95, Pheterogeneity = .10). In subgroup analyses, similar association was found in the group whose control size was lower than 300.

CONCLUSION

In conclusion, our results suggested that there was no sufficient evidence to support the association between IL-1RN 86-bp VNTR polymorphism and IS. Further large epidemiologic studies need to be done to confirm these findings.

Investigation of the role of interleukin-1 receptor antagonist VNTR variant on the Behçet's disease

Objective

Behçet's disease (BD), a chronic multisystem inflammatory disorder, is mainly characterized by relapsing periods of a wide range of clinical symptoms. Several cytokine genes may play important roles in the pathogenesis of BD. Therefore, interleukin-1 receptor antagonist (IL-1Ra) gene 86bp variable number tandem repeat (VNTR) variant was investigated in patients with BD in a Turkish population.

Methods

One hundred nine patients (60 females, 49 males; the mean age±standard deviation [SD] was 36.56±9.571 years) with BD and one hundred healthy individuals (54 females, 46 males; the mean age±SD was 36.64±2.294 years) were examined in the study. For genotyping, polymerase chain reaction-restriction fragment length polymorphism analysis was employed. Data were analyzed using Statistical Package for Social Sciences (SPSS) 22.0 (IBM Corp.; Armonk, NY, USA) (p<0.05).

Results

The genotype distribution and allele frequencies of the IL-1Ra VNTR variant did not differ significantly between the patients and the controls (p>0.05). The frequency of the a1/a1, a1/a2 genotypes and a1, a2 alleles were the most common both in patients and healthy controls (p=0.37, p=0.26, and p=0.53, respectively). Also, no statistically significant difference was found between the IL-1Ra VNTR variant genotypes and clinical characteristics (p>0.05).

Conclusion

The results of this study do not support an association between the IL-1Ra VNTR variant and the risk of BD in a Turkish population. However, further studies of this variant with larger sample sizes and different ethnicities are required for confirmation.

Association of TNFRSF1B Promoter Polymorphisms with Human Disease: Further Studies Examining T-Regulatory Cells Are Required

The TNFR2 receptor is expressed by highly active regulatory T cells, and thus constitutes an important therapeutic target for the treatment of autoimmune disease and cancer. Disease susceptibility as well as the potential response to therapies directed at TNFR2 could be significantly impacted by genetic variation in the promoter of the TNFRSF1B gene that codes for the TNFR2 protein. To date, only a few studies have examined the association of TNFRSF1B promoter variation with disease, and the potential impact on T-regulatory cell (Treg) number and function has not been examined. We propose that copy number variation of a key transcription factor binding site has a significant effect on TNFRSF1B promoter activity, and should be considered in studies of disease susceptibility and especially with regard to variation in the level of TNFR2 expression on Tregs.

Characterization of a REST-Regulated Internal Promoter in the Schizophrenia Genome-Wide Associated Gene MIR137

MIR137 has been identified as a candidate gene for schizophrenia from genome-wide association studies via association with an intronic single nucleotide polymorphism (SNP), rs1625579. The location of the SNP suggests one mechanism in which transcriptional or posttranscriptional regulation of miR-137 expression could underlie schizophrenia. We identified and validated a novel promoter of the MIR137 gene adjacent to miR-137 itself which can direct the expression of distinct mRNA isoforms encoding miR-137. Analysis of both endogenous gene expression and reporter gene assays determined that this internal promoter is regulated by repressor element-1 silencing transcription factor (REST), which has previously been associated with pathways linked to schizophrenia. Distinct isoforms of REST mediate differential expression at this locus, suggesting the relative levels of these isoforms are important for miR-137 expression profiles. The internal promoter contains a variable number tandem repeat (VNTR) domain adjacent to the pre-miR-137 sequence. The reporter gene activity directed by this promoter was modified by the genotype of the VNTR. Differential expression was also observed in response to cocaine, which is known to regulate the REST pathway in SH-SY5Y cells. Our data support the hypothesis that a "gene × environment" interaction could modify the level of miR-137 expression via this internal promoter and that the genotype of the VNTR could modulate transcriptional responses. We demonstrate that this promoter region is not in disequilibrium with rs1625579 and therefore would supply a distinct pathway to potentially alter miR-137 levels in response to environmental cues.

Case-control study of the PERIOD3 clock gene length polymorphism and colorectal adenoma formation

Clock genes are expressed in a self-perpetuating, circadian pattern in virtually every tissue including the human gastrointestinal tract. They coordinate cellular processes critical for tumor development, including cell proliferation, DNA damage response and apoptosis. Circadian rhythm disturbances have been associated with an increased risk for colon cancer and other cancers. This mechanism has not been elucidated, yet may involve dysregulation of the 'period' (PER) clock genes, which have tumor suppressor properties. A variable number tandem repeat (VNTR) in the PERIOD3 (PER3) gene has been associated with sleep disorders, differences in diurnal hormone secretion, and increased premenopausal breast cancer risk. Susceptibility related to PER3 has not been examined in conjunction with adenomatous polyps. This exploratory case-control study was the first to test the hypothesis that the 5-repeat PER3 VNTR sequence is associated with increased odds of adenoma formation. Information on demographics, medical history, occupation and lifestyle was collected prior to colonoscopy. Cases (n=49) were individuals with at least one histopathologically confirmed adenoma. Controls (n=97) included patients with normal findings or hyperplastic polyps not requiring enhanced surveillance. Unconditional multiple logistic regression was used to calculate odds ratios (ORs) with 95% confidence intervals (CIs), after adjusting for potential confounding. Adenomas were detected in 34% of participants. Cases were more likely to possess the 5-repeat PER3 genotype relative to controls (4/5 OR, 2.1; 95% CI, 0.9-4.8; 5/5 OR, 5.1; 95% CI, 1.4-18.1; 4/5+5/5 OR, 2.5; 95% CI, 1.7-5.4). Examination of the Oncomine microarray database indicated lower PERIOD gene expression in adenomas relative to adjacent normal tissue. Results suggest a need for follow-up in a larger sample.

Isolation of Leptospira from a phocid: acute renal failure and mortality from Leptospirosis in rehabilitated northern elephant seals (Mirounga angustirostris), California, USA

During rehabilitation, acute renal failure due to leptospirosis occurred in eight male northern elephant seals (Mirounga angustirostris) that stranded along the central California coast in 2011. Characteristic histologic lesions including renal tubular degeneration, necrosis, and mineralization, and mild lymphoplasmacytic interstitial nephritis were noted in the six animals examined. Immunohistochemistry, bacterial culture, and PCR were positive in 2/3, 2/3, and 3/4 seals, respectively, and 6/8 had high serum antibody titers to Leptospira interrogans serovar pomona. Pulsed-field gel electrophoresis confirmed one isolate as serovar pomona. Variable number tandem repeat (VNTR) analysis showed both elephant seal isolates were identical to each other but distinct from those isolated from California sea lions (Zalophus californianus). The time from stranding to onset of azotemia was 1 to 38 (median=24) days, suggesting some seals were infected at the rehabilitation facility. Based on temporal and spatial incidence of infection, transmission among elephant seals likely occurred during rehabilitation. Molecular (VNTR) analysis of the two isolates indicates there is a unique L. interrogans serovar pomona genotype in elephant seals, and sea lions were not the source of infection prior to or during rehabilitation. This study confirms the susceptibility of northern elephant seals to leptospirosis, indicates intraspecies transmission during rehabilitation, and reports the first isolation and preliminary characterization of leptospires from elephant seals.

Significant association of interleukin-4 gene intron 3 VNTR polymorphism with susceptibility to knee osteoarthritis

OBJECTIVE

Interleukin-4 (IL-4) is a strong chondroprotective cytokine and polymorphisms within this gene may be a risk factor for osteoarthritis (OA). We aimed to investigate genotype and allele frequencies of IL-4 gene intron 3 variable number of tandem repeats (VNTR) polymorphism in patients with knee OA in a Turkish population.

METHODS

The study included 202 patients with knee OA and 180 healthy controls. Genomic DNA was isolated and IL-4 gene 70 bp VNTR polymorphism determined by using polymerase chain reaction (PCR) with specific primers followed by restriction fragment length polymorphism (RFLP) analysis.

RESULTS

Our result show that there was statistically significant difference between knee OA patients and control group with respect to IL-4 genotype distribution and allele frequencies (p=0.000, OR: 0.20, 95% CI: 0.10-0.41, OR: 0.22, 95% CI: 0.12-0.42, respectively).

CONCLUSIONS

Our findings suggest that there is an association of IL-4 gene intron 3 VNTR polymorphism with susceptibility of a person for development of knee OA. As a result, IL-4 gene intron 3 VNTR polymorphism could be a genetic marker in OA in a Turkish study population. This is the first association study that evaluates the associations between IL-4 gene VNTR polymorphism and knee OA.

Do dopamine gene variants and prenatal smoking interactively predict youth externalizing behavior?

Externalizing behaviors (encompassing antisocial, impulsive, and substance use behaviors) are pervasive and impairing across a multitude of settings and developmental contexts. These behaviors, though often investigated separately, are highly comorbid. Prenatal tobacco exposure in interaction with various genetic influences has predicted later externalizing behavior, and recent evidence supports investigating sex differences in these patterns. In the current study, we extend this work by (a) examining two functional genetic markers in the dopamine system: the transporter gene (DAT1) and the dopamine receptor D4 gene (DRD4) in interaction with prenatal tobacco exposure to predict a latent composite of externalizing behavior and (b) testing whether these patterns differ by sex of youth in a community sample of adolescents (n=176). The relatively small sample is partially offset by high quality, multi-method prospective measurement. We assessed prenatal tobacco exposure using prospective repeated cotinine-corrected reports and externalizing behaviors were assessed utilizing multiple measures across three waves. The interaction between DAT1 (but not DRD4) and prenatal tobacco exposure was statistically significant in boys, and patterns appeared to differ by sex. Risk for externalizing behaviors for exposed boys increased linearly as a function of the 10r DAT1 allele. For exposed girls, there was a trend such that DAT1 heterozygotes had a marginally higher risk than homozygotes. This pattern was not explained by passive gene-environment correlation. Elucidating sex-specific pathways through which early adverse exposures and genetic susceptibilities contribute to externalizing behavior can inform early targeted prevention efforts for those children at highest risk.

Identity tests: determination of cell line cross-contamination

Cell line cross-contamination is a phenomenon that arises as a result of the continuous cell line culture. It has been estimated that around 20% of the cell lines are misidentified, therefore it is necessary to carry out quality control tests for the detection of this issue. Since cell line cross-contamination discovery, different methods have been applied, such as isoenzyme analysis for inter-species cross-contamination; HLA typing, and DNA fingerprinting using short tandem repeat and a variable number of tandem repeat for intra-species cross-contamination. The cell banks in this sense represent the organizations responsible for guaranteeing the authenticity of cell lines for future research and clinical uses.