Haplotypes and the systematic analysis of genetic variation in genes and genomes

Analysis of 1276 Haplotype-Resolved Genomes Allows Characterization of Cis- and Trans-Abundant Genes

Many methods for haplotyping have materialized, but their application on a significant scale has been rare to date. Here we summarize analyses that were carried out in 1092 genomes from the 1000 Genomes Consortium and validated in an unprecedented number of 184 PGP genomes that have been experimentally haplotype-resolved by application of the Long-Fragment Read (LFR) technology. These analyses provided first insights into the diplotypic nature of human genomes and its potential functional implications. Thus, protein-changing variants were not randomly distributed between the two homologues of 18,121 autosomal protein-coding genes but occurred significantly more frequently in cis than in trans configurations in virtually each of the 1276 phased genomes. This resulted in global cis/trans ratios of ~60:40, establishing "cis abundance" as a universal characteristic of diploid human genomes. This phenomenon was based on two different classes of genes, a larger one exhibiting cis configurations of protein-changing variants in excess, so-called "cis-abundant" genes, and a smaller one of "trans-abundant" genes. These two gene classes, which together constitute a common diplotypic exome, were further functionally distinguished by means of gene ontology (GO) and pathway enrichment analysis. Moreover, they were distinguishable in terms of their effects on the human interactome, where they constitute distinct cis and trans modules, as shown with network propagation on a large integrated protein-protein interaction network. These analyses, recently performed with updated database and analysis tools, further consolidated the characterization of cis- and trans-abundant genes while expanding previous results. In this chapter, we present the key results along with the materials and methods to motivate readers to investigate these findings independently and gain further insights into the diplotypic nature of genes and genomes.

Genetic Fine Mapping and Genomic Annotation Defines Causal Mechanisms at A Novel Colorectal Cancer Susceptibility Locus in Han Chinese

Genome-wide association studies of colorectal cancer (CRC) have identified two risk SNPs. The characterization of these risk regions in diverse racial groups with different linkage disequilibrium structure would aid in localizing the causal variants. Herein, fine mapping of the established CRC loci was carried out in 1,508 cases and 1,482 controls obtained from the Han Chinese population. One distinct association signal was identified at these loci, where fine mapping implicated rs1010208 as a functional locus. Next, the candidate target genes of functional SNP rs1010208 were analyzed using data from TCGA databases by expression quantitative trait loci analysis method; the data from Peking University People's Hospital were utilized for verification. The dual-luciferase reporter system analysis confirmed that rs1010208 is a regulatory region that can be mutated to decrease the expression of HINT1, resulting in proliferation and invasiveness of CRC.

A strategy for building and using a human reference pangenome

In March 2019, 45 scientists and software engineers from around the world converged at the University of California, Santa Cruz for the first pangenomics codeathon. The purpose of the meeting was to propose technical specifications and standards for a usable human pangenome as well as to build relevant tools for genome graph infrastructures. During the meeting, the group held several intense and productive discussions covering a diverse set of topics, including advantages of graph genomes over a linear reference representation, design of new methods that can leverage graph-based data structures, and novel visualization and annotation approaches for pangenomes. Additionally, the participants self-organized themselves into teams that worked intensely over a three-day period to build a set of pipelines and tools for specific pangenomic applications. A summary of the questions raised and the tools developed are reported in this manuscript.

A strategy for building and using a human reference pangenome

In March 2019, 45 scientists and software engineers from around the world converged at the University of California, Santa Cruz for the first pangenomics codeathon. The purpose of the meeting was to propose technical specifications and standards for a usable human pangenome as well as to build relevant tools for genome graph infrastructures. During the meeting, the group held several intense and productive discussions covering a diverse set of topics, including advantages of graph genomes over a linear reference representation, design of new methods that can leverage graph-based data structures, and novel visualization and annotation approaches for pangenomes. Additionally, the participants self-organized themselves into teams that worked intensely over a three-day period to build a set of pipelines and tools for specific pangenomic applications. A summary of the questions raised and the tools developed are reported in this manuscript.

Significant abundance of cis configurations of coding variants in diploid human genomes

To fully understand human genetic variation and its functional consequences, the specific distribution of variants between the two chromosomal homologues of genes must be known. The 'phase' of variants can significantly impact gene function and phenotype. To assess patterns of phase at large scale, we have analyzed 18 121 autosomal genes in 1092 statistically phased genomes from the 1000 Genomes Project and 184 experimentally phased genomes from the Personal Genome Project. Here we show that genes with cis-configurations of coding variants are more frequent than genes with trans-configurations in a genome, with global cis/trans ratios of ∼60:40. Significant cis-abundance was observed in virtually all genomes in all populations. Moreover, we identified a large group of genes exhibiting cis-configurations of protein-changing variants in excess, so-called 'cis-abundant genes', and a smaller group of 'trans-abundant genes'. These two gene categories were functionally distinguishable, and exhibited strikingly different distributional patterns of protein-changing variants. Underlying these phenomena was a shared set of phase-sensitive genes of importance for adaptation and evolution. This work establishes common patterns of phase as key characteristics of diploid human exomes and provides evidence for their functional significance, highlighting the importance of phase for the interpretation of protein-coding genetic variation and gene function.

Genetic Diversity of Echinococcus granulosus Genotype G1 in Xinjiang, Northwest of China

Cystic echinococcosis (CE) caused by E. granulosus is a serious helminthic zoonosis in humans, livestock and wildlife. Xinjiang is one of high endemic province for CE in China. A total of 55 sheep and cattle livers containing echinococcal cysts were collected from slaughterhouses in Changji and Yining City, northern region of Xinjiang. PCR was employed for cloning 2 gene fragments, 12S rRNA and CO1 for analysis of phylogenetic diversity of E. granulosus. The results showed that all the samples collected were identified as G1 genotype of E. granulosus. Interestingly, YL5 and CJ75 strains were the older branches compared to those strains from France, Argentina, Australia. CO1 gene fragment showed 20 new genotype haploids and 5 new genotype haplogroups (H1-H5) by the analysis of Network 5.0 software, and the YLY17 strain was identified as the most ancestral haplotype. The major haplotypes, such as CJ75 and YL5 strains, showed identical to the isolates from Middle East. The international and domestic trade of livestock might contribute to the dispersal of different haplotypes for E. granulosus evolution.

The possible role of EZH2 and DNMT1 polymorphisms in sporadic triple-negative breast carcinoma in southern Chinese females

Triple-negative breast cancer (TNBC) has a more invasive and metastatic potential than the other types of breast cancer and hence is associated with poor prognosis. Zeste homolog 2 (EZH2) and DNA methyltransferase 1 (DNMT1) could lead to tumorigenesis by separately methylating histone H3K27 and CpG islands in tumor suppressor genes. In order to investigate the association between oncogenesis and the distribution of single nucleotide polymorphisms (SNPs) of EZH2, DNMT1, a case-control study on SNPs in TNBC cases from south China was conducted. A total of 13 SNPs were genotyped from 234 cases of TNBC tissues, and 300 normal blood samples from age-matched control group were analyzed using Snapshot technology. The expressions of EZH2 and DNMT1 were examined in the 234 cases of TNBC tissues by immunohistochemistry (IHC). The T allele of rs2288349 and the C allele of rs16999593 increase the risk of TNBC, with relative risk coefficients of 1.76 and 1.69, respectively (p < 0.001). The TC genotypes of rs2288349 and rs16999593 were higher in TNBC compared with the control group; the cancer risk increased to 5.27 and 4.13, respectively (p < 0.001). There were no significant differences between the frequencies of the other 10 SNPs and the risk of TNBC (p > 0.05). Five common haplotypes (>8 % frequency) were identified with a cumulative frequency of 96 % in the controls, while the haplotypes of AAGTAG, GGGTGA, and GACCAG were significantly increased in the control group compared to that in patients (p < 0.05). The G allele of rs10274701 significantly increased the EZH2 expression level in TNBC (p = 0.01). This is the first study to demonstrate a significant association between TNBC risk and the polymorphisms of EZH2 and DNMT1, and our researches indicate that the SNPs of EZH2 and DNMT1 are risk predictors for TNBC.

Bioinformatics and evolution of vertebrate nociceptin and opioid receptors

G protein-coupled receptors (GPCRs) are ancestrally related membrane proteins on cells that mediate the pharmacological effect of most drugs and neurotransmitters. GPCRs are the largest group of membrane receptor proteins encoded in the human genome. One of the most famous types of GPCRs is the opioid receptors. Opioid family receptors consist of four closely related proteins expressed in all vertebrate brains and spinal cords examined to date. The three classical types of opioid receptors shown unequivocally to mediate analgesia in animal models and in humans are the mu- (MOR), delta- (DOR), and kappa-(KOR) opioid receptor proteins. The fourth and most recent member of the opioid receptor family discovered is the nociceptin or orphanin FQ receptor (ORL). The role of ORL and its ligands in producing analgesia is not as clear, with both analgesic and hyperalgesic effects reported. All four opioid family receptor genes were cloned from expressed mRNA in a number of vertebrate species, and there are enough sequences presently available to carry out bioinformatic analysis. This chapter presents the results of a comparative analysis of vertebrate opioid receptors using pharmacological studies, bioinformatics, and the latest data from human whole-genome studies. Results confirm our initial hypotheses that the four opioid receptor genes most likely arose by whole-genome duplication, that there is an evolutionary vector of opioid receptor type divergence in sequence and function, and that the hMOR gene shows evidence of positive selection or adaptive evolution in Homo sapiens.

Multiple haplotype-resolved genomes reveal population patterns of gene and protein diplotypes

To fully understand human biology and link genotype to phenotype, the phase of DNA variants must be known. Here we present a comprehensive analysis of haplotype-resolved genomes to assess the nature and variation of haplotypes and their pairs, diplotypes, in European population samples. We use a set of 14 haplotype-resolved genomes generated by fosmid clone-based sequencing, complemented and expanded by up to 372 statistically resolved genomes from the 1000 Genomes Project. We find immense diversity of both haploid and diploid gene forms, up to 4.1 and 3.9 million corresponding to 249 and 235 per gene on average. Less than 15% of autosomal genes have a predominant form. We describe a 'common diplotypic proteome', a set of 4,269 genes encoding two different proteins in over 30% of genomes. We show moreover an abundance of cis configurations of mutations in the 386 genomes with an average cis/trans ratio of 60:40, and distinguishable classes of cis- versus trans-abundant genes. This work identifies key features characterizing the diplotypic nature of human genomes and provides a conceptual and analytical framework, rich resources and novel hypotheses on the functional importance of diploidy.

Whole-genome haplotyping approaches and genomic medicine

Genomic information reported as haplotypes rather than genotypes will be increasingly important for personalized medicine. Current technologies generate diploid sequence data that is rarely resolved into its constituent haplotypes. Furthermore, paradigms for thinking about genomic information are based on interpreting genotypes rather than haplotypes. Nevertheless, haplotypes have historically been useful in contexts ranging from population genetics to disease-gene mapping efforts. The main approaches for phasing genomic sequence data are molecular haplotyping, genetic haplotyping, and population-based inference. Long-read sequencing technologies are enabling longer molecular haplotypes, and decreases in the cost of whole-genome sequencing are enabling the sequencing of whole-chromosome genetic haplotypes. Hybrid approaches combining high-throughput short-read assembly with strategic approaches that enable physical or virtual binning of reads into haplotypes are enabling multi-gene haplotypes to be generated from single individuals. These techniques can be further combined with genetic and population approaches. Here, we review advances in whole-genome haplotyping approaches and discuss the importance of haplotypes for genomic medicine. Clinical applications include diagnosis by recognition of compound heterozygosity and by phasing regulatory variation to coding variation. Haplotypes, which are more specific than less complex variants such as single nucleotide variants, also have applications in prognostics and diagnostics, in the analysis of tumors, and in typing tissue for transplantation. Future advances will include technological innovations, the application of standard metrics for evaluating haplotype quality, and the development of databases that link haplotypes to disease.

Common genetic variants of the human steroid 21-hydroxylase gene (CYP21A2) are related to differences in circulating hormone levels

Purpose

Systematic evaluation of the potential relationship between the common genetic variants of CYP21A2 and hormone levels.

Methods

The relationships of CYP21A2 intron 2 polymorphisms and haplotypes with diverse baseline and stimulated blood hormone levels were studied in 106 subjects with non-functioning adrenal incidentaloma (NFAI). The rationale for using NFAI subjects is dual: i) their baseline hormone profiles do not differ from those of healthy subjects and ii) hormone levels after stimulation tests are available.

Results

The carriers (N = 27) of a well-defined CYP21A2 haplotype cluster (c5) had significantly elevated levels of cortisol (p = 0.0110), and 17-hydroxyprogesterone (p = 0.0001) after ACTH stimulation, and 11-deoxycortisol after metyrapone administration (p = 0.0017), but the hormone values were in normal ranges. In addition, the carriers (N = 33) of the C allele of the rs6462 polymorphism had a higher baseline aldosterone level (p = 0.0006). The prevalence of these genetic variants of CYP21A2 did not differ between NFAI and healthy subjects.

Conclusions

The common CYP21A2 variants presumably exert the same effect on hormone levels in the healthy and disease-affected populations. Therefore, they may contribute to complex diseases such as some cardiovascular diseases, and may influence the genotype-phenotype correlation in patients with congenital adrenal hyperplasia (CAH) including the individual need for hormone substitution.

Genetic variations of MMP9 gene and intracerebral hemorrhage outcome: a cohort study in Chinese Han population

OBJECTIVE

To investigate the association between genetic variations of matrix metalloproteinase 9 (MMP9) gene and intracerebral hemorrhage (ICH) outcome in Chinese Han population.

METHODS

The clinical data and peripheral blood samples from the patients with ICH were collected. The patients were followed up for 3 months, and poor outcome was defined as death or dependency (modified Rankin scale score of 3-6). MassARRAY Analyser was used to genotype the tagger single nucleotide polymorphisms (SNPs) of MMP9 gene. Construction of haplotypes and genetic comparisons were performed by employing PLINK 1.0.7 software.

RESULTS

181 patients with ICH were recruited between September 2009 and October 2010. Information on genetic variations and follow-up assessments were available for 169 (93.4%) patients. Independent patients at 90 days were younger than those who died or dependent (57.82 ± 13.47 vs. 66.99 ± 11.49 years, p<0.01). In addition, independent patients had lower National Institutes of Health Stroke Scale (NIHSS) score (4[1-6] vs. 13[9-19], p<0.01). There were no significant associations (all p>0.05) between MMP9 genetic variations (alleles, genotypes and haplotypes) and ICH outcome after adjustment for conventional risk factors.

CONCLUSIONS

The genetic variations of MMP9 gene were not significantly associated with ICH outcome at 90 days in Chinese Han population.

Intraspecific evolution of human RCCX copy number variation traced by haplotypes of the CYP21A2 gene

The RCCX region is a complex, multiallelic, tandem copy number variation (CNV). Two complete genes, complement component 4 (C4) and steroid 21-hydroxylase (CYP21A2, formerly CYP21B), reside in its variable region. RCCX is prone to nonallelic homologous recombination (NAHR) such as unequal crossover, generating duplications and deletions of RCCX modules, and gene conversion. A series of allele-specific long-range polymerase chain reaction coupled to the whole-gene sequencing of CYP21A2 was developed for molecular haplotyping. By means of the developed techniques, 35 different kinds of CYP21A2 haplotype variant were experimentally determined from 112 unrelated European subjects. The number of the resolved CYP21A2 haplotype variants was increased to 61 by bioinformatic haplotype reconstruction. The CYP21A2 haplotype variants could be assigned to the haplotypic RCCX CNV structures (the copy number of RCCX modules) in most cases. The genealogy network constructed from the CYP21A2 haplotype variants delineated the origin of RCCX structures. The different RCCX structures were located in tight groups. The minority of groups with identical RCCX structure occurred once in the network, implying monophyletic origin, but the majority of groups occurred several times and in different locations, indicating polyphyletic origin. The monophyletic groups were often created by single unequal crossover, whereas recurrent unequal crossover events generated some of the polyphyletic groups. As a result of recurrent NAHR events, more CYP21A2 haplotype variants with different allele patterns belonged to the same RCCX structure. The intraspecific evolution of RCCX CNV described here has provided a reasonable expectation for that of complex, multiallelic, tandem CNVs in humans.

Determination of haplotypes at structurally complex regions using emulsion haplotype fusion PCR

Background

Genotyping and massively-parallel sequencing projects result in a vast amount of diploid data that is only rarely resolved into its constituent haplotypes. It is nevertheless this phased information that is transmitted from one generation to the next and is most directly associated with biological function and the genetic causes of biological effects. Despite progress made in genome-wide sequencing and phasing algorithms and methods, problems assembling (and reconstructing linear haplotypes in) regions of repetitive DNA and structural variation remain. These dynamic and structurally complex regions are often poorly understood from a sequence point of view. Regions such as these that are highly similar in their sequence tend to be collapsed onto the genome assembly. This is turn means downstream determination of the true sequence haplotype in these regions poses a particular challenge. For structurally complex regions, a more focussed approach to assembling haplotypes may be required.

Results

In order to investigate reconstruction of spatial information at structurally complex regions, we have used an emulsion haplotype fusion PCR approach to reproducibly link sequences of up to 1kb in length to allow phasing of multiple variants from neighbouring loci, using allele-specific PCR and sequencing to detect the phase. By using emulsion systems linking flanking regions to amplicons within the CNV, this led to the reconstruction of a 59kb haplotype across the DEFA1A3 CNV in HapMap individuals.

Conclusion

This study has demonstrated a novel use for emulsion haplotype fusion PCR in addressing the issue of reconstructing structural haplotypes at multiallelic copy variable regions, using the DEFA1A3 locus as an example.

Investigation of allelic heterogeneity of the CCK-A receptor gene in paranoid schizophrenia

The cholecystokinin type A receptor (CCKAR) gene has been found to be associated with positive symptoms in patients with schizophrenia but the results reported to date are inconsistent. Considering the involvement of allelic heterogeneity in poor replication of the CCKAR finding, we genotyped five single nucleotide polymorphisms (SNPs) located in the 5' putative regulatory region of the CCKAR gene in a Chinese case-control sample and then applied the 5-SNP haplotype analysis to extract allelic heterogeneity information. The results showed that three individual haplotypes were strongly associated with increased risk of schizophrenia (corrected P = 2.9 × 10(-4), P = 2.5 × 10(-5), and P = 1.4 × 10(-5), respectively) and their combination gave an odds ratio (OR) of 6.12 with 95% CI 3.67-10.21 (P = 6.7 × 10(-15)). The haplotypes were also associated with some clinical symptoms including hallucination, suspiciousness, and hostility. Our work provided further evidence in support of the CCKAR hypothesis of schizophrenia and also suggested that haplotype-based association analysis may be a powerful approach for identification of allelic heterogeneity of a disease-underlying gene, which is very likely to be attributable to poor replication of an initial finding due to the reduction of sample power and the complexity of genetic architectures.

THE INCOMPLETE PERFECT PHYLOGENY HAPLOTYPE PROBLEM

The problem of resolving genotypes into haplotypes, under the perfect phylogeny model, has been under intensive study recently. All studies so far handled missing data entries in a heuristic manner. We prove that the perfect phylogeny haplotype problem is NP-complete when some of the data entries are missing, even when the phylogeny is rooted. We define a biologically motivated probabilistic model for genotype generation and for the way missing data occur. Under this model, we provide an algorithm, which takes an expected polynomial time. In tests on simulated data, our algorithm quickly resolves the genotypes under high rates of missing entries.

Fosmid-based whole genome haplotyping of a HapMap trio child: evaluation of Single Individual Haplotyping techniques

Determining the underlying haplotypes of individual human genomes is an essential, but currently difficult, step toward a complete understanding of genome function. Fosmid pool-based next-generation sequencing allows genome-wide generation of 40-kb haploid DNA segments, which can be phased into contiguous molecular haplotypes computationally by Single Individual Haplotyping (SIH). Many SIH algorithms have been proposed, but the accuracy of such methods has been difficult to assess due to the lack of real benchmark data. To address this problem, we generated whole genome fosmid sequence data from a HapMap trio child, NA12878, for which reliable haplotypes have already been produced. We assembled haplotypes using eight algorithms for SIH and carried out direct comparisons of their accuracy, completeness and efficiency. Our comparisons indicate that fosmid-based haplotyping can deliver highly accurate results even at low coverage and that our SIH algorithm, ReFHap, is able to efficiently produce high-quality haplotypes. We expanded the haplotypes for NA12878 by combining the current haplotypes with our fosmid-based haplotypes, producing near-to-complete new gold-standard haplotypes containing almost 98% of heterozygous SNPs. This improvement includes notable fractions of disease-related and GWA SNPs. Integrated with other molecular biological data sets, this phase information will advance the emerging field of diploid genomics.

A comprehensively molecular haplotype-resolved genome of a European individual

Independent determination of both haplotype sequences of an individual genome is essential to relate genetic variation to genome function, phenotype, and disease. To address the importance of phase, we have generated the most complete haplotype-resolved genome to date, "Max Planck One" (MP1), by fosmid pool-based next generation sequencing. Virtually all SNPs (>99%) and 80,000 indels were phased into haploid sequences of up to 6.3 Mb (N50 ~1 Mb). The completeness of phasing allowed determination of the concrete molecular haplotype pairs for the vast majority of genes (81%) including potential regulatory sequences, of which >90% were found to be constituted by two different molecular forms. A subset of 159 genes with potentially severe mutations in either cis or trans configurations exemplified in particular the role of phase for gene function, disease, and clinical interpretation of personal genomes (e.g., BRCA1). Extended genomic regions harboring manifold combinations of physically and/or functionally related genes and regulatory elements were resolved into their underlying "haploid landscapes," which may define the functional genome. Moreover, the majority of genes and functional sequences were found to contain individual or rare SNPs, which cannot be phased from population data alone, emphasizing the importance of molecular phasing for characterizing a genome in its molecular individuality. Our work provides the foundation to understand that the distinction of molecular haplotypes is essential to resolve the (inherently individual) biology of genes, genomes, and disease, establishing a reference point for "phase-sensitive" personal genomics. MP1's annotated haploid genomes are available as a public resource.

A recombination hotspot leads to sequence variability within a novel gene (AK005651) and contributes to type 1 diabetes susceptibility

More than 25 loci have been linked to type 1 diabetes (T1D) in the nonobese diabetic (NOD) mouse, but identification of the underlying genes remains challenging. We describe here the positional cloning of a T1D susceptibility locus, Idd11, located on mouse chromosome 4. Sequence analysis of a series of congenic NOD mouse strains over a critical 6.9-kb interval in these mice and in 25 inbred strains identified several haplotypes, including a unique NOD haplotype, associated with varying levels of T1D susceptibility. Haplotype diversity within this interval between congenic NOD mouse strains was due to a recombination hotspot that generated four crossover breakpoints, including one with a complex conversion tract. The Idd11 haplotype and recombination hotspot are located within a predicted gene of unknown function, which exhibits decreased expression in relevant tissues of NOD mice. Notably, it was the recombination hotspot that aided our mapping of Idd11 and confirms that recombination hotspots can create genetic variation affecting a common polygenic disease. This finding has implications for human genetic association studies, which may be affected by the approximately 33,000 estimated hotspots in the genome.

The evolution of vertebrate opioid receptors

The proteins that mediate the analgesic and other effects of opioid drugs and endogenous opioid peptides are known as opioid receptors. Opioid receptors consist of a family of four closely-related proteins belonging to the large superfamily of G-protein coupled receptors. The three types of opioid receptors shown unequivocally to mediate analgesia in animal models are the mu (MOR), delta (DOR), and kappa (KOR) opioid receptor proteins. The role of the fourth member of the opioid receptor family, the nociceptin or orphanin FQ receptor (ORL), is not as clear as hyperalgesia, analgesia, and no effect was reported after administration of ORL agonists. There are now cDNA sequences for all four types of opioid receptors that are expressed in the brain of six species from three different classes of vertebrates. This review presents a comparative analysis of vertebrate opioid receptors using bioinformatics and data from recent human genome studies. Results indicate that opioid receptors arose by gene duplication, that there is a vector of opioid receptor divergence, and that MOR shows evidence of rapid evolution.

Polymorphisms of the serotonin transporter gene (5-HTTLPR, A/G SNP in 5-HTTLPR, and STin2 VNTR) and their relation to personality traits in healthy individuals from Russia

OBJECTIVE

Numerous studies have reported association of the serotonin transporter gene (5-HTT) polymorphisms and neuroticism and traits characterizing sociability and activity. This study aimed to define a single genotype effect of three polymorphic markers in the 5-HTT gene (5-HTTLPR, A/G SNP in 5-HTTLPR and STin2 VNTR) and to check possible association of the 5-HTT haplotypes and personality traits [assessed with Eysenck Personality Inventory (EPI) and Temperament and Character Inventory (TCI) questionnaires] in 301 healthy young individuals.

METHODS

To investigate single genotype and haplotype effects of all polymorphic markers, multivariate analysis of variance and haplotype trend regression analyses were conducted correspondingly.

RESULTS

Individuals with STin2.10 allele scored significantly lower on Neuroticism (EPI) (P=0.007) and Harm Avoidance (P=0.005) in the overall sample. The same pattern of association was reported in women: carriers of STin2.10 allele scored lower on Harm Avoidance (TCI) (P=0.008). Haplotype trend regression analyses revealed that carriers of S12 haplotype had lower sociability-related traits such as Extraversion (EPI) and Novelty Seeking (TCI), whereas Harm Avoidance (TCI) (anxiety-related trait) was higher. Opposite association was observed for S10 haplotype: Extraversion (EPI) score was higher, whereas Harm Avoidance (TCI) score was lower in carriers of this haplotype.

CONCLUSION

As single polymorphism effect of STin2 was observed in relation to anxiety-related traits, opposite S10 and S12 haplotype effects on Neuroticism and Harm Avoidance could be explained by the larger impact of STin2 polymorphism. Controversially, we consider that the variance in sociability-related traits is related to specific haplotypes of 5-HTT gene.

Haplotype analysis of the endothelial nitric oxide synthase gene in asthma

Nitric oxide (NO) is an important mediator of physiologic processes in the airways. Evidence exists that genetic factors affect NO formation and contribute to the pathophysiology of asthma. The aims of this study were to determine the endothelial NO synthase (eNOS) haplotypes in Czech asthmatics and control subjects and examine their relation to asthma. We analyzed a total of six polymorphisms. Two SNPs in the promoter (C-786T and C-691T), two variants in the introns (27-bp repeat in intron 4 and G11T in intron 23), and two others in the exons (C774T in exon 6 and G894T in exon 7) were genotyped in 610 subjects (asthma, n = 294; healthy controls, n = 316), and a case-control association study was conducted. No significant differences in allele or genotype frequencies for individual polymorphisms were observed between patients with asthma and controls after correction for multiple comparisons. Nevertheless, a G to T exchange in intron 23 was related with specific sensitization for feather (p = 0.008, p(corr) < 0.05). However, the common haplotype -786T/-691C/27-bp 5 repeat variant/774C/894G/11T was associated with lower risk of asthma (p = 0.001, p(corr) < 0.05, odds ratio = 0.58, 95% confidence interval = 0.46-0.73). These findings suggest that endothelial NOS variants may be one of the factors participating in protection or susceptibility to asthma in our population.

Risk for premenstrual dysphoric disorder is associated with genetic variation in ESR1, the estrogen receptor alpha gene

BACKGROUND

Premenstrual dysphoric disorder (PMDD) is a heritable mood disorder that is triggered by gonadal steroids during the luteal phase in susceptible women.

METHODS

We performed haplotype analyses of estrogen receptors alpha and beta (ESR1 and ESR2) in 91 women with prospectively confirmed PMDD and 56 control subjects to investigate possible sources of the genetic susceptibility to affective dysregulation induced by normal levels of gonadal steroids. We also examined associations with the valine (Val)158methionine (Met) single nucleotide polymorphism (SNP) of the gene for catechol-O-methyltransferase (COMT), an enzyme involved in estrogen metabolism and prefrontal cortical activation.

RESULTS

Four SNPs in intron 4 of ESR1 showed significantly different genotype and allele distributions between patients and control subjects. Significant case-control differences were seen in sliding-window analyses of two-, three-, and four-marker haplotypes but only in those haplotypes containing SNPs in intron 4 that were positive in the single-locus analysis. No significant associations were observed with ESR2 or with the COMT Val158Met polymorphism, although the significant associations with ESR1 were observed only in those with the Val/Val genotype.

CONCLUSIONS

These are the first positive (albeit preliminary) genetic findings in this reproductive endocrine-related mood disorder and involve the receptor for a hormone that is pathogenically relevant.

Mapping translational research in personalized therapeutics: from molecular markers to health policy

Translational research is frequently used in the bioscience literature to refer to the translation of basic science into practical applications at the point of patient care. With the introduction of theragnostics, a new medical subspecialty that fuses therapeutics and diagnostic medicine with the goal of providing individualized pharmacotherapy, we suggest that the focus of translational research is shifting. We identify two bottlenecks or gaps in translational research for theragnostics: GAP1 translation from basic science to first-in-human proof-of-concept; and GAP2 translation from clinical proof-of-concept to development of evidence-based personalized treatment guidelines. GAP1 translational research in theragnostics is usually performed in traditional craft-based studies with small sample sizes and led by independent academic or industry researchers. In contrast, GAP2 translational investigations typically rely on large research consortiums and population-based biobanks that couple biomarker information with longitudinal ‘real-life’ observational data on a broad range of pharmacological phenotypes. Despite an abundance of research on the use of biobanks in disease gene discovery, there has been little conceptual work on whether and to what extent population biobanks can be utilized for translating genomics discoveries to practical treatment guidelines for theragnostic tests.

Tryptophan hydroxylase-1 gene variants associated with schizophrenia

BACKGROUND

Serotonin (5-HT) has been implicated in the pathophysiology of schizophrenia. Tryptophan hydroxylase (TPH) is the rate-limiting enzyme in the biosynthesis of serotonin (5-HT), and as such it might be related to the pathogenesis of schizophrenia. Two isoforms are known, TPH-1 and TPH-2. TPH-1 association with schizophrenia is debated.

METHODS

A case-control design was employed for gene-disease association in 155 schizophrenic psychosis patients and 253 healthy controls, all North European Caucasians. Six single nucleotide polymorphisms (SNPs) with a haplotype block structure spanning over 23 kb of the total TPH-1 29 kb were analyzed. Linkage disequilibrium and haplotype analyses were performed. Bonferroni correction was used for multiple testing.

RESULTS

Single marker association analyses showed two SNPs significantly associated with schizophrenia. Several haplotypes were associated with the disease. A "sliding window" analysis attributed the strongest disease association to a haplotype configuration localized between the promoter region and intron 3.

CONCLUSIONS

Our data indicate that TPH-1 associates with schizophrenia. It appears that specific combinations of promoter variants vis-à-vis gene transcript variants contribute to genetic predisposition to the disease.

Tryptophan hydroxylase-1 gene variants associate with a group of suicidal borderline women

Alterations in the serotonin (5-HT) system have been related to impulsive aggression and suicidal behavior, common features of the borderline personality disorder (BPD). Tryptophan hydroxylase (TPH) is the rate-limiting enzyme in 5-HT biosynthesis. Two isoforms are known, TPH-1 and TPH-2. TPH-1 has been correlated to various psychiatric and behavioral disorders by gene polymorphism association studies. We aimed to determine whether specific TPH-1 haplotypes associate with BPD. A case-control design was employed. The control group included 98 women without psychiatric history. In all, 95 patients were included, all Caucasian women with a BPD diagnosis who had attempted suicide at least twice during their lifetime. Exclusion criteria were: (i) substance dependence; (ii) dementia or other irreversible organic brain syndromes; (iii) psychotic disorders or major depressive illness with melancholic features; (iv) life-threatening eating disorders. Six single-nucleotide polymorphisms (SNPs) were found at significant linkage disequilibrium across 23 kb of the TPH-1 gene in both patients and controls, suggesting a haplotype block structure. While no individual SNP showed association, several haplotypes associated with the BPD group. In particular, one six-SNP haplotype was absent from the control group while representing about one-quarter of all haplotypes in the BPD group (corrected P<<10(-5)). A 'sliding window' analysis attributed the strongest disease association to haplotype configurations located between the gene promoter and intron 3. We conclude that TPH-1 associates with BPD in suicidal women. Our data support the expectation that haplotype analysis is superior to single locus analysis in gene-disease, case-control association studies.

Haplotype analysis reveals tryptophan hydroxylase (TPH) 1 gene variants associated with major depression

BACKGROUND

Tryptophan hydroxylase (TPH) is the rate-limiting enzyme in the biosynthesis of serotonin (5-HT) and might be related to the pathogenesis of major depression (MD). Two isoforms are known, TPH-1 and TPH-2. Tryptophan hydroxylase-1 association with MD is still debated.

METHODS

A single nucleotide polymorphism (SNP) screening strategy was used to define TPH-1 haplotypes spanning over 23 kilobase (kb) of the 29 kb gene length. Genotyping was performed in 228 MD patients and 253 healthy control subjects.

RESULTS

Six SNPs were found at linkage disequilibrium in both patients and control subjects, suggesting a haplotype block structure. Single marker association analyses showed only one SNP significantly associated with MD. Several haplotypes were associated with MD. When all six locus haplotypes were divided into two groups, above or below a 5% threshold, the compound haplotype group below a 5% frequency resulted as associated with the disease (31.6% vs. 18.0% in control subjects, p < 10(-5)). A "sliding window" analysis attributed the strongest disease association to a haplotype configuration localized between introns 7 and 8 (p < 10(-5)).

CONCLUSIONS

Haplotype analysis indicates that TPH-1 associates with MD. The most common TPH-1 variants appear to carry no risk, while some of the less frequent variants might contribute to genetic predisposition to MD.

Common Genetic Variation at the Endothelial Nitric Oxide Synthase Locus and Relations to Brachial Artery Vasodilator Function in the Community

Background— Sequence variants at the endothelial nitric oxide synthase ( NOS3 ) locus have been associated with endothelial function measures, but replication has been limited.

Methods and Results— In reference pedigrees, we characterized linkage disequilibrium structure at the NOS3 locus using 33 common single nucleotide polymorphisms (SNPs). Eighteen SNPs that capture underlying common variation were genotyped in unrelated Framingham Heart Study participants (49.5% women; mean age, 62 years) with measured brachial artery flow-mediated dilation (n=1446) or hyperemic flow velocity (n=1043). Within 3 defined blocks of strong linkage disequilibrium that spanned NOS3 , 11 SNPs captured >80% of common haplotypic variation. Among men, there were nominally significant associations between 8 NOS3 SNPs (minimum P =0.002) and between haplotypes (minimum P =0.002) and either flow-mediated dilation or hyperemic flow velocity. In women, we did not observe significant associations between NOS3 SNPs or haplotypes and endothelial function measures. To correct for multiple testing, we constructed 1000 bootstrapped null data sets and found that empirical probability values exceeded 0.05 for both phenotypes.

Conclusions— A parsimonious set of SNPs captures common genetic variation at the NOS3 locus. A conservative interpretation of our results is that, accounting for multiple testing, we did not observe statistically significant relations between NOS3 sequence variants and endothelial function measures in either sex. The nominal associations of select NOS3 variants with endothelial function in men (unadjusted for multiple testing) should be viewed as hypothesis-generating observations and may merit testing in other cohorts and experimental designs.

The genome Assembly Archive: a new public resource

With the genome assembly archive, it is possible to examine the raw data that underlies the DNA sequence in any sequenced genome

Oestrogen receptor alpha gene haplotype and postmenopausal breast cancer risk: a case control study

INTRODUCTION

Oestrogen receptor alpha, which mediates the effect of oestrogen in target tissues, is genetically polymorphic. Because breast cancer development is dependent on oestrogenic influence, we have investigated whether polymorphisms in the oestrogen receptor alpha gene (ESR1) are associated with breast cancer risk.

METHODS

We genotyped breast cancer cases and age-matched population controls for one microsatellite marker and four single-nucleotide polymorphisms (SNPs) in ESR1. The numbers of genotyped cases and controls for each marker were as follows: TAn, 1514 cases and 1514 controls; c.454-397C --> T, 1557 cases and 1512 controls; c.454-351A --> G, 1556 cases and 1512 controls; c.729C --> T, 1562 cases and 1513 controls; c.975C --> G, 1562 cases and 1513 controls. Using logistic regression models, we calculated odds ratios (ORs) and 95% confidence intervals (CIs). Haplotype effects were estimated in an exploratory analysis, using expectation-maximisation algorithms for case-control study data.

RESULTS

There were no compelling associations between single polymorphic loci and breast cancer risk. In haplotype analyses, a common haplotype of the c.454-351A --> G or c.454-397C --> T and c.975C --> G SNPs appeared to be associated with an increased risk for ductal breast cancer: one copy of the c.454-351A --> G and c.975C --> G haplotype entailed an OR of 1.19 (95% CI 1.06-1.33) and two copies with an OR of 1.42 (95% CI 1.15-1.77), compared with no copies, under a model of multiplicative penetrance. The association with the c.454-397C --> T and c.975C --> G haplotypes was similar. Our data indicated that these haplotypes were more influential in women with a high body mass index. Adjustment for multiple comparisons rendered the associations statistically non-significant.

CONCLUSION

We found suggestions of an association between common haplotypes in ESR1 and the risk for ductal breast cancer that is stronger in heavy women.

Genetic variation and pharmacogenomics: concepts, facts, and challenges

The analysis of genetic variation in candidate genes is an issue of central importance in pharmacogenomics. The specific approaches taken will have a critical impact on the successful identification of disease genes, the molecular correlates of drug response, and the establishment of meaningful relationships between genetic variants and phenotypes of biomedical and pharmaceutical importance in general. Against a historical background, this article distinguishes different approaches to candidate gene analysis, reflecting different stages in human genome research. Only recently has it become feasible to analyze genetic variation systematically at the ultimate level of resolution, ie, the DNA sequence. In this context, the importance of haplotype-based approaches to candidate gene analysis has at last been recognized; the determination of the specific combinations of variants for each of the two sequences of a gene defined as a haplotype is essential. An up-to-date summary of such maximum resolution data on the amount, nature, and structure of genetic variation in candidate genes will be given. These data demonstrate abundant gene sequence and haplotype diversity. Numerous individually different forms of a gene may exist. This presents major challenges to the analysis of relationships between genetic variation, gene function, and phenotype. First solutions seem within reach. The implications of naturally occurring variation for pharmacogenomics and “personalized” medicine are now evident. Future approaches to the identification, evaluation, and prioritization of drug targets, the optimization of clinical trials, and the development of efficient therapies must be based on in-depth knowledge of candidate gene variation as an essential prerequisite.

Haplotype analysis and identification of genes for a complex trait: examples from schizophrenia

For more than a decade there has been intensive research into the genetic etiology of schizophrenia, yet it is only recently that the first findings of promising genes associating with the disorder have been reported. Linkage analyses in families collected from different populations have provided relatively well defined genomic loci. These have been typically followed by fine mapping studies using single nucleotide polymorphisms (SNPs). A number of analysis programs have been produced to test SNPs and their haplotypes for association. Typically association has been established to specific haplotypes representing an allelic variant of the corresponding gene. The inherent problem of multiple testing in the analysis of haplotypes needs to be addressed fully, to determine if any of these recent findings can be considered as confirmed susceptibility genes for schizophrenia. However, informative haplotypes have provided a way to define allelic variants of genes associated with schizophrenia in numerous study samples, and are a useful tool in characterizing the extent of allelic diversity of putative schizophrenia susceptibility genes within different populations.