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.

Short-rib polydactyly and jeune syndromes are caused by mutations in WDR60

Short-rib polydactyly syndromes (SRPS I-V) are a group of lethal congenital disorders characterized by shortening of the ribs and long bones, polydactyly, and a range of extraskeletal phenotypes. A number of other disorders in this grouping, including Jeune and Ellis-van Creveld syndromes, have an overlapping but generally milder phenotype. Collectively, these short-rib dysplasias (with or without polydactyly) share a common underlying defect in primary cilium function and form a subset of the ciliopathy disease spectrum. By using whole-exome capture and massive parallel sequencing of DNA from an affected Australian individual with SRPS type III, we detected two novel heterozygous mutations in WDR60, a relatively uncharacterized gene. These mutations segregated appropriately in the unaffected parents and another affected family member, confirming compound heterozygosity, and both were predicted to have a damaging effect on the protein. Analysis of an additional 54 skeletal ciliopathy exomes identified compound heterozygous mutations in WDR60 in a Spanish individual with Jeune syndrome of relatively mild presentation. Of note, these two families share one novel WDR60 missense mutation, although haplotype analysis suggested no shared ancestry. We further show that WDR60 localizes at the base of the primary cilium in wild-type human chondrocytes, and analysis of fibroblasts from affected individuals revealed a defect in ciliogenesis and aberrant accumulation of the GLI2 transcription factor at the centrosome or basal body in the absence of an obvious axoneme. These findings show that WDR60 mutations can cause skeletal ciliopathies and suggest a role for WDR60 in ciliogenesis.

Matters Arising: Zhu et al, "A novel gene variation of TNFα associated with ankylosing spondylitis: A reconfirmed study"

In a recently published study1 involving 79 ankylosing spondylitis (AS) patients and 132 unrelated healthy blood donors, Zhu et al report association of a tumour necrosis factor α (TNFα) single nucleotide polymorphism (SNP), rs1799724, with AS.1 I have many concerns with this paper and its conclusions...

Finnish HLA studies confirm the increased risk conferred by HLA-B27 homozygosity in ankylosing spondylitis

Objective: To determine the influence of HLA-B27 homozygosity and HLA-DRB1 alleles in the susceptibility to, and severity of, ankylosing spondylitis in a Finnish population. Methods: 673 individuals from 261 families with ankylosing spondylitis were genotyped for HLA-DRB1 alleles and HLA-B27 heterozygosity/ homozygosity. The frequencies of HLA-B27 homozygotes in probands from these families were compared with the expected number of HLA-B27 homozygotes in controls under Hardy-Weinberg equilibrium (HWE). The effect of HLA-DRB1 alleles was assessed using a logistic regression procedure conditioned on HLA-B27 and case-control analysis. Results: HLA-B27 was detected in 93% of cases of ankylosing spondylitis. An overrepresentation of HLA-B27 homozygotes was noted in ankylosing spondylitis (11%) compared with the expected number of HLA-B27 homozygotes under HWE (4%) (odds ratio (OR) = 3.3 (95% confidence interval, 1.6 to 6.8), p = 0.002). HLA-B27 homozygosity was marginally associated with reduced BASDAI (HLA-B27 homozygotes, 4.5 (1.6); HLA-B27 heterozygotes, 5.4 (1.8) (mean (SD)), p = 0.05). Acute anterior uveitis (AAU) was present in significantly more HLA-B27 positive cases (50%) than HLA-B27 negative cases (16%) (OR = 5.4 (1.7 to 17), p<0.004). HLA-B27 positive cases had a lower average age of symptom onset (26.7 (8.0) years) compared with HLA-B27 negative cases (35.7 (11.2) years) (p<0.0001). Conclusions: HLA-627 homozygosity is associated with a moderately increased risk of ankylosing spondylitis compared with HLA-β27 heterozygosity. HLA-B27 positive cases had an earlier age of onset of ankylosing spondylitis than HLA-B27 negative cases and were more likely to develop AAU. HLA-DRB1 alleles may influence the age of symptom onset of ankylosing spondylitis.

Genetic testing for haemochromatosis in patients with chondrocalcinosis

Hereditary haemochromatosis (HH) is the most common lethal monogenic human disease, affecting roughly 1 in 300 white northern Europeans. Homozygosity for the C282Y polymorphism within the HFE gene causes more than 80% of cases, with compound heterozygosity of the C282Y and H63D polymorphism also increasing susceptibility to disease. The aim of this study was to determine the frequency of the C282Y and H63D polymorphisms in the disease, and to assess the risk of HH in heterozygotes for the C282Y polymorphism. 128 patients were recruited because of either radiographic chondrocalcinosis (at least bicompartmental knee disease or joints other than the knee involved) or CPPD pseudogout. Genotyping of the HFE C282Y and H63D mutations was performed using PCR/SSP and genotypes for the C282Y polymorphism confirmed by PCR/RFLP. Historical white European control data were used for comparison. Two previously undiagnosed C282Y homozygotes (1.6%), and 16 C282Y heterozygotes (12.5%), including four (3.1%) C282Y/ H63D compound heterozygotes were identified. This represents a significant overrepresentation of C282Y homozygotes (relative risk 3.4, p-0.037), but the number of heterozygotes was not significantly increased. At a cost per test of £1 for each subject, screening all patients with chondrocalcinosis using the above ascertainment criteria costs only £64 for each case of haemochromatosis identified, clearly a highly cost effective test given the early mortality associated with untreated haemochromatosis. Routine screening for haemochromatosis in patients with appreciable chondrocatcinosis is recommended.

A novel ACVR1 mutation in the glycine/serine-rich domain found in the most benign case of a fibrodysplasia ossificans progressiva variant reported to date

Fibrodysplasia Ossificans Progressiva (FOP) is a rare, autosomal dominant condition, classically characterised by heterotopic ossification beginning in childhood and congenital great toe malformations; occurring in response to a c.617 G>A ACVR1 mutation in the functionally important glycine/serine-rich domain of exon 6. Here we describe a novel c.587 T>C mutation in the glycine/serine-rich domain of ACVR1, associated with delayed onset of heterotopic ossification and an exceptionally mild clinical course. Absence of great toe malformations, the presence of early ossification of the cervical spine facets joints, plus mild bilateral camptodactyly of the 5th fingers, together with a novel ACVR1 mutation, are consistent with the 'FOP-variant' syndrome. The c.587 T>C mutation replaces a conserved leucine with proline at residue 196. Modelling of the mutant protein reveals a steric clash with the kinase domain that will weaken interactions with FKBP12 and induce exposure of the glycine/serine-rich repeat. The mutant receptor is predicted to be hypersensitive to ligand stimulation rather than being constitutively active, consistent with the mild clinical phenotype. This case extends our understanding of the 'FOP-variant' syndrome.

The effect of LRP5 polymorphisms on bone mineral density is apparent in childhood

Bone mass acquired during childhood is the primary determinant of adult bone mineral density (BMD) and osteoporosis risk. Bone accrual is subject to genetic influences. Activating and inactivating LRP5 gene mutations elicit extreme bone phenotypes, while more common LRP5 polymorphisms are associated with normal variation of BMD. Our aim was to test the hypothesis that LRP5 gene polymorphisms influence bone mass acquisition during childhood. The association between LRP5 gene polymorphisms and bone size and mineralization was examined in 819 unrelated British Caucasian children (n = 429 boys) aged 9 years. Height, weight, pubertal status (where available), total-body and spinal bone area, bone mineral content (BMC), BMD, and area-adjusted BMC (aBMC) were assessed. Dual-energy X-ray absorptiometry (DXA)-gene associations were assessed by linear regression, with adjustment for age, gender, pubertal status, and body size parameters. There were 140, 79, 12, and 2 girls who achieved Tanner stages I-IV, respectively, and 179 and 32 boys who achieved Tanner stages I and II, respectively. The rs2306862 (N740N) coding polymorphism in exon 10 of the LRP5 gene was associated with spinal BMD and aBMC (each P = 0.01) and total-body BMD and aBMC (P = 0.04 and 0.03, respectively). Adjusting for pubertal stage strengthened associations between this polymorphism and spinal BMD and aBMC (P = 0.01 and 0.002, respectively). Individuals homozygous for the T allele had greater spinal BMD and aBMC scores than those homozygous for the C allele. A dose effect was apparent as the mean spinal BMD and aBMC of heterozygous TC individuals were intermediate between those of their TT and CC counterparts. The N740N polymorphism in exon 10 of LRP5 was associated with spinal BMD and aBMC in pre- and early pubertal children. These results indicate that LRP5 influences volumetric bone density in childhood, possibly through effects on trabecular bone formation.

Novel ANKH amino terminus mutation (Pro5Ser) associated with early-onset calcium pyrophosphate disease with associated phosphaturia

This report describes a 32-year-old woman presenting since childhood with progressive calcium pyrophosphate disease (CPPD), characterized by severe arthropathy and chondrocalcinosis involving multiple peripheral joints and intervertebral disks. Because ANKH mutations have been previously described in familial CPPD, the proband's DNA was assessed at this locus by direct sequencing of promoter and coding regions and revealed 3 sequence variants in ANKH. Sequences of exon 1 revealed a novel isolated nonsynonymous mutation (c.13 C>T), altering amino acid in codon 5 from proline to serine (CCG>TCG). Sequencing of parental DNA revealed an identical mutation in the proband's father but not the mother. Subsequent clinical evaluation demonstrated extensive chondrocalcinosis and degenerative arthropathy in the proband's father. In summary, we report a novel mutation, not previously described, in ANKH exon 1, wherein serine replaces proline, in a case of early-onset severe CPPD associated with metabolic abnormalities, with similar findings in the proband's father.

Genetic, functional, and histopathological evaluation of two C-terminal BRCA1 missense variants

Background: The vast majority of BRCA1 missense sequence variants remain uncharacterised for their possible effect on protein expression and function, and therefore are unclassified in terms of their pathogenicity. BRCA1 plays diverse cellular roles and it is unlikely that any single functional assay will accurately reflect the total cellular implications of missense mutations in this gene. Objective: To elucidate the effect of two BRCA1 variants, 5236G>C (G1706A) and 5242C>A (A1708E) on BRCA1 function, and to survey the relative usefulness of several assays to direct the characterisation of other unclassified variants in BRCA genes. Methods and Results: Data from a range of bioinformatic, genetic, and histopathological analyses, and in vitro functional assays indicated that the 1708E variant was associated with the disruption of different cellular functions of BRCA1. In transient transfection experiments in T47D and 293T cells, the 1708E product was mislocalised to the cytoplasm and induced centrosome amplification in 293T cells. The 1708E variant also failed to transactivate transcription of reporter constructs in mammalian transcriptional transactivation assays. In contrast, the 1706A variant displayed a phenotype comparable to wildtype BRCA1 in these assays. Consistent with functional data, tumours from 1708E carriers showed typical BRCA1 pathology, while tumour material from 1706A carriers displayed few histopathological features associated with BRCA1 related tumours. Conclusions: A comprehensive range of genetic, bioinformatic, and functional analyses have been combined for the characterisation of BRCA1 unclassified sequence variants. Consistent with the functional analyses, the combined odds of causality calculated for the 1706A variant after multifactorial likelihood analysis (1:142) indicates a definitive classification of this variant as "benign". In contrast, functional assays of the 1708E variant indicate that it is pathogenic, possibly through subcellular mislocalisation. However, the combined odds of 262:1 in favour of causality of this variant does not meet the minimal ratio of 1000:1 for classification as pathogenic, and A1708E remains formally designated as unclassified. Our findings highlight the importance of comprehensive genetic information, together with detailed functional analysis for the definitive categorisation of unclassified sequence variants. This combination of analyses may have direct application to the characterisation of other unclassified variants in BRCA1 and BRCA2.

Polymorphisms of the CYP2D6 gene increase susceptibility to ankylosing spondylitis

Ankylosing spondylitis (AS) is a common and highly familial rheumatic disorder. The sibling recurrence risk ratio for the disease is 63 and heritability assessed in twins > 90%. Although MHC genes, including HLA-B27, contribute only 20-50% of the genetic risk for the disease, no non-MHC gene has yet been convincingly demonstrated to influence either susceptibility to the disease or its phenotypic expression. Previous linkage and association studies have suggested the presence of a susceptibility gene for AS close to, or within, the cytochrome P450 2D6 gene (CYP2D6, debrisoquine hydroxylase) located at chromosome 22q13.1. We performed a linkage study of chromosome 22 in 200 families with AS affected sibling-pairs. Association of alleles of the CYP2D6 gene was examined by both case-control and within-family means. For case-control studies, 617 unrelated individuals with AS (361 probands from sibling-pair and parent-case trio families and 256 unrelated non-familial sporadic cases) and 402 healthy ethnically matched controls were employed. For within-family association studies, 361 families including 161 parent-case trios and 200 affected sibling-pair families were employed. Homozygosity for poor metabolizer alleles was found to be associated with AS. Heterozygosity for the most frequent poor metabolizer allele (CYP2D6*4) was not associated with increased susceptibility to AS. Significant within-family association of CYP2D6*4 alleles and AS was demonstrated. Weak linkage was also demonstrated between CYP2D6 and AS. We postulate that altered metabolism of a natural toxin or antigen by the CYP2D6 gene may increase susceptibility to AS.

Characterization of 26 new microsatellite loci in the dugong (Dugong dugon)

Twenty-six microsatellite loci have been isolated from a dugong (Dugong dugon). The average heterozygosity was 0.52 with two to 10 alleles per locus surveyed from 50 individuals. The markers are suitable for genetic mark-recapture (PID = 5 × 10-16) in dugongs and they could also be used to quantify physical tag loss, estimate relatedness, assign paternity, elucidate population structure and identify migrants.

New microsatellite loci for Carcharhinid sharks (Carcharhinus tilstoni and C. sorrah) and their cross-amplification in other shark species

Twelve microsatellite DNA markers were isolated in the spot-tail shark (Carcharhinus sorrah) and nine were isolated in Australian black-tip shark (Carcharhinus tilstoni). These loci plus 18 others developed for sharks from the genera Negaprion, Ginglymostoma, Carcharodon and Isurus were tested for amplification success on four species of Carcharhinus (including C. sorrah and C. tilstoni) and four other species representing three diverse families. Cross-amplification was most common within families. Five loci were subsequently tested for polymorphism on 50 C. sorrah and 60 C. tilstoni. The number of alleles per locus was two to 24 and the average heterozygosity was 0.54 (range 0.16-0.87) for C. sorrah and 0.64 (range 0.44-0.78) for C. tilstoni. These loci may be useful tools for genetic analyses of the Carcharhinidae.

Rapid Global Expansion of the Fungal Disease Chytridiomycosis into Declining and Healthy Amphibian Populations.

The fungal disease chytridiomycosis, caused by Batrachochytrium dendrobatidis, is enigmatic because it occurs globally in both declining and apparently healthy (non-declining) amphibian populations. This distribution has fueled debate concerning whether, in sites where it has recently been found, the pathogen was introduced or is endemic. In this study, we addressed the molecular population genetics of a global collection of fungal strains from both declining and healthy amphibian populations using DNA sequence variation from 17 nuclear loci and a large fragment from the mitochondrial genome. We found a low rate of DNA polymorphism, with only two sequence alleles detected at each locus, but a high diversity of diploid genotypes. Half of the loci displayed an excess of heterozygous genotypes, consistent with a primarily clonal mode of reproduction. Despite the absence of obvious sex, genotypic diversity was high (44 unique genotypes out of 59 strains). We provide evidence that the observed genotypic variation can be generated by loss of heterozygosity through mitotic recombination. One strain isolated from a bullfrog possessed as much allelic diversity as the entire global sample, suggesting the current epidemic can be traced back to the outbreak of a single clonal lineage. These data are consistent with the current chytridiomycosis epidemic resulting from a novel pathogen undergoing a rapid and recent range expansion. The widespread occurrence of the same lineage in both healthy and declining populations suggests that the outcome of the disease is contingent on environmental factors and host resistance.