A genomewide search for type 2 diabetes-susceptibility genes in indigenous Australians

The prevalence of type 2 diabetes among Australian residents is 7.5%; however, prevalence rates up to six times higher have been reported for indigenous Australian communities. Epidemiological evidence implicates genetic factors in the susceptibility of indigenous Australians to type 2 diabetes and supports the hypothesis of the thrifty genotype, but, to date, the nature of the genetic predisposition is unknown. We have ascertained clinical details from a community of indigenous Australian descent in North Stradbroke Island, Queensland. In this population, the phenotype is characterized by severe insulin resistance. We have conducted a genomewide scan, at an average resolution of 10 cM, for type 2 diabetes-susceptibility genes in a large multigeneration pedigree from this community. Parametric linkage analysis undertaken using FASTLINK version 4.1p yielded a maximum two-point LOD score of +2.97 at marker D2S2345. Multipoint analysis yielded a peak LOD score of +3.9

A genomewide survey of developmentally relevant genes in Ciona intestinalis - IV. Genes for HMG transcriptional regulators, bZip and GATA/Gli/Zic/Snail

Many kinds of transcription factors and regulators play key roles in a variety of developmental processes. In the present survey, genes encoding proteins with conserved HMG-box, bZip domains, and some types of zinc finger motifs were surveyed in the completely sequenced genome of Ciona intestinalis. In the present analysis, 21 HMG-box-containing genes and 26 bZip genes were identified as well as four small groups of zinc finger genes in the Ciona genome. The results also showed that a less redundant set of genes is present in the Ciona genome compared with vertebrate genomes. In addition, cDNA clones for almost all genes identified have been cloned and distributed as a Ciona intestinalis Gene Collection Release I. The present comprehensive analysis therefore provides a means to study the role of these transcription factors in developmental processes of basal chordates.

A genomewide survey of developmentally relevant genes in Ciona intestinalis - V. Genes for receptor tyrosine kinase pathway and Notch signaling pathway

In the present survey, we identified most of the genes involved in the receptor tyrosine kinase (RTK), mitogen activated protein kinase (MAPK) and Notch signaling pathways in the draft genome sequence of Ciona intestinalis, a basal chordate. Compared to vertebrates, most of the genes found in the Ciona genome had fewer paralogues, although several genes including ephrin, Eph and fringe appeared to have multiplied or duplicated independently in the ascidian genome. In contrast, some genes including kit/flt, PDGF and Trk receptor tyrosine kinases were not found in the present survey, suggesting that these genes are innovations in the vertebrate lineage or lost in the ascidian lineage. The gene set identified in the present analysis provides an insight into genes for the RTK, MAPK and Notch signaling pathways in the ancient chordate genome and thereby how chordates evolved these signaling pathway.

A genomewide association study of smoking relapse in four European population-based samples.

OBJECTIVES: Genomewide association studies (GWAS) have identified clear evidence of genetic markers for nicotine dependence. Other smoking phenotypes have been tested, but the results are less consistent. The tendency to relapse versus the ability to maintain long-term abstinence has received little attention in genetic studies; thus, our aim was to provide a better biological understanding of this phenotype through the identification of genetic loci associated with smoking relapse. METHODS: We carried out a GWAS on data from two European population-based collections, including a total of 835 cases (relapsers) and 990 controls (abstainers). Top-ranked findings from the discovery phase were tested for replication in two additional independent European population-based cohorts. RESULTS: Of the seven top markers from the discovery phase, none were consistently associated with smoking relapse across all samples and none reached genomewide significance. A single-nucleotide polymorphism rs1008509, within the Xylosyltransferase II (XYLT2) gene, was suggestively associated with smoking relapse in the discovery phase (β=-0.504; P=5.6E-06) and in the first replication sample (ALSPAC) (β=-0.27; P=0.004; n=1932), but not in the second sample (KORA) (β=0.19; P=0.138; n=912). We failed to identify an association between loci implicated previously in other smoking phenotypes and smoking relapse. CONCLUSION: Although no genomewide significant findings emerged from this study, we found that loci implicated in other smoking phenotypes were not associated with smoking relapse, which suggests that the neurobiology of smoking relapse and long-term abstinence may be distinct from biological mechanisms implicated in the development of nicotine dependence.

Genomewide association study using a high-density single nucleotide polymorphism array and case-control design identifies a novel essential hypertension susceptibility locus in the promoter region of endothelial NO synthase.

Essential hypertension is a multifactorial disorder and is the main risk factor for renal and cardiovascular complications. The research on the genetics of hypertension has been frustrated by the small predictive value of the discovered genetic variants. The HYPERGENES Project investigated associations between genetic variants and essential hypertension pursuing a 2-stage study by recruiting cases and controls from extensively characterized cohorts recruited over many years in different European regions. The discovery phase consisted of 1865 cases and 1750 controls genotyped with 1M Illumina array. Best hits were followed up in a validation panel of 1385 cases and 1246 controls that were genotyped with a custom array of 14 055 markers. We identified a new hypertension susceptibility locus (rs3918226) in the promoter region of the endothelial NO synthase gene (odds ratio: 1.54 [95% CI: 1.37-1.73]; combined P=2.58 · 10(-13)). A meta-analysis, using other in silico/de novo genotyping data for a total of 21 714 subjects, resulted in an overall odds ratio of 1.34 (95% CI: 1.25-1.44; P=1.032 · 10(-14)). The quantitative analysis on a population-based sample revealed an effect size of 1.91 (95% CI: 0.16-3.66) for systolic and 1.40 (95% CI: 0.25-2.55) for diastolic blood pressure. We identified in silico a potential binding site for ETS transcription factors directly next to rs3918226, suggesting a potential modulation of endothelial NO synthase expression. Biological evidence links endothelial NO synthase with hypertension, because it is a critical mediator of cardiovascular homeostasis and blood pressure control via vascular tone regulation. This finding supports the hypothesis that there may be a causal genetic variation at this locus.

A genomewide linkage study on suicidality in major depressive disorder confirms evidence for linkage to 2p12.

Suicidal behavior is commonly associated with depression. Twin studies indicate that both suicidality and major depressive disorder (MDD) are heritable. However, epidemiological evidence suggests that the inheritance of suicidality is likely to be independent of the underlying psychiatric disorder, implying a distinct genetic contribution to suicidality. We conducted a genomewide linkage search aiming to detect genomic loci that may harbor susceptibility genes contributing to risk for suicidality in recurrent MDD. Affected sibling pair (ASP) variance components analysis was performed using the Depression Network cohort of 971 ASPs. The quantitative trait measuring suicidality as a broad phenotype, encompassing ideation and suicide attempts, was established from Schedules for Clinical Assessment in Neuropsychiatry interview items. We examined 1,060 genotyped microsatellite markers with an average spacing of 3.3 cM. Empirical thresholds for linkage evidence were set by whole-genome simulations (LOD = 2.71 for genomewide significance, 1.71 for suggestive linkage). No genomewide significant findings were found. Marker D3S1234 on 3p14 achieved suggestive linkage and yielded a maximum LOD of 1.853 (P = 0.0017), loci 9p24.3 and 18q22-q23 achieved LOD scores >1.5. We found some support for linkage to 2p12 (LOD = 1.2, P = 0.0087) which was previously implicated in linkage studies of suicidality. Our follow-up meta-analysis of five studies showed strong linkage to this region (P = 2 × 10(-6) ). In conclusion, this study analyzed suicidality as a continuous trait in MDD. We found modest evidence for linkage on 3p14. Our meta-analysis supports previous evidence of linkage to suicidality on 2p12. Some candidate genes in these regions may plausibly be implicated in suicidality.

Genomewide association study of atazanavir pharmacokinetics and hyperbilirubinemia in AIDS Clinical Trials Group protocol A5202.

BACKGROUND: Atazanavir-associated hyperbilirubinemia can cause premature discontinuation of atazanavir and avoidance of its initial prescription. We used genomewide genotyping and clinical data to characterize determinants of atazanavir pharmacokinetics and hyperbilirubinemia in AIDS Clinical Trials Group protocol A5202. METHODS: Plasma atazanavir pharmacokinetics and indirect bilirubin concentrations were characterized in HIV-1-infected patients randomized to atazanavir/ritonavir-containing regimens. A subset had genomewide genotype data available. RESULTS: Genomewide assay data were available from 542 participants, of whom 475 also had data on estimated atazanavir clearance and relevant covariates available. Peak bilirubin concentration and relevant covariates were available for 443 participants. By multivariate analysis, higher peak on-treatment bilirubin levels were found to be associated with the UGT1A1 rs887829 T allele (P=6.4×10), higher baseline hemoglobin levels (P=4.9×10), higher baseline bilirubin levels (P=6.7×10), and slower plasma atazanavir clearance (P=8.6×10). For peak bilirubin levels greater than 3.0 mg/dl, the positive predictive value of a baseline bilirubin level of 0.5 mg/dl or higher with hemoglobin concentrations of 14 g/dl or higher was 0.51, which increased to 0.85 with rs887829 TT homozygosity. For peak bilirubin levels of 3.0 mg/dl or lower, the positive predictive value of a baseline bilirubin level less than 0.5 mg/dl with a hemoglobin concentration less than 14 g/dl was 0.91, which increased to 0.96 with rs887829 CC homozygosity. No polymorphism predicted atazanavir pharmacokinetics at genomewide significance. CONCLUSION: Atazanavir-associated hyperbilirubinemia is best predicted by considering UGT1A1 genotype, baseline bilirubin level, and baseline hemoglobin level in combination. Use of ritonavir as a pharmacokinetic enhancer may have abrogated genetic associations with atazanavir pharmacokinetics.

A genomewide scan for early-onset coronary artery disease in 438 families: the GENECARD Study.

A family history of coronary artery disease (CAD), especially when the disease occurs at a young age, is a potent risk factor for CAD. DNA collection in families in which two or more siblings are affected at an early age allows identification of genetic factors for CAD by linkage analysis. We performed a genomewide scan in 1,168 individuals from 438 families, including 493 affected sibling pairs with documented onset of CAD before 51 years of age in men and before 56 years of age in women. We prospectively defined three phenotypic subsets of families: (1) acute coronary syndrome in two or more siblings; (2) absence of type 2 diabetes in all affected siblings; and (3) atherogenic dyslipidemia in any one sibling. Genotypes were analyzed for 395 microsatellite markers. Regions were defined as providing evidence for linkage if they provided parametric two-point LOD scores >1.5, together with nonparametric multipoint LOD scores >1.0. Regions on chromosomes 3q13 (multipoint LOD = 3.3; empirical P value <.001) and 5q31 (multipoint LOD = 1.4; empirical P value <.081) met these criteria in the entire data set, and regions on chromosomes 1q25, 3q13, 7p14, and 19p13 met these criteria in one or more of the subsets. Two regions, 3q13 and 1q25, met the criteria for genomewide significance. We have identified a region on chromosome 3q13 that is linked to early-onset CAD, as well as additional regions of interest that will require further analysis. These data provide initial areas of the human genome where further investigation may reveal susceptibility genes for early-onset CAD.

Sequentially testing for a gene-drug interaction in a genomewide analysis

Assaying a large number of genetic markers from patients in clinical trials is now possible in order to tailor drugs with respect to efficacy. The statistical methodology for analysing such massive data sets is challenging. The most popular type of statistical analysis is to use a univariate test for each genetic marker, once all the data from a clinical study have been collected. This paper presents a sequential method for conducting an omnibus test for detecting gene-drug interactions across the genome, thus allowing informed decisions at the earliest opportunity and overcoming the multiple testing problems from conducting many univariate tests. We first propose an omnibus test for a fixed sample size. This test is based on combining F-statistics that test for an interaction between treatment and the individual single nucleotide polymorphism (SNP). As SNPs tend to be correlated, we use permutations to calculate a global p-value. We extend our omnibus test to the sequential case. In order to control the type I error rate, we propose a sequential method that uses permutations to obtain the stopping boundaries. The results of a simulation study show that the sequential permutation method is more powerful than alternative sequential methods that control the type I error rate, such as the inverse-normal method. The proposed method is flexible as we do not need to assume a mode of inheritance and can also adjust for confounding factors. An application to real clinical data illustrates that the method is computationally feasible for a large number of SNPs. Copyright (c) 2007 John Wiley & Sons, Ltd.

A genomewide survey of developmentally relevant genes in Ciona intestinalis - III. Genes for Fox, ETS, nuclear receptors and NF kappa B

A survey against the draft genome sequence and the cDNA/EST database of Ciona intestinalis identified a number of genes encoding transcription factors regulating a variety of processes including development. In the present study, we describe almost complete sets of genes for Fox, ETS-domain transcription factors, nuclear receptors, and NFkappaB as well as other factors regulating NFkappaB activity, with their phylogenetic nature. Vertebrate Fox transcription factors are currently delineated into 17 subfamilies: FoxA to FoxQ. The present survey yielded 29 genes of this family in the Ciona genome, 24 of which were Ciona orthologues of known Fox genes. In addition, we found 15 ETS aenes, 17 nuclear receptor genes, and several NFkappaB signaling pathway genes in the Ciona genome. The number of Ciona genes in each family is much smaller than that of vertebrates, which represents a simplified feature of the ascidian genome. For example, humans have two NFkappaB genes, three Rel genes, and five NFAT genes, while Ciona has one gene for each family. The Ciona genome also contains smaller numbers of genes for the NFkappaB regulatory system, i.e. after the split of ascidians/vertebrates, vertebrates evolved a more complex NFkappaB system. The present results therefore provide molecular information for the investigation of complex developmental processes, and an insight into chordate evolution.

A genomewide association mapping study using ultrasound-scanned information identifies potential genomic regions and candidate genes affecting carcass traits in Nellore cattle

The aim of this study was to identify candidate genes and genomic regions associated with ultrasound-derived measurements of the rib-eye area (REA), backfat thickness (BFT) and rumpfat thickness (RFT) in Nellore cattle. Data from 640 Nellore steers and young bulls with genotypes for 290 863 single nucleotide polymorphisms (SNPs) were used for genomewide association mapping. Significant SNP associations were explored to find possible candidate genes related to physiological processes. Several of the significant markers detected were mapped onto functional candidate genes including ARFGAP3, CLSTN2 and DPYD for REA; OSBPL3 and SUDS3 for BFT; and RARRES1 and VEPH1 for RFT. The physiological pathway related to lipid metabolism (CLSTN2, OSBPL3, RARRES1 and VEPH1) was identified. The significant markers within previously reported QTLs reinforce the importance of the genomic regions, and the other loci offer candidate genes that have not been related to carcass traits in previous investigations.

Genomewide association for a dominant pigmentation gene in sheep

Most published genomewide association studies (GWAS) in sheep have investigated recessively inherited monogenic traits. The objective here was to assess the feasibility of performing GWAS for a dominant trait for which the genetic basis was already known. A total of 42 Manchega and Rasa Aragonesa sheep that segregate solid black or white coat pigmentation were genotyped using the SNP50 BeadChip. Previous analysis in Manchegas demonstrated a complete association between the pigmentation trait and alleles of the MC1R gene, setting an a priori expectation for GWAS. Multiple methods were used to identify and quantify the strength of population substructure between black and white animals, before allelic association testing was performed for 49 034 SNPs. Following correction for substructure, GWAS identified the most strongly associated SNP (s26449) was also the closest to the MC1R gene. The finding was strongly supported by the permutation tree-based random forest (RF) analysis. Importantly, GWAS identified unlinked SNP with only slightly lower p-values than for s26449. Random forest analysis indicated these were false positives, suggesting interpretation based on both approaches was beneficial. The results indicate that a combined analytical approach can be successful in studies where a modest number of animals are available and substantial population stratification exists.

A genomewide survey of basic helix–loop–helix factors in Drosophila

The basic helix–loop–helix (bHLH) transcription factors play important roles in the specification of tissue type during the development of animals. We have used the information contained in the recently published genomic sequence of Drosophila melanogaster to identify 12 additional bHLH proteins. By sequence analysis we have assigned these proteins to families defined by Atonal, Hairy-Enhancer of Split, Hand, p48, Mesp, MYC/USF, and the bHLH-Per, Arnt, Sim (PAS) domain. In addition, one single protein represents a unique family of bHLH proteins. mRNA in situ analysis demonstrates that the genes encoding these proteins are expressed in several tissue types but are particularly concentrated in the developing nervous system and mesoderm.

Genomewide significant linkage to migrainous Hheadache on chromosome 5q21

Familial typical migraine is a common, complex disorder that shows strong familial aggregation. Using latent-class analysis (LCA), we identified subgroups of people with migraine/severe headache in a community sample of 12,245 Australian twins (60% female), drawn from two cohorts of individuals aged 23-90 years who completed an interview based on International Headache Society criteria. We report results from genomewide linkage analyses involving 756 twin families containing a total of 790 independent sib pairs ( 130 affected concordant, 324 discordant, and 336 unaffected concordant for LCA-derived migraine). Quantitative-trait linkage analysis produced evidence of significant linkage on chromosome 5q21 and suggestive linkage on chromosomes 8, 10, and 13. In addition, we replicated previously reported typical-migraine susceptibility loci on chromosomes 6p12.2-p21.1 and 1q21-q23, the latter being within 3 cM of the rare autosomal dominant familial hemiplegic migraine gene (ATP1A2), a finding which potentially implicates ATP1A2 in familial typical migraine for the first time. Linkage analyses of individual migraine symptoms for our six most interesting chromosomes provide tantalizing hints of the phenotypic and genetic complexity of migraine. Specifically, the chromosome 1 locus is most associated with phonophobia; the chromosome 5 peak is predominantly associated with pulsating headache; the chromosome 6 locus is associated with activity-prohibiting headache and photophobia; the chromosome 8 locus is associated with nausea/vomiting and moderate/severe headache; the chromosome 10 peak is most associated with phonophobia and photophobia; and the chromosome 13 peak is completely due to association with photophobia. These results will prove to be invaluable in the design and analysis of future linkage and linkage disequilibrium studies of migraine.