Virus genome dynamics under different propagation pressures: reconstruction of whole genome haplotypes of west nile viruses from NGS data

Background

Extensive focus is placed on the comparative analyses of consensus genotypes in the study of West Nile virus (WNV) emergence. Few studies account for genetic change in the underlying WNV quasispecies population variants. These variants are not discernable in the consensus genome at the time of emergence, and the maintenance of mutation-selection equilibria of population variants is greatly underestimated. The emergence of lineage 1 WNV strains has been studied extensively, but recent epidemics caused by lineage 2 WNV strains in Hungary, Austria, Greece and Italy emphasizes the increasing importance of this lineage to public health. In this study we explored the quasispecies dynamics of minority variants that contribute to cell-tropism and host determination, i.e. the ability to infect different cell types or cells from different species from Next Generation Sequencing (NGS) data of a historic lineage 2 WNV strain.

Results

Minority variants contributing to host cell membrane association persist in the viral population without contributing to the genetic change in the consensus genome. Minority variants are shown to maintain a stable mutation-selection equilibrium under positive selection, particularly in the capsid gene region.

Conclusions

This study is the first to infer positive selection and the persistence of WNV haplotype variants that contribute to viral fitness without accompanying genetic change in the consensus genotype, documented solely from NGS sequence data. The approach used in this study streamlines the experimental design seeking viral minority variants accurately from NGS data whilst minimizing the influence of associated sequence error.     

Common variants in the type 2 diabetes KCNQ1 gene are associated with impairments in insulin secretion during hyperglycaemic glucose clamp

Background

Genome-wide association studies in Japanese populations recently identified common variants in the KCNQ1 gene to be associated with type 2 diabetes. We examined the association of these variants within KCNQ1 with type 2 diabetes in a Dutch population, investigated their effects on insulin secretion and metabolic traits and on the risk of developing complications in type 2 diabetes patients.

Methodology

The KCNQ1 variants rs151290, rs2237892, and rs2237895 were genotyped in a total of 4620 type 2 diabetes patients and 5285 healthy controls from the Netherlands. Data on macrovascular complications, nephropathy and retinopathy were available in a subset of diabetic patients. Association between genotype and insulin secretion/action was assessed in the additional sample of 335 individuals who underwent a hyperglycaemic clamp.

Principal Findings

We found that all the genotyped KCNQ1 variants were significantly associated with type 2 diabetes in our Dutch population, and the association of rs151290 was the strongest (OR 1.20, 95% CI 1.07–1.35, p = 0.002). The risk C-allele of rs151290 was nominally associated with reduced first-phase glucose-stimulated insulin secretion, while the non-risk T-allele of rs2237892 was significantly correlated with increased second-phase glucose-stimulated insulin secretion (p = 0.025 and 0.0016, respectively). In addition, the risk C-allele of rs2237892 was associated with higher LDL and total cholesterol levels (p = 0.015 and 0.003, respectively). We found no evidence for an association of KCNQ1 with diabetic complications.

Conclusions

Common variants in the KCNQ1 gene are associated with type 2 diabetes in a Dutch population, which can be explained at least in part by an effect on insulin secretion. Furthermore, our data suggest that KCNQ1 is also associated with lipid metabolism.     

Next Generation Sequencing Uncovers Unexpected Bacterial Pathogens in Ticks in Western Europe

Background and Aims

Ticks are highly susceptible to global environmental and socio-economical changes. Several tick-borne pathogens have been reported in new geographical regions while new species, strains or genetic variants of tick-borne microorganisms are continually being detected. However, tick-borne pathogens are still poorly understood, and it is estimated that half of all human tick-borne disease has an unknown origin. Therefore in order to prevent these diseases, more effort is required to identify unknown or unexpected tick-borne pathogens. Ixodes ricinus is the vector for a broad range of bacterial pathogens and the most prevalent tick in Europe. The aim of the present study was to evaluate the capability of Next Generation Sequencing (NGS) to extend the inventory of pathogenic bacteria carried by this species of tick in France.

Methods

RNA and DNA were extracted from 1450 I. ricinus questing nymphs collected by flagging in Alsace, France. RNA was pooled and used for NGS. Following de novo assembly, bacterial contigs were assigned to the closest known taxonomy. DNA was used for real time PCR to confirm taxonomic species assignment of NGS-derived contigs for the doubtful cases, and for determination of prevalence.

Results

We have generated a global in-depth picture of tick-borne bacteria. We identified RNA from the main pathogenic bacterial species known to be transmitted by I. ricinus. In addition we also identified unanticipated bacterial species for which we have estimated the prevalence within those ticks inhabiting the studied areas.

Conclusions

The data obtained from this study has proven that NGS has an enormous potential to detect the unexpected and provides the means to monitor pathogen occurrence.     

Contribution of FKBP5 Genetic Variation to Gemcitabine Treatment and Survival in Pancreatic Adenocarcinoma

Purpose

FKBP51, (FKBP5), is a negative regulator of Akt. Variability in FKBP5 expression level is a major factor contributing to variation in response to chemotherapeutic agents including gemcitabine, a first line treatment for pancreatic cancer. Genetic variation in FKBP5 could influence its function and, ultimately, treatment response of pancreatic cancer.

Experimental Design

We set out to comprehensively study the role of genetic variation in FKBP5 identified by Next Generation DNA resequencing on response to gemcitabine treatment of pancreatic cancer by utilizing both tumor and germline DNA samples from 43 pancreatic cancer patients, including 19 paired normal-tumor samples. Next, genotype-phenotype association studies were performed with overall survival as well as with FKBP5 gene expression in tumor using the same samples in which resequencing had been performed, followed by functional genomics studies.

Results

In-depth resequencing identified 404 FKBP5 single nucleotide polymorphisms (SNPs) in normal and tumor DNA. SNPs with the strongest associations with survival or FKBP5 expression were subjected to functional genomic study. Electromobility shift assay showed that the rs73748206 “A(T)” SNP altered DNA-protein binding patterns, consistent with significantly increased reporter gene activity, possibly through its increased binding to Glucocorticoid Receptor (GR). The effect of rs73748206 was confirmed on the basis of its association with FKBP5 expression by affecting the binding to GR in lymphoblastoid cell lines derived from the same patients for whom DNA was used for resequencing.

Conclusion

This comprehensive FKBP5 resequencing study provides insights into the role of genetic variation in variation of gemcitabine response.     

Genetic Variation of the Fc Gamma Receptor 3B Gene and Association with Rheumatoid Arthritis

Background

Fc gamma receptors (FcγRs) play a crucial role in immunity by linking IgG antibody-mediated responses with cellular effector and regulatory functions. Genetic variants in these receptors have been previously identified as risk factors for several chronic inflammatory conditions. The present study aimed to investigate the presence of copy number variations (CNVs) in the FCGR3B gene and its potential association with the autoimmune disease rheumatoid arthritis (RA).

Methodology/Principal Findings

CNV of the FCGR3B gene was studied using Multiplex Ligation Dependent Probe Amplification (MLPA) in 518 Dutch RA patients and 304 healthy controls. Surprisingly, three independent MLPA probes targeting the FCGR3B promoter measured different CNV frequencies, with probe#1 and #2 measuring 0 to 5 gene copies and probe#3 showing little evidence of CNV. Quantitative-PCR correlated with the copy number results from MLPA probe#2, which detected low copy number (1 copy) in 6.7% and high copy number (≥3 copies) in 9.4% of the control population. No significant difference was observed between RA patients and the healthy controls, neither in the low copy nor the high copy number groups (p-values = 0.36 and 0.71, respectively). Sequencing of the FCGR3B promoter region revealed an insertion/deletion (indel) that explained the disparate CNV results of MLPA probe#1. Finally, a non-significant trend was found between the novel -256A>TG indel and RA (40.7% in healthy controls versus 35.9% in RA patients; P = 0.08).

Conclusions/Significance

The current study highlights the complexity and poor characterization of the FCGR3B gene sequence, indicating that the design and interpretation of genotyping assays based on specific probe sequences must be performed with caution. Nonetheless, we confirmed the presence of CNV and identified novel polymorphisms in the FCGR3B gene in the Dutch population. Although no association was found between RA and FCGR3B CNV, the possible protective effect of the -256A>TG indel polymorphism must be addressed in larger studies.     

Genome-Wide Association Study Identified CNP12587 Region Underlying Height Variation in Chinese Females

Introduction

Human height is a highly heritable trait considered as an important factor for health. There has been limited success in identifying the genetic factors underlying height variation. We aim to identify sequence variants associated with adult height by a genome-wide association study of copy number variants (CNVs) in Chinese.

Methods

Genome-wide CNV association analyses were conducted in 1,625 unrelated Chinese adults and sex specific subgroup for height variation, respectively. Height was measured with a stadiometer. Affymetrix SNP6.0 genotyping platform was used to identify copy number polymorphisms (CNPs). We constructed a genomic map containing 1,009 CNPs in Chinese individuals and performed a genome-wide association study of CNPs with height.

Results

We detected 10 significant association signals for height (p<0.05) in the whole population, 9 and 11 association signals for Chinese female and male population, respectively. A copy number polymorphism ({"type":"entrez-protein","attrs":{"text":"CNP12587","term_id":"893597589","term_text":"CNP12587"}}CNP12587, chr18:54081842-54086942, p = 2.41×10⁻⁴) was found to be significantly associated with height variation in Chinese females even after strict Bonferroni correction (p = 0.048). Confirmatory real time PCR experiments lent further support for CNV validation. Compared to female subjects with two copies of the CNP, carriers of three copies had an average of 8.1% decrease in height. An important candidate gene, ubiquitin-protein ligase NEDD4-like (NEDD4L), was detected at this region, which plays important roles in bone metabolism by binding to bone formation regulators.

Conclusions

Our findings suggest the important genetic variants underlying height variation in Chinese.     

Integrative Analysis Reveals Relationships of Genetic and Epigenetic Alterations in Osteosarcoma

Background

Osteosarcomas are the most common non-haematological primary malignant tumours of bone, and all conventional osteosarcomas are high-grade tumours showing complex genomic aberrations. We have integrated genome-wide genetic and epigenetic profiles from the EuroBoNeT panel of 19 human osteosarcoma cell lines based on microarray technologies.

Principal Findings

The cell lines showed complex patterns of DNA copy number changes, where genomic copy number gains were significantly associated with gene-rich regions and losses with gene-poor regions. By integrating the datasets, 350 genes were identified as having two types of aberrations (gain/over-expression, hypo-methylation/over-expression, loss/under-expression or hyper-methylation/under-expression) using a recurrence threshold of 6/19 (>30%) cell lines. The genes showed in general alterations in either DNA copy number or DNA methylation, both within individual samples and across the sample panel. These 350 genes are involved in embryonic skeletal system development and morphogenesis, as well as remodelling of extracellular matrix. The aberrations of three selected genes, CXCL5, DLX5 and RUNX2, were validated in five cell lines and five tumour samples using PCR techniques. Several genes were hyper-methylated and under-expressed compared to normal osteoblasts, and expression could be reactivated by demethylation using 5-Aza-2′-deoxycytidine treatment for four genes tested; AKAP12, CXCL5, EFEMP1 and IL11RA. Globally, there was as expected a significant positive association between gain and over-expression, loss and under-expression as well as hyper-methylation and under-expression, but gain was also associated with hyper-methylation and under-expression, suggesting that hyper-methylation may oppose the effects of increased copy number for detrimental genes.

Conclusions

Integrative analysis of genome-wide genetic and epigenetic alterations identified dependencies and relationships between DNA copy number, DNA methylation and mRNA expression in osteosarcomas, contributing to better understanding of osteosarcoma biology.     

Genetic variants of diabetes risk and incident cardiovascular events in chronic coronary artery disease

Objective

To determine whether information from genetic risk variants for diabetes is associated with cardiovascular events incidence.

Methods

From the about 30 known genes associated with diabetes, we genotyped single-nucleotide polymorphisms at the 10 loci most associated with type-2 diabetes in 425 subjects from the MASS-II Study, a randomized study in patients with multi-vessel coronary artery disease. The combined genetic information was evaluated by number of risk alleles for diabetes. Performance of genetic models relative to major cardiovascular events incidence was analyzed through Kaplan-Meier curve comparison and Cox Hazard Models and the discriminatory ability of models was assessed for cardiovascular events by calculating the area under the ROC curve.

Results

Genetic information was able to predict 5-year incidence of major cardiovascular events and overall-mortality in non-diabetic individuals, even after adjustment for potential confounders including fasting glycemia. Non-diabetic individuals with high genetic risk had a similar incidence of events then diabetic individuals (cumulative hazard of 33.0 versus 35.1% of diabetic subjects). The addition of combined genetic information to clinical predictors significantly improved the AUC for cardiovascular events incidence (AUC = 0.641 versus 0.610).

Conclusions

Combined information of genetic variants for diabetes risk is associated to major cardiovascular events incidence, including overall mortality, in non-diabetic individuals with coronary artery disease.

Clinical Trial Registration Information

Medicine, Angioplasty, or Surgery Study (MASS II). Unique identifier: ISRCTN66068876 URL.     

Structural,Genetic, and Functional Signatures of Disordered Neuro-Immunological Development in Autism Spectrum Disorder

Background

Numerous linkage studies have been performed in pedigrees of Autism Spectrum Disorders, and these studies point to diverse loci and etiologies of autism in different pedigrees. The underlying pattern may be identified by an integrative approach, especially since ASD is a complex disorder manifested through many loci.

Method

Autism spectrum disorder (ASD) was studied through two different and independent genome-scale measurement modalities. We analyzed the results of copy number variation in autism and triangulated these with linkage studies.

Results

Consistently across both genome-scale measurements, the same two molecular themes emerged: immune/chemokine pathways and developmental pathways.

Conclusion

Linkage studies in aggregate do indeed share a thematic consistency, one which structural analyses recapitulate with high significance. These results also show for the first time that genomic profiling of pathways using a recombination distance metric can capture pathways that are consistent with those obtained from copy number variations (CNV).     

Evaluation of Five Candidate Genes from GWAS for Association with Oligozoospermia in a Han Chinese Population

Background

Oligozoospermia is one of the severe forms of idiopathic male infertility. However, its pathology is largely unknown, and few genetic factors have been defined. Our previous genome-wide association study (GWAS) has identified four risk loci for non-obstructive azoospermia (NOA).

Objective

To investigate the potentially functional genetic variants (including not only common variants, but also less-common and rare variants) of these loci on spermatogenic impairment, especially oligozoospermia.

Design, Setting, and Participants

A total of 784 individuals with oligozoospermia and 592 healthy controls were recruited to this study from March 2004 and January 2011.

Measurements

We conducted a two-stage study to explore the association between oligozoospermia and new makers near NOA risk loci. In the first stage, we used next generation sequencing (NGS) in 96 oligozoospermia cases and 96 healthy controls to screen oligozoospermia-susceptible genetic variants. Next, we validated these variants in a large cohort containing 688 cases and 496 controls by SNPscan for high-throughput Single Nucleotide Polymorphism (SNP) genotyping.

Results and Limitations

Totally, we observed seven oligozoospermia associated variants (rs3791185 and rs2232015 in PRMT6, rs146039840 and rs11046992 in Sox5, rs1129332 in PEX10, rs3197744 in SIRPA, rs1048055 in SIRPG) in the first stage. In the validation stage, rs3197744 in SIRPA and rs11046992 in Sox5 were associated with increased risk of oligozoospermia with an odds ratio (OR) of 4.62 (P  =  0.005, 95%CI 1.58-13.4) and 1.82 (P  =  0.005, 95%CI 1.01-1.64), respectively. Further investigation in larger populations and functional characterizations are needed to validate our findings.

Conclusions

Our study provides evidence of independent oligozoospermia risk alleles driven by variants in the potentially functional regions of genes discovered by GWAS. Our findings suggest that integrating sequence data with large-scale genotyping will serve as an effective strategy for discovering risk alleles in the future.     

Genetic variants in toll-like receptors are not associated with rheumatoid arthritis susceptibility or anti-tumour necrosis factor treatment outcome

Background

Several studies point to a role of Toll-like receptors (TLRs) in the development of rheumatoid arthritis (RA). We investigated if genetic variants in TLR genes are associated with RA and response to tumour necrosis factor blocking (anti-TNF) medication.

Methodology and Principal Findings

22 single nucleotide polymorphisms (SNPs) in seven TLR genes were genotyped in a Dutch cohort consisting of 378 RA patients and 294 controls. Significantly associated variants were investigated in replication cohorts from The Netherlands, United Kingdom and Sweden (2877 RA patients and 2025 controls). 182 of the Dutch patients were treated with anti-TNF medication. Using these patients and a replication cohort (269 Swedish patients) we analysed if genetic variants in TLR genes were associated with anti-TNF outcome. In the discovery phase of the study we found a significant association of SNPs rs2072493 in TLR5 and rs3853839 in TLR7 with RA disease susceptibility. Meta-analysis of discovery and replication cohorts did not confirm these findings. SNP rs2072493 in TLR5 was associated with anti-TNF outcome in the Dutch but not in the Swedish population.

Conclusion

We conclude that genetic variants in TLRs do not play a major role in susceptibility for developing RA nor in anti-TNF treatment outcome in a Caucasian population.     

Reference genome of wild goat (capra aegagrus) and sequencing of goat breeds provide insight into genic basis of goat domestication

Background

Domestic goats (Capra hircus) have been selected to play an essential role in agricultural production systems, since being domesticated from their wild progenitor, bezoar (Capra aegagrus). A detailed understanding of the genetic consequences imparted by the domestication process remains a key goal of evolutionary genomics.

Results

We constructed the reference genome of bezoar and sequenced representative breeds of domestic goats to search for genomic changes that likely have accompanied goat domestication and breed formation. Thirteen copy number variation genes associated with coat color were identified in domestic goats, among which ASIP gene duplication contributes to the generation of light coat-color phenotype in domestic goats. Analysis of rapidly evolving genes identified genic changes underlying behavior-related traits, immune response and production-related traits.

Conclusion

Based on the comparison studies of copy number variation genes and rapidly evolving genes between wild and domestic goat, our findings and methodology shed light on the genetic mechanism of animal domestication and will facilitate future goat breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1606-1) contains supplementary material, which is available to authorized users.     

Genetic Variants in Meiotic Program Initiation Pathway Genes Are Associated with Spermatogenic Impairment in a Han Chinese Population

Background

The meiotic program initiation pathway genes (CYP26B1, NANOS1 and STRA8) have been proposed to play key roles in spermatogenesis.

Objective

To elucidate the exact role of the genetic variants of the meiosis initiation genes in spermatogenesis, we genotyped the potential functional genetic variants of CYP26B1, NANOS1 and STRA8 genes, and evaluated their effects on spermatogenesis in our study population.

Design, Setting, and Participants

In this study, all subjects were volunteers from the affiliated hospitals of Nanjing Medical University between March 2004 and July 2009 (NJMU Infertile Study). Total 719 idiopathic infertile cases were recruited and divided into three groups according to WHO semen parameters: 201 azoospermia patients (no sperm in the ejaculate even after centrifugation), 155 oligozoospermia patients (sperm counts <20×10⁶/ml) and 363 infertility/normozoospermia subjects (sperm counts >20×10⁶/ml). The control group consisted of 383 subjects with normal semen parameters, all of which had fathered at least one child without assisted reproductive technologies.

Measurements

Eight single nucleotide polymorphisms (SNPs) in CYP26B1, NANOS1 and STRA8 genes were determined by TaqMan allelic discrimination assay in 719 idiopathic infertile men and 383 healthy controls.

Results and Limitations

The genetic variant rs10269148 of STRA8 gene showed higher risk of spermatogenic impairment in the groups of abnormospermia (including azoospermia subgroup and oligozoospermia subgroup) and azoospermia than the controls with odds ratios and 95% confidence intervals of 2.52 (1.29–4.94) and 2.92 (1.41–6.06), respectively (P = 0.006, 0.002 respective). Notably, larger sample size studies and in vivo or in vitro functional studies are needed to substantiate the biological roles of these variants.

Conclusions

Our results provided epidemiological evidence supporting the involvement of genetic polymorphisms of the meiotic program initiation genes in modifying the risk of azoospermia and oligozoospermia in a Han-Chinese population.     

Validating Discovered Cis-Acting Regulatory Genetic Variants: Application of an Allele Specific Expression Approach to HapMap Populations

Background

Localising regulatory variants that control gene expression is a challenge for genome research. Several studies have recently identified non-coding polymorphisms associated with inter-individual differences in gene expression. These approaches rely on the identification of signals of association against a background of variation due to other genetic and environmental factors. A complementary approach is to use an Allele-Specific Expression (ASE) assay, which is more robust to the effects of environmental variation and trans-acting genetic factors.

Methodology/Principal Findings

Here we apply an ASE method which utilises heterozygosity within an individual to compare expression of the two alleles of a gene in a single cell. We used individuals from three HapMap population groups and analysed the allelic expression of genes with cis-regulatory regions previously identified using total gene expression studies. We were able to replicate the results in five of the six genes tested, and refined the cis- associated regions to a small number of variants. We also showed that by using multi-populations it is possible to refine the associated cis-effect DNA regions.

Conclusions/Significance

We discuss the efficacy and drawbacks of both total gene expression and ASE approaches in the discovery of cis-acting variants. We show that the ASE approach has significant advantages as it is a cleaner representation of cis-acting effects. We also discuss the implication of using different populations to map cis-acting regions and the importance of finding regulatory variants which contribute to human phenotypic variation.     

Genetic variants of the FADS gene cluster and ELOVL gene family, colostrums LC-PUFA levels, breastfeeding, and child cognition

Introduction

Breastfeeding effects on cognition are attributed to long-chain polyunsaturated fatty acids (LC-PUFAs), but controversy persists. Genetic variation in fatty acid desaturase (FADS) and elongase (ELOVL) enzymes has been overlooked when studying the effects of LC-PUFAs supply on cognition. We aimed to: 1) to determine whether maternal genetic variants in the FADS cluster and ELOVL genes contribute to differences in LC-PUFA levels in colostrum; 2) to analyze whether these maternal variants are related to child cognition; and 3) to assess whether children''s variants modify breastfeeding effects on cognition.

Methods

Data come from two population-based birth cohorts (n = 400 mother-child pairs from INMA-Sabadell; and n = 340 children from INMA-Menorca). LC-PUFAs were measured in 270 colostrum samples from INMA-Sabadell. Tag SNPs were genotyped both in mothers and children (13 in the FADS cluster, 6 in ELOVL2, and 7 in ELOVL5). Child cognition was assessed at 14 mo and 4 y using the Bayley Scales of Infant Development and the McCarthy Scales of Children''s Abilities, respectively.

Results

Children of mothers carrying genetic variants associated with lower FADS1 activity (regulating AA and EPA synthesis), higher FADS2 activity (regulating DHA synthesis), and with higher EPA/AA and DHA/AA ratios in colostrum showed a significant advantage in cognition at 14 mo (3.5 to 5.3 points). Not being breastfed conferred an 8- to 9-point disadvantage in cognition among children GG homozygote for rs174468 (low FADS1 activity) but not among those with the A allele. Moreover, not being breastfed resulted in a disadvantage in cognition (5 to 8 points) among children CC homozygote for rs2397142 (low ELOVL5 activity), but not among those carrying the G allele.

Conclusion

Genetically determined maternal supplies of LC-PUFAs during pregnancy and lactation appear to be crucial for child cognition. Breastfeeding effects on cognition are modified by child genetic variation in fatty acid desaturase and elongase enzymes.     

Genomic copy number variation in Mus musculus

Background

Copy number variation is an important dimension of genetic diversity and has implications in development and disease. As an important model organism, the mouse is a prime candidate for copy number variant (CNV) characterization, but this has yet to be completed for a large sample size. Here we report CNV analysis of publicly available, high-density microarray data files for 351 mouse tail samples, including 290 mice that had not been characterized for CNVs previously.

Results

We found 9634 putative autosomal CNVs across the samples affecting 6.87 % of the mouse reference genome. We find significant differences in the degree of CNV uniqueness (single sample occurrence) and the nature of CNV-gene overlap between wild-caught mice and classical laboratory strains. CNV-gene overlap was associated with lipid metabolism, pheromone response and olfaction compared to immunity, carbohydrate metabolism and amino-acid metabolism for wild-caught mice and classical laboratory strains, respectively. Using two subspecies of wild-caught Mus musculus, we identified putative CNVs unique to those subspecies and show this diversity is better captured by wild-derived laboratory strains than by the classical laboratory strains. A total of 9 genic copy number variable regions (CNVRs) were selected for experimental confirmation by droplet digital PCR (ddPCR).

Conclusion

The analysis we present is a comprehensive, genome-wide analysis of CNVs in Mus musculus, which increases the number of known variants in the species and will accelerate the identification of novel variants in future studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1713-z) contains supplementary material, which is available to authorized users.     

seqCNA: an R package for DNA copy number analysis in cancer using high-throughput sequencing

Background

Deviations in the amount of genomic content that arise during tumorigenesis, called copy number alterations, are structural rearrangements that can critically affect gene expression patterns. Additionally, copy number alteration profiles allow insight into cancer discrimination, progression and complexity. On data obtained from high-throughput sequencing, improving quality through GC bias correction and keeping false positives to a minimum help build reliable copy number alteration profiles.

Results

We introduce seqCNA, a parallelized R package for an integral copy number analysis of high-throughput sequencing cancer data. The package includes novel methodology on (i) filtering, reducing false positives, and (ii) GC content correction, improving copy number profile quality, especially under great read coverage and high correlation between GC content and copy number. Adequate analysis steps are automatically chosen based on availability of paired-end mapping, matched normal samples and genome annotation.

Conclusions

seqCNA, available through Bioconductor, provides accurate copy number predictions in tumoural data, thanks to the extensive filtering and better GC bias correction, while providing an integrated and parallelized workflow.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-178) contains supplementary material, which is available to authorized users.     

Genetic Variants and Early Cigarette Smoking and Nicotine Dependence Phenotypes in Adolescents

Background

While the heritability of cigarette smoking and nicotine dependence (ND) is well-documented, the contribution of specific genetic variants to specific phenotypes has not been closely examined. The objectives of this study were to test the associations between 321 tagging single-nucleotide polymorphisms (SNPs) that capture common genetic variation in 24 genes, and early smoking and ND phenotypes in novice adolescent smokers, and to assess if genetic predictors differ across these phenotypes.

Methods

In a prospective study of 1294 adolescents aged 12–13 years recruited from ten Montreal-area secondary schools, 544 participants who had smoked at least once during the 7–8 year follow-up provided DNA. 321 single-nucleotide polymorphisms (SNPs) in 24 candidate genes were tested for an association with number of cigarettes smoked in the past 3 months, and with five ND phenotypes (a modified version of the Fagerstrom Tolerance Questionnaire, the ICD-10 and three clusters of ND symptoms representing withdrawal symptoms, use of nicotine for self-medication, and a general ND/craving symptom indicator).

Results

The pattern of SNP-gene associations differed across phenotypes. Sixteen SNPs in seven genes (ANKK1, CHRNA7, DDC, DRD2, COMT, OPRM1, SLC6A3 (also known as DAT1)) were associated with at least one phenotype with a p-value <0.01 using linear mixed models. After permutation and FDR adjustment, none of the associations remained statistically significant, although the p-values for the association between rs557748 in OPRM1 and the ND/craving and self-medication phenotypes were both 0.076.

Conclusions

Because the genetic predictors differ, specific cigarette smoking and ND phenotypes should be distinguished in genetic studies in adolescents. Fifteen of the 16 top-ranked SNPs identified in this study were from loci involved in dopaminergic pathways (ANKK1/DRD2, DDC, COMT, OPRM1, and SLC6A3).

Impact

Dopaminergic pathways may be salient during early smoking and the development of ND.