Identification of CAD candidate genes in GWAS loci and their expression in vascular cells

Recent genome-wide association studies (GWAS) have identified 35 loci that significantly associate with coronary artery disease (CAD) susceptibility. The majority of the genes represented in these loci have not previously been studied in the context of atherosclerosis. To characterize the roles of these candidate genes in the vessel wall, we determined their expression levels in endothelial, smooth muscle, and macrophage cells isolated from healthy, prelesioned, and lesioned mouse aortas. We also performed expression quantitative locus (eQTL) mapping of these genes in human endothelial cells under control and proatherogenic conditions. Of the 57 genes studied, 31 were differentially expressed in one or more cell types in disease state in mice, and the expression levels of 8 were significantly associated with the CAD SNPs in human cells, 7 of which were also differentially expressed in mice. By integrating human and mouse results, we predict that PPAP2B, GALNT4, MAPKAPK5, TCTN1, SRR, SNF8, and ICAM1 play a causal role in the susceptibility to atherosclerosis through a role in the vasculature. Additionally, we highlight the genetic complexity of a subset of CAD loci through the differential expression of multiple candidate genes per locus and the involvement of genes that lie outside linkage disequilibrium blocks.     

A Parent-of-Origin Effect Determines the Susceptibility of a Non-Informative F1 Population to Trypanosoma cruzi Infection In Vivo

The development of Chagas disease is determined by a complex interaction between the genetic traits of both the protozoan parasite, T. cruzi, and the infected host. This process is regulated by multiple genes that control different aspects of the host-parasite interaction. While determination of the relevant genes in humans is extremely difficult, it is feasible to use inbred mouse strains to determine the genes and loci responsible for host resistance to infection. In this study, we investigated the susceptibility of several inbred mouse strains to infection with the highly virulent Y strain of T. cruzi and found a considerable difference in susceptibility between A/J and C57BL/6 mice. We explored the differences between these two mouse strains and found that the A/J strain presented higher mortality, exacerbated and uncontrolled parasitemia and distinct histopathology in the target organs, which were associated with a higher parasite burden and more extensive tissue lesions. We then employed a genetic approach to assess the pattern of inheritance of the resistance phenotype in an F1 population and detected a strong parent-of-origin effect determining the susceptibility of the F1 male mice. This effect is unlikely to result from imprinted genes because the inheritance of this susceptibility was affected by the direction of the parental crossing. Collectively, our genetic approach of using the F1 population suggests that genes contained in the murine chromosome X contribute to the natural resistance against T. cruzi infection. Future linkage studies may reveal the locus and genes participating on the host resistance process reported herein.     

Integration of Mouse and Human Genome-Wide Association Data Identifies KCNIP4 as an Asthma Gene

Asthma is a common chronic respiratory disease characterized by airway hyperresponsiveness (AHR). The genetics of asthma have been widely studied in mouse and human, and homologous genomic regions have been associated with mouse AHR and human asthma-related phenotypes. Our goal was to identify asthma-related genes by integrating AHR associations in mouse with human genome-wide association study (GWAS) data. We used Efficient Mixed Model Association (EMMA) analysis to conduct a GWAS of baseline AHR measures from males and females of 31 mouse strains. Genes near or containing SNPs with EMMA p-values <0.001 were selected for further study in human GWAS. The results of the previously reported EVE consortium asthma GWAS meta-analysis consisting of 12,958 diverse North American subjects from 9 study centers were used to select a subset of homologous genes with evidence of association with asthma in humans. Following validation attempts in three human asthma GWAS (i.e., Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG) and two human AHR GWAS (i.e., SHARP, DAG), the Kv channel interacting protein 4 (KCNIP4) gene was identified as nominally associated with both asthma and AHR at a gene- and SNP-level. In EVE, the smallest KCNIP4 association was at rs6833065 (P-value 2.9e-04), while the strongest associations for Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG were 1.5e-03, 1.0e-03, 3.1e-03 at rs7664617, rs4697177, rs4696975, respectively. At a SNP level, the strongest association across all asthma GWAS was at rs4697177 (P-value 1.1e-04). The smallest P-values for association with AHR were 2.3e-03 at rs11947661 in SHARP and 2.1e-03 at rs402802 in DAG. Functional studies are required to validate the potential involvement of KCNIP4 in modulating asthma susceptibility and/or AHR. Our results suggest that a useful approach to identify genes associated with human asthma is to leverage mouse AHR association data.     

Copy Number Variation of the Beta Defensin Gene Cluster on Chromosome 8p Influences the Bacterial Microbiota within the Nasopharynx of Otitis-Prone Children

As there is increasing evidence that aberrant defensin expression is related to susceptibility for infectious disease and inflammatory disorders, we sought to determine if copy number of the beta-defensin gene cluster located on chromosome 8p23.1 (DEFB107, 106, 105, 104, 103, DEFB4 and SPAG11), that shows copy number variation as a block, was associated with susceptibility to otitis media (OM). The gene DEFB103 within this complex encodes human beta defensin-3 (hBD-3), an antimicrobial peptide (AP) expressed by epithelial cells that line the mammalian airway, important for defense of mucosal surfaces and previously shown to have bactericidal activity in vitro against multiple human pathogens, including the three that predominate in OM. To this end, we conducted a retrospective case-control study of 113 OM prone children and 267 controls aged five to sixty months. We identified the copy number of the above defined beta-defensin gene cluster (DEFB-CN) in each study subject by paralogue ratio assays. The mean DEFB-CN was indistinguishable between subjects classified as OM prone based on a recent history of multiple episodes of OM and control subjects who had no history of OM (4.4±0.96 versus 4.4±1.08, respectively: Odds Ratio [OR]: 1.16 (95% CI: 0.61, 2.20). Despite a lack of direct association, we observed a statistically significant correlation between DEFB-CN and nasopharyngeal bacterial colonization patterns. Collectively, our findings suggested that susceptibility to OM might be mediated by genetic variation among individuals, wherein a DEFB-CN less than 4 exerts a marked influence on the microbiota of the nasopharynx, specifically with regard to colonization by the three predominant bacterial pathogens of OM.     

Quantitative trait loci in a bacterially induced model of inflammatory bowel disease

Inflammatory bowel diseases (IBDs) are complex disorders caused by a combination of environmental, microbial, and genetic factors. Genome-wide association studies in humans have successfully identified multiple genes and loci associated with disease susceptibility, but the mechanisms by which these loci interact with each other and/or with environmental factors (i.e., intestinal microbiota) to cause disease are poorly understood. Helicobacter hepaticus-induced intestinal inflammation in mice is an ideal model system for elucidating the genetic basis of IBD susceptibility in a bacterially induced system, as there are significant differences in H. hepaticus-induced disease susceptibility among inbred mouse strains. Infected A/J mice develop acute overexpression of proinflammatory cytokines followed 2?C3?months later by chronic cecal inflammation, whereas infected C57BL/6 mice fail to develop cecal inflammation or increased cytokine expression. The goal of this project was to use quantitative trait locus (QTL) mapping to evaluate genetic factors that contribute to the differential disease susceptibility between these two mouse strains. Using acute cecal IL-12/23p40 expression as a biomarker for disease susceptibility, QTL analysis of H. hepaticus-infected F₂ mice revealed involvement of multiple loci. The loci with the strongest association were located on Chromosome 3 and Chromosome 17, with logarithm of odds (LOD) scores of 6.89 and 3.09, respectively. Cecal expression of IL-12/23p40 in H. hepaticus-infected C57BL/6J-Chr3^A/J/NaJ chromosome substitution mice had an intermediate phenotype, significantly higher than in resistant C57BL/6 but lower than in susceptible A/J mice, confirming the importance of this locus to the immune response to H. hepaticus infection.     

Refining Susceptibility Loci of Chronic Obstructive Pulmonary Disease with Lung eqtls

Chronic obstructive pulmonary disease (COPD) is the fourth leading cause of mortality worldwide. Recent genome-wide association studies (GWAS) have identified robust susceptibility loci associated with COPD. However, the mechanisms mediating the risk conferred by these loci remain to be found. The goal of this study was to identify causal genes/variants within susceptibility loci associated with COPD. In the discovery cohort, genome-wide gene expression profiles of 500 non-tumor lung specimens were obtained from patients undergoing lung surgery. Blood-DNA from the same patients were genotyped for 1,2 million SNPs. Following genotyping and gene expression quality control filters, 409 samples were analyzed. Lung expression quantitative trait loci (eQTLs) were identified and overlaid onto three COPD susceptibility loci derived from GWAS; 4q31 (HHIP), 4q22 (FAM13A), and 19q13 (RAB4B, EGLN2, MIA, CYP2A6). Significant eQTLs were replicated in two independent datasets (n = 363 and 339). SNPs previously associated with COPD and lung function on 4q31 (rs1828591, rs13118928) were associated with the mRNA expression of HHIP. An association between mRNA expression level of FAM13A and SNP rs2045517 was detected at 4q22, but did not reach statistical significance. At 19q13, significant eQTLs were detected with EGLN2. In summary, this study supports HHIP, FAM13A, and EGLN2 as the most likely causal COPD genes on 4q31, 4q22, and 19q13, respectively. Strong lung eQTL SNPs identified in this study will need to be tested for association with COPD in case-control studies. Further functional studies will also be needed to understand the role of genes regulated by disease-related variants in COPD.     

Genetic dissection of testis weight in mice: quantitative trait locus analysis using F2 intercrosses between strains with extreme testis weight, and association study using Y-consomic strains

In the present study, dissection of genetic bases of testis weight in mice was performed. Autosomes and the X chromosome were searched using traditional quantitative trait locus (QTL) scans, and the Y chromosome was searched by association studies of Y-consomic strains. QTL analysis was performed in ??DDD?×???CBA F₂ mice; the inbred mouse DDD has the heaviest testes, whereas the inbred mouse CBA has the lightest testes. Two significant testis weight QTLs were identified on chromosomes 1 and X. A DDD allele was associated with increased and decreased testis weight at the locus on chromosomes 1 and X, respectively. In the reciprocal cross ??CBA?×???DDD F₂ mice, QTL on chromosome 1, and not on chromosome X, had a significant effect on testis weight. The DDD allele at the X-linked locus could not sustain testis weight in combination with the Y chromosome of the CBA strain. The Y chromosome per se had a significant effect on testis weight, i.e., DH-Chr Y^DDD had significantly heavier testes than DH-Chr Y^CBA. On the basis of the results of Y-chromosome-wide association studies using 17 Y-consomic strains, variations in Uty, Usp9y, and Sry were significantly associated with testis weight. Thus, testis weight is a complex quantitative phenotype controlled by multiple genes on autosomes and sex chromosomes and their interactions.     

Identification of KIF3A as a novel candidate gene for childhood asthma using RNA expression and population allelic frequencies differences

Background

Asthma is a chronic inflammatory disease with a strong genetic predisposition. A major challenge for candidate gene association studies in asthma is the selection of biologically relevant genes.

Methodology/Principal Findings

Using epithelial RNA expression arrays, HapMap allele frequency variation, and the literature, we identified six possible candidate susceptibility genes for childhood asthma including ADCY2, DNAH5, KIF3A, PDE4B, PLAU, SPRR2B. To evaluate these genes, we compared the genotypes of 194 predominantly tagging SNPs in 790 asthmatic, allergic and non-allergic children. We found that SNPs in all six genes were nominally associated with asthma (p<0.05) in our discovery cohort and in three independent cohorts at either the SNP or gene level (p<0.05). Further, we determined that our selection approach was superior to random selection of genes either differentially expressed in asthmatics compared to controls (p = 0.0049) or selected based on the literature alone (p = 0.0049), substantiating the validity of our gene selection approach. Importantly, we observed that 7 of 9 SNPs in the KIF3A gene more than doubled the odds of asthma (OR = 2.3, p<0.0001) and increased the odds of allergic disease (OR = 1.8, p<0.008). Our data indicate that KIF3A rs7737031 (T-allele) has an asthma population attributable risk of 18.5%. The association between KIF3A rs7737031 and asthma was validated in 3 independent populations, further substantiating the validity of our gene selection approach.

Conclusions/Significance

Our study demonstrates that KIF3A, a member of the kinesin superfamily of microtubule associated motors that are important in the transport of protein complexes within cilia, is a novel candidate gene for childhood asthma. Polymorphisms in KIF3A may in part be responsible for poor mucus and/or allergen clearance from the airways. Furthermore, our study provides a promising framework for the identification and evaluation of novel candidate susceptibility genes.     

Whole Exome Sequencing Suggests Much of Non-BRCA1/BRCA2 Familial Breast Cancer Is Due to Moderate and Low Penetrance Susceptibility Alleles

The identification of the two most prevalent susceptibility genes in breast cancer, BRCA1 and BRCA2, was the beginning of a sustained effort to uncover new genes explaining the missing heritability in this disease. Today, additional high, moderate and low penetrance genes have been identified in breast cancer, such as P53, PTEN, STK11, PALB2 or ATM, globally accounting for around 35 percent of the familial cases. In the present study we used massively parallel sequencing to analyze 7 BRCA1/BRCA2 negative families, each having at least 6 affected women with breast cancer (between 6 and 10) diagnosed under the age of 60 across generations. After extensive filtering, Sanger sequencing validation and co-segregation studies, variants were prioritized through either control-population studies, including up to 750 healthy individuals, or case-control assays comprising approximately 5300 samples. As a result, a known moderate susceptibility indel variant (CHEK2 1100delC) and a catalogue of 11 rare variants presenting signs of association with breast cancer were identified. All the affected genes are involved in important cellular mechanisms like DNA repair, cell proliferation and survival or cell cycle regulation. This study highlights the need to investigate the role of rare variants in familial cancer development by means of novel high throughput analysis strategies optimized for genetically heterogeneous scenarios. Even considering the intrinsic limitations of exome resequencing studies, our findings support the hypothesis that the majority of non-BRCA1/BRCA2 breast cancer families might be explained by the action of moderate and/or low penetrance susceptibility alleles.     

Identification of candidate genes for dyslexia susceptibility on chromosome 18

Background

Six independent studies have identified linkage to chromosome 18 for developmental dyslexia or general reading ability. Until now, no candidate genes have been identified to explain this linkage. Here, we set out to identify the gene(s) conferring susceptibility by a two stage strategy of linkage and association analysis.

Methodology/Principal Findings

Linkage analysis: 264 UK families and 155 US families each containing at least one child diagnosed with dyslexia were genotyped with a dense set of microsatellite markers on chromosome 18. Association analysis: Using a discovery sample of 187 UK families, nearly 3000 SNPs were genotyped across the chromosome 18 dyslexia susceptibility candidate region. Following association analysis, the top ranking SNPs were then genotyped in the remaining samples. The linkage analysis revealed a broad signal that spans approximately 40 Mb from 18p11.2 to 18q12.2. Following the association analysis and subsequent replication attempts, we observed consistent association with the same SNPs in three genes; melanocortin 5 receptor (MC5R), dymeclin (DYM) and neural precursor cell expressed, developmentally down-regulated 4-like (NEDD4L).

Conclusions

Along with already published biological evidence, MC5R, DYM and NEDD4L make attractive candidates for dyslexia susceptibility genes. However, further replication and functional studies are still required.     

A GWA study reveals genetic loci for body conformation traits in Chinese Laiwu pigs and its implications for human BMI

Pigs share numerous physiological and phenotypic similarities with human and thus have been considered as a good model in nonrodent mammals for the study of genetic basis of human obesity. Researches on candidate genes for obesity traits have successfully identified some common genes between humans and pigs. However, few studies have assessed how many similarities exist between the genetic architecture of obesity in pigs and humans by large-scale comparative genomics. Here, we performed a genome-wide association study (GWAS) using the porcine 60 K SNP Beadchip for BMI and other four conformation traits at three different ages in a Chinese Laiwu pig population, which shows a large variability in fat deposition. In total, 35 SNPs were found to be significant at Bonferroni-corrected 5 % chromosome-wise level (P = 2.13 × 10^?5) and 88 SNPs had suggestive (P < 10^?4) association with the conformation traits. Some SNPs showed age-dependent association. Intriguingly, out of 32 regions associated with BMI in pigs, 18 were homologous with the loci for BMI in humans. Furthermore, five closest genes to GWAS peaks including HIF1AN, SMYD3, COX10, SLMAP, and GBE1 have been already associated with BMI in humans, which makes them very promising candidates for these QTLs. The result of GO analysis provided strong support to the fact that mitochondria and synapse play important roles in obesity susceptibility, which is consistent with previous findings on human obesity, and it also implicated new gene sets related to chromatin modification and Ig-like C2-type 5 domain. Therefore, these results not only provide new insights into the genetic architecture of BMI in pigs but also highlight that humans and pigs share the significant overlap of obesity-related genes.     

Genetic association study of adaptor protein complex 4 with cerebral palsy in a Han Chinese population

Adaptor protein complex 4 (AP-4) plays a key role in vesicle formation, trafficking, and sorting processes that are critical for brain development and function. AP-4 consists of four subunits encoded by the AP4E1, AP4B1, AP4M1, and AP4S1 genes. A number of studies have pointed to the involvement of AP-4–mediated vesicular trafficking pathways in the etiology of cerebral palsy (CP), the most notable of which are the causative mutations that have recently been identified in each of the AP-4 genes in different CP families. We postulated, therefore, that variations in AP-4 genes might influence an indivual’s susceptibility to CP. In the present study, 16 SNPs were genotyped among 517 CP patients and 502 healthy controls from the Han Chinese population. We systematically analyzed the association of the AP4E1, AP4B1, AP4M1, and AP4S1 genes with CP on the basis of clinical characteristics. No significant associations were found between these variants and the overall risk of CP. Subgroup analysis showed that rs1217401 of AP4B1 was significantly associated with CP as a sequela of hypoxic-ischemic encephalopathy (HIE) (CP + HIE) (allele: p = 0.042151; genotype: p = 4.46 × 10^?6). Our results indicate that the 16 variants studied in the genes of the four subunits of AP-4 have no detectable effects on the overall susceptibility to CP, but AP4B1 appears to be a susceptibility gene for CP + HIE in the Han Chinese population.     

Gene-gene interaction between tuberculosis candidate genes in a South African population

In a complex disease such as tuberculosis (TB) it is increasingly evident that gene-gene interactions play a far more important role in an individual??s susceptibility to develop the disease than single polymorphisms on their own, as one gene can enhance or hinder the expression of another gene. Gene-gene interaction analysis is a new approach to elucidate susceptibility to TB. The possibility of gene-gene interactions was assessed, focusing on 11 polymorphisms in nine genes (DC-SIGN, IFN-??, IFNGR1, IL-8, IL-1Ra, MBL, NRAMP1, RANTES, and SP-D) that have been associated with TB, some repeatedly. An optimal model, which best describes and predicts TB case?Ccontrol status, was constructed. Significant interactions were detected between eight pairs of variants. The models fitted the observed data extremely well, with p?<?0.0001 for all eight models. A highly significant interaction was detected between INFGR1 and NRAMP1, which is not surprising because macrophage activation is greatly enhanced by IFN-?? and IFN-?? response elements that are present in the human NRAMP1 promoter region, providing further evidence for their interaction. This study enabled us to test the theory that disease outcome may be due to interaction of several gene effects. With eight instances of statistically significant gene-gene interactions, the importance of epistasis is clearly identifiable in this study. Methods for studying gene-gene interactions are based on a multilocus and multigene approach, consistent with the nature of complex-trait diseases, and may provide the paradigm for future genetic studies of TB.     

Genetic requirement for pneumococcal ear infection

Background

Ear infection or otitis media (OM) accounts for most bacterial respiratory infections in children in both developed and developing nations. Streptococcus pneumoniae, nontypeable Haemophilus influenzae, and Moraxella catarrhalis are the major OM pathogens. However, little is known about the genetic basis of bacterial OM largely due to practical difficulties in conducting research in ear infection models and genetically manipulating clinical isolates. Here, we report the first genome-scale in vivo screen for bacterial genes required for ear infection in a chinchilla model by signature tagged mutagenesis (STM), a high throughput mutant screen technique.

Methodology/Principal Findings

STM strains were constructed with a multi-drug resistant OM isolate ST556 (serotype 19F) and screened in a chinchilla OM model. Out of 5,280 mutants tested, 248 mutants were substantially underrepresented in the mutant pools recovered from the middle ear fluids of the infected chinchillas, indicating the impaired ability to survive and replicate in the middle ears due to genetic disruptions in the chromosome of strain ST556. Further DNA sequencing analysis mapped the mutations to 169 pneumococcal genes. Surprisingly, only 52 of these genes were required for pneumococcal nasopharyngeal colonization in a murine model. This infection site-specific gene requirement was verified by targeted mutagenesis in the selected genes.

Conclusions/Significance

These findings suggest that there are a subset of pneumococcal genes required for ear infection and that these may be distinct from those required for nasal colonization. Our data thus provide comprehensive gene targets for mechanistic understanding of pneumococcal ear infection. Finally, this study has also developed a model for future genome-scale search for virulence determinants in other pathogens associated with ear infections.     

Insights on ADAMTS proteases and ADAMTS-like proteins from mammalian genetics

The mammalian ADAMTS superfamily comprises 19 secreted metalloproteinases and 7 ADAMTS-like proteins, each the product of a distinct gene. Thus far, all appear to be relevant to extracellular matrix function or to cell–matrix interactions. Most ADAMTS functions first emerged from analysis of spontaneous human and animal mutations and genetically engineered animals. The clinical manifestations of Mendelian disorders resulting from mutations in ADAMTS2, ADAMTS10, ADAMTS13, ADAMTS17, ADAMTSL2 and ADAMTSL4 identified essential roles for each gene, but also suggested potential cooperative functions of ADAMTS proteins. These observations were extended by analysis of spontaneous animal mutations, such as in bovine ADAMTS2, canine ADAMTS10, ADAMTS17 and ADAMTSL2 and mouse ADAMTS20. These human and animal disorders are recessive and their manifestations appear to result from a loss-of-function mechanism. Genome-wide analyses have determined an association of some ADAMTS loci such as ADAMTS9 and ADAMTS7, with specific traits and acquired disorders. Analysis of genetically engineered rodent mutations, now achieved for over half the superfamily, has provided novel biological insights and animal models for the respective human genetic disorders and suggested potential candidate genes for related human phenotypes. Engineered mouse mutants have been interbred to generate combinatorial mutants, uncovering cooperative functions of ADAMTS proteins in morphogenesis. Specific genetic models have provided crucial insights on mechanisms of osteoarthritis (OA), a common adult-onset degenerative condition. Engineered mutants will facilitate interpretation of exome variants identified in isolated birth defects and rare genetic conditions, as well as in genome-wide screens for trait and disease associations. Mammalian forward and reverse genetics, together with genome-wide analysis, together constitute a powerful force for revealing the functions of ADAMTS proteins in physiological pathways and health disorders. Their continuing use, together with genome-editing technology and the ability to generate stem cells from mutants, presents numerous opportunities for advancing basic knowledge, human disease pathways and therapy.