ParseCNV integrative copy number variation association software with quality tracking

A number of copy number variation (CNV) calling algorithms exist; however, comprehensive software tools for CNV association studies are lacking. We describe ParseCNV, unique software that takes CNV calls and creates probe-based statistics for CNV occurrence in both case–control design and in family based studies addressing both de novo and inheritance events, which are then summarized based on CNV regions (CNVRs). CNVRs are defined in a dynamic manner to allow for a complex CNV overlap while maintaining precise association region. Using this approach, we avoid failure to converge and non-monotonic curve fitting weaknesses of programs, such as CNVtools and CNVassoc, and although Plink is easy to use, it only provides combined CNV state probe-based statistics, not state-specific CNVRs. Existing CNV association methods do not provide any quality tracking information to filter confident associations, a key issue which is fully addressed by ParseCNV. In addition, uncertainty in CNV calls underlying CNV associations is evaluated to verify significant results, including CNV overlap profiles, genomic context, number of probes supporting the CNV and single-probe intensities. When optimal quality control parameters are followed using ParseCNV, 90% of CNVs validate by polymerase chain reaction, an often problematic stage because of inadequate significant association review. ParseCNV is freely available at http://parsecnv.sourceforge.net.     

Introduction to galectins

Good evidence suggest roles of galectins in cancer, immunity and inflammation, and development, but a unifying picture of their biological function is lacking. Instead galectins appear to have a particularly diverse, bewildering but intriguing array of activities both inside and outside cells--"clear truths and mysteries are inextricably twined". Fortunately this has not discouraged but rather enthused a large number of good galectin researchers, some of which have contributed to this special issue of Glycoconjugate Journal to provide a personal, critical status of the field. Here we will give a brief introduction to the galectins as a protein family with some comments on nomenclature.     

Activation of the neutrophil nicotinamide adenine dinucleotide phosphate oxidase by galectin-1

Galectins are a group of lactose-binding proteins widely distributed in nature. Twelve mammalian galectins have so far been identified, but their functions are to a large extent unknown. In this work we study galectin-1 in its interaction with human neutrophils, with regard to both cell surface binding and activation of the superoxide-producing NADPH-oxidase. We show that galectin-1 is able to activate the neutrophil NADPH-oxidase, provided that the cells have been primed by extravasation from the blood into the tissue, an activation pattern that is similar to that of galectin-3. Using in vitro priming protocols, the galectin-1 responsiveness was found to correlate to granule mobilization and galectin-1 binding to the cells, suggesting the presence of granule-localized receptors that are up-regulated to the cell surface upon priming. By galectin-1 overlay of fractionated neutrophils we identified potential galectin-1 receptor candidates localized in the membranes of the secretory vesicle and gelatinase granules. The binding of galectin-1 and galectin-3 to neutrophil proteins was compared, as were the dose dependencies for activation by the two lectins. The results suggest that, although similarities are found between the two galectins, they appear to activate the NADPH-oxidase using different receptors. In conclusion, galectin-1 appears to have proinflammatory functions, mediated through activation of the neutrophil respiratory burst.     

The Role of ARF6 in Biliary Atresia

Background & Aims

Altered extrahepatic bile ducts, gut, and cardiovascular anomalies constitute the variable phenotype of biliary atresia (BA).

Methods

To identify potential susceptibility loci, Caucasian children, normal (controls) and with BA (cases) at two US centers were compared at >550000 SNP loci. Systems biology analysis was carried out on the data. In order to validate a key gene identified in the analysis, biliary morphogenesis was evaluated in 2-5-day post-fertilization zebrafish embryos after morpholino-antisense oligonucleotide knockdown of the candidate gene ADP ribosylation factor-6 (ARF6, Mo-arf6).

Results

Among 39 and 24 cases at centers 1 and 2, respectively, and 1907 controls, which clustered together on principal component analysis, the SNPs rs3126184 and rs10140366 in a 3’ flanking enhancer region for ARF6 demonstrated higher minor allele frequencies (MAF) in each cohort, and 63 combined cases, compared with controls (0.286 vs. 0.131, P = 5.94x10^-7, OR 2.66; 0.286 vs. 0.13, P = 5.57x10^-7, OR 2.66). Significance was enhanced in 77 total cases, which included 14 additional BA genotyped at rs3126184 only (p = 1.58x10^-2, OR = 2.66). Pathway analysis of the 1000 top-ranked SNPs in CHP cases revealed enrichment of genes for EGF regulators (p<1 x10^-7), ERK/MAPK and CREB canonical pathways (p<1 x10^-34), and functional networks for cellular development and proliferation (p<1 x10^-45), further supporting the role of EGFR-ARF6 signaling in BA. In zebrafish embryos, Mo-arf6 injection resulted in a sparse intrahepatic biliary network, several biliary epithelial cell defects, and poor bile excretion to the gall bladder compared with uninjected embryos. Biliary defects were reproduced with the EGFR-blocker AG1478 alone or with Mo-arf6 at lower doses of each agent and rescued with arf6 mRNA.

Conclusions

The BA-associated SNPs identify a chromosome 14q21.3 susceptibility locus encompassing the ARF6 gene. arf6 knockdown in zebrafish implicates early biliary dysgenesis as a basis for BA, and also suggests a role for EGFR signaling in BA pathogenesis.     

Pathway-Based Genome-Wide Association Studies for Plasma Triglycerides in Obese Females and Normal-Weight Controls

Pathway-based analysis as an alternative approach can provide complementary information to single-marker genome-wide association studies (GWASs), which always ignore the epistasis and does not have sufficient power to find rare variants. In this study, using genotypes from a genome-wide association study (GWAS), pathway-based association studies were carried out by a modified Gene Set Enrichment Algorithm (GSEA) method (GenGen) for triglyceride in 1028 unrelated European-American extremely obese females (BMI≥35kg/m²) and normal-weight controls (BMI<25kg/m²), and another pathway association analysis (ICSNPathway) was also used to verify the GenGen result in the same data. The GO0009110 pathway (vitamin anabolism) was among the strongest associations with triglyceride (empirical P<0.001); the result remained significant after FDR correction (P = 0.022). MMAB, an obesity-related locus, included in this pathway. The ABCG1 and BCL6 gene was found in several triglyceride-related pathways (empirical P<0.05), which were also replicated by ICSNPathway (empirical P<0.05, FDR<0.05). We also performed single-marked GWAS using PLINK for TG levels (log-transformed). Significant associations were found between ASTN2 gene SNPs and plasma triglyceride levels (rs7035794, P = 2.24×10⁻¹⁰). Our study suggested that vitamin anabolism pathway, BCL6 gene pathways and ASTN2 gene may contribute to the genetic variation of plasma triglyceride concentrations.     

Phenotype restricted genome-wide association study using a gene-centric approach identifies three low-risk neuroblastoma susceptibility Loci

Neuroblastoma is a malignant neoplasm of the developing sympathetic nervous system that is notable for its phenotypic diversity. High-risk patients typically have widely disseminated disease at diagnosis and a poor survival probability, but low-risk patients frequently have localized tumors that are almost always cured with little or no chemotherapy. Our genome-wide association study (GWAS) has identified common variants within FLJ22536, BARD1, and LMO1 as significantly associated with neuroblastoma and more robustly associated with high-risk disease. Here we show that a GWAS focused on low-risk cases identified SNPs within DUSP12 at 1q23.3 (P = 2.07×10⁻⁶), DDX4 and IL31RA both at 5q11.2 (P = 2.94×10⁻⁶ and 6.54×10⁻⁷ respectively), and HSD17B12 at 11p11.2 (P = 4.20×10⁻⁷) as being associated with the less aggressive form of the disease. These data demonstrate the importance of robust phenotypic data in GWAS analyses and identify additional susceptibility variants for neuroblastoma.     

Role of BMI-associated loci identified in GWAS meta-analyses in the context of common childhood obesity in European Americans

Obesity is a serious health concern for children and adolescents, particularly in Western societies, where its incidence is now considered to have reached epidemic proportions. A number of genetic determinants of adult BMI have already been established through genome wide association studies (GWAS), most recently from the GIANT meta-analysis of such datasets combined. In this current study of European Americans, we examined the 32 loci detected in that GIANT study in the context of common childhood obesity within a cohort of 1,097 cases (defined as BMI ≥95th percentile), together with 2,760 lean controls (defined as BMI <50th percentile), aged between 2 and 18 years old. Nine of these single-nucleotide polymorphims (SNPs) yielded at least nominal evidence for association with common childhood obesity, namely at the FTO, TMEM18, NRXN3, MC4R, SEC16B, GNPDA2, TNNI3K, QPCTL, and BDNF loci. However, overall 28 of the 32 loci showed directionally consistent effects to that of the adult BMI meta-analysis. We conclude that among the 32 loci that have been reported to associate with adult BMI in the largest meta-analysis of BMI to date, at least nine also contribute to the determination of common obesity in childhood in European Americans, as demonstrated by their associations in our pediatric cohort.     

Phasing of many thousands of genotyped samples

Haplotypes are an important resource for a large number of applications in human genetics, but computationally inferred haplotypes are subject to switch errors that decrease their utility. The accuracy of computationally inferred haplotypes increases with sample size, and although ever larger genotypic data sets are being generated, the fact that existing methods require substantial computational resources limits their applicability to data sets containing tens or hundreds of thousands of samples. Here, we present HAPI-UR (haplotype inference for unrelated samples), an algorithm that is designed to handle unrelated and/or trio and duo family data, that has accuracy comparable to or greater than existing methods, and that is computationally efficient and can be applied to 100,000 samples or more. We use HAPI-UR to phase a data set with 58,207 samples and show that it achieves practical runtime and that switch errors decrease with sample size even with the use of samples from multiple ethnicities. Using a data set with 16,353 samples, we compare HAPI-UR to Beagle, MaCH, IMPUTE2, and SHAPEIT and show that HAPI-UR runs 18× faster than all methods and has a lower switch-error rate than do other methods except for Beagle; with the use of consensus phasing, running HAPI-UR three times gives a slightly lower switch-error rate than Beagle does and is more than six times faster. We demonstrate results similar to those from Beagle on another data set with a higher marker density. Lastly, we show that HAPI-UR has better runtime scaling properties than does Beagle so that for larger data sets, HAPI-UR will be practical and will have an even larger runtime advantage. HAPI-UR is available online (see Web Resources).     

Results of Genome-Wide Analyses on Neurodevelopmental Phenotypes at Four-Year Follow-Up following Cardiac Surgery in Infancy

Background

Adverse neurodevelopmental sequelae are reported among children who undergo early cardiac surgery to repair congenital heart defects (CHD). APOE genotype has previously been determined to contribute to the prediction of these outcomes. Understanding further genetic causes for the development of poor neurobehavioral outcomes should enhance patient risk stratification and improve both prevention and treatment strategies.

Methods

We performed a prospective observational study of children who underwent cardiac surgery before six months of age; this included a neurodevelopmental evaluation between their fourth and fifth birthdays. Attention and behavioral skills were assessed through parental report utilizing the Attention Deficit-Hyperactivity Disorder-IV scale preschool edition (ADHD-IV), and Child Behavior Checklist (CBCL/1.5-5), respectively. Of the seven investigated, three neurodevelopmental phenotypes met genomic quality control criteria. Linear regression was performed to determine the effect of genome-wide genetic variation on these three neurodevelopmental measures in 316 subjects.

Results

This genome-wide association study identified single nucleotide polymorphisms (SNPs) associated with three neurobehavioral phenotypes in the postoperative children ADHD-IV Impulsivity/Hyperactivity, CBCL/1.5-5 PDPs, and CBCL/1.5-5 Total Problems. The most predictive SNPs for each phenotype were: a LGALS8 intronic SNP, rs4659682, associated with ADHD-IV Impulsivity (P = 1.03×10⁻⁶); a PCSK5 intronic SNP, rs2261722, associated with CBCL/1.5-5 PDPs (P = 1.11×10⁻⁶); and an intergenic SNP, rs11617488, 50 kb from FGF9, associated with CBCL/1.5-5 Total Problems (P = 3.47×10⁻⁷). 10 SNPs (3 for ADHD-IV Impulsivity, 5 for CBCL/1.5-5 PDPs, and 2 for CBCL/1.5-5 Total Problems) had p<10⁻⁵.

Conclusions

No SNPs met genome-wide significance for our three neurobehavioral phenotypes; however, 10 SNPs reached a threshold for suggestive significance (p<10⁻⁵). Given the unique nature of this cohort, larger studies and/or replication are not possible. Studies to further investigate the mechanisms through which these newly identified genes may influence neurodevelopment dysfunction are warranted.