Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa.

PURPOSE: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. METHODS: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. RESULTS: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease. CONCLUSIONS: The MDH1 gene is not the cause of RP28-linked arRP. Our experimental strategy shows that long-range genomic PCR followed by UHTs provides an excellent system to perform a thorough screening of candidate genes for hereditary retinal degeneration.

Identity by descent and candidate gene mapping of Richieri-Costa and Pereira syndrome

The Richieri-Costa-Pereira syndrome is a rare autosomal recessive disorder characterized by short stature, Robin sequence, cleft mandible, pre/postaxial anomalies and clubfoot. of 15 families reported with this disorder 14 are from Brazil suggesting a founder effect. We studied 15 families using identity-by-descent as a hypothesis to attempt gene localization We have examined through linkage analysis 497 polymorphicmarkers and also performed direct sequencing of exons for 10 candidate genes selected on the basis of their expression in the developing mandible and limb. No evidence for allele sharing at any locus tested or mutations in candidate genes was found. Additional higher resolution mapping, new families and other candidate genes might improve future chances of gene identification. (C) 2003 Wiley-Liss, Inc.

Molecular analysis of carbohydrate N-acetylgalactosamine 4-O sulfotransferase 8 (CHST8) as a candidate gene for bovine spongiform encephalopathy susceptibility

Endogenous prion proteins (PrP) play the central role in the pathogenesis of transmissible spongiform encephalopathies. The carbohydrate N-acetylgalactosamine 4-O sulfotransferase 8 (CHST8) promotes the conversion of the cellular PrP(C) into the pathogenic PrP(d). Six sequence variants within the CHST8 gene were identified by comparative sequencing and genotyped for a sample of 623 animals comprising bovine spongiform encephalopathy (BSE)-affected and healthy control cows representing German Fleckvieh (German Simmental), German Holstein (Holstein-Friesian) and Brown Swiss. Significant differences in the allele, genotype and haplotype frequencies between BSE-affected and healthy cows indicate an association of sequence variant g.37254017G>T with the development of the disease in Brown Swiss cattle.

The bovine aristaless-like homeobox 4 (ALX4) as a candidate gene for syndactyly

The ALX4 (aristaless-like homeobox 4) gene encodes a paired-type homeodomain transcriptional activator and plays a major role in anterior-posterior pattern formation during limb development. Here, the cloning, genomic structure and expression of the bovine ortholog of the ALX4 gene are reported. The bovine ALX4 gene consists of four exons and is located on BTA15q28-->q29 in a region syntenic to HSA11p11.2. The transcribed ALX4 mRNA encodes a 397-amino-acid protein showing a paired-type homeodomain and a C-terminal stretch of amino acids known as the OAR- or aristaless domain. The predicted protein shares 92.5% identity to human and mouse ALX4 proteins and all three species share almost complete identity in the conserved domains. ALX4 expression was detected by reverse transcriptase polymerase chain reaction in bovine fetal limb bones. The ALX4 gene was evaluated as a candidate gene for bovine syndactyly which has been mapped on the telomeric region of cattle chromosome 15. Sequencing of the four exons with flanking sequences of the bovine ALX4 gene from a panel of 14 affected animals belonging to German Holstein, German Fleckvieh and crossbreds, and 27 unaffected individuals from German Holstein revealed five silent SNPs within the coding region out of eleven SNPs in total. Four SNPs were polymorphic in the affected animals, but in comparison to the genotyped unaffected individuals the genotype distribution showed no evidence for an association to the phenotype. Therefore our data indicate that the ALX4 gene can probably be excluded as candidate gene for bovine syndactyly in the examined animals.

Evaluation of the CTSL2 gene as a candidate gene for alopecia X in Pomeranians and Keeshonden

Alopecia X is a noninflammatory, progressive, bilateral symmetric alopecia in dogs. The disease is mainly found in Nordic breeds. The breed predisposition and a strong familial accumulation suggest a hereditary background. We analyzed the cathepsin L2 gene (CTSL2) as a candidate for alopecia X. The comparative sequencing of 14 affected and 18 control animals revealed ten polymorphisms; however, none of these polymorphisms affected the coding sequence. Haplotype analysis did not reveal an association of one particular CTSL2 haplotype with the disease phenotype; therefore, we conclude that the CTSL2 gene is probably not the causative gene for alopecia X.

Fine-mapping and mutation analysis of TRPM1: a candidate gene for leopard complex (LP) spotting and congenital stationary night blindness in horses.

Leopard Complex spotting occurs in several breeds of horses and is caused by an incompletely dominant allele (LP). Homozygosity for LP is also associated with congenital stationary night blindness (CSNB) in Appaloosa horses. Previously, LP was mapped to a 6 cm region on ECA1 containing the candidate gene TRPM1 (Transient Receptor Potential Cation Channel, Subfamily M, Member 1) and decreased expression of this gene, measured by qRT-PCR, was identified as the likely cause of both spotting and ocular phenotypes. This study describes investigations for a mutation causing or associated with the Leopard Complex and CSNB phenotype in horses. Re-sequencing of the gene and associated splice sites within the 105 624 bp genomic region of TRPM1 led to the discovery of 18 SNPs. Most of the SNPs did not have a predictive value for the presence of LP. However, one SNP (ECA1:108,249,293 C>T) found within intron 11 had a strong (P < 0.0005), but not complete, association with LP and CSNB and thus is a good marker but unlikely to be causative. To further localize the association, 70 SNPs spanning over two Mb including the TRPM1 gene were genotyped in 192 horses from three different breeds segregating for LP. A single 173 kb haplotype associated with LP and CSNB (ECA1: 108,197,355- 108,370,150) was identified. Illumina sequencing of 300 kb surrounding this haplotype revealed 57 SNP variants. Based on their localization within expressed sequences or regions of high sequence conservation across mammals, six of these SNPs were considered to be the most likely candidate mutations. While the precise function of TRPM1 remains to be elucidated, this work solidifies its functional role in both pigmentation and night vision. Further, this work has identified several potential regulatory elements of the TRPM1 gene that should be investigated further in this and other species.

Candidate-Gene Approach in Posttraumatic Stress Disorder After Urban Violence: Association Analysis of the Genes Encoding Serotonin Transporter, Dopamine Transporter, and BDNF

Posttraumatic stress disorder (PTSD) is a prevalent, disabling anxiety disorder marked by behavioral and physiologic alterations which commonly follows a chronic course. Exposure to a traumatic event constitutes a necessary, but not sufficient, factor. There is evidence from twin studies supporting a significant genetic predisposition to PTSD. However, the precise genetic loci still remain unclear. The objective of the present study was to identify, in a case-control study, whether the brain-derived neurotrophic factor (BDNF) val66met polymorphism (rs6265), the dopamine transporter (DAT1) three prime untranslated region (3`UTR) variable number of tandem repeats (VNTR), and the serotonin transporter (5-HTTPRL) short/long variants are associated with the development of PTSD in a group of victims of urban violence. All polymorphisms were genotyped in 65 PTSD patients as well as in 34 victims of violence without PTSD and in a community control group (n = 335). We did not find a statistical significant difference between the BDNF val66met and 5-HTTPRL polymorphism and the traumatic phenotype. However, a statistical association was found between DAT1 3`UTR VNTR nine repeats and PTSD (OR = 1.82; 95% CI, 1.20-2.76). This preliminary result confirms previous reports supporting a susceptibility role for allele 9 and PTSD.

Molecular variation at a candidate gene implicated in the regulation of fire ant social behavior

The fire ant Solenopsis invicta and its close relatives display an important social polymorphism involving differences in colony queen number. Colonies are headed by either a single reproductive queen (monogyne form) or multiple queens (polygyne form). This variation in social organization is associated with variation at the gene Gp-9, with monogyne colonies harboring only B-like allelic variants and polygyne colonies always containing b-like variants as well. We describe naturally occurring variation at Gp-9 in fire ants based on 185 full-length sequences, 136 of which were obtained from S. invicta collected over much of its native range. While there is little overall differentiation between most of the numerous alleles observed, a surprising amount is found in the coding regions of the gene, with such substitutions usually causing amino acid replacements. This elevated coding-region variation may result from a lack of negative selection acting to constrain amino acid replacements over much of the protein, different mutation rates or biases in coding and non-coding sequences, negative selection acting with greater strength on non-coding than coding regions, and/or positive selection acting on the protein. Formal selection analyses provide evidence that the latter force played an important role in the basal b-like lineages coincident with the emergence of polygyny. While our data set reveals considerable paraphyly and polyphyly of S. invicta sequences with respect to those of other fire ant species, the b-like alleles of the socially polymorphic species are monophyletic. An expanded analysis of colonies containing alleles of this clade confirmed the invariant link between their presence and expression of polygyny. Finally, our discovery of several unique alleles bearing various combinations of b-like and B-like codons allows us to conclude that no single b-like residue is completely predictive of polygyne behavior and, thus, potentially causally involved in its expression. Rather, all three typical b-like residues appear to be necessary.

Post-genomic update on a classical candidate gene for coronary artery disease: ESR1.

BACKGROUND: After age, sex is the most important risk factor for coronary artery disease (CAD). The mechanism through which women are protected from CAD is still largely unknown, but the observed sex difference suggests the involvement of the reproductive steroid hormone signaling system. Genetic association studies of the gene-encoding Estrogen Receptor α (ESR1) have shown conflicting results, although only a limited range of variation in the gene has been investigated. METHODS AND RESULTS: We exploited information made available by advanced new methods and resources in complex disease genetics to revisit the question of ESR1's role in risk of CAD. We performed a meta-analysis of 14 genome-wide association studies (CARDIoGRAM discovery analysis, N=≈87,000) to search for population-wide and sex-specific associations between CAD risk and common genetic variants throughout the coding, noncoding, and flanking regions of ESR1. In addition to samples from the MIGen (N=≈6000), WTCCC (N=≈7400), and Framingham (N=≈3700) studies, we extended this search to a larger number of common and uncommon variants by imputation into a panel of haplotypes constructed using data from the 1000 Genomes Project. Despite the widespread expression of ERα in vascular tissues, we found no evidence for involvement of common or low-frequency genetic variation throughout the ESR1 gene in modifying risk of CAD, either in the general population or as a function of sex. CONCLUSIONS: We suggest that future research on the genetic basis of sex-related differences in CAD risk should initially prioritize other genes in the reproductive steroid hormone biosynthesis system.

ELOVL6 Genetic Variation Is Related to Insulin Sensitivity: A New Candidate Gene in Energy Metabolism

BACKGROUND: The elongase of long chain fatty acids family 6 (ELOVL6) is an enzyme that specifically catalyzes the elongation of saturated and monounsaturated fatty acids with 12, 14 and 16 carbons. ELOVL6 is expressed in lipogenic tissues and it is regulated by sterol regulatory element binding protein 1 (SREBP-1). OBJECTIVE: We investigated whether ELOVL6 genetic variation is associated with insulin sensitivity in a population from southern Spain. DESIGN: We undertook a prospective, population-based study collecting phenotypic, metabolic, nutritional and genetic information. Measurements were made of weight and height and the body mass index (BMI) was calculated. Insulin resistance was measured by homeostasis model assessment. The type of dietary fat was assessed from samples of cooking oil taken from the participants' kitchens and analyzed by gas chromatography. Five SNPs of the ELOVL6 gene were analyzed by SNPlex. RESULTS: Carriers of the minor alleles of the SNPs rs9997926 and rs6824447 had a lower risk of having high HOMA_IR, whereas carriers of the minor allele rs17041272 had a higher risk of being insulin resistant. An interaction was detected between the rs6824447 polymorphism and the intake of oil in relation with insulin resistance, such that carriers of this minor allele who consumed sunflower oil had lower HOMA_IR than those who did not have this allele (P = 0.001). CONCLUSIONS: Genetic variations in the ELOVL6 gene were associated with insulin sensitivity in this population-based study.

Candidate gene analysis of ocular toxoplasmosis in Brazil: evidence for a role for toll-like receptor 9 (TLR9)

Toxoplasma gondii infection is an important mediator of ocular disease in Brazil more frequently than reported from elsewhere. Infection and pathology are characterized by a strong proinflammatory response which in mice is triggered by interaction of the parasite with the toll-like receptor (TLR)/MyD88 pathway. A powerful way to identify the role of TLRs in humans is to determine whether polymorphisms at these loci influence susceptibility to T. gondii-mediated pathologies. Here we report on a small family-based study (60 families; 68 affected offspring) undertaken in Brazil which was powered for large effect sizes using single nucleotide polymorphisms with minor alleles frequencies > 0.3. Of markers in TLR2, TLR5 and TLR9 that met these criteria, we found an association Family Based Association Tests [(FBAT) Z score = 4.232; p = 1.5 x 10-5; p corrected = 1.2 x 10-4] between the C allele (frequency = 0.424; odds ratio = 7; 95% confidence interval 1.6-30.8) of rs352140 at TLR9 and toxoplasmic retinochoroiditis in Brazil. This supports the hypothesis that direct interaction between T. gondii and TLR9 may trigger proinflammatory responses that lead to severe pathologies such as the ocular disease that is associated with this infection in Brazil.

Comparative analyses of classical phenotypic method and ribosomal RNA gene sequencing for identification of medically relevant Candida species

As the distribution of Candida species and their susceptibility to antifungal agents have changed, a new means of accurately and rapidly identifying these species is necessary for the successful early resolution of infection and the subsequent reduction of morbidity and mortality. The current work aimed to evaluate ribosomal RNA gene sequencing for the identification of medically relevant Candida species in comparison with a standard phenotypic method. Eighteen reference strains (RSs), 69 phenotypically identified isolates and 20 inconclusively identified isolates were examined. Internal transcribed spaces (ITSs) and D1/D2 of the 26S ribosomal RNA gene regions were used as targets for sequencing. Additionally, the sequences of the ITS regions were used to establish evolutionary relationships. The sequencing of the ITS regions was successful for 88% (94/107) of the RS and isolates, whereas 100% of the remaining 12% (13/107) of the samples were successfully analysed by sequencing the D1/D2 region. Similarly, genotypic analysis identified all of the RS and isolates, including the 20 isolates that were not phenotypically identified. Phenotypic analysis, however, misidentified 10% (7/69) of the isolates. Phylogenetic analysis allowed the confirmation of the relationships between evolutionarily close species. Currently, the use of genotypic methods is necessary for the correct identification of Candida species.

Candidate gene study of TRAIL and TRAIL receptors: association with response to interferon beta therapy in multiple sclerosis patients.

TRAIL and TRAIL Receptor genes have been implicated in Multiple Sclerosis pathology as well as in the response to IFN beta therapy. The objective of our study was to evaluate the association of these genes in relation to the age at disease onset (AAO) and to the clinical response upon IFN beta treatment in Spanish MS patients. We carried out a candidate gene study of TRAIL, TRAILR-1, TRAILR-2, TRAILR-3 and TRAILR-4 genes. A total of 54 SNPs were analysed in 509 MS patients under IFN beta treatment, and an additional cohort of 226 MS patients was used to validate the results. Associations of rs1047275 in TRAILR-2 and rs7011559 in TRAILR-4 genes with AAO under an additive model did not withstand Bonferroni correction. In contrast, patients with the TRAILR-1 rs20576-CC genotype showed a better clinical response to IFN beta therapy compared with patients carrying the A-allele (recessive model: p = 8.88×10(-4), pc = 0.048, OR = 0.30). This SNP resulted in a non synonymous substitution of Glutamic acid to Alanine in position 228 (E228A), a change previously associated with susceptibility to different cancer types and risk of metastases, suggesting a lack of functionality of TRAILR-1. In order to unravel how this amino acid change in TRAILR-1 would affect to death signal, we performed a molecular modelling with both alleles. Neither TRAIL binding sites in the receptor nor the expression levels of TRAILR-1 in peripheral blood mononuclear cell subsets (monocytes, CD4+ and CD8+ T cells) were modified, suggesting that this SNP may be altering the death signal by some other mechanism. These findings show a role for TRAILR-1 gene variations in the clinical outcome of IFN beta therapy that might have relevance as a biomarker to predict the response to IFN beta in MS.

Examining ERBB2 as a candidate gene for susceptibility to leprosy (Hansen’s disease) in Brazil

Leprosy remains prevalent in Brazil. ErbB2 is a receptor for leprosy bacilli entering Schwann cells, which mediates Mycobacterium leprae-induced demyelination and the ERBB2 gene lies within a leprosy susceptibility locus on chromosome 17q11-q21. To determine whether polymorphisms at the ERBB2 locus contribute to this linkage peak, three haplotype tagging single nucleotide polymorphisms (tag-SNPs) (rs2517956, rs2952156, rs1058808) were genotyped in 72 families (208 cases; 372 individuals) from the state of Pará (PA). All three tag-SNPs were associated with leprosy per se [best SNP rs2517959 odds ratio (OR) = 2.22; 95% confidence interval (CI) 1.37-3.59; p = 0.001]. Lepromatous (LL) (OR = 3.25; 95% CI 1.37-7.70; p = 0.007) and tuberculoid (TT) (OR = 1.79; 95% CI 1.04-3.05; p = 0.034) leprosy both contributed to the association, which is consistent with the previous linkage to chromosome 17q11-q21 in the population from PA and supports the functional role of ErbB2 in disease pathogenesis. To attempt to replicate these findings, six SNPs (rs2517955, rs2517956, rs1810132, rs2952156, rs1801200, rs1058808) were genotyped in a population-based sample of 570 leprosy cases and 370 controls from the state of Rio Grande do Norte (RN) and the results were analysed using logistic regression analysis. However, none of the associations were replicated in the RN sample, whether analysed for leprosy per se, LL leprosy, TT leprosy, erythema nodosum leprosum or reversal reaction conditions. The role of polymorphisms at ERBB2 in controlling susceptibility to leprosy in Brazil therefore remains unclear.

Association of genetic variation with systolic and diastolic blood pressure among African Americans: the Candidate Gene Association Resource study.

The prevalence of hypertension in African Americans (AAs) is higher than in other US groups; yet, few have performed genome-wide association studies (GWASs) in AA. Among people of European descent, GWASs have identified genetic variants at 13 loci that are associated with blood pressure. It is unknown if these variants confer susceptibility in people of African ancestry. Here, we examined genome-wide and candidate gene associations with systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the Candidate Gene Association Resource (CARe) consortium consisting of 8591 AAs. Genotypes included genome-wide single-nucleotide polymorphism (SNP) data utilizing the Affymetrix 6.0 array with imputation to 2.5 million HapMap SNPs and candidate gene SNP data utilizing a 50K cardiovascular gene-centric array (ITMAT-Broad-CARe [IBC] array). For Affymetrix data, the strongest signal for DBP was rs10474346 (P= 3.6 × 10(-8)) located near GPR98 and ARRDC3. For SBP, the strongest signal was rs2258119 in C21orf91 (P= 4.7 × 10(-8)). The top IBC association for SBP was rs2012318 (P= 6.4 × 10(-6)) near SLC25A42 and for DBP was rs2523586 (P= 1.3 × 10(-6)) near HLA-B. None of the top variants replicated in additional AA (n = 11 882) or European-American (n = 69 899) cohorts. We replicated previously reported European-American blood pressure SNPs in our AA samples (SH2B3, P= 0.009; TBX3-TBX5, P= 0.03; and CSK-ULK3, P= 0.0004). These genetic loci represent the best evidence of genetic influences on SBP and DBP in AAs to date. More broadly, this work supports that notion that blood pressure among AAs is a trait with genetic underpinnings but also with significant complexity.