Genome-wide pathway analysis of a genome-wide association study on psoriasis and Behcet’s disease

The aim of this study was to identify candidate causal single nucleotide polymorphisms (SNPs) and candidate causal mechanisms of psoriasis and Behcets’s disease (BD) and to generate an SNP → gene → pathway hypothesis. A psoriasis genome-wide association study (GWAS) dataset that included 436,192 SNPs in 1,409 psoriasis cases and 1,436 controls of European descent and a BD GWAS dataset that contained 310,324 SNPs in 1,215 BD cases and 1,278 controls were used in this study. Identify candidate causal SNPs and pathways (ICSNPathway) analysis was applied to the GWAS datasets. ICSNPathway analysis identified 15 candidate causal SNPs and 28 candidate causal pathways. The top five candidate causal SNPs were rs1063478 (P = 1.45E−10), rs8084 (P = 2.20E−08), rs7192 (P = 5.18E−08), rs20541 (P = 5.30E−06), and rs1130838 (P = 5.65E−06), which with the exception of rs20541 [interleukin (IL)-13] are at human leukocyte antigen (HLA) loci. These candidate causal SNPs and pathways provided ten hypothetical biological mechanisms. The most strongly associated pathway concerned HLA. When HLA loci were excluded, ICSNPathway analysis provided one hypothetical biological mechanism. rs20541 (non_synonymous_coding) → IL-13 → dendritic cell involvement in the regulation of Th1 and Th2 development, and the GATA3 pathway. ICSNPathway analysis identified four candidate causal SNPs, eleven candidate causal pathways, and three hypothetical biological mechanisms. One of them was as follows: rs2072895 (non_synonymous_coding & splice-site) and rs2735059 (non_synonymous_coding) → HLA-F → type I diabetes mellitus, antigen processing and presentation, and autoimmune thyroid disease. The application of ICSNPathway analysis to GWAS dataset of psoriasis and BD resulted in the identification of candidate causal SNPs and candidate pathways that might contribute to psoriasis susceptibility.     

Evaluating the effective numbers of independent tests and significant <Emphasis Type="Italic">p</Emphasis>-value thresholds in commercial genotyping arrays and public imputation reference datasets

Current genome-wide association studies (GWAS) use commercial genotyping microarrays that can assay over a million single nucleotide polymorphisms (SNPs). The number of SNPs is further boosted by advanced statistical genotype-imputation algorithms and large SNP databases for reference human populations. The testing of a huge number of SNPs needs to be taken into account in the interpretation of statistical significance in such genome-wide studies, but this is complicated by the non-independence of SNPs because of linkage disequilibrium (LD). Several previous groups have proposed the use of the effective number of independent markers (M _e) for the adjustment of multiple testing, but current methods of calculation for M _e are limited in accuracy or computational speed. Here, we report a more robust and fast method to calculate M _e. Applying this efficient method [implemented in a free software tool named Genetic type 1 error calculator (GEC)], we systematically examined the M _e, and the corresponding p-value thresholds required to control the genome-wide type 1 error rate at 0.05, for 13 Illumina or Affymetrix genotyping arrays, as well as for HapMap Project and 1000 Genomes Project datasets which are widely used in genotype imputation as reference panels. Our results suggested the use of a p-value threshold of ~10⁻⁷ as the criterion for genome-wide significance for early commercial genotyping arrays, but slightly more stringent p-value thresholds ~5 × 10⁻⁸ for current or merged commercial genotyping arrays, ~10⁻⁸ for all common SNPs in the 1000 Genomes Project dataset and ~5 × 10⁻⁸ for the common SNPs only within genes.     

MHC region and risk of systemic lupus erythematosus in African American women

The major histocompatibility complex (MHC) on chromosome 6p21 is a key contributor to the genetic basis of systemic lupus erythematosus (SLE). Although SLE affects African Americans disproportionately compared to European Americans, there has been no comprehensive analysis of the MHC region in relationship to SLE in African Americans. We conducted a screening of the MHC region for 1,536 single nucleotide polymorphisms (SNPs) and the deletion of the C4A gene in a SLE case–control study (380 cases, 765 age-matched controls) nested within the prospective Black Women’s Health Study. We also genotyped 1,509 ancestral informative markers throughout the genome to estimate European ancestry to control for population stratification due to population admixture. The most strongly associated SNP with SLE was the rs9271366 (odds ratio, OR = 1.70, p = 5.6 × 10⁻⁵) near the HLA-DRB1 gene. Conditional haplotype analysis revealed three other SNPs, rs204890 (OR = 1.86, p = 1.2 × 10⁻⁴), rs2071349 (OR = 1.53, p = 1.0 × 10⁻³), and rs2844580 (OR = 1.43, p = 1.3 × 10⁻³), to be associated with SLE independent of the rs9271366 SNP. In univariate analysis, the OR for the C4A deletion was 1.38, p = 0.075, but after simultaneous adjustment for the other four SNPs the odds ratio was 1.01, p = 0.98. A genotype score combining the four newly identified SNPs showed an additive risk according to the number of high-risk alleles (OR = 1.67 per high-risk allele, p < 0.0001). Our strongest signal, the rs9271366 SNP, was also associated with higher risk of SLE in a previous Chinese genome-wide association study (GWAS). In addition, two SNPs found in a GWAS of European ancestry women were confirmed in our study, indicating that African Americans share some genetic risk factors for SLE with European and Chinese subjects. In summary, we found four independent signals in the MHC region associated with risk of SLE in African American women.     

Genomewide association study for susceptibility genes contributing to familial Parkinson disease

Five genes have been identified that contribute to Mendelian forms of Parkinson disease (PD); however, mutations have been found in fewer than 5% of patients, suggesting that additional genes contribute to disease risk. Unlike previous studies that focused primarily on sporadic PD, we have performed the first genomewide association study (GWAS) in familial PD. Genotyping was performed with the Illumina HumanCNV370Duo array in 857 familial PD cases and 867 controls. A logistic model was employed to test for association under additive and recessive modes of inheritance after adjusting for gender and age. No result met genomewide significance based on a conservative Bonferroni correction. The strongest association result was with SNPs in the GAK/DGKQ region on chromosome 4 (additive model: p = 3.4 × 10⁻⁶; OR = 1.69). Consistent evidence of association was also observed to the chromosomal regions containing SNCA (additive model: p = 5.5 × 10⁻⁵; OR = 1.35) and MAPT (recessive model: p = 2.0 × 10⁻⁵; OR = 0.56). Both of these genes have been implicated previously in PD susceptibility; however, neither was identified in previous GWAS studies of PD. Meta-analysis was performed using data from a previous case–control GWAS, and yielded improved p values for several regions, including GAK/DGKQ (additive model: p = 2.5 × 10⁻⁷) and the MAPT region (recessive model: p = 9.8 × 10⁻⁶; additive model: p = 4.8 × 10⁻⁵). These data suggest the identification of new susceptibility alleles for PD in the GAK/DGKQ region, and also provide further support for the role of SNCA and MAPT in PD susceptibility. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. N. Pankratz and J. B. Wilk are joint first authors.     

Candidate genes for non-diabetic ESRD in African Americans: a genome-wide association study using pooled DNA

African Americans have increased susceptibility to non-diabetic (non-DM) forms of end-stage renal disease (ESRD) and extensive evidence supports a genetic contribution. A genome-wide association study (GWAS) using pooled DNA was performed in 1,000 African Americans to detect associated genes. DNA from 500 non-DM ESRD cases and 500 non-nephropathy controls was quantified using gel electrophoresis and spectrophotometric analysis and pools of 50 case and 50 control DNA samples were created. DNA pools were genotyped in duplicate on the Illumina HumanHap550-Duo BeadChip. Normalization methods were developed and applied to array intensity values to reduce inter-array variance. Allele frequencies were calculated from normalized channel intensities and compared between case and control pools. Three SNPs had p values of <1.0E−6: rs4462445 (ch 13), rs4821469 (ch 22) and rs8077346 (ch 17). After normalization, top scoring SNPs (n = 65) were genotyped individually in 464 of the original cases and 478 of the controls, with replication in 336 non-DM ESRD cases and 363 non-nephropathy controls. Sixteen SNPs were associated with non-DM ESRD (p < 7.7E−4, Bonferroni corrected). Twelve of these SNPs are in or near the MYH9 gene. The four non-MYH9 SNPs that were associated with non-DM ESRD in the pooled samples were not associated in the replication set. Five SNPs that were modestly associated in the pooled samples were more strongly associated in the replication and/or combined samples. This is the first GWAS for non-DM ESRD in African Americans using pooled DNA. We demonstrate strong association between non-DM ESRD in African Americans with MYH9, and have identified additional candidate loci.     

Identification, characterisation and application of single nucleotide polymorphisms for diversity assessment in cassava (Manihot esculenta Crantz)

To monitor genetic diversity in the field it is important that it is measured accurately. Here, we elucidate the potential of single nucleotide polymorphisms (SNPs) for measuring genetic diversity in cassava. The nature and frequency of SNPs was characterised and their utility in genetic diversity assessment compared to that of simple sequence repeats (SSRs). This was achieved by direct sequencing of amplicons in diverse cassava varieties. A total of 26 SNPs were identified from quality sequences of nine genes, giving an estimated frequency of one SNP every 121 nucleotides. Nucleotide diversity ranged from 7.8 × 10⁻⁴ to 5.6 × 10⁻³. Average haplotype-based polymorphic information content (PIC = 0.414) was higher than for individual SNPs (PIC = 0.228). The Mantel test indicated interdependence (r = 0.219; P < 0.001) between SNP and SSR genotypic data. Individual SNPs had lower PIC values than SSRs. For this reason larger numbers of SNPs may be necessary to achieve the same level of discrimination among genotypes provided by SSRs.     

Replication study of novel risk variants in six genes with type 2 diabetes and related quantitative traits in the Han Chinese lean individuals

To replicating the associations of type 2 diabetes (T2D) and six novel reported variants in Han Chinese lean individuals of first episode T2D, a total of six high risk single nucleotide polymorphisms (SNPs) from the BCL11A, DUSP9, IRS1, CENTD2, ADRA2A, and CDKAL1 genes were examined. Candidate six SNPs were genotyped in 761 T2D patients and 433 control subjects, and associations between the six SNPs and Body Mass Index (BMI), Fasting Plasma Glucose (FPG) and Two Hours Oral Glucose Tolerance Test (2hOGTT) were also investigated. CDKAL1 provided the strongest evidence for replication, where rs7754840 was associated with T2D (odds ratio = 1.54, per copy of the risk C allele, P = 8.10 × 10⁻⁷). SNP rs5945326 at DUSP9 showed modest significance (odds ratio = 0.81, per copy of the protective G allele, P = 0.02). After adjusting the confounders of age, gender and BMI, the above results remain significant for both rs7754840 (P < 1.0 × 10⁻⁴) and rs5945326 (P = 0.043) respectively. After correcting for multiple testing, however, only the association between T2D and rs7754840 at CDKAL1 (P < 1×10⁻⁴) remains significant. In addition, the risk C allele of CDKAL1 rs7754840 was significantly associated with increased FPG levels (P = 3.8 × 10⁻⁴). The association between genetic variant in CDKAL1 gene was detected in the Han Chinese lean individuals. The correlation between rs7754840-C allele and increased FPG levels is consistent with the potential function of CDKAL1 gene in pancreatic islets.     

A genome‐wide association study suggests new candidate genes for milk production traits in Chinese Holstein cattle

A genome‐wide association study (GWAS) was conducted on 15 milk production traits in Chinese Holstein. The experimental population consisted of 445 cattle, each genotyped by the GGP (GeneSeek genomic profiling)‐BovineLD V3 SNP chip, which had 26 151 public SNPs in its manifest file. After data cleaning, 20 326 SNPs were retained for the GWAS. The phenotypes were estimated breeding values of traits, provided by a public dairy herd improvement program center that had been collected once a month for 3 years. Two statistical models, a fixed‐effect linear regression model and a mixed‐effect linear model, were used to estimate the association effects of SNPs on each of the phenotypes. Genome‐wide significant and suggestive thresholds were set at 2.46E‐06 and 4.95E‐05 respectively. The two statistical models concurrently identified two genome‐wide significant (P < 0.05) SNPs on milk production traits in this Chinese Holstein population. The positional candidate genes, which were the ones closest to these two identified SNPs, were EEF2K (eukaryotic elongation factor 2 kinase) and KLHL1 (kelch like family member 1). These two genes could serve as new candidate genes for milk yield and lactation persistence, yet their roles need to be verified in further function studies.     

A genome‐wide association study suggests several novel candidate genes for carcass traits in Chinese Simmental beef cattle

A genome‐wide association study (GWAS) was conducted for two carcass traits in Chinese Simmental beef cattle. The experimental population consisted of 1301 individuals genotyped with the Illumina BovineHD SNP BeadChip (770K). After quality control, 671 990 SNPs and 1217 individuals were retained for the GWAS. The phenotypic traits included carcass weight and bone weight, which were measured after the cattle were slaughtered at 16 to 18 months of age. Three statistical models—a fixed polygene model, a random polygene model and a composite interval mapping polygene model—were used for the GWAS. The genome‐wide significance threshold after Bonferroni correction was 7.44E‐08 (= 0.05/671 990). In this study, we detected eight and seven SNPs significantly associated with carcass weight and bone weight respectively. In total, 11 candidate genes were identified within or close to these significant SNPs. Of these, we found several novel candidate genes, including PBX1, GCNT4, ALDH1A2, LCORL and WDFY3, to be associated with carcass weight and bone weight in Chinese Simmental beef cattle, and their functional roles need to be verified in further studies.     

A risk-associated single nucleotide polymorphism of <Emphasis Type="Italic">SMAD7</Emphasis> is common to colorectal,gastric, and lung cancers in a Han Chinese population

SMAD7 has been demonstrated to antagonize TGF-β-mediated fibrosis, carcinogenesis, and inflammation. Two previous genome-wide association studies identified three single nucleotide polymorphisms (SNPs) (rs4939827, rs12953717 and rs4464148) in SMAD7 to be associated with colorectal cancer in a Western population. We conducted the first case–control study in a Han Chinese population to explore the associations between these three SNPs and colorectal, gastric, and lung cancers. Of the three SNPs, only rs12953717 was strongly associated with the three types of cancer, fitting the overdominant model. Compared with the CC/TT (CC combined with TT) genotype, the adjusted odds ratios for the CT genotype were 2.002 (95% CI, 1.250–3.207, P = 0.004), 1.678 (95% CI, 1.048–2.689, P = 0.031), 3.825 (95% CI, 2.310–6.335, P < 1 × 10⁻⁴), and 2.294 (95% CI, 1.537–3.343, P < 1 × 10⁻⁴), respectively, for colorectal, gastric, lung, and combined cancers. These outcomes suggest that rs12953717 is a common risk marker of these three types of cancer in the Han Chinese.     

EST-derived single nucleotide polymorphism markers for assembling genetic and physical maps of the barley genome

In a panel of seven genotypes, 437 expressed sequence tag (EST)-derived DNA fragments were sequenced. Single nucleotide polymorphisms (SNPs) that were polymorphic between the parents of three mapping populations were mapped by heteroduplex analysis and a genome-wide consensus map comprising 216 EST-derived SNPs and 4 InDel (insertion/deletion) markers was constructed. The average frequency of SNPs amounted to 1/130 bp and 1/107.8 bp for a set of randomly selected and a set of mapped ESTs, respectively. The calculated nucleotide diversities (π) ranged from 0 to 40.0 × 10⁻³ (average 3.1 × 10⁻³) and 0.52 × 10⁻³ to 39.51 × 10^–3 (average 4.37 × 10⁻³) for random and mapped ESTs, respectively. The polymorphism information content value for mapped SNPs ranged from 0.24 to 0.50 with an average of 0.34. As expected, combination of SNPs present in an amplicon (haplotype) exhibited a higher information content ranging from 0.24 to 0.85 with an average of 0.50. Cleaved amplified polymorphic sequence assays (including InDels) were designed for a total of 87 (39.5%) SNP markers. The high abundance of SNPs in the barley genome provides avenues for the systematic development of saturated genetic maps and their integration with physical maps. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Both R. Kota and R.K. Varshney contributed equally to this work.     

Genome-wide association study for growth and fatness traits in Chinese Sujiang pigs

Growth and fatness traits are complex and economically important traits in the pig industry. The molecular basis underlying porcine growth and fatness traits remains largely unknown. To uncover genetic loci and candidate genes for these traits, we explored the GeneSeek GGP Porcine 80K SNP chip to perform a GWAS for seven growth and fatness traits in 365 individuals from the Sujiang pig, a recently developed breed in China. We identified two, 17, one and 11 SNPs surpassing the suggestively significant threshold (P < 1.86 × 10⁻⁵) for body weight, chest circumference, chest width and backfat thickness respectively. Of these SNPs, 20 represent novel genetic loci, and five and four SNPs were respectively associated with chest circumference and backfat thickness at a genome-wide significant threshold (P < 9.31 × 10⁻⁷). Eight SNPs had a pleiotropic effect on both chest circumference and backfat thickness. The most remarkable locus resided in a region between 72.95 and 76.27 Mb on pig chromosome 4, harboring a number of previously reported quantitative trait loci related to backfat deposition. In addition to two reported genes (PLAG1 and TAS2R38), we identified four genes including GABRB3, ZNF106, XKR4 and MGAM as novel candidates for body weight and backfat thickness at the mapped loci. Our findings provide insights into the genetic architecture of porcine growth and fatness traits and potential markers for selective breeding of Chinese Sujiang pigs.     

The environmental signals of stable carbon isotope in various tree-ring components of Pinus tabulaeformis

Chinese pine (Pinus tabulaeformis) trees were sampled in the Helan Mountain, northwest China. The stable carbon isotope (δ¹³C) values of whole wood, holocellulose and alpha-cellulose in tree rings over 30 years (1968–1997) were measured to study the δ¹³C response of different tree-ring components to past environmental change. There were obvious differences in the δ¹³C values of the three components. The Pearson correlation coefficient between the δ¹³C of alpha-cellulose and that of holocellulose was 0.547 (ρ < 0.01); between alpha-cellulose and whole wood, the coefficient was −0.126 (ρ > 0.10); between holocellulose and whole wood, the coefficient was −0.056. Correlation function analyses indicated that the δ¹³C content of tree-ring alpha-cellulose correlated strongly with the average temperature from June to August (r = 0.427, ρ < 0.05), more than that of holocellulose (0.324, ρ < 0.10) or total wood (−0.245, ρ > 0.10). Significant correlations were observed between δ¹³C of tree-ring alpha-cellulose and the precipitation from the current year’s February to July (r = −0.514, ρ < 0.01) that were much higher than that of holocellulose (−0.481, ρ < 0.05) or total wood (−0.249, ρ > 0.10). A significant correlation (−0.545, ρ < 0.01) was also found between the ring width and the δ¹³C residual chronologies. These results suggest that more past environmental information is retained in the δ¹³C of tree-ring alpha-cellulose. Thus, the δ¹³C of alpha-cellulose of tree rings is the most suitable among the studied parameters for reconstructing the past climatic conditions during the growing season. The δ¹³C values of other organic compounds in Pinus tabulaeformis xylem were affected by the external environment after carbon was fixed from the atmosphere.     

Genome‐wide association study for the level of serum electrolytes in Italian Large White pigs

Calcium, magnesium and phosphorus are essential electrolytes involved in a large number of biological processes. Imbalance of these minerals in blood may indicate clinically relevant conditions and are important in inferring acute or chronic pathologies in humans and animals. In this work, we carried out a genome‐wide association study (GWAS) for the level of these three electrolytes in the serum of 843 performance‐tested Italian Large White pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip, and GWAS was carried out using genome‐wide efficient mixed‐model association. For the level of Ca²⁺, eight single nucleotide polymorphisms (SNPs) were significant, considering a false discovery rate (FDR) < 0.05, and another eight were above the moderate association threshold (P_{nominal value} < 5.00E‐05). These SNPs are distributed in four porcine chromosomes (SSC): SSC8, SSC11, SSC12 and SSC13. In particular, a few putative different signals of association detected on SSC13 and one on SSC12 were in genes or close to genes involved in calcium metabolism (P2RY1, RAP2B, SLC9A9, C3orf58, TSC22D2, PLCH1 and CACNB1). Only one SNP (on SSC7) and six SNPs (on SSC2 and SSC7) showed moderate association with the level of magnesium and phosphorus respectively. The association signals for these two latter minerals might identify genes not known thus far for playing a role in their biological functions and regulations. In conclusion, our GWAS contributed to increased knowledge on the role that calcium, magnesium and phosphorus may play in the genetically determined physiological mechanisms affecting the natural variability of mineral levels in mammalian blood.     

Continuous cultures of <Emphasis Type="Italic">Pseudomonas putida</Emphasis> mt-2 overcome catabolic function loss under real case operating conditions

The long-term performance and stability of Pseudomonas putida mt-2 cultures, a toluene-sensitive strain harboring the genes responsible for toluene biodegradation in the archetypal plasmid pWW0, was investigated in a chemostat bioreactor functioning under real case operating conditions. The process was operated at a dilution rate of 0.1 h⁻¹ under toluene loading rates of 259 ± 23 and 801 ± 78 g m⁻³ h⁻¹ (inlet toluene concentrations of 3.5 and 10.9 g m⁻³, respectively). Despite the deleterious effects of toluene and its degradation intermediates, the phenotype of this sensitive P. putida culture rapidly recovered from a 95% Tol⁻ population at day 4 to approx. 100% Tol⁺ cells from day 13 onward, sustaining elimination capacities of 232 ± 10 g m⁻³ h⁻¹ at 3.5 g Tol m⁻³ and 377 ± 13 g m⁻³ h⁻¹ at 10.9 g Tol m⁻³, which were comparable to those achieved by highly tolerant strains such as P. putida DOT T1E and P. putida F1 under identical experimental conditions. Only one type of Tol⁻ variant, harboring a TOL-like plasmid with a 38.5 kb deletion (containing the upper and meta operons for toluene biodegradation), was identified.     

Genotype frequency and F ST analysis of polymorphisms in immunoregulatory genes in Chinese and Caucasian populations

Selection and genetic drift can create genetic differences between populations. Cytokines and chemokines play an important role in both hematopoietic development and the inflammatory response. We compared the genotype frequencies of 45 SNPs in 30 cytokine and chemokine genes in two healthy Chinese populations and one Caucasian population. Several SNPs in IL4 had substantial genetic differentiation between the Chinese and Caucasian populations (F _ST ~0.40), and displayed a strikingly different haplotype distribution. To further characterize common genetic variation in worldwide populations at the IL4 locus, we genotyped 9 SNPs at the IL4 gene in the Human Diversity Panel’s (N = 1056) individuals from 52 world geographic regions. We observed low haplotype diversity, yet strikingly different haplotype frequencies between non-African populations, which may indicate different selective pressures on the IL4 gene in different parts of the world. SNPs in CSF2, IL6, IL10, CTLA4, and CX3CR1 showed moderate genetic differentiation between the Chinese and Caucasian populations (0.15 < F _ST < 0.25). These results suggest that there is substantial genetic diversity in immune genes and exploration of SNP associations with immune-related diseases that vary in incidence across these two populations may be warranted.     

Genome‐wide association study for pigmentation traits in Chinese Holstein population

With the Illumina BovineSNP50K BeadChip, we performed a genome‐wide association study (GWAS) for two pigmentation traits in a Chinese Holstein population: proportion of black (PB) and teat colour (TC). A case–control design was used. Cases were the cows with PB <0.30 (n = 129) and TC <2 points (n = 140); controls were those with PB >0.90 (n = 58) and TC >4 points (n = 281). The RM test of roadtrips (version 1.2) was applied to detect SNPs for the two traits with 42 883 and 42 741 SNPs respectively. A total of nine and 12 genome‐wide significant (P < 0.05) SNPs associated with PB and TC respectively were identified. Of these, two SNPs for PB were located within the KIT and IGFBP7 genes, and the other four SNPs were 23~212 kb away from the PDGFRA gene on BTA6; nine SNPs associated with TC were located within or 21~78.8 kb away from known genes on chromosomes 4, 11, 22, 23 and 24. By combing through our GWAS results and the biological functions of the genes, we suggest that the KIT, IGFBP7, PDGFRA, MITF, ING3 and WNT16 genes are promising candidates for PB and TC in Holstein cattle, providing a basis for further investigation on the genetic mechanism of pigmentation formation.     

Genome-Wide Association Analyses Identify SPOCK as a Key Novel Gene Underlying Age at Menarche

For females, menarche is a most significant physiological event. Age at menarche (AAM) is a trait with high genetic determination and is associated with major complex diseases in women. However, specific genes for AAM variation are largely unknown. To identify genetic factors underlying AAM variation, a genome-wide association study (GWAS) examining about 380,000 SNPs was conducted in 477 Caucasian women. A follow-up replication study was performed to validate our major GWAS findings using two independent Caucasian cohorts with 854 siblings and 762 unrelated subjects, respectively, and one Chinese cohort of 1,387 unrelated subjects—all females. Our GWAS identified a novel gene, SPOCK (Sparc/Osteonectin, CWCV, and Kazal-like domains proteoglycan), which had seven SNPs associated with AAM with genome-wide false discovery rate (FDR) q<0.05. Six most significant SNPs of the gene were selected for validation in three independent replication cohorts. All of the six SNPs were replicated in at least one cohort. In particular, SNPs rs13357391 and rs1859345 were replicated both within and across different ethnic groups in all three cohorts, with p values of 5.09×10⁻³ and 4.37×10⁻³, respectively, in the Chinese cohort and combined p values (obtained by Fisher's method) of 5.19×10⁻⁵ and 1.02×10⁻⁴, respectively, in all three replication cohorts. Interestingly, SPOCK can inhibit activation of MMP-2 (matrix metalloproteinase-2), a key factor promoting endometrial menstrual breakdown and onset of menstrual bleeding. Our findings, together with the functional relevance, strongly supported that the SPOCK gene underlies variation of AAM.     

Calbindin 1, fibroblast growth factor 20, and α-synuclein in sporadic Parkinson’s disease

Parkinson’s disease (PD), one of the most common human neurodegenerative disorders, is characterized by the loss of dopaminergic neurons in the substantia nigra of the midbrain. Our recent case-control association study of 268 SNPs in 121 candidate genes identified α-synuclein (SNCA) as a susceptibility gene for sporadic PD (P = 1.7 × 10⁻¹¹). We also replicated the association of fibroblast growth factor 20 (FGF20) with PD (P = 0.0089). To find other susceptibility genes, we added 34 SNPs to the previous screen. Of 302 SNPs in a total 137 genes, but excluding SNCA, SNPs in NDUFV2, FGF2, CALB1 and B2M showed significant association (P < 0.01; 882 cases and 938 control subjects). We replicated the association analysis for these SNPs in a second independent sample set (521 cases and 1,003 control subjects). One SNP, rs1805874 in calbindin 1 (CALB1), showed significance in both analyses (P = 7.1 × 10⁻⁵; recessive model). When the analysis was stratified relative to the SNCA genotype, the odds ratio of CALB1 tended to increase according to the number of protective alleles in SNCA. In contrast, FGF20 was significant only in the subgroup of SNCA homozygote of risk allele. CALB1 is a calcium-binding protein that widely is expressed in neurons. A relative sparing of CALB1-positive dopaminergic neurons is observed in PD brains, compared with CALB1-negative neurons. Our genetic analysis suggests that CALB1 is associated with PD independently of SNCA, and that FGF20 is associated with PD synergistically with SNCA. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Association between variants in the 5′-untranslated region of the bovine <Emphasis Type="Italic">MC4R</Emphasis> gene and two growth traits in Nanyang cattle

Melanocortin-4 receptor (MC4R) is one of five G-protein-coupled receptors binding melanocortins that is implicated in the control of feeding behavior and energy homeostasis. Six cattle populations (n = 594), including four Chinese indigenous breeds, Chinese Holstein, and a meat type breed (Angus), were used to detect single nucleotide polymorphisms in 5′-untranslated region of MC4R gene by means of PCR–SSCP and DNA sequencing. Four linked SNPs (g.[−293C>G; −193A>T; −192T>G; −129A>G]) were identified. The g.−293C>G and g.−129A>G could be genotyped with a PCR–RFLP using TaiI in three combined genotypes (AA, AB and BB). The two linked SNPs were associated with body weight and daily gain in Nanyang aged 6 months (P < 0.05), but they had no significant effect on body weight and daily gain in Nanyang aged 24 months (P > 0.05).