Investigating the genetic association between ERAP1 and ankylosing spondylitis

A strong association between ERAP1 and ankylosing spondylitis (AS) was recently identified by the Wellcome Trust Case Control Consortium and the Australo-Anglo-American Spondylitis Consortium (WTCCC-TASC) study. ERAP1 is highly polymorphic with strong linkage disequilibrium evident across the gene. We therefore conducted a series of experiments to try to identify the primary genetic association(s) with ERAP1. We replicated the original associations in an independent set of 730 patients and 1021 controls, resequenced ERAP1 to define the full extent of coding polymorphisms and tested all variants in additional association studies. The genetic association with ERAP1 was independently confirmed; the strongest association was with rs30187 in the replication set (P = 3.4 × 103). When the data were combined with the original WTCCC-TASC study the strongest association was with rs27044 (P = 1.1 × 10-9). We identified 33 sequence polymorphisms in ERAP1, including three novel and eight known non-synonymous polymorphisms. We report several new associations between AS and polymorphisms distributed across ERAP1 from the extended case-control study, the most significant of which was with rs27434 (P = 4.7 × 10-7). Regression analysis failed to identify a primary association clearly; we therefore used data from HapMap to impute genotypes for an additional 205 non-coding SNPs located within and adjacent to ERAP1. A number of highly significant associations (P < 5 × 10-9) were identified in regulatory sequences which are good candidates for causing susceptibility to AS, possibly by regulating ERAP1 expression. © 2009 The Author(s).

Genetic analyses in a sample of individuals with high or low BMD shows association with multiple Wnt pathway genes

INTRODUCTION Although the high heritability of BMD variation has long been established, few genes have been conclusively shown to affect the variation of BMD in the general population. Extreme truncate selection has been proposed as a more powerful alternative to unselected cohort designs in quantitative trait association studies. We sought to test these theoretical predictions in studies of the bone densitometry measures BMD, BMC, and femoral neck area, by investigating their association with members of the Wnt pathway, some of which have previously been shown to be associated with BMD in much larger cohorts, in a moderate-sized extreme truncate selected cohort (absolute value BMD Z-scores = 1.5-4.0; n = 344). MATERIALS AND METHODS Ninety-six tag-single nucleotide polymorphism (SNPs) lying in 13 Wnt signaling pathway genes were selected to tag common genetic variation (minor allele frequency [MAF] > 5% with an r(2) > 0.8) within 5 kb of all exons of 13 Wnt signaling pathway genes. The genes studied included LRP1, LRP5, LRP6, Wnt3a, Wnt7b, Wnt10b, SFRP1, SFRP2, DKK1, DKK2, FZD7, WISP3, and SOST. Three hundred forty-four cases with either high or low BMD were genotyped by Illumina Goldengate microarray SNP genotyping methods. Association was tested either by Cochrane-Armitage test for dichotomous variables or by linear regression for quantitative traits. RESULTS Strong association was shown with LRP5, polymorphisms of which have previously been shown to influence total hip BMD (minimum p = 0.0006). In addition, polymorphisms of the Wnt antagonist, SFRP1, were significantly associated with BMD and BMC (minimum p = 0.00042). Previously reported associations of LRP1, LRP6, and SOST with BMD were confirmed. Two other Wnt pathway genes, Wnt3a and DKK2, also showed nominal association with BMD. CONCLUSIONS This study shows that polymorphisms of multiple members of the Wnt pathway are associated with BMD variation. Furthermore, this study shows in a practical trial that study designs involving extreme truncate selection and moderate sample sizes can robustly identify genes of relevant effect sizes involved in BMD variation in the general population. This has implications for the design of future genome-wide studies of quantitative bone phenotypes relevant to osteoporosis.

Breakthroughs in genetic studies of ankylosing spondylitis

Ankylosing spondylitis (AS), the prototypic seronegative arthropathy, is known to be highly heritable, with >90% of the risk of developing the disease determined genetically. As with most common heritable diseases, progress in identifying the genes involved using family-based or candidate gene approaches has been slow. The recent development of the genome-wide association study approach has revolutionized genetic studies of such diseases. Early studies in ankylosing spondylitis have produced two major breakthroughs in the identification of genes contributing roughly one third of the population attributable risk of the disease, and pointing directly to a potential therapy. These exciting findings highlight the potential of future more comprehensive genetic studies of determinants of disease risk and clinical manifestations, and are the biggest advance in our understanding of the causation of the disease since the discovery of the association with HLA-B27.

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.

Genetic and genomic studies of PADI4 in rheumatoid arthritis

Objectives. Strong genetic association of rheumatoid arthritis (RA) with PADI4 (peptidyl arginine deiminase) has previously been described in Japanese, although this was not confirmed in a subsequent study in the UK. We therefore undertook a further study of genetic association between PADI4 and RA in UK Caucasians and also studied expression of PADI4 in the peripheral blood of patients with RA. Methods. Seven single-nucleotide polymorphisms (SNP) were genotyped using polymerase chain reaction (PCR)-restriction fragment length polymorphism in 111 RA cases and controls. A marker significantly associated with RA (PADI4_100, rs#2240339) in this first data set (P = 0.03) was then tested for association in a larger group of 439 RA patients and 428 controls. PADI4 transcription was also assessed by real-time quantitative PCR using RNA extracted from peripheral blood mononuclear cells from 13 RA patients and 11 healthy controls. Results. A single SNP was weakly associated with RA (P = 0.03) in the initial case-control study, a single SNP (PADI4_100) and a two marker haplotype of that SNP and the neighbouring SNP (PADI4_04) were significantly associated with RA (P = 0.02 and P = 0.03 respectively). PADI4_100 was not associated with RA in a second sample set. PADI4 expression was four times greater in cases than controls (P = 0.004), but expression levels did not correlate with the levels of markers of inflammation. Conclusion. PADI4 is significantly overexpressed in the blood of RA patients but genetic variation within PADI4 is not a major risk factor for RA in Caucasians.

Genome-wide association study identifies novel genetic variants contributing to variation in blood metabolite levels

Metabolites are small molecules involved in cellular metabolism, which can be detected in biological samples using metabolomic techniques. Here we present the results of genome-wide association and meta-analyses for variation in the blood serum levels of 129 metabolites as measured by the Biocrates metabolomic platform. In a discovery sample of 7,478 individuals of European descent, we find 4,068 genome- and metabolome-wide significant (Z-test, P<1.09 × 10−9) associations between single-nucleotide polymorphisms (SNPs) and metabolites, involving 59 independent SNPs and 85 metabolites. Five of the fifty-nine independent SNPs are new for serum metabolite levels, and were followed-up for replication in an independent sample (N=1,182). The novel SNPs are located in or near genes encoding metabolite transporter proteins or enzymes (SLC22A16, ARG1, AGPS and ACSL1) that have demonstrated biomedical or pharmaceutical importance. The further characterization of genetic influences on metabolic phenotypes is important for progress in biological and medical research.

Pharmacogenetic Variation at CYP2D6, CYP2C9, and CYP2C19: Population Genetic and Forensic Aspects

Pharmacogenetics deals with genetically determined variation in drug response. In this context, three phase I drug-metabolizing enzymes, CYP2D6, CYP2C9, and CYP2C19, have a central role, affecting the metabolism of about 20-30% of clinically used drugs. Since genes coding for these enzymes in human populations exhibit high genetic polymorphism, they are of major pharmacogenetic importance. The aims of this study were to develop new genotyping methods for CYP2D6, CYP2C9, and CYP2C19 that would cover the most important genetic variants altering the enzyme activity, and, for the first time, to describe the distribution of genetic variation at these loci on global and microgeographic scales. In addition, pharmacogenetics was applied to a postmortem forensic setting to elucidate the role of genetic variation in drug intoxications, focusing mainly on cases related to tricyclic antidepressants, which are commonly involved in fatal drug poisonings in Finland. Genetic variability data were obtained by genotyping new population samples by the methods developed based on PCR and multiplex single-nucleotide primer extension reaction, as well as by collecting data from the literature. Data consisted of 138, 129, and 146 population samples for CYP2D6, CYP2C9, and CYP2C19, respectively. In addition, over 200 postmortem forensic cases were examined with respect to drug and metabolite concentrations and genotypic variation at CYP2D6 and CYP2C19. The distribution of genetic variation within and among human populations was analyzed by descriptive statistics and variance analysis and by correlating the genetic and geographic distances using Mantel tests and spatial autocorrelation. The correlation between phenotypic and genotypic variation in drug metabolism observed in postmortem cases was also analyzed statistically. The genotyping methods developed proved to be informative, technically feasible, and cost-effective. Detailed molecular analysis of CYP2D6 genetic variation in a global survey of human populations revealed that the pattern of variation was similar to those of neutral genomic markers. Most of the CYP2D6 diversity was observed within populations, and the spatial pattern of variation was best described as clinal. On the other hand, genetic variants of CYP2D6, CYP2C9, and CYP2C19 associated with altered enzymatic activity could reach extremely high frequencies in certain geographic regions. Pharmacogenetic variation may also be significantly affected by population-specific demographic histories, as seen within the Finnish population. When pharmacogenetics was applied to a postmortem forensic setting, a correlation between amitriptyline metabolic ratios and genetic variation at CYP2D6 and CYP2C19 was observed in the sample material, even in the presence of confounding factors typical for these cases. In addition, a case of doxepin-related fatal poisoning was shown to be associated with a genetic defect at CYP2D6. Each of the genes studied showed a distinct variation pattern in human populations and high frequencies of altered activity variants, which may reflect the neutral evolution and/or selective pressures caused by dietary or environmental exposure. The results are relevant also from the clinical point of view since the genetic variation at CYP2D6, CYP2C9, and CYP2C19 already has a range of clinical applications, e.g. in cancer treatment and oral anticoagulation therapy. This study revealed that pharmacogenetics may also contribute valuable information to the medicolegal investigation of sudden, unexpected deaths.

Genetic differentiation analysis for the identification of complementary parental pools for sorghum hybrid breeding in Ethiopia

Key message The potential for exploiting heterosis for sorghum hybrid production in Ethiopia with improved local adaptation and farmers preferences has been investigated and populations suitable for initial hybrid development have been identified. Abstract Hybrids in sorghum have demonstrated increased productivity and stability of performance in the developed world. In Ethiopia, the uptake of hybrid sorghum has been limited to date, primarily due to poor adaptation and absence of farmer’s preferred traits in existing hybrids. This study aimed to identify complementary parental pools to develop locally adapted hybrids, through an analysis of whole genome variability of 184 locally adapted genotypes and introduced hybrid parents (R and B). Genetic variability was assessed using genetic distance, model-based STRUCTURE analysis and pair-wise comparison of groups. We observed a high degree of genetic similarity between the Ethiopian improved inbred genotypes and a subset of landraces adapted to lowland agro-ecology with the introduced R lines. This coupled with the genetic differentiation from existing B lines, indicated that these locally adapted genotype groups are expected to have similar patterns of heterotic expression as observed between introduced R and B line pools. Additionally, the hybrids derived from these locally adapted genotypes will have the benefit of containing farmers preferred traits. The groups most divergent from introduced B lines were the Ethiopian landraces adapted to highland and intermediate agro-ecologies and a subset of lowland-adapted genotypes, indicating the potential for increased heterotic response of their hybrids. However, these groups were also differentiated from the R lines, and hence are different from the existing complementary heterotic pools. This suggests that although these groups could provide highly divergent parental pools, further research is required to investigate the extent of heterosis and their hybrid performance.

The study of genetic diversity and population structure of the Korean fleshy shrimp, Fenneropenaeus chinensis, using newly developed microsatellite markers

The fleshy shrimp, Fenneropenaeus chinensis, is the family of Penaeidae and one of the most economically important marine culture species in Korea. However, its genetic characteristics have never been studied. In this study, a total of 240 wild F. chinensis individuals were collected from four locations as follows: Narodo (NRD, n = 60), Beopseongpo (BSP, n = 60), Chaesukpo (CSP, n = 60), and Cheonsuman (CSM, n = 60). Genetic variability and the relationships among four wild F. chinensis populations were analyzed using 13 newly developed microsatellite loci. Relatively high levels of genetic variability (mean allelic richness = 16.87; mean heterozygosity = 0.845) were found among localities. Among the 52 population loci, 13 showed significant deviation from the Hardy–Weinberg equilibrium. Neighbor-joining, principal coordinate, and molecular variance analyses revealed the presence of three subpopulations (NRD, CSM, BSP and CSP), which was consistent with clustering based on genetic distance. The mean observed heterozygosity values of the NRD, CSM, BSP, and CSP populations were 0.724, 0.821, 0.814, and 0.785 over all loci, respectively. These genetic variability and differentiation results of the four wild populations can be applied for future genetic improvement using selective breeding and to design suitable management guidelines for Korean F. chinensis culture.

Use of allozyme markers to determine the genetic structure of hatchery population of Thai pangas, Pangasius hypophthalmus of Jessore, Bangladesh

Genetic structure of hatchery population of Thai pangas (Pangasius hypophthalmus) of Jessore region, Bangladesh has been investigated from 1 January 2004 to 31 December 2004. Samples for this study were collected from five fish hatcheries viz. Asrom, Banchte Shekha, Chowdhury, Maola and Rezaul Haque. The enzymes were encoded by 15 gene loci: Adh-1*, Est-1*, G3pdh-2*, Gpi-1*, Gpi-2*, Idhp-1*, Idhp-2*, Ldh-1*, Ldh-2*, Mdh-1*, Mdh-2*, Pgm*, Sdh-1*, Sdh-2* and Sod*. Among them four (Est-1*, G3pdh-2*, Gpi-2*and Pgm*) were found to be polymorphic in different populations but only Gpi-2* was polymorphic in all the sampled populations. The mean proportion of polymorphic loci per population was the highest (26.7%) in Banchte Shekha hatchery while the mean proportion of heterozygous loci was 13.33% per individual in Banchte Shekha and Maola hatcheries. The UPGMA dendrogram of Nei's (1972) genetic distances indicated a relationship between the genetic distance and geographical difference. High genetic variability in stocks of Thai pangas was observed in the Banchte Shekha and Maola hatcheries and less variability was found in the other three hatcheries.

Genetic diversity of Chinese native pigs inferred from protein electrophoresis

We examined protein polymorphism of 20 native pig breeds in China and 3 introduced pig breeds. Thirty loci have been investigated, among which six loci were found to be polymorphic. Especially, the polymorphism of malate dehydrogenase (MDH), adenylate kinase (AK), and two new alleles of adenosine deaminase (ADA) had not been reported in domestic pigs and wild pigs. The percentage of polymorphic loci (P), the mean heterozygosity (H), and the mean number of alleles (A) are 0.200, 0.065, and 1.300, respectively. The degree of genetic variability of Chinese pigs as a whole was higher than that of goats, lower than that of cattle and horses, and similar to that of sheep. Using the gene frequencies of the 30 loci, Nei's genetic distance among the 20 native breeds in China and 3 introduced pig breeds was calculated by the formula of Nei. The program NEIGHBOR in PHYLIP 3.5c was chosen to construct an UPGMA tree and a NJ tree. Our results show that, of the total genetic variation found in the native pig breeds in China, 31% (0.31) is ascribable to genetic differences among breeds. About 69% of the total genetic variation is found within breeds. Most breeds are in linkage disequilibrium. The patterns of genetic similarities between the Chinese native pig breeds were not in agreement with the proposed pig type classification.

A note on recent advances in the genetic characterization of tilapia stocks in Lake Victoria region

Oreochromis esculenta, the original "ngege" is virtually extinct in Lake Victoria, and is limited to satellite lakes and reservoirs in the greater Lake Victoria region. Oreochromis variabilis can still be found in Lake Victoria and some satellite lakes in the Kyoga System, but in small numbers and only at a few localities (WANDERA and KAUFMAN, unpub. data). Little is known about the influence that species translocations have had on the genetic structure of these crucial fishery species, and even the source of the parent stocks for the introductions remain obscure. Genetic variability was examined within and among allopatric populations of three species in the tilapiine genus Oreochromis: O. esculentus (endemic to Lakes Victoria and Kyoga), and two exotic species introduced to Lake Victoria in the late 1950's to supplement the failing fisheries for native tilapiines, O. niloticus and O. leucostictus.

Population genetic structure of the Yangtze finless porpoise (Neophocaena phocaenoides asiaeorientalis): Implications for management and conservation

Understanding the population genetic structure is a prerequisite for conservation of a species. The degree of genetic variability characteristic of the mitochondrial DNA control region has been widely exploited in studies of population genetic structure and can be useful in identifying meaningful population subdivisions. To estimate the genetic profile of the Yangtze finless porpoise (Neophocaena phocaenoides asiaeorientalis), an endangered freshwater population endemic to China, the complete mtDNA control region was examined in 39 individuals belonging to seven different stocks inhabiting the middle and lower reaches of the Yangtze River. Very low genetic diversity was found (nucleotide diversity 0.0011 +/- 0.0002 and haplotypic diversity 0.65 +/- 0.05). The mtDNA genetic pattern of the Yangtze population appears to indicate a founder event in its evolutionary history and to support the marine origin for this population. Analyses by F-st and Phi(st) yielded statistically significant population genetic structure (F-st = 0.44, P < 0.05; phi(st) = 0.36, P < 0.05). These results may have significant implications for the management and conservation of the Yangtze finless porpoise in the future.

Population genetic structure of Sargassum thunbergii (Fucales, Phaeophyta) detected by RAPD and ISSR markers

Genetic variation of four populations of Sargassum thunbergii (Mert.) O. Kuntze and one outgroup of S. fusiforme (Harv.) Setchell from Shandong peninsula of China was studied with random amplified polymorphic DNA (RAPD) and inter-simple sequence repeat (ISSR) markers. A total of 28 RAPD primers and 19 ISSR primers were amplified, showing 174 loci and 125 loci, respectively. Calculation of genetic diversity with different indicators (P%, percentage of polymorphic loci; H, the expected heterozygosity; I, Shannon's information index) revealed low or moderate levels of genetic variations within each S. thunbergii population. High genetic differentiations were determined with pairwise Nei's unbiased genetic distance (D) and fixation index (F-ST) between the populations. The Mantel test showed that two types of matrices of D and FST were highly correlated, whether from RAPD or ISSR data, r=0.9310 (P = 0.008) and 0.9313 (P=0.009) respectively. Analysis of molecular variance (AMOVA) was used to apportion the variations between and within the S. thunbergii populations. It indicated that the variations among populations were higher than those within populations, being 57.57% versus 42.43% by RAPD and 59.52% versus 40.08% by ISSR, respectively. Furthermore, the Mantel test suggested that the genetic differentiations between the four populations were related to the geographical distances (r > 0.5), i.e., they conformed to the IBD (isolation by distance) model, as expected from UPGMA (unweighted pair group method with arithmetic averages) cluster analysis. As a whole, the high genetic structuring between the four S. thunbergii populations along distant locations was clearly indicated in the RAPD and ISSR analyses (r > 0.8) in our study.

Genetic variation within and among populations of Rhodiola alsia (Crassulaceae) native to the Tibetan Plateau as detected by ISSR markers

Genetic variation of 10 Rhodiola alsia ( Crassulaceae) populations from the Qinghai - Tibet Plateau of China was investigated using intersimple sequence repeat (ISSR) markers. R. alsia is an endemic species of the Qinghai - Tibet Plateau. Of the 100 primers screened, 13 were highly polymorphic. Using these primers, 140 discernible DNA fragments were generated with 112 (80%) being polymorphic, indicating pronounced genetic variation at the species level. Also there were high levels of polymorphism at the population level with the percentage of polymorphic bands (PPB) ranging from 63.4 to 88.6%. Analysis of molecular variance (AMOVA) showed that the genetic variation was mainly found among populations (70.3%) and variance within populations was 29.7%. The main factors responsible for the high level of differentiation among populations are probably the isolation from other populations and clonal propagation of this species. Occasional sexual reproduction might occur in order to maintain high levels of variation within populations. Environmental conditions could also influence population genetic structure as they occur in severe habitats. The strong genetic differentiation among populations in our study indicates that the conservation of genetic variability in R. alsia requires maintenance of as many populations as possible.