Use of molecular variation in the NCBI dbSNP database

Endothelial MAP2K1 mutations in arteriovenous malformation activate the RAS/MAPK pathway

Arteriovenous malformation (AVM) is a locally destructive congenital vascular anomaly caused by somatic mutations in MAP2K1. The mutation is isolated to endothelial cells (ECs). The purpose of this study was to determine the effects of mutant MAP2K1 on EC signaling and vascular network formation. Pathway effects were studied using both mutant MAP2K1 (K57N) human AVM tissue and human umbilical vein endothelial cells (HUVECs) engineered to overexpress the MAP2K1 (K57N) mutation. Western blot was used to determine cell signaling along the RAS/MAPK pathway. Geltrex tube formation assays were performed to assess EC vascular network formation. Cells were treated with a MAP2K1 inhibitor (Trametinib) to determine its effect on signaling and vascular tube formation. Human mutant MAP2K1-AVM ECs had similar baseline MEK1 and ERK1/2 expression with controls; however, mutant MAP2K1-AVM ECs produced significantly more phosphorylated ERK1/2 than wild-type ECs. Mutant MAP2K1 HUVECs demonstrated significantly more phosphorylated ERK1/2 than control HUVECs. Trametinib reduced the phosphorylation of ERK1/2 in mutant cells and prevented the ability of ECs to form vascular networks. AVM MAP2K1 mutations activate RAS/MAPK signaling in ECs. ERK activation and vascular network formation are reduced with Trametinib. Pharmacotherapy using MAP2K1 inhibitors may prevent the formation or progression of AVMs.

ELMOD3, a novel causative gene, associated with human autosomal dominant nonsyndromic and progressive hearing loss

Autosomal dominant nonsyndromic hearing loss (ADNSHL) is a highly genetically heterogeneous disorder. Up to date only approximately 37 ADNSHL-causing genes have been identified. The goal of this study was to determine the causative gene in a five-generation Chinese family with ADNSHL. A Chinese family was ascertained. Simultaneously, two affected individuals and one normal hearing control from the family were analyzed by whole exome capture sequencing. To assess the functional effect of the identified variant, in-vitro studies were performed. novel missense variant, c.512A>G (p.His171Arg) in exon 8 of the ELMO domain-containing 3 (ELMOD3) gene, was identified as a causative variant in this family affected by late-onset and progressive ADNSHL. The variant was validated by Sanger sequencing and found to co-segregate with the phenotype within the pedigree and was absent in 500 ethnically matched unrelated normal hearing control subjects. To our knowledge, this is the first report of a family with ADNSHL caused by ELMOD3 mutation. Western blots and immunofluorescence staining demonstrated that p.His171Arg resulted in abnormal expression levels of ELMOD3 and abnormal subcellular localization. Furthermore, the analysis of the stability of the wild-type (WT) and mutant ELMOD3 protein shows that the decay of p.His171Arg is faster than that of the WT, suggesting a shorter halflife of the c.512A > G variant. A novel variant in the ELMOD3 gene, encoding a member of the engulfment and cell motility (ELMO) family of GTPase-activating proteins, was identified for the first time as responsible for ADNSHL.

Somatic MAP2K1 Mutations Are Associated with Extracranial Arteriovenous Malformation

Arteriovenous malformation (AVM) is a fast-flow, congenital vascular anomaly that may arise anywhere in the body. AVMs typically progress, causing destruction of surrounding tissue and, sometimes, cardiac overload. AVMs are difficult to control; they often re-expand after embolization or resection, and pharmacologic therapy is unavailable. We studied extracranial AVMs in order to identify their biological basis. We performed whole-exome sequencing (WES) and whole-genome sequencing (WGS) on AVM tissue from affected individuals. Endothelial cells were separated from non-endothelial cells by immune-affinity purification. We used droplet digital PCR (ddPCR) to confirm mutations found by WES and WGS, to determine whether mutant alleles were enriched in endothelial or non-endothelial cells, and to screen additional AVM specimens. In seven of ten specimens, WES and WGS detected and ddPCR confirmed somatic mutations in mitogen activated protein kinase kinase 1 (MAP2K1), the gene that encodes MAP-extracellular signal-regulated kinase 1 (MEK1). Mutant alleles were enriched in endothelial cells and were not present in blood or saliva. 9 of 15 additional AVM specimens contained mutant MAP2K1 alleles. Mutations were missense or small in-frame deletions that affect amino acid residues within or adjacent to the protein's negative regulatory domain. Several of these mutations have been found in cancers and shown to increase MEK1 activity. In summary, somatic mutations in MAP2K1 are a common cause of extracranial AVM. The likely mechanism is endothelial cell dysfunction due to increased MEK1 activity. MEK1 inhibitors, which are approved to treat several forms of cancer, are potential therapeutic agents for individuals with extracranial AVM.

Loci associated with adult stature also affect calf birth survival in cattle

Background

Understanding the underlying pleiotropic relationships among quantitative traits is necessary in order to predict correlated responses to artificial selection. The availability of large-scale next-generation sequence data in cattle has provided an opportunity to examine whether pleiotropy is responsible for overlapping QTL in multiple economic traits. In the present study, we examined QTL affecting cattle stillbirth, calf size, and adult stature located in the same genomic region.

Results

A genome scan using imputed whole genome sequence variants revealed one QTL with large effects on the service sire calving index (SCI), and body conformation index (BCI) at the same location (~39 Mb) on chromosome 6 in Nordic Red cattle. The targeted region was analyzed for SCI and BCI component traits. The QTL peak included LCORL and NCAPG genes, which had been reported to influence fetal growth and adult stature in several species. The QTL exhibited large effects on calf size and stature in Nordic Red cattle. Two deviant haplotypes (HAP1 and HAP2) were resolved which increased calf size at birth, and affected adult body conformation. However, the haplotypes also resulted in increased calving difficulties and calf mortality due to increased calf size at birth. Haplotype locations overlapped, however linkage disequilibrium (LD) between the sites was low, suggesting that two independent mutations were responsible for similar effects. The difference in prevalence between the two haplotypes in Nordic Red subpopulations suggested independent origins in different populations.

Conclusions

Results of our study identified QTL with large effects on body conformation and service sire calving traits on chromosome 6 in cattle. We present robust evidence that variation at the LCORL and NCAPG locus affects calf size at birth and adult stature. We suggest the two deviant haplotypes within the QTL were due to two independent mutations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0202-3) contains supplementary material, which is available to authorized users.

Genome-wide association study using high-density single nucleotide polymorphism arrays and whole-genome sequences for clinical mastitis traits in dairy cattle

Mastitis is a mammary disease that frequently affects dairy cattle. Despite considerable research on the development of effective prevention and treatment strategies, mastitis continues to be a significant issue in bovine veterinary medicine. To identify major genes that affect mastitis in dairy cattle, 6 chromosomal regions on Bos taurus autosome (BTA) 6, 13, 16, 19, and 20 were selected from a genome scan for 9 mastitis phenotypes using imputed high-density single nucleotide polymorphism arrays. Association analyses using sequence-level variants for the 6 targeted regions were carried out to map causal variants using whole-genome sequence data from 3 breeds. The quantitative trait loci (QTL) discovery population comprised 4,992 progeny-tested Holstein bulls, and QTL were confirmed in 4,442 Nordic Red and 1,126 Jersey cattle. The targeted regions were imputed to the sequence level. The highest association signal for clinical mastitis was observed on BTA 6 at 88.97 Mb in Holstein cattle and was confirmed in Nordic Red cattle. The peak association region on BTA 6 contained 2 genes: vitamin D-binding protein precursor (GC) and neuropeptide FF receptor 2 (NPFFR2), which, based on known biological functions, are good candidates for affecting mastitis. However, strong linkage disequilibrium in this region prevented conclusive determination of the causal gene. A different QTL on BTA 6 located at 88.32 Mb in Holstein cattle affected mastitis. In addition, QTL on BTA 13 and 19 were confirmed to segregate in Nordic Red cattle and QTL on BTA 16 and 20 were confirmed in Jersey cattle. Although several candidate genes were identified in these targeted regions, it was not possible to identify a gene or polymorphism as the causal factor for any of these regions.

Ring finger protein 39 genetic variants associate with HIV-1 plasma viral loads and its replication in cell culture

BACKGROUND

The human immunodeficiency virus (HIV-1) exploits host proteins to complete its life cycle. Genome-wide siRNA approaches suggested that host proteins affect HIV-1 replication. However, the results barely overlapped. RING finger protein 39 (RNF39) has been identified from genome-wide association studies. However, its function during HIV-1 replication remains unclear.

METHODS AND RESULTS

We investigated the relationship between common RNF39 genetic variants and HIV-1 viral loads. The effect of RNF39 protein knockdown or overexpression on HIV-1 replication was then investigated in different cell lines. Two genetic variants were associated with HIV-1 viral loads. Patients with the ht1-GG/GG haplotype presented lower RNF39 expression levels and lower HIV-1 viral load. RNF39 knockdown inhibited HIV-1 expression.

CONCLUSIONS

RNF39 protein may be involved in HIV-1 replication as observed in genetic studies on patients with HIV-1 and in in vitro cell cultures.

Coronary artery aneurysms occurrence risk analysis between Kawasaki disease and LRP1B gene in Taiwanese children

Background: Kawasaki disease (KD) is an acute and systemic vasculitis. Its complications in coronary artery aneurysms (CAA) make KD one of the leading causes of acquired cardiovascular diseases in childhood. Low density lipoprotein receptor-related protein 1B (LRP1B) is abundantly expressed in the medial layer of coronary arteries and involved in endothelium inflammations.

Purpose: We aimed to identify the role of LRP1B in CAA formation during KD progression.

Methods: we investigated genetic variations in LRP1B in a Taiwanese cohort of 258 KD patients (83 with CAA and 175 without CAA complications). We used univariate and multivariate regression analyses to identify the associations between LRP1B genetic variations and KD patients.

Results: CAA formation in KD was significantly associated with the LRP1B (rs6707826) genetic variant (p = 0.007). By using multivariate regression analysis, significant correlations were observed between KD with CAA complications and the presence of the TT+TG genotypes for the LRP1B rs6707826 single-nucleotide polymorphism (full model: odds ratio = 2.82; 95% CI = 1.33–5.78).

Conclusion: Our results suggest that genetic polymorphism of LRP1B gene may be used as a genetic marker for the diagnosis and prognosis of the CAA formation in KD and contribute to genetic profiling studies for personalized medicine.

A 660-Kb deletion with antagonistic effects on fertility and milk production segregates at high frequency in Nordic Red cattle: additional evidence for the common occurrence of balancing selection in livestock

In dairy cattle, the widespread use of artificial insemination has resulted in increased selection intensity, which has led to spectacular increase in productivity. However, cow fertility has concomitantly severely declined. It is generally assumed that this reduction is primarily due to the negative energy balance of high-producing cows at the peak of lactation. We herein describe the fine-mapping of a major fertility QTL in Nordic Red cattle, and identify a 660-kb deletion encompassing four genes as the causative variant. We show that the deletion is a recessive embryonically lethal mutation. This probably results from the loss of RNASEH2B, which is known to cause embryonic death in mice. Despite its dramatic effect on fertility, 13%, 23% and 32% of the animals carry the deletion in Danish, Swedish and Finnish Red Cattle, respectively. To explain this, we searched for favorable effects on other traits and found that the deletion has strong positive effects on milk yield. This study demonstrates that embryonic lethal mutations account for a non-negligible fraction of the decline in fertility of domestic cattle, and that associated positive effects on milk yield may account for part of the negative genetic correlation. Our study adds to the evidence that structural variants contribute to animal phenotypic variation, and that balancing selection might be more common in livestock species than previously appreciated.

We report the identification of a large deletion encompassing four genes and the demonstration of its negative effect on fertility in Nordic Red dairy cattle. We show that this deletion is recessively lethal (homozygous embryos die) and therefore, when carrier cows are mated to carrier bulls, there is a high risk of embryonic mortality. As a result, chances of insemination failure are higher for such matings. Surprisingly, despite its negative effect, the deletion is frequent in Nordic Red cattle. We show that this high frequency may be a consequence of the fact that the deletion is associated with increased milk production and therefore selected for. Due to increased levels of inbreeding resulting from the widespread use of artificial insemination, such recessive lethal alleles may account for a non-negligible fraction of the reduction in fertility observed in cattle.

Association between GRIN3A gene polymorphism in Kawasaki disease and coronary artery aneurysms in Taiwanese children

Kawasaki disease (KD) is pediatric systemic vasculitis with the classic complication of coronary artery aneurysm (CAA). It is the leading cause of acquired cardiovascular diseases in children. Some severe cases present with multi-organ involvement or neurological dysfunction. To identify the role of the glutamate receptor, ionotropic, N-methyl-d-aspartate 3A (GRIN3A) in KD, we investigated genetic variations in GRIN3A in a Taiwanese cohort of 262 KD patients (76 with and 186 without CAA complications). We used univariate and multivariate regression analyses to identify the associations between clinical characteristics and GRIN3A genetic variations in KD. According to univariate regression analysis, CAA formation in KD was significantly associated with fever duration (p < 0.0001), first Intravenous immunoglobulin (IVIG) used (days after day one of fever) (p < 0.0001), and the GRIN3A (rs7849782) genetic variant (p < 0.001). KD patients with GG+GC genotype showed a lower rate of developing CAA (GG+GC genotype: odds ratio = 0.26; 95% CI = 0.14–0.46). Significant associations were identified between KD with CAA complication and the GRIN3A (rs7849782) genetic variant by using multivariate regression analysis. Specifically, significant correlations were observed between KD with CAA complications and the presence of GG+GC genotypes for the GRIN3A rs7849782 single-nucleotide polymorphism (full model: odds ratio = 0.25; 95% CI = 0.14–0.46). Our results suggest that a polymorphism of the GRIN3A gene may play a role in KD pathogenesis.

Sorting nexin 24 genetic variation associates with coronary artery aneurysm severity in Kawasaki disease patients

Background

The sorting nexin (SNX) family is involved in endocytosis and protein trafficking and plays multiple roles in various diseases. The role of SNX proteins in Kawasaki disease (KD) is not known. We attempted to test whether genetic SNX variation associates with the risk of coronary artery aneurysm (CAA) formation in KD.

Methods and results

Chi-square tests were used to identify SNX24 genetic variants associated with KD susceptibility and CAA formation in KD; models were adjusted for fever duration and time of first administration of intravenous immunoglobulin. We obtained clinical characteristics and genotypes from KD patients (76 with CAA and 186 without CAA) in a population-based retrospective KD cohort study (n = 262). Clinical and genetic factors were associated with CAA formation in KD. In addition, endothelial cell inflammation was evaluated. Significant correlation was observed between KD with CAA complications and the rs28891 single-nucleotide polymorphism in SNX24. Patients with CC + CT genotypes had lesser CAA complications. In lipopolysaccharide-treated human umbilical vein endothelial cells, siRNA knockdown of SNX24 significantly decreased gene expression of the proinflammatory cytokines IL-1 beta, IL-6, and IL-8.

Conclusions

Polymorphisms in SNX24 may be used as genetic markers for the diagnosis and prognosis of CAA formation in KD.

Variants in ZNRD1 gene predict HIV-1/AIDS disease progression in a Han Chinese population in Taiwan

Patients demonstrate notable variations in disease progression following human immunodeficiency virus (HIV) infection. We aimed to identify ZNRD1 and RNF39 genetic variants linked to AIDS progression. We conducted a genetic association study in HIV-1-infected Han Chinese patients residing in Taiwan. The clinical characteristics of 143 HIV-1-infected patients were measured, and patients were split into 2 groups: AIDS progression and AIDS non-progression. Genotyping of ZNRD1 and RNF39 was performed in all participants. We found that patients in the AIDS progression group had higher HIV-1 viral loads and lower CD4 cell counts than did patients in the AIDS non-progression group. The frequency of the AA genotype of ZNRD1 (rs16896970) was lower in the AIDS progression group than in the AIDS non-progression group. Patients with AA genotypes had lower levels of HIV-1 viral loads and higher levels of CD4 cell counts than did patients with AG+GG genotypes. AIDS progression in patients with the AA group is significantly different from that in patients with the AG and GG groups by using Kaplan-Meier survival analysis. The hazard ratio for progression was lower in the AA group than in the AG and GG groups. We identified a SNP that contributes to AIDS progression in HIV-1-infected patients in this population. This SNP had a significant protective influence on AIDS progression, and polymorphisms of the ZNRD1 gene may play a role in the pathogenesis of HIV-1 infection.

BTNL2 gene polymorphisms may be associated with susceptibility to Kawasaki disease and formation of coronary artery lesions in Taiwanese children

The butyrophilin-like 2 (BTNL2) gene is a member of the B7 receptor family that probably functions as a T cell costimulatory molecule. Because altered T cell functions are implicated in dysregulation of the immune response seen in Kawasaki disease (KD), it is reasonable to speculate that BTNL2 gene is involved in the pathophysiology of KD. The purpose of this study was to investigate whether polymorphisms of the BTNL2 gene are associated with KD and the development of coronary artery lesions (CALs) in Taiwanese children. Nine-three patients with KD and 669 ethnically matched healthy controls were genotyped for BTNL2 gene rs1555115 C/G and rs2395158 A/G polymorphisms. The frequency of GG genotype of rs 1555115 was significantly higher in KD patients compared with controls (2.2% vs 0.2%, P = 0.012). The odds ratio for developing KD in individuals with rs 1555115 GG genotype was 14.7 (95% confidence interval, 2.04-105.5, P = 0.003) compared with individuals with rs 1555115 CG and CC genotypes. No significant difference was observed in the genotype and allelic frequencies of rs 2395158 polymorphism between KD patients and controls. However, the frequency of the G allele of rs 2395158 was significantly higher in KD patients with CALs than in those without CALs (P = 0.001). No significant difference was observed in the genotype and allelic frequencies of rs 1555115 polymorphism between KD patients with and without CALs. In conclusion, our results suggest that BTNL2 gene polymorphisms might be genetic markers of KD susceptibility and risk of coronary artery complication in Taiwanese children.

dbSNP in the detail and copy number complexities

dbSNP is a general catalog of genetic polymorphism maintained by NCBI, mainly collating information for single nucleotide variations, many of which will be single nucleotide polymorphisms (SNPs), but also including small indels. It takes submissions from many sources, now also including large numbers of sequence variants identified by next-generation sequencing. A number of differently designed studies have attempted to estimate the error rates in data archived in dbSNP. Most recently, a study added to earlier studies identifying specific issues for duplicons and copy number variations (CNVs); earlier analyses have focused on stop codons, splice sites, and the general content of dbSNP. This article overviews dbSNP itself, these studies, and their implications.

HLA-E gene polymorphism associated with susceptibility to Kawasaki disease and formation of coronary artery aneurysms

OBJECTIVE

Kawasaki disease (KD) is a pediatric systemic vasculitis of unknown cause for which a genetic influence is supposed. The purpose of this study was to identify possible genetic variants in the major histocompatibility complex (MHC) region that are associated with KD and the development of coronary artery aneurysms (CAAs) in a Taiwanese population.

METHODS

The 168 genetic variants covering the MHC locus were analyzed in an association study of a Taiwanese cohort of 93 KD patients and 680 unrelated healthy children matched for sex and age with the study patients.

RESULTS

Eleven single-nucleotide polymorphisms (SNPs) were associated with the occurrence of KD. The SNP located at the 3'-untranslated region of HLA-E (rs2844724) was highly associated (P < 1 x 10(-7)). In addition, the frequency of the C allele was higher in KD patients without CAAs than in controls (P < 0.001) due to a significantly increased frequency of the CC and CT genotypes. Plasma levels of soluble HLA-E were significantly higher in KD patients than in controls regardless of the presence of CAAs. Furthermore, there was a trend toward higher plasma levels of soluble HLA-E in KD patients with the CT and TT genotypes of the HLA-E gene polymorphism.

CONCLUSION

Our results suggest that the HLA-E gene polymorphism may play a role in the pathogenesis of KD.

Relative effects of mutability and selection on single nucleotide polymorphisms in transcribed regions of the human genome

MOTIVATION

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation in humans. However, the factors that affect SNP density are poorly understood. The goal of this study was to estimate the relative effects of mutability and selection on SNP density in transcribed regions of human genes. It is important for prediction of the regions that harbor functional polymorphisms.

RESULTS

We used frequency-validated SNPs resulting from single-nucleotide substitutions. SNPs were subdivided into five functional categories: (i) 5' untranslated region (UTR) SNPs, (ii) 3' UTR SNPs, (iii) synonymous SNPs, (iv) SNPs producing conservative missense mutations, and (v) SNPs producing radical missense mutations. Each of these categories was further subdivided into nine mutational categories on the basis of the single-nucleotide substitution type. Thus, 45 functional/mutational categories were analyzed. The relative mutation rate in each mutational category was estimated on the basis of published data. The proportion of segregating sites (PSSs) for each functional/mutational category was estimated by dividing the observed number of SNPs by the number of potential sites in the genome for a given functional/mutational category. By analyzing each functional group separately, we found significant positive correlations between PSSs and relative mutation rates (Spearman's correlation coefficient, at least r = 0.96, df = 9, P < 0.001). We adjusted the PSSs for the mutation rate and found that the functional category had a significant effect on SNP density (F = 5.9, df = 4, P = 0.001), suggesting that selection affects SNP density in transcribed regions of the genome. We used analyses of variance and covariance to estimate the relative effects of selection (functional category) and mutability (relative mutation rate) on the PSSs and found that approximately 87% of variation in PSS was due to variation in the mutation rate and approximately 13% was due to selection, suggesting that the probability that a site located in a transcribed region of a gene is polymorphic mostly depends on the mutability of the site.

A simple and accurate SNP scoring strategy based on typeIIS restriction endonuclease cleavage and matrix-assisted laser desorption/ionization mass spectrometry

BACKGROUND

We describe the development of a novel matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF)-based single nucleotide polymorphism (SNP) scoring strategy, termed Restriction Fragment Mass Polymorphism (RFMP) that is suitable for genotyping variations in a simple, accurate, and high-throughput manner. The assay is based on polymerase chain reaction (PCR) amplification and mass measurement of oligonucleotides containing a polymorphic base, to which a typeIIS restriction endonuclease recognition was introduced by PCR amplification. Enzymatic cleavage of the products leads to excision of oligonucleotide fragments representing base variation of the polymorphic site whose masses were determined by MALDI-TOF MS.

RESULTS

The assay represents an improvement over previous methods because it relies on the direct mass determination of PCR products rather than on an indirect analysis, where a base-extended or fluorescent report tag is interpreted. The RFMP strategy is simple and straightforward, requiring one restriction digestion reaction following target amplification in a single vessel. With this technology, genotypes are generated with a high call rate (99.6%) and high accuracy (99.8%) as determined by independent sequencing.

CONCLUSION

The simplicity, accuracy and amenability to high-throughput screening analysis should make the RFMP assay suitable for large-scale genotype association study as well as clinical genotyping in laboratories.

Shifting paradigm of association studies: value of rare single-nucleotide polymorphisms

Currently, single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) of >5% are preferentially used in case-control association studies of common human diseases. Recent technological developments enable inexpensive and accurate genotyping of a large number of SNPs in thousands of cases and controls, which can provide adequate statistical power to analyze SNPs with MAF <5%. Our purpose was to determine whether evaluating rare SNPs in case-control association studies could help identify causal SNPs for common diseases. We suggest that slightly deleterious SNPs (sdSNPs) subjected to weak purifying selection are major players in genetic control of susceptibility to common diseases. We compared the distribution of MAFs of synonymous SNPs with that of nonsynonymous SNPs (1) predicted to be benign, (2) predicted to be possibly damaging, and (3) predicted to be probably damaging by PolyPhen. Our sources of data were the International HapMap Project, ENCODE, and the SeattleSNPs project. We found that the MAF distribution of possibly and probably damaging SNPs was shifted toward rare SNPs compared with the MAF distribution of benign and synonymous SNPs that are not likely to be functional. We also found an inverse relationship between MAF and the proportion of nsSNPs predicted to be protein disturbing. On the basis of this relationship, we estimated the joint probability that a SNP is functional and would be detected as significant in a case-control study. Our analysis suggests that including rare SNPs in genotyping platforms will advance identification of causal SNPs in case-control association studies, particularly as sample sizes increase.

A Hidden Markov Model to estimate population mixture and allelic copy-numbers in cancers using Affymetrix SNP arrays

Background

Affymetrix SNP arrays can interrogate thousands of SNPs at the same time. This allows us to look at the genomic content of cancer cells and to investigate the underlying events leading to cancer. Genomic copy-numbers are today routinely derived from SNP array data, but the proposed algorithms for this task most often disregard the genotype information available from germline cells in paired germline-tumour samples. Including this information may deepen our understanding of the "true" biological situation e.g. by enabling analysis of allele specific copy-numbers. Here we rely on matched germline-tumour samples and have developed a Hidden Markov Model (HMM) to estimate allelic copy-number changes in tumour cells. Further with this approach we are able to estimate the proportion of normal cells in the tumour (mixture proportion).

Results

We show that our method is able to recover the underlying copy-number changes in simulated data sets with high accuracy (above 97.71%). Moreover, although the known copy-numbers could be well recovered in simulated cancer samples with more than 70% cancer cells (and less than 30% normal cells), we demonstrate that including the mixture proportion in the HMM increases the accuracy of the method. Finally, the method is tested on HapMap samples and on bladder and prostate cancer samples.

Conclusion

The HMM method developed here uses the genotype calls of germline DNA and the allelic SNP intensities from the tumour DNA to estimate allelic copy-numbers (including changes) in the tumour. It differentiates between different events like uniparental disomy and allelic imbalances. Moreover, the HMM can estimate the mixture proportion, and thus inform about the purity of the tumour sample.

Genome resources and comparative analysis tools for cardiovascular research

Disorders of the cardiovascular system are often caused by the interaction of genetic and environmental factors that jointly contribute to individual susceptibility. Genomic data and bioinformatics tools generated from genome projects, coupled with functional verification, offer novel approaches to study both rare single-gene and complex multigenic cardiovascular diseases. These approaches include gene mapping using genome variation, especially single-nucleotide polymorphisms and comparative genomics within and between species. This chapter illustrates the major genome resources, associated bioinformatics tools, and their potential application in cardiovascular research.

Genetic polymorphisms in DPF3 associated with risk of breast cancer and lymph node metastases

BACKGROUND

Several studies have identified rare genetic variations responsible for many cases of familial breast cancer but their contribution to total breast cancer incidence is relatively small. More common genetic variations with low penetrance have been postulated to account for a higher proportion of the population risk of breast cancer.

METHODS AND RESULTS

In an effort to identify genes that influence non-familial breast cancer risk, we tested over 25,000 single nucleotide polymorphisms (SNPs) located within approximately 14,000 genes in a large-scale case-control study in 254 German women with breast cancer and 268 age-matched women without malignant disease. We identified a marker on chromosome 14q24.3-q31.1 that was marginally associated with breast cancer status (OR = 1.5, P = 0.07). Genotypes for this SNP were also significantly associated with indicators of breast cancer severity, including presence of lymph node metastases (P = 0.006) and earlier age of onset (P = 0.01). The association with breast cancer status was replicated in two independent samples (OR = 1.35, P = 0.05). High-density association fine mapping showed that the association spanned about 80 kb of the zinc-finger gene DPF3 (also known as CERD4). One SNP in intron 1 was found to be more strongly associated with breast cancer status in all three sample collections (OR = 1.6, P = 0.003) as well as with increased lymph node metastases (P = 0.01) and tumor size (P = 0.01).

CONCLUSION

Polymorphisms in the 5' region of DPF3 were associated with increased risk of breast cancer development, lymph node metastases, age of onset, and tumor size in women of European ancestry. This large-scale association study suggests that genetic variation in DPF3 contributes to breast cancer susceptibility and severity.

Different applications of polymerases with and without proofreading activity in single-nucleotide polymorphism analysis

With the completion of the human genome project, single-nucleotide polymorphisms (SNPs) have become the focus of intense study in biomedical research. Polymerase-mediated primer extension has been employed in a variety of SNP assays. However, these SNP assays using polymerase without proofreading function are compromised by their low reliability. Using a newly developed short amplicon harboring restriction enzyme site, EcoR-I, we were able to compare the single-base discrimination abilities of polymerases with and without proofreading function in primer extension in a broad range of annealing temperatures. Thermodynamic analysis demonstrated a striking single-nucleotide discrimination ability of polymerases with proofreading function. Using unmodified 3'-end allele-specific primers, only template-dependent products were generated by polymerase with proofreading activity. This powerful single-base discrimination ability of exo(+) polymerases was further evaluated in primer extension using three types of 3' terminally modified allele-specific primers. As compared with the poor fidelity in primer extension of polymerases lacking 3' exonuclease activity, this study provides convincing evidence that the use of proofreading polymerases in combination with 3'-end modified allele-specific primers can be a powerful new strategy for the development of SNP assays.

Polymorphisms in human soluble epoxide hydrolase

Human soluble epoxide hydrolase (hsEH) metabolizes a variety of epoxides to the corresponding vicinal diols. Arachidonic and linoleic acid epoxides are thought to be endogenous substrates for hsEH. Enzyme activity in humans shows high interindividual variation (e.g., 500-fold in liver) suggesting the existence of regulatory and/or structural gene polymorphisms. We resequenced each of the 19 exons of the hsEH gene (EPHX2) from 72 persons representing black, Asian, and white populations. A variety of polymorphisms was found, six of which result in amino acid substitutions. Amino acid variants were localized on the crystal structure of the mouse sEH, resulting in the prediction that at least two of these (Arg287Gln and Arg103Cys) might significantly affect enzyme function. The six variants of the hsEH cDNA corresponding to each single polymorphism and one corresponding to a double polymorphism were then constructed by site-directed mutagenesis and expressed in insect cells. As predicted, Arg287Gln and the double mutant Arg287Gln/Arg103Cys showed decreased enzyme activity using trans-stilbene oxide, trans-diphenylpropene oxide, and 14,15-epoxyeicosatrienoic acid as substrates. Lys55Arg and Cys154Tyr mutants had elevated activity for all three substrates. Detailed kinetic studies revealed that the double mutant Arg287Gln/Arg103Cys showed significant differences in Km and Vmax. In addition, stability studies showed that the double mutant was less stable than wild-type protein when incubated at 37 degrees C. These results suggest that at least six hsEH variants exist in the human population and that at least four of these may influence hsEH-mediated metabolism of exogenous and endogenous epoxide substrates in vivo.

Human diallelic insertion/deletion polymorphisms

We report the identification and characterization of 2,000 human diallelic insertion/deletion polymorphisms (indels) distributed throughout the human genome. Candidate indels were identified by comparison of overlapping genomic or cDNA sequences. Average confirmation rate for indels with a > or =2-nt allele-length difference was 58%, but the confirmation rate for indels with a 1-nt length difference was only 14%. The vast majority of the human diallelic indels were monomorphic in chimpanzees and gorillas. The ratio of deletionrcolon;insertion mutations was 4.1. Allele frequencies for the indels were measured in Europeans, Africans, Japanese, and Native Americans. New alleles were generally lower in frequency than old alleles. This tendency was most pronounced for the Africans, who are likely to be closest among the four groups to the original modern human population. Diallelic indels comprise approximately 8% of all human polymorphisms. Their abundance and ease of analysis make them useful for many applications.

The matrix coalescent and an application to human single-nucleotide polymorphisms

The "matrix coalescent" is a reformulation of the familiar coalescent process of population genetics. It ignores the topology of the gene tree and treats the coalescent as a Markov process describing the decay in the number of ancestors of a sample of genes as one proceeds backward in time. The matrix formulation of this process is convenient when the population changes in size, because such changes affect only the eigenvalues of the transition matrix, not the eigenvectors. The model is used here to calculate the expectation of the site frequency spectrum under various assumptions about population history. To illustrate how this method can be used with data, we then use it in conjunction with a set of SNPs to test hypotheses about the history of human population size.

Candidate-gene approaches for studying complex genetic traits: practical considerations

Association studies with candidate genes have been widely used for the study of complex diseases. However, this approach has been criticized because of non-replication of results and limits on its ability to include all possible causative genes and polymorphisms. These challenges have led to pessimism about the candidate-gene approach and about the genetic analysis of complex diseases in general. We believe that these criticisms can be usefully countered with an appeal to the principles of epidemiological investigation.

In silico approaches to mechanistic and predictive toxicology: an introduction to bioinformatics for toxicologists

Bioinformatics, or in silico biology, is a rapidly growing field that encompasses the theory and application of computational approaches to model, predict, and explain biological function at the molecular level. This information rich field requires new skills and new understanding of genome-scale studies in order to take advantage of the rapidly increasing amount of sequence, expression, and structure information in public and private databases. Toxicologists are poised to take advantage of the large public databases in an effort to decipher the molecular basis of toxicity. With the advent of high-throughput sequencing and computational methodologies, expressed sequences can be rapidly detected and quantitated in target tissues by database searching. Novel genes can also be isolated in silico, while their function can be predicted and characterized by virtue of sequence homology to other known proteins. Genomic DNA sequence data can be exploited to predict target genes and their modes of regulation, as well as identify susceptible genotypes based on single nucleotide polymorphism data. In addition, highly parallel gene expression profiling technologies will allow toxicologists to mine large databases of gene expression data to discover molecular biomarkers and other diagnostic and prognostic genes or expression profiles. This review serves to introduce to toxicologists the concepts of in silico biology most relevant to mechanistic and predictive toxicology, while highlighting the applicability of in silico methods using select examples.

Differential occurrence of mutations causative of eye diseases in the Chinese population

Ethnic differences and geographic variations affect the frequencies and nature of human mutations. In the literature, descriptions of causative mutations of eye diseases in the Chinese population are few. In this paper we attempt to reveal molecular information on genetic eye diseases involving Chinese patients from published and unpublished works by us and other groups. Our studies on candidate genes of eye diseases in the Chinese population in Hong Kong include MYOC and TISR for primary open angle glaucoma, RHO and RP1 for retinitis pigmentosa, ABCA4 and APOE for age-related macular degeneration, RB1 for retinoblastoma, APC for familial adenomatous polyposis with congenital hypertrophy of retinal pigment epithelium, BIGH3/TGFBI for corneal dystrophies, PAX6 for aniridia and Reiger syndrome, CRYAA and CRYBB2 for cataracts, and mtDNA for Leber hereditary optic neuropathy. We have revealed novel mutations in most of these genes, and in RHO, RP1, RB1, BIGH3, and PAX6 we have reported mutations that contribute to better understanding of the functions and properties of the respective gene products. We showed absence of MYOC does not necessarily cause glaucoma. No disease causative mutations have been identified in MYOC or ABCA4. There are similarities in the patterns of sequence alterations and phenotype-genotype associations in comparison with other ethnic groups, while the MYOC, RB1, APC, and PAX6 genes have more Chinese-specific sequence alterations. Establishment of a mutation database specific for the Chinese is essential for identification of genetic markers with diagnostic, prognostic, or pharmacological values.

Reading between the LINEs: human genomic variation induced by LINE-1 retrotransposition

The insertion of mobile elements into the genome represents a new class of genetic markers for the study of human evolution. Long interspersed elements (LINEs) have amplified to a copy number of about 100,000 over the last 100 million years of mammalian evolution and comprise approximately 15% of the human genome. The majority of LINE-1 (L1) elements within the human genome are 5' truncated copies of a few active L1 elements that are capable of retrotransposition. Some of the young L1 elements have inserted into the human genome so recently that populations are polymorphic for the presence of an L1 element at a particular chromosomal location. L1 insertion polymorphisms offer several advantages over other types of polymorphisms for human evolution studies. First, they are typed by rapid, simple, polymerase chain reaction (PCR)-based assays. Second, they are stable polymorphisms that rarely undergo deletion. Third, the presence of an L1 element represents identity by descent, because the probability is negligible that two different young L1 repeats would integrate independently between the exact same two nucleotides. Fourth, the ancestral state of L1 insertion polymorphisms is known to be the absence of the L1 element, which can be used to root plots/trees of population relationships. Here we report the development of a PCR-based display for the direct identification of dimorphic L1 elements from the human genome. We have also developed PCR-based assays for the characterization of six polymorphic L1 elements within the human genome. PCR analysis of human/rodent hybrid cell line DNA samples showed that the polymorphic L1 elements were located on several different chromosomes. Phylogenetic analysis of nonhuman primate DNA samples showed that all of the recently integrated "young" L1 elements were restricted to the human genome and absent from the genomes of nonhuman primates. Analysis of a diverse array of human populations showed that the allele frequencies and level of heterozygosity for each of the L1 elements was variable. Polymorphic L1 elements represent a new source of identical-by-descent variation for the study of human evolution. [The sequence data described in this paper have been submitted to the GenBank data library under accession nos. AF242435-AF242451.]

Central mutation databases--a review

The Internet has been a key component in the coordination of the diverse group of scientists involved in the Human Genome Project. Nowhere has this contribution been more critical than in the maintenance and exchange of information about genetic variation and mutation. Whereas the majority of DNA sequence is generated and stored by a relatively few sites, a far greater number of researchers investigate the variations in that sequence from sites scattered worldwide. It falls to central databases to utilize the Internet to assemble data from these sites and make them available to the greater human genomic community.

Sequence variation database project at the European Bioinformatics Institute

The sequence variation project at EBI aims to create a unified resource for browsing and searching sequence differences. Technical advances in reading in new data types and in validating and cross-referencing entries are reported. It is suggested that the hardest problems in unifying mutation databases are related to intellectual property rights. The concept of copylefting is introduced as a potential solution to these.

MuStaR and other software for locus-specific mutation databases

As the human genome sequencing project nears completion, there has been a vast increase in the rate at which disease and nondisease associated variant sequences are being sought and detected. This has heightened the need for software with which to accumulate allelic variant (mutation) data, and with which to make the data accessible to the scientific community. Many ad hoc solutions have been developed by those interested in specific genes and diseases, and the creation of central databases which hold data for all genes has provided an alternative repository for some of the locus data. Despite this, few specialised software tools exist for researchers to create their own locus-specific allelic variant databases. This article describes methods available to potential curators, including software systems developed with the sole purpose of generating locus-specific mutation databases. In particular, the authors' own software, MuStaRtrade mark, is described. MuStaRtrade mark allows curators to maintain a database on a laptop computer if desired, while being able to export the data to an automatically generated Website which will run on any cgi compliant Web server. Searching the database and the submission of new mutations are made possible through fill-in Web forms. A number of other software tools which may be of use to curators are also described.