Genome-wide gene-gene interaction analysis for cardiovascular disease

Numerous studies have been carried out to try to better understand the genetic predisposition for cardiovascular disease. Although it is widely believed that multifactorial diseases such as cardiovascular disease is the result from effects of many genes which working alone or interact with other genes, most genetic studies have been focused on identifying of cardiovascular disease susceptibility genes and usually ignore the effects of gene-gene interactions in the analysis. The current study applies a novel linkage disequilibrium based statistic for testing interactions between two linked loci using data from a genome-wide study of cardiovascular disease. A total of 53,394 single nucleotide polymorphisms (SNPs) are tested for pair-wise interactions, and 8,644 interactions are found to be significant with p-values less than 3.5×10-11. Results indicate that known cardiovascular disease susceptibility genes tend not to have many significantly interactions. One SNP in the CACNG1 (calcium channel, voltage-dependent, gamma subunit 1) gene and one SNP in the IL3RA (interleukin 3 receptor, alpha) gene are found to have the most significant pair-wise interactions. Findings from the current study should be replicated in other independent cohort to eliminate potential false positive results.^

Association of folate metabolic pathway genes with human spina bifida meningomyelocele

Neural tube defects including spina bifida meningomyelocele (SBMM) are common malformations of the brain and spinal cord, and include all abnormalities resulting from lack of closure of the developing neural tube during embryological development.^ The specific aims of this study were to determine if single nucleotide polymorphic variants (SNPs) in the folate/homocysteine metabolic pathway genes confer a risk for NTD susceptibility within this SBMM population.^ In completion of the first specific aim, two novel SNPs were identified in the FOLR1 gene in Chromosome 11of patients including one in non-coding exon 1 with a C → T transition at nucleotide position 71578317 and another in non-coding exon 3 with a T → G transversion at nucleotide position 71579123. It will be important to determine if these variants are present in the respective parents of these individuals. If they are in fact de novo variants, then these SNPs may be more likely to contribute to the birth defect.^ The second project aim was to analyze genotypes associated with SBMM risk by transmission disequilibrium tests (TDT) and association was detected on several SNPs across the folate metabolic pathway genes in this population. SNPs with significant RC-TDT values were found within the DHFR gene (rs1650723), the MTRR gene (rs327592), the FOLR2 gene (rs13908), four tightly linked variants in the FOLR3 gene (rs7925545, rs7926875, rs7926987, rs7926360) and a variant in the SLC19A1 gene (rs1888530). The product of each of these genes performs a vital function in the folate metabolic pathway. It is conceivable, therefore, that if the individual SNP or SNPs can be proven to perturb the function in some way that they may be involved in the disruption of folate metabolism and in the resulting birth defect. Validating the results of this study in other independent populations will further strengthen the evidence that dysfunction of folate enzymes and receptors may confer SBMM risk in humans. ^

Genome-wide algorithm for detecting CNV associations with diseases

SNP genotyping arrays have been developed to characterize single-nucleotide polymorphisms (SNPs) and DNA copy number variations (CNVs). The quality of the inferences about copy number can be affected by many factors including batch effects, DNA sample preparation, signal processing, and analytical approach. Nonparametric and model-based statistical algorithms have been developed to detect CNVs from SNP genotyping data. However, these algorithms lack specificity to detect small CNVs due to the high false positive rate when calling CNVs based on the intensity values. Association tests based on detected CNVs therefore lack power even if the CNVs affecting disease risk are common. In this research, by combining an existing Hidden Markov Model (HMM) and the logistic regression model, a new genome-wide logistic regression algorithm was developed to detect CNV associations with diseases. We showed that the new algorithm is more sensitive and can be more powerful in detecting CNV associations with diseases than an existing popular algorithm, especially when the CNV association signal is weak and a limited number of SNPs are located in the CNV.^

Genome-wide association study of LDL-cholesterol over ApoB ratio

Genome-Wide Association Study analytical (GWAS) methods were applied in a large biracial sample of individuals to investigate variation across the genome for its association with a surrogate low-density lipoprotein (LDL) particle size phenotype, the ratio of LDL-cholesterol level over ApoB level. Genotyping was performed on the Affymetrix 6.0 GeneChip with approximately one million single nucleotide polymorphisms (SNPs). The ratio of LDL cholesterol to ApoB was calculated, and association tests used multivariable linear regression analysis with an additive genetic model after adjustment for the covariates sex, age and BMI. Association tests were performed separately in African Americans and Caucasians. There were 9,562 qualified individuals in the Caucasian group and 3,015 qualified individuals in the African American group. Overall, in Caucasians two statistically significant loci were identified as being associated with the ratio of LDL-cholesterol over ApoB: rs10488699 (p<5 x10-8, 11q23.3 near BUD13) and the SNP rs964184 (p<5 x10-8 11q23.3 near ZNF259). We also found rs12286037 ((p<4x10-7) (11q23.3) near APOA5/A4/C3/A1 with suggestive associate in the Caucasian sample. In exploratory analyses, a difference in the pattern of association between individuals taking and not taking LDL-cholesterol lowering medications was observed. Individuals who were not taking medications had smaller p-value than those taking medication. In the African-American group, there were no significant (p<5x10-8) or suggestive associations (p<4x10-7) with the ratio of LDL-cholesterol over ApoB after adjusting for age, BMI, and sex and comparing individuals with and without LDL-cholesterol lowering medication. Conclusions: There were significant and suggestive associations between SNP genotype and the ratio of LDL-cholesterol to ApoB in Caucasians, but these associations may be modified by medication treatment.^

Genetic variants in the mammalian target of rapamycin (mTOR) signaling pathway as predictors of survival and clinical response in women with ovarian cancer

Background. The mTOR pathway is commonly altered in human tumors and promotes cell survival and proliferation. Preliminary evidence suggests this pathway's involvement in chemoresistance to platinum and taxanes, first line therapy for epithelial ovarian cancer. A pathway-based approach was used to identify individual germline single nucleotide polymorphisms (SNPs) and cumulative effects of multiple genetic variants in mTOR pathway genes and their association with clinical outcome in women with ovarian cancer. ^ Methods. The case-series was restricted to 319 non-Hispanic white women with high grade ovarian cancer treated with surgery and platinum-based chemotherapy. 135 SNPs in 20 representative genes in the mTOR pathway were genotyped. Hazard ratios (HRs) for death and Odds ratios (ORs) for failure to respond to primary therapy were estimated for each SNP using the multivariate Cox proportional hazards model and multivariate logistic regression model, respectively, while adjusting for age, stage, histology and treatment sequence. A survival tree analysis of SNPs with a statistically significant association (p<0.05) was performed to identify higher order gene-gene interactions and their association with overall survival. ^ Results. There was no statistically significant difference in survival by tumor histology or treatment regimen. The median survival for the cohort was 48.3 months. Seven SNPs were significantly associated with decreased survival. Compared to those with no unfavorable genotypes, the HR for death increased significantly with the increasing number of unfavorable genotypes and women in the highest risk category had HR of 4.06 (95% CI 2.29–7.21). The survival tree analysis also identified patients with different survival patterns based on their genetic profiles. 13 SNPs on five different genes were found to be significantly associated with a treatment response, defined as no evidence of disease after completion of primary therapy. Rare homozygous genotype of SNP rs6973428 showed a 5.5-fold increased risk compared to the wild type carrying genotypes. In the cumulative effect analysis, the highest risk group (individuals with ≥8 unfavorable genotypes) was significantly less likely to respond to chemotherapy (OR=8.40, 95% CI 3.10–22.75) compared to the low risk group (≤4 unfavorable genotypes). ^ Conclusions. A pathway-based approach can demonstrate cumulative effects of multiple genetic variants on clinical response to chemotherapy and survival. Therapy targeting the mTOR pathway may modify outcome in select patients.^

Molecular epidemiological studies on genetic predispositions to squamous cell carcinoma of the head and neck

Two molecular epidemiological studies were conducted to examine associations between genetic variation and risk of squamous cell carcinoma of the head and neck (SCCHN). In the first study, we hypothesized that genetic variation in p53 response elements (REs) may play roles in the etiology of SCCHN. We selected and genotyped five polymorphic p53 REs as well as a most frequently studied p53 codon 72 (Arg72Pro, rs1042522) polymorphism in 1,100 non-Hispanic White SCCHN patients and 1,122 age-and sex-matched cancer-free controls recruited at The University of Texas M. D. Anderson Cancer Center. In multivariate logistic regression analysis with adjustment for age, sex, smoking and drinking status, marital status and education level, we observed that the EOMES rs3806624 CC genotype had a significant effect of protection against SCCHN risk (adjusted odds ratio= 0.79, 95% confidence interval =0.64–0.98), compared with the -838TT+CT genotypes. Moreover, a significantly increased risk associated with the combined genotypes of p53 codon 72CC and EOMES -838TT+CT was observed, especially in the subgroup of non-oropharyneal cancer patients. The values of false-positive report probability were also calculated for significant findings. In the second study, we assessed the association between SCCHN risk and four potential regulatory single nucleotide polymorphisms (SNPs) of DEC1 (deleted in esophageal cancer 1) gene, a candidate tumor suppressor gene for esophageal cancer. After adjustment for age, sex, and smoking and drinking status, the variant -606CC (i.e., -249CC) homozygotes had a significantly reduced SCCHN risk (adjusted odds ratio = 0.71, 95% confidence interval = 0.52–0.99), compared with the -606TT homozygotes. Stratification analyses showed that a reduced risk associated with the -606CC genotype was more pronounced in subgroups of non-smokers, non-drinkers, younger subjects (defined as ≤ 57 years), carriers of TP53 Arg/Arg (rs1042522) genotype, patients with oropharyngeal cancer or late-stage SCCHN. Further in silico analysis revealed that the -249 T-to-C change led to a gain of a transcription factor binding site. Additional functional analysis showed that the -249T-to-C change significantly enhanced transcriptional activity of the DEC1 promoter and the DNA-protein binding activity. We conclude that the DEC1 promoter -249 T>C (rs2012775) polymorphism is functional, modulating susceptibility to SCCHN among non-Hispanic Whites. Additional large-scale, preferably population-based studies are needed to validate our findings.^

Genetic variants associated with metastatic prostate cancer

Introduction. Distant metastasis remains the leading cause of death among prostate cancer patients. Several genetic susceptibility loci associated with Prostate cancer have been identified by the Genome Wide Association Studies (GWAS). To date, few studies have explored the ability of these SNPs to identify metastatic prostate cancer. Based on the identification of genetic variants as predictors of aggressive disease, a case comparison study of prostate cancer patients was designed to explore the association of 96 GWAS single nucleotide polymorphisms (SNPs) with metastatic disease. ^ Method. 1242 histologically confirmed prostate cancer patients, with and without metastatic disease, were enrolled into the study. Data were collected from personal interviews, hospital database and abstraction of medical records. Ninety six SNPs identified from GWAS studies based on their associations with prostate cancer risk were genotyped in the study population. Univariate and multivariate logistic regression analyses were used to explore the relationships of the variants with metastatic prostate cancer in Whites and African American men. ^ Results. Four SNPs showed independent associations with metastatic prostate cancer (rs721048 in EHBP1 (2p15), rs3025039 in VEGF (6p12), rs11228565 in Intergenic(11q13.2) and rs2735839 in KLK3(19q13.33)) in the White population. For SNP rs2735839 in KLK3, genotype GA was 1.71 times as likely to be associated with metastatic prostate cancer diagnosis as genotype AA after adjusting for other significant SNPs and covariates (95% CI, 1.12-2.60; p=0.012). In men of African descent, three SNPs: rs1512268 in NKX3-1(8p21.2), rs12155172 in intergenic (7p15.3) & rs10486567 in JAZF1 (7p15.2) were positively associated with metastatic disease in the multivariate analysis. The strongest SNP was rs1512268 heterozygous genotype AG in NKX3-1(8p21.2) which was associated with 3.97-fold increased risk of metastatic prostate cancer diagnosis (95% CI, 1.69-9.34; p =0.002). ^ Conclusion. Genetic variants associated with metastatic prostate cancer were different in Whites and African American men. Given the high mortality rate recorded in men diagnosed with metastatic prostate tumor, further studies are needed to validate associations and establish their clinical application.^

Populations genomics analysis of legume host preference for specific rhizobial genotypes in the Rhizobium leguminosarum bv. viciae symbioses

Rhizobium leguminosarum bv. viciae establishes root nodule symbioses with several legume genera. Although most isolates are equally effective in establishing symbioses with all host genera, previous evidence suggests that hosts select specific rhizobial genotypes among those present in the soil. We have used population genomics to further investigate this observation. P. sativum, L. culinaris, V. sativa, and V. faba plants were used to trap rhizobia from a well-characterized soil, and pooled genomic DNAs from one-hundred isolates from each plant were sequenced. Sequence reads were aligned to the R. leguminosarum bv. viciae 3841 reference genome. High overall conservation of sequences was observed in all subpopulations, although several multigenic regions were absent from the soil population. A large fraction (16-22%) of sequence reads could not be recruited to the reference genome, suggesting that they represent sequences specific to that particular soil population. Although highly conserved, the 16S-23S rRNA gene region presented single nucleotide polymorphisms (SNPs) regarding the reference genome, but no striking differences could be found among plant-selected subpopulations. Plant-specific SNP patterns were, however, clearly observed within the nod gene cluster, supporting the existence of a plant preference for specific rhizobial genotypes. This was also shown after genome-wide analysis of SNP patterns.

Estimación del coste de almacenamiento de una solución de preservación de la privacidad en archivos BAM reales

A menudo los científicos secuencian el ADN de un gran número de personas con el objetivo de determinar qué genes se asocian con determinadas enfermedades. Esto permite meóon del genoma humano. El precio de un perfil genómico completo se ha posicionado por debajo de los 200 dólares y este servicio lo ofrecen muchas compañías, la mayor parte localizadas en EEUU. Como consecuencia, en unos pocos a~nos la mayoría de las personas procedentes de los países desarrollados tendrán los medios para tener su ADN secuenciado. Alrededor del 0.5% del ADN de cada persona (que corresponde a varios millones de nucleótidos) es diferente del genoma de referencia debido a variaciones genéticas. Así que el genoma contiene información altamente sensible y personal y representa la identidad biológica óon sobre el entorno o estilo de vida de uno (a menudo facilmente obtenible de las redes sociales), sería posible inferir el fenotipo del individuo. Multiples GWAS (Genome Wide Association Studies) realizados en los últimos a~nos muestran que la susceptibilidad de un paciente a tener una enfermedad en particular, como el Alzheimer, cáncer o esquizofrenia, puede ser predicha parcialmente a partir de conjuntos de sus SNP (Single Nucleotide Polimorphism). Estos resultados pueden ser usados para medicina genómica personalizada (facilitando los tratamientos preventivos y diagnósticos), tests de paternidad genéticos y tests de compatibilidad genética para averiguar a qué enfermedades pueden ser susceptibles los descendientes. Estos son algunos de los beneficios que podemos obtener usando la información genética, pero si esta información no es protegida puede ser usada para investigaciones criminales y por compañías aseguradoras. Este hecho podría llevar a discriminaci ón genética. Por lo que podemos concluir que la privacidad genómica es fundamental por el hecho de que contiene información sobre nuestra herencia étnica, nuestra predisposición a múltiples condiciones físicas y mentales, al igual que otras características fenotópicas, ancestros, hermanos y progenitores, pues los genomas de cualquier par de individuos relacionados son idénticos al 99.9%, contrastando con el 99.5% de dos personas aleatorias. La legislación actual no proporciona suficiente información técnica sobre como almacenar y procesar de forma segura los genomas digitalizados, por lo tanto, es necesaria una legislación mas restrictiva ---ABSTRACT---Scientists typically sequence DNA from large numbers of people in order to determine genes associated with particular diseases. This allows to improve the modern healthcare and to provide a better understanding of the human genome. The price of a complete genome profile has plummeted below $200 and this service is ofered by a number of companies, most of them located in the USA. Therefore, in a few years, most individuals in developed countries will have the means of having their genomes sequenced. Around 0.5% of each person's DNA (which corresponds to several millions of nucleotides) is diferent from the reference genome, owing to genetic variations. Thus, the genome contains highly personal and sensitive information, and it represents our ultimate biological identity. By combining genomic data with information about one's environment or lifestyle (often easily obtainable from social networks), could make it possible to infer the individual's phenotype. Multiple Genome Wide Association Studies (GWAS) performed in recent years have shown that a patient's susceptibility to particular diseases, such as Alzheimer's, cancer, or schizophrenia, can be partially predicted from sets of his SNPs. This results can be used for personalized genomic medicine (facilitating preventive treatment and diagnosis), genetic paternity tests, ancestry and genealogical testing, and genetic compatibility tests in order to have knowledge about which deseases would the descendant be susceptible to. These are some of the betefts we can obtain using genoma information, but if this information is not protected it can be used for criminal investigations and insurance purposes. Such issues could lead to genetic discrimination. So we can conclude that genomic privacy is fundamental due to the fact that genome contains information about our ethnic heritage, predisposition to numerous physical and mental health conditions, as well as other phenotypic traits, and ancestors, siblings, and progeny, since genomes of any two closely related individuals are 99.9% identical, in contrast with 99.5%, for two random people. The current legislation does not ofer suficient technical information about safe and secure ways of storing and processing digitized genomes, therefore, there is need for more restrictive legislation.

Direct genetic analysis by matrix-assisted laser desorption/ionization mass spectrometry

An approach to analyzing single-nucleotide polymorphisms (SNPs) found in the human genome has been developed that couples a recently developed invasive cleavage assay for nucleic acids with detection by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). The invasive cleavage assay is a signal amplification method that enables the analysis of SNPs by MALDI-TOF MS directly from human genomic DNA without the need for initial target amplification by PCR. The results presented here show the successful genotyping by this approach of twelve SNPs located randomly throughout the human genome. Conventional Sanger sequencing of these SNP positions confirmed the accuracy of the MALDI-TOF MS analysis results. The ability to unambiguously detect both homozygous and heterozygous genotypes is clearly demonstrated. The elimination of the need for target amplification by PCR, combined with the inherently rapid and accurate nature of detection by MALDI-TOF MS, gives this approach unique and significant advantages in the high-throughput genotyping of large numbers of SNPs, useful for locating, identifying, and characterizing the function of specific genes.

Complex promoter and coding region β2-adrenergic receptor haplotypes alter receptor expression and predict in vivo responsiveness

The human β2-adrenergic receptor gene has multiple single-nucleotide polymorphisms (SNPs), but the relevance of chromosomally phased SNPs (haplotypes) is not known. The phylogeny and the in vitro and in vivo consequences of variations in the 5′ upstream and ORF were delineated in a multiethnic reference population and an asthmatic cohort. Thirteen SNPs were found organized into 12 haplotypes out of the theoretically possible 8,192 combinations. Deep divergence in the distribution of some haplotypes was noted in Caucasian, African-American, Asian, and Hispanic-Latino ethnic groups with >20-fold differences among the frequencies of the four major haplotypes. The relevance of the five most common β2-adrenergic receptor haplotype pairs was determined in vivo by assessing the bronchodilator response to β agonist in asthmatics. Mean responses by haplotype pair varied by >2-fold, and response was significantly related to the haplotype pair (P = 0.007) but not to individual SNPs. Expression vectors representing two of the haplotypes differing at eight of the SNP loci and associated with divergent in vivo responsiveness to agonist were used to transfect HEK293 cells. β2-adrenergic receptor mRNA levels and receptor density in cells transfected with the haplotype associated with the greater physiologic response were ≈50% greater than those transfected with the lower response haplotype. The results indicate that the unique interactions of multiple SNPs within a haplotype ultimately can affect biologic and therapeutic phenotype and that individual SNPs may have poor predictive power as pharmacogenetic loci.

Population structure and recent evolution of Plasmodium falciparum

Plasmodium falciparum is the agent of malignant malaria, one of mankind's most severe maladies. The parasite exhibits antigenic polymorphisms that have been postulated to be ancient. We have proposed that the extant world populations of P. falciparum have derived from one single parasite, a cenancestor, within the last 5,000–50,000 years. This inference derives from the virtual or complete absence of synonymous nucleotide polymorphisms at genes not involved in immune or drug responses. Seeking to conciliate this claim with extensive antigenic polymorphism, we first note that allele substitutions or polymorphisms can arise very rapidly, even in a single generation, in large populations subject to strong natural selection. Second, new alleles can arise not only by single-nucleotide mutations, but also by duplication/deletion of short simple-repeat DNA sequences, a process several orders of magnitude faster than single-nucleotide mutation. We analyze three antigenic genes known to be extremely polymorphic: Csp, Msp-1, and Msp-2. We identify regions consisting of tandem or proximally repetitive short DNA sequences, including some previously unnoticed. We conclude that the antigenic polymorphisms are consistent with the recent origin of the world populations of P. falciparum inferred from the analysis of nonantigenic genes.

DNA sequences of Alu elements indicate a recent replacement of the human autosomal genetic complement.

DNA sequences of neutral nuclear autosomal loci, compared across diverse human populations, provide a previously untapped perspective into the mode and tempo of the emergence of modern humans and a critical comparison with published clonally inherited mitochondrial DNA and Y chromosome measurements of human diversity. We obtained over 55 kilobases of sequence from three autosomal loci encompassing Alu repeats for representatives of diverse human populations as well as orthologous sequences for other hominoid species at one of these loci. Nucleotide diversity was exceedingly low. Most individuals and populations were identical. Only a single nucleotide difference distinguished presumed ancestral alleles from descendants. These results differ from those expected if alleles from divergent archaic populations were maintained through multiregional continuity. The observed virtual lack of sequence polymorphism is the signature of a recent single origin for modern humans, with general replacement of archaic populations.

On signal sequence polymorphisms and diseases of distribution.

We report a previously unappreciated property of the signals that target organelle-specific proteins to their subcellular sites of action. Such targeting sequences are shown to be polymorphic. We discovered this polymorphism when we cloned the mitochondrial manganese-containing superoxide dismutase from cell lines of normal individuals and patients with genetic diseases of premature aging and compared their sequences to each other and to those previously reported. The polymorphism consists of a single nucleotide change in the region of the DNA that encodes the signal sequence such that either an alanine or valine is present. Subsequently, eight cell lines were analyzed and all three possible combinations of the two signal sequences were observed. Such signal sequence polymorphisms could result in diseases of distribution, where essential proteins are not properly targeted, thereby leading to absolute or relative deficiencies of critical enzymes within specific cellular compartments. Progeria and related syndromes may be diseases of distribution.

Mapping nucleosome position at single base-pair resolution by using site-directed hydroxyl radicals.

A base-pair resolution method for determining nucleosome position in vitro has been developed to com- plement existing, less accurate methods. Cysteaminyl EDTA was tethered to a recombinant histone octamer via a mutant histone H4 with serine 47 replaced by cysteine. When assembled into nucleosome core particles, the DNA could be cut site specifically by hydroxyl radical-catalyzed chain scission by using the Fenton reaction. Strand cleavage occurs mainly at a single nucleotide close to the dyad axis of the core particle, and assignment of this location via the symmetry of the nucleosome allows base-pair resolution mapping of the histone octamer position on the DNA. The positions of the histone octamer and H3H4 tetramer were mapped on a 146-bp Lytechinus variegatus 5S rRNA sequence and a twofold-symmetric derivative. The weakness of translational determinants of nucleosome positioning relative to the overall affinity of the histone proteins for this DNA is clearly demonstrated. The predominant location of both histone octamer and H3H4 tetramer assembled on the 5S rDNA is off center. Shifting the nucleosome core particle position along DNA within a conserved rotational phase could be induced under physiologically relevant conditions. Since nucleosome shifting has important consequences for chromatin structure and gene regulation, an approach to the thermodynamic characterization of this movement is proposed. This mapping method is potentially adaptable for determining nucleosome position in chromatin in vivo.