Characterization of novel microsatellite markers derived from Korean rose bitterling (Rhodeus uyekii) genomic library

Korean rose bitterling (Rhodeus uyekii) is a freshwater fish endemic to Korea. Natural populations of this species have experienced severe declines as a result of habitat fragmentation and water pollution. To conserve and restore R. uyekii, the genetic diversity of this species needs to be assessed at the population level. Eighteen novel polymorphic microsatellite loci for R. uyekii were developed using an enriched partial genomic library. Polymorphisms at these loci were studied in 150 individuals collected from three populations. The number of alleles at each locus ranged from 3 to 47 (mean = 17.1). Within the populations, the observed heterozygosity ranged from 0.032 to 1.000, expected heterozygosity from 0.082 to 0.967, and polymorphism information content from 0.078 to 0.950. Six loci showed significant deviation from Hardy-Weinberg equilibrium after Bonferroni’s correction, and no significant linkage disequilibrium was detected between most locus pairs, except in three cases. These highly informative microsatellite markers should be useful for genetic population structure analyses of R. uyekii.

Development of polymorphic microsatellite markers suitable for genetic linkage mapping of olive flounder Paralichthys olivaceus

Microsatellite markers are important for gene mapping and for marker-assisted selection. Sixty-five polymorphic microsatellite markers were developed with an enriched partial genomic library from olive flounder Paralichthys olivaceus an important commercial fish species in Korea. The variability of these markers was tested in 30 individuals collected from the East Sea (Korea). The number of alleles for each locus ranged from 2 to 33 (mean, 17.1). Observed and expected heterozygosity as well as polymorphism information content varied from 0.313 to 1.000 (mean, 0.788), from 0.323 to 0.977 (mean, 0.820), and from 0.277 to 0.960 (mean, 0.787), respectively. Nine loci showed significant deviation from the Hardy-Weinberg equilibrium after sequential Bonferroni correction. Analysis with MICROCHECKER suggested the presence of null alleles at five of these loci with estimated null allele frequencies of 0.126-0.285. These new microsatellite markers from genomic libraries will be useful for constructing a P. olivaceus linkage map.

Genome-wide association study identifies a variant in HDAC9 associated with large vessel ischemic stroke

Genetic factors have been implicated in stroke risk, but few replicated associations have been reported. We conducted a genome-wide association study (GWAS) for ischemic stroke and its subtypes in 3,548 affected individuals and 5,972 controls, all of European ancestry. Replication of potential signals was performed in 5,859 affected individuals and 6,281 controls. We replicated previous associations for cardioembolic stroke near PITX2 and ZFHX3 and for large vessel stroke at a 9p21 locus. We identified a new association for large vessel stroke within HDAC9 (encoding histone deacetylase 9) on chromosome 7p21.1 (including further replication in an additional 735 affected individuals and 28,583 controls) (rs11984041; combined P = 1.87 × 10 -11; odds ratio (OR) = 1.42, 95% confidence interval (CI) = 1.28-1.57). All four loci exhibited evidence for heterogeneity of effect across the stroke subtypes, with some and possibly all affecting risk for only one subtype. This suggests distinct genetic architectures for different stroke subtypes.

Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn's disease

Ankylosing spondylitis (AS) is a common inflammatory arthritic condition. Overt inflammatory bowel disease (IBD) occurs in about 10% of AS patients, and in addition 70% of AS cases may have subclinical terminal ileitis. Spondyloarthritis is also common in IBD patients. We therefore tested Crohn's disease susceptibility genes for association with AS, aiming to identify pleiotropic genetic associations with both diseases. Genotyping was carried out using Sequenom and Applied Biosystems TaqMan and OpenArray technologies on 53 markers selected from 30 Crohn's disease associated genomic regions. We tested genotypes in a population of unrelated individual cases (n = 2,773) and controls (n = 2,215) of white European ancestry for association with AS. Statistical analysis was carried out using a Cochran-Armitage test for trend in PLINK. Strong association was detected at chr1q32 near KIF21B (rs11584383, P = 1.66 x 10-10, odds ratio (OR) = 0.74, 95% CI:0.68-0.82). Association with disease was also detected for 2 variants within STAT3 (rs6503695, P = 4.6×10-4. OR = 0.86 (95% CI:0.79-0.93); rs744166, P = 2.6×10-5, OR = 0.84 (95% CI:0.77-0.91)). Association was confirmed for IL23R (rs11465804, P = 1.2×10-5, OR = 0.65 (95% CI:0.54-0.79)), and further associations were detected for IL12B (rs10045431, P = 5.261025, OR = 0.83 (95% CI:0.76-0.91)), CDKAL1 (rs6908425, P = 1.1×10-4, OR = 0.82 (95% CI:0.74-0.91)), LRRK2/MUC19 (rs11175593, P = 9.9×10-5, OR = 1.92 (95% CI: 1.38-2.67)), and chr13q14 (rs3764147, P = 5.9×10-4, OR = 1.19 (95% CI: 1.08-1.31)). Excluding cases with clinical IBD did not significantly affect these findings. This study identifies chr1q32 and STAT3 as ankylosing spondylitis susceptibility loci. It also further confirms association for IL23R and detects suggestive association with another 4 loci. STAT3 is a key signaling molecule within the Th17 lymphocyte differentiation pathway and further enhances the case for a major role of this T-lymphocyte subset in ankylosing spondylitis. Finally these findings suggest common aetiopathogenic pathways for AS and Crohn's disease and further highlight the involvement of common risk variants across multiple diseases.

N-ethyl-N-nitrosourea (ENU) induced mutations within the Klotho gene lead to ectopic calcification and reduced lifespan in mouse models

Ectopic calcification (EC), which is the pathological deposition of calcium and phosphate in extra-skeletal tissues, may be associated with hypercalcaemic and hyperphosphataemic disorders, or it may occur in the absence of metabolic abnormalities. In addition, EC may be inherited as part of several monogenic disorders and studies of these have provided valuable insights into the metabolic pathways regulating mineral metabolism. For example, studies of tumoural calcinosis, a disorder characterised by hyperphosphataemia and progressive EC, have revealed mutations of fibroblast growth factor 23 (FGF23), polypeptide N-acetyl galactosaminyltransferase 3 (GALNT3) and klotho (KL), which are all part of a phosphate-regulating pathway. However, such studies in humans are limited by the lack of available large families with EC, and to facilitate such studies we assessed the progeny of mice treated with the chemical mutagen N-ethyl-N-nitrosourea (ENU) for EC. This identified two mutants with autosomal recessive forms of EC, and reduced lifespan, designated Ecalc1 and Ecalc2. Genetic mapping localized the Ecalc1 and Ecalc2 loci to a 11.0 Mb region on chromosome 5 that contained the klotho gene (Kl), and DNA sequence analysis identified nonsense (Gln203Stop) and missense (Ile604Asn) Kl mutations in Ecalc1 and Ecalc2 mice, respectively. The Gln203Stop mutation, located in KL1 domain, was severely hypomorphic and led to a 17-fold reduction of renal Kl expression. The Ile604Asn mutation, located in KL2 domain, was predicted to impair klotho protein stability and in vitro expression studies in COS-7 cells revealed endoplasmic reticulum retention of the Ile604Asn mutant. Further phenotype studies undertaken in Ecalc1 (kl203X/203X) mice demonstrated elevations in plasma concentrations of phosphate, FGF23 and 1,25-dihydroxyvitamin D. Thus, two allelic variants of Kl that develop EC and represent mouse models for tumoural calcinosis have been established. © 2015 Esapa et al.

Identity-by-Descent Mapping to Detect Rare Variants Conferring Susceptibility to Multiple Sclerosis

Genome-wide association studies (GWAS) have identified around 60 common variants associated with multiple sclerosis (MS), but these loci only explain a fraction of the heritability of MS. Some missing heritability may be caused by rare variants that have been suggested to play an important role in the aetiology of complex diseases such as MS. However current genetic and statistical methods for detecting rare variants are expensive and time consuming. 'Population-based linkage analysis' (PBLA) or so called identity-by-descent (IBD) mapping is a novel way to detect rare variants in extant GWAS datasets. We employed BEAGLE fastIBD to search for rare MS variants utilising IBD mapping in a large GWAS dataset of 3,543 cases and 5,898 controls. We identified a genome-wide significant linkage signal on chromosome 19 (LOD = 4.65; p = 1.9×10-6). Network analysis of cases and controls sharing haplotypes on chromosome 19 further strengthened the association as there are more large networks of cases sharing haplotypes than controls. This linkage region includes a cluster of zinc finger genes of unknown function. Analysis of genome wide transcriptome data suggests that genes in this zinc finger cluster may be involved in very early developmental regulation of the CNS. Our study also indicates that BEAGLE fastIBD allowed identification of rare variants in large unrelated population with moderate computational intensity. Even with the development of whole-genome sequencing, IBD mapping still may be a promising way to narrow down the region of interest for sequencing priority. © 2013 Lin et al.

Five years of GWAS discovery

The past five years have seen many scientific and biological discoveries made through the experimental design of genome-wide association studies (GWASs). These studies were aimed at detecting variants at genomic loci that are associated with complex traits in the population and, in particular, at detecting associations between common single-nucleotide polymorphisms (SNPs) and common diseases such as heart disease, diabetes, auto-immune diseases, and psychiatric disorders. We start by giving a number of quotes from scientists and journalists about perceived problems with GWASs. We will then briefly give the history of GWASs and focus on the discoveries made through this experimental design, what those discoveries tell us and do not tell us about the genetics and biology of complex traits, and what immediate utility has come out of these studies. Rather than giving an exhaustive review of all reported findings for all diseases and other complex traits, we focus on the results for auto-immune diseases and metabolic diseases. We return to the perceived failure or disappointment about GWASs in the concluding section. © 2012 The American Society of Human Genetics.

Recurrence risk modelling of the genetic susceptibility to ankylosing spondylitis

Objectives - It has long been suspected that susceptibility to ankylosing spondylitis (AS) is influenced by genes lying distant to the major histocompatibility complex. This study compares genetic models of AS to assess the most likely mode of inheritance, using recurrence risk ratios in relatives of affected subjects. Methods - Recurrence risk ratios in different degrees of relatives were determined using published data from studies specifically designed to address the question. The methods of Risch were used to determine the expected recurrence risk ratios in different degrees of relatives, assuming equal first degree relative recurrence risk between models. Goodness of fit was determined by χ2 comparison of the expected number of affected subjects with the observed number, given equal numbers of each type of relative studied. Results - The recurrence risks in different degrees of relatives were: monozygotic (MZ) twins 63% (17/27), first degree relatives 8.2% (441/5390), second degree relatives 1.0% (8/834), and third degree relatives 0.7% (7/997). Parent-child recurrence risk (7.9%, 37/466) was not significantly different from the sibling recurrence risk (8.2%, 404/4924), excluding a significant dominance genetic component to susceptibility. Poor fitting models included single gene, genetic heterogeneity, additive, two locus multiplicative, and one locus and residual polygenes (χ2 > 32 (two degrees of freedom), p < 10-6 for all models). The best fitting model studied was a five locus model with multiplicative interaction between loci (χ2 = 1.4 (two degrees of freedom), p = 0.5). Oligogenic multiplicative models were the best fitting over a range of population prevalences and first degree recurrence risk rates. Conclusions - This study suggests that of the genetic models tested, the most likely model operating in AS is an oligogenic model with predominantly multiplicative interaction between loci.

Association of variants in MMEL1 and CTLA4 with rheumatoid arthritis in the Han Chinese population

Background: The genome-wide association study era has made great progress in identifying susceptibility genes and genetic loci for rheumatoid arthritis (RA) in populations of White European ancestry. However, few studies have tried to dissect disease aetiopathogenesis in other ethnic populations. Objective: To investigate these associations in the Han Chinese population. Methods: Haplotypes from the HapMap database Chinese population were used to select tag-single-nucleotide polymorphisms (SNPs) (r2 =0.8) across 19 distinct RA genomic regions. A two phase case-control association study was performed, with 169 SNPs genotyped in phase I (n=571 cases, n=880 controls), and 64 SNPs achieving p<0.2 in the first phase being genotyped in phase II (n=464 cases, n=822 controls). Association statistics were calculated using permutation tests both unadjusted and adjusted for the number of markers studied. Results: Robust association was detected for MMEL1 and CTLA4 , and modest association was identified for another six loci: PADI4 , STAT4 , PRDM1 , CDK6 , TRAF1-C5 and KIF5A-PIP4K2C. All three markers genotyped in MMEL1 demonstrated association, with peak signal for rs3890745 (p=2.6×10 -5unadjusted, p=0.003 adjusted, OR=0.79). For CTLA4 , significance was detected for three of five variants showing association, with peak association for marker rs12992492 (p=4.3×10-5 unadjusted, p=0.0021 adjusted, OR=0.77). Lack of association of common variants in PTPN22 with RA in Han Chinese was confirmed. Conclusion: This study identifies MMEL1 and CTLA4 as RA susceptibility genes, provides suggestive evidence of association for a further six loci in the Han Chinese population and confirms lack of PTPN22 association in Asian populations. It also confirms the value of multiethnic population studies to help dissect disease aetiopathogenesis.

Dissection of class III major histocompatibility complex haplotypes associated with rheumatoid arthritis

Objective. We have previously identified a single-nucleotide polymorphism (SNP) haplotype involving the lymphotoxin α (LTA) and tumor necrosis factor (TNF) loci (termed haplotype LTA-TNF2) on chromosome 6 that shows differential association with rheumatoid arthritis (RA) on HLA-DRB1*0404 and *0401 haplotypes, suggesting the presence of additional non-HLA-DRB1 RA susceptibility genes on these haplotypes. To refine this association, we performed a case-control association study using both SNPs and microsatellite markers in haplotypes matched either for HLA-DRB1*0404 or for HLA-DRB1*0401. Methods. Fourteen SNPs lying between HLA-DRB1 and LTA were genotyped in 87 DRB1*04-positive families. High-density microsatellite typing was performed using 24 markers spanning 2,500 kb centered around the TNF gene in 305 DRB1*0401 or *0404 cases and 400 DRB1*0401 or *0404 controls. Single-marker, 2-marker, and 3-marker minihaplotypes were constructed and their frequencies compared between the DRB1*0401 and DRB1*0404 matched case and control haplotypes. Results. Marked preservation of major histocompatibility complex haplotypes was seen, with chromosomes carrying LTA-TNF2 and either DRB1*0401 or DRB1*0404 both carrying an identical SNP haplotype across the 1-Mb region between TNF and HLA-DRB1. Using microsatellite markers, we observed two 3-marker minihaplotypes that were significantly overrepresented in the DRB1*0404 case haplotypes (P = 0.00024 and P = 0.00097). Conclusion. The presence of a single extended SNP haplotype between LTA-TNF2 and both DRB1*0401 and DRB1*0404 is evidence against this region harboring the genetic effects in linkage disequillbrium with LTA-TNF2. Two RA-associated haplotypes on the background of DRB1*0404 were identified in a 126-kb region surrounding and centromeric to the TNF locus.

Promise and pitfalls of the Immunochip

Genomewide association studies (GWAS) have proven a powerful hypothesis-free method to identify common disease-associated variants. Even quite large GWAS, however, have only at best identified moderate proportions of the genetic variants contributing to disease heritability. To provide cost-effective genotyping of common and rare variants to map the remaining heritability and to fine-map established loci, the Immunochip Consortium has developed a 200,000 SNP chip that has been produced in very large numbers for a fraction of the cost of GWAS chips. This chip provides a powerful tool for immunogenetics gene mapping.

Genome-wide association study for radiographic vertebral fractures: A potential role for the 16q24 BMD locus

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged >55years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey-Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey-Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han-Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p<5×10-8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p=4.6×10-8. However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% CI: 0.98-1.14; p=0.17), displaying high degree of heterogeneity (I2=57%; Qhet p=0.0006). Under Han-Eskin alternative random effects model the summary effect was significant (p=0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size >1.25) may still be consistent with an effect size <1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.

The role of germline polymorphisms in the T-cell receptor in susceptibility to ankylosing spondylitis

The role of germline polymorphisms of the T-cell receptor A/D and B loci in susceptibility to ankylosing spondylitis was investigated by linkage studies using microsatellite markers in 215 affected sibling pairs. The presence of a significant susceptibility gene (lambda ≤ 1.6) at the TCRA/D locus was excluded (LOD score < -2.0). At the TCRB locus, there was weak evidence of the presence of a susceptibility gene (P = 0.01, LOD score 1.1). Further family studies will be required to determine whether this is a true or false-positive finding. It is unlikely that either the TCRA/D or TCRB loci contain genes responsible for more than a moderate proportion of the non-MHC genetic susceptibility to ankylosing spondylitis.

Suggestive linkage of 2p22-25 and 11q12-13 with low bone mineral density at the lumbar spine in the Irish population

Osteoporosis is a disease characterized by low bone mineral density (BMD) and poor bone quality. Peak bone density is achieved by the third decade of life, after which bone is maintained by a balanced cycle of bone resorption and synthesis. Age-related bone loss occurs as the bone resorption phase outweighs the bone synthesis phase of bone metabolism. Heritability accounts for up to 90% of the variability in BMD. Chromosomal loci including 1p36, 2p22-25, 11q12-13, parathyroid hormone receptor type 1 (PTHR1), interleukin-6 (IL-6), interleukin 1 alpha (IL-1α) and type II collagen A1/vitamin D receptor (COL11A1/VDR) have been linked or shown suggestive linkage with BMD in other populations. To determine whether these loci predispose to low BMD in the Irish population, we investigated 24 microsatellite markers at 7 chromosomal loci by linkage studies in 175 Irish families of probands with primary low BMD (T-score ≤ -1.5). Nonparametric analysis was performed using the maximum likelihood variance estimation and traditional Haseman-Elston tests on the Mapmaker/Sibs program. Suggestive evidence of linkage was observed with lumbar spine BMD at 2p22-25 (maximum LOD score 2.76) and 11q12-13 (MLS 2.55). One region, 1p36, approached suggestive linkage with femoral neck BMD (MLS 2.17). In addition, seven markers achieved LOD scores > 1.0, D2S149, D11S1313, D11S987, D11S1314 including those encompassing the PTHR1 (D3S3559, D3S1289) for lumbar spine BMD and D2S149 for femoral neck BMD. Our data suggest that genes within a these chromosomal regions are contributing to a predisposition to low BMD in the Irish population.

Polymorphisms in a desaturase 2 ortholog associate with cuticular hydrocarbon and male mating success variation in a natural population of Drosophila serrata

Elucidating the nature of genetic variation underlying both sexually selected traits and the fitness components of sexual selection is essential to understanding the broader consequences of sexual selection as an evolutionary process. To date, there have been relatively few attempts to connect the genetic variance in sexually selected traits with segregating DNA sequence polymorphisms. We set out to address this in a well-characterized sexual selection system - the cuticular hydrocarbons (CHCs) of Drosophila serrata - using an indirect association study design that allowed simultaneous estimation of the genetic variance in CHCs, sexual fitness and single nucleotide polymorphism (SNP) effects in an outbred population. We cloned and sequenced an ortholog of the D. melanogaster desaturase 2 gene, previously shown to affect CHC biosynthesis in D. melanogaster, and associated 36 SNPs with minor allele frequencies > 0.02 with variance in CHCs and sexual fitness. Three SNPs had significant multivariate associations with CHC phenotype (q-value < 0.05). At these loci, minor alleles had multivariate effects on CHCs that were weakly associated with the multivariate direction of sexual selection operating on these traits. Two of these SNPs had pleiotropic associations with male mating success, suggesting these variants may underlie responses to sexual selection due to this locus. There were 15 significant male mating success associations (q-value < 0.1), and interestingly, we detected a nonrandom pattern in the relationship between allele frequency and direction of effect on male mating success. The minor-frequency allele usually reduced male mating success, suggesting a positive association between male mating success and total fitness at this locus.