Common sequence variants of the macrophage scavenger receptor 1 gene are associated with prostate cancer risk

Rare germline mutations of macrophage scavenger receptor 1 (MSR1) gene were reported to be associated with prostate cancer risk in families with hereditary prostate cancer (HPC) and in patients with non-HPC (Xu et al. 2002). To further evaluate the role of MSR1 in prostate cancer susceptibility, at Johns Hopkins Hospital, we studied five common variants of MSR1 in 301 patients with non-HPC who underwent prostate cancer treatment and in 250 control subjects who participated in prostate cancer-screening programs and had normal digital rectal examination and PSA levels (<4 ng/ml). Significantly different allele frequencies between case subjects and control subjects were observed for each of the five variants (P value range.01-.04). Haplotype analyses provided consistent findings, with a significant difference in the haplotype frequencies from a global score test (P=.01). Because the haplotype that is associated with the increased risk for prostate cancer did not harbor any of the known rare mutations, it appears that the observed association of common variants and prostate cancer risk are independent of the effect of the known rare mutations. These results consistently suggest that MSR1 may play an important role in prostate carcinogenesis.     

Nucleotide variation,haplotype structure,and association with end-stage renal disease of the human interleukin-1 gene cluster

A dense gene-based SNP map was constructed across a 360-kb region containing the interleukin-1 gene cluster (IL1A, IL1B, and IL1RN), focusing on IL1RN. In total, 95 polymorphisms were confirmed or identified primarily by direct sequencing. Polymorphisms were precisely mapped to completed BAC and genomic sequences spanning this region. The polymorphisms were typed in 443 case-control subjects from Caucasian and African American groups. Consecutive pair-wise marker linkage disequilibrium was not strictly correlated with distance and ranged from D'=0.0079 to 1.000 and D'=0.0521 to 1.0000 in Caucasians and African Americans, respectively. Single markers and haplotypes in IL1 cluster genes were evaluated for association with end-stage renal disease (ESRD). Eleven SNPs show some evidence of association with ESRD, with the strongest associations in two IL1A variants, one SNP, rs1516792-3, in intron 5 (p=0.0015) and a 4-bp insertion/deletion within the 3'UTR, rs16347-2 (p=0.0024), among African Americans with non-T2DM-associated ESRD.     

Pseudomonas aeruginosa flagellin and alginate elicit very distinct gene expression patterns in airway epithelial cells: implications for cystic fibrosis disease

Infection with the opportunistic pathogen Pseudomonas aeruginosa remains a major health concern. Two P. aeruginosa phenotypes relevant in human disease include motility and mucoidy. Motility is characterized by the presence of flagella and is essential in the establishment of acute infections, while mucoidy, defined by the production of the exopolysaccharide alginate, is critical in the development of chronic infections, such as the infections seen in cystic fibrosis patients. Indeed, chronic infection of the lung by mucoid P. aeruginosa is a major cause of morbidity and mortality in cystic fibrosis patients. We have used Calu-3 human airway epithelial cells to investigate global responses to infection with motile and mucoid P. aeruginosa. The response of airway epithelial cells to exposure to P. aeruginosa motile strains is characterized by a specific increase in gene expression in pathways controlling inflammation and host defense. By contrast, the response of airway epithelia to the stimuli presented by mucoid P. aeruginosa is not proinflammatory and, hence, may not be conducive to the effective elimination of the pathogen. The pattern of gene expression directed by flagellin, but not alginate, includes innate host defense genes, proinflammatory cytokines, and chemokines. By contrast, infection with alginate-producing P. aeruginosa results in an overall attenuation of host responses and an antiapoptotic effect.     

Population structure of Hispanics in the United States: the multi-ethnic study of atherosclerosis

Using ∼60,000 SNPs selected for minimal linkage disequilibrium, we perform population structure analysis of 1,374 unrelated Hispanic individuals from the Multi-Ethnic Study of Atherosclerosis (MESA), with self-identification corresponding to Central America (n = 93), Cuba (n = 50), the Dominican Republic (n = 203), Mexico (n = 708), Puerto Rico (n = 192), and South America (n = 111). By projection of principal components (PCs) of ancestry to samples from the HapMap phase III and the Human Genome Diversity Panel (HGDP), we show the first two PCs quantify the Caucasian, African, and Native American origins, while the third and fourth PCs bring out an axis that aligns with known South-to-North geographic location of HGDP Native American samples and further separates MESA Mexican versus Central/South American samples along the same axis. Using k-means clustering computed from the first four PCs, we define four subgroups of the MESA Hispanic cohort that show close agreement with self-identification, labeling the clusters as primarily Dominican/Cuban, Mexican, Central/South American, and Puerto Rican. To demonstrate our recommendations for genetic analysis in the MESA Hispanic cohort, we present pooled and stratified association analysis of triglycerides for selected SNPs in the LPL and TRIB1 gene regions, previously reported in GWAS of triglycerides in Caucasians but as yet unconfirmed in Hispanic populations. We report statistically significant evidence for genetic association in both genes, and we further demonstrate the importance of considering population substructure and genetic heterogeneity in genetic association studies performed in the United States Hispanic population.     

Expression profiling reveals differential gene induction underlying specific and non-specific memory for pheromones in mice

Memory for the mating male’s pheromones in female mice is thought to require synaptic changes in the accessory olfactory bulb (AOB). Induction of this memory depends on release of glutamate in response to pheromonal exposure coincident with release of norepinephrine (NE) in the AOB following mating. A similar memory for pheromones can also be induced artificially by local infusion of the GABA_A receptor antagonist bicuculline into the AOB. The natural memory formed by exposure to pheromones during mating is specific to the pheromones sensed by the female during mating. In contrast, the artificial memory induced by bicuculline is non-specific and results in the female mice recognizing all pheromones as if they were from the mating male. Although protein synthesis has been shown to be essential for development of pheromone memory, the gene expression cascades critical for memory formation are not known. We investigated changes in gene expression in the AOB using oligonucleotide microarrays during mating-induced pheromone memory (MIPM) as well as bicuculline-induced pheromone memory (BIPM). We found the set of genes induced during MIPM and BIPM are largely non-overlapping and Ingenuity Pathway Analysis revealed that the signaling pathways in MIPM and BIPM also differ. The products of genes induced during MIPM are associated with synaptic function, indicating the possibility of modification at specific synapses, while those induced during BIPM appear to possess neuron-wide functions, which would be consistent with global cellular changes. Thus, these results begin to provide a mechanistic explanation for specific and non-specific memories induced by pheromones and bicuculline infusion respectively.     

The intracellular form of human MAGP1 elicits a complex and specific transcriptional response

Microfibril-associated glycoprotein-1 (MAGP1) is found associated with microfibrils in the extracellular matrix (ECM). In humans, MAGP1 is expressed as two alternatively spliced isoforms: MAGP1A, the extracellular microfibril-associated form; and MAGP1B, an exclusively intracellular isoform derived from the skipping of exon 3. The biological function of MAGP1B is unknown. We performed gene expression profiling to study the cellular response to MAGP1B using whole-genome genechips. We found that MAGP1B specifically induces the expression of genes linked to cell adhesion, motility, metabolism, gene expression, development and signal transduction. Versican, a gene product involved in the structure and functional regulation of the ECM, showed the highest up-regulation in response to MAGP1B. These studies suggest a dual role for MAGP1, with extracellular MAGP1A involved in ECM function, and intracellular MAGP1B modulating the expression of genes that function in cell adhesion, migration and control of ECM deposition.     

Family-Based Association Analysis Confirms the Role of the Chromosome 9q21.32 Locus in the Susceptibility of Diabetic Nephropathy

A genome-wide association scan of type 1 diabetic patients from the GoKinD collections previously identified four novel diabetic nephropathy susceptibility loci that have subsequently been shown to be associated with diabetic nephropathy in unrelated patients with type 2 diabetes. To expand these findings, we examined whether single nucleotide polymorphisms (SNPs) at these susceptibility loci were associated with diabetic nephropathy in patients from the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection. Six SNPs across the four loci identified in the GoKinD collections and 7 haplotype tagging SNPs, were genotyped in 66 extended families of European ancestry. Pedigrees from this collection contained an average of 18.5 members, including 2 to 14 members with type 2 diabetes. Among diabetic family members, the 9q21.32 locus approached statistical significance with advanced diabetic nephropathy (P = 0.037 [adjusted P = 0.222]). When we expanded our definition of diabetic nephropathy to include individuals with high microalbuminuria, the strength of this association improved significantly (P = 1.42×10⁻³ [adjusted P = 0.009]). This same locus also trended toward statistical significance with variation in urinary albumin excretion in family members with type 2 diabetes (P = 0.032 [adjusted P = 0.192]) and in analyses expanded to include all relatives (P = 0.019 [adjusted P = 0.114]). These data increase support that SNPs identified in the GoKinD collections on chromosome 9q21.32 are true diabetic nephropathy susceptibility loci.     

Analysis of family- and population-based samples in cohort genome-wide association studies

Cohort studies typically sample unrelated individuals from a population, although family members of index cases may also be recruited to investigate shared familial risk factors. Recruitment of family members may be incomplete or ancillary to the main cohort, resulting in a mixed sample of independent family units, including unrelated singletons and multiplex families. Multiple methods are available to perform genome-wide association (GWA) analysis of binary or continuous traits in families, but it is unclear whether methods known to perform well on ascertained pedigrees, sibships, or trios are appropriate in analysis of a mixed unrelated cohort and family sample. We present simulation studies based on Multi-Ethnic Study of Atherosclerosis (MESA) pedigree structures to compare the performance of several popular methods of GWA analysis for both quantitative and dichotomous traits in cohort studies. We evaluate approaches suitable for analysis of families, and combined the best performing methods with population-based samples either by meta-analysis, or by pooled analysis of family- and population-based samples (mega-analysis), comparing type 1 error and power. We further assess practical considerations, such as availability of software and ability to incorporate covariates in statistical modeling, and demonstrate our recommended approaches through quantitative and binary trait analysis of HDL cholesterol (HDL-C) in 2,553 MESA family- and population-based African-American samples. Our results suggest linear modeling approaches that accommodate family-induced phenotypic correlation (e.g., variance-component model for quantitative traits or generalized estimating equations for dichotomous traits) perform best in the context of combined family- and population-based cohort GWAS.     

Evaluation of a SNP map of 6q24–27 confirms diabetic nephropathy loci and identifies novel associations in type 2 diabetes patients with nephropathy from an African-American population

Previously, we performed a genome scan for type 2 diabetes (T2DM) using 638 African-American (AA) affected sibling pairs from 247 families; non-parametric linkage analysis suggested evidence of linkage at 6q24–27 (LOD 2.26). To comprehensively evaluate this region, we performed a two-stage association study by first constructing a SNP map of 754 SNPs selected from HapMap on the basis of linkage disequilibrium (LD) in 300 AAT2DM end-stage renal disease (ESRD) subjects, 311 AA controls, 43 European American controls and 45 Yoruba Nigerian samples (Set 1). Replication analyses were conducted in an independent population of 283 AA T2DM-ESRD subjects and 282 AA controls (Set 2). In addition, we adjusted for the impact of admixture on association results by using ancestry informative markers (AIMs). In Stage 1, 137 (18.2%) SNPs showed nominal evidence of association (P < 0.05) in one or more of tests of association: allelic (n = 33), dominant (n = 36), additive (n = 29), or recessive (n = 34) genotypic models, and 2- (n = 47) and 3-SNP (n = 43) haplotypic analyses. These SNPs were selected for follow-up genotyping. Stage 2 analyses confirmed association with a predicted 2-SNP “risk” haplotype in the PARK2 gene. Also, two intergenic SNPs showed consistent genotypic association with T2DM-ESRD: rs12197043 and rs4897081. Combined analysis of all subjects from both stages revealed nominal associations with 17 SNPs within genes, including suggestive associations in ESR1 and PARK2. This study confirms known diabetic nephropathy loci and identifies potentially novel susceptibility variants located within 6q24–27 in AA. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.