SNPs for parentage testing and traceability in globally diverse breeds of sheep

DNA-based parentage determination accelerates genetic improvement in sheep by increasing pedigree accuracy. Single nucleotide polymorphism (SNP) markers can be used for determining parentage and to provide unique molecular identifiers for tracing sheep products to their source. However, the utility of a particular "parentage SNP" varies by breed depending on its minor allele frequency (MAF) and its sequence context. Our aims were to identify parentage SNPs with exceptional qualities for use in globally diverse breeds and to develop a subset for use in North American sheep. Starting with genotypes from 2,915 sheep and 74 breed groups provided by the International Sheep Genomics Consortium (ISGC), we analyzed 47,693 autosomal SNPs by multiple criteria and selected 163 with desirable properties for parentage testing. On average, each of the 163 SNPs was highly informative (MAF≥0.3) in 48±5 breed groups. Nearby polymorphisms that could otherwise confound genetic testing were identified by whole genome and Sanger sequencing of 166 sheep from 54 breed groups. A genetic test with 109 of the 163 parentage SNPs was developed for matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry. The scoring rates and accuracies for these 109 SNPs were greater than 99% in a panel of North American sheep. In a blinded set of 96 families (sire, dam, and non-identical twin lambs), each parent of every lamb was identified without using the other parent's genotype. In 74 ISGC breed groups, the median estimates for probability of a coincidental match between two animals (PI), and the fraction of potential adults excluded from parentage (PE) were 1.1×10(-39) and 0.999987, respectively, for the 109 SNPs combined. The availability of a well-characterized set of 163 parentage SNPs facilitates the development of high-throughput genetic technologies for implementing accurate and economical parentage testing and traceability in many of the world's sheep breeds.

Genetic underpinnings of variation in brain ventricular volume: A genome-wide association study

The ventricular system is a critical component of the central nervous system (CNS) that is formed early in the developmental stages and remains functional through the lifetime. Changes in the ventricular system can be easily discerned via neuroimaging procedures and most of the time it reflects changes in the physiology of the CNS. In this study we attempted to identify specific genes associated with variation in ventricular volume in humans. Methods. We conducted a genome wide association (GWA) analysis of the volume of the lateral ventricles among 1605 individuals of European ancestry from two community based cohorts, the Genetics of Microangiopathic Brain Injury (GMBI; N=814) and Atherosclerosis Risk in Communities (ARIC; N=791). Significant findings from the analysis were tested for replication in both the cohorts and then meta-analyzed to get an estimate of overall significance. Results. In our GWA analyses, no single nucleotide polymorphism (SNP) reached a genome-wide significance of p<10−8. There were 25 SNPs in GMBI and 9 SNPs in ARIC that reached a threshold of p<10 −5. However, none of the top SNPs from each cohort were replicated in the other. In the meta-analysis, no SNP reached the genome-wide threshold of 5×10−8, but we identified five novel SNPs associated with variation in ventricular volume at the p<10 −5 level. Strongest association was for rs2112536 in an intergenic region on chromosome 5q33 (Pmeta= 8.46×10−7 ). The remaining four SNPs were located on chromosome 3q23 encompassing the gene for Calsyntenin-2 (CLSTN2). The SNPs with strongest association in this region were rs17338555 (Pmeta= 5.28×10 −6), rs9812091 (Pmeta= 5.89×10−6 ), rs9812283 (Pmeta= 5.97×10−6) and rs9833213 (Pmeta= 6.96×10−6). Conclusions. This GWA study of ventricular volumes in the community-based cohorts of European descent identifies potential locus on chromosomes 3 and 5. Further characterization of these loci may provide insights into pathophysiology of ventricular involvement in various neurological diseases.^

Instrumental variable method for the causal relationship between soluble receptor for advanced glycation end-products (sRAGE) and risk of pancreatic cancer

This thesis project is motivated by the potential problem of using observational data to draw inferences about a causal relationship in observational epidemiology research when controlled randomization is not applicable. Instrumental variable (IV) method is one of the statistical tools to overcome this problem. Mendelian randomization study uses genetic variants as IVs in genetic association study. In this thesis, the IV method, as well as standard logistic and linear regression models, is used to investigate the causal association between risk of pancreatic cancer and the circulating levels of soluble receptor for advanced glycation end-products (sRAGE). Higher levels of serum sRAGE were found to be associated with a lower risk of pancreatic cancer in a previous observational study (255 cases and 485 controls). However, such a novel association may be biased by unknown confounding factors. In a case-control study, we aimed to use the IV approach to confirm or refute this observation in a subset of study subjects for whom the genotyping data were available (178 cases and 177 controls). Two-stage IV method using generalized method of moments-structural mean models (GMM-SMM) was conducted and the relative risk (RR) was calculated. In the first stage analysis, we found that the single nucleotide polymorphism (SNP) rs2070600 of the receptor for advanced glycation end-products (AGER) gene meets all three general assumptions for a genetic IV in examining the causal association between sRAGE and risk of pancreatic cancer. The variant allele of SNP rs2070600 of the AGER gene was associated with lower levels of sRAGE, and it was neither associated with risk of pancreatic cancer, nor with the confounding factors. It was a potential strong IV (F statistic = 29.2). However, in the second stage analysis, the GMM-SMM model failed to converge due to non- concaveness probably because of the small sample size. Therefore, the IV analysis could not support the causality of the association between serum sRAGE levels and risk of pancreatic cancer. Nevertheless, these analyses suggest that rs2070600 was a potentially good genetic IV for testing the causality between the risk of pancreatic cancer and sRAGE levels. A larger sample size is required to conduct a credible IV analysis.^

Metagenomic Anlaysis of microsymbiont selection by the legume plant host

Rhizobium leguminosarum bv.viciae is able to establish nitrogen-fixing symbioses with legumes of the genera Pisum, Lens, Lathyrus and Vicia. Classic studies using trap plants (Laguerre et al., Young et al.) provided evidence that different plant hosts are able to select different rhizobial genotypes among those available in a given soil. However, these studies were necessarily limited by the paucity of relevant biodiversity markers. We have now reappraised this problem with the help of genomic tools. A well-characterized agricultural soil (INRA Bretennieres) was used as source of rhizobia. Plants of Pisum sativum, Lens culinaris, Vicia sativa and V. faba were used as traps. Isolates from 100 nodules were pooled, and DNA from each pool was sequenced (BGI-Hong Kong; Illumina Hiseq 2000, 500 bp PE libraries, 100 bp reads, 12 Mreads). Reads were quality filtered (FastQC, Trimmomatic), mapped against reference R. leguminosarum genomes (Bowtie2, Samtools), and visualized (IGV). An important fraction of the filtered reads were not recruited by reference genomes, suggesting that plant isolates contain genes that are not present in the reference genomes. For this study, we focused on three conserved genomic regions: 16S-23S rDNA, atpD and nodDABC, and a Single Nucleotide Polymorphism (SNP) analysis was carried out with meta / multigenomes from each plant. Although the level of polymorphism varied (lowest in the rRNA region), polymorphic sites could be identified that define the specific soil population vs. reference genomes. More importantly, a plant-specific SNP distribution was observed. This could be confirmed with many other regions extracted from the reference genomes (data not shown). Our results confirm at the genomic level previous observations regarding plant selection of specific genotypes. We expect that further, ongoing comparative studies on differential meta / multigenomic sequences will identify specific gene components of the plant-selected genotypes

Effect of the Legume Plant Host on Rhizobial Soil Population Dynamics

Rhizobium leguminosarum bv viciae (Rlv) is a soil bacterium able to establish specific root-nodule symbioses with legumes of four different genera: Pisum, Vicia, Lens and Lathyrus. Rlv isolates from nodules of any of these legumes can nodulate any of them; however, it has been shown that plants select specific rhizobial genotypes from those present in the soil (1,2). We have previously shown this at the genomic level by following a population genomics approach. Pool genomic sequences from 100 isolates from each of four plant species: P. sativum, L. culinaris, V. faba and V. sativa, show different, specific profiles at the single nucleotide polymorphism (SNP) level for relevant genes. In this work, the extent of Rlv selection from a well-characterized soil population by different legume plant hosts: P. sativum, L. culinaris, V. faba and V. sativa, after a medium-term mesocosm study is described. Direct soil isolates from each of these mesocosm studies have been tested for specific rhizobial genes (glnII and fnrN) and symbiotic genes (nodC and nifH). Different populations were characterized further by Sanger sequencing of both the rpoB phylogenetic marker gene and the symbiotic genes nodC and nifH. The distribution and size of the rhizobial population for each legume host showed changes during the medium-term mesocosm study. Particularly, a non-symbiotic group of rhizobia was enriched by all four hosts, in contrast to the symbiotic rhizobia profile, which was specific for each legume plant host.

The Pangenome of the Plant-Selected Rhizobium leguminosarum bv viciae Sub-popuations

Rhizobium leguminosarum bv viciae (Rlv) is a bacterium able to establish effective symbioses with four different legume genera: Pisum, Lens, Lathyrus and Vicia. Classic studies using trap plants have previously shown that, given a choice, different plants prefer specific genotypes of rhizobia, which are adapted to the host (1, 2). In previous work we have performed a Pool-Seq analysis bases on pooled DNA samples from Rlv nodule isolates obtained from Pisum sativum, Lens culinaris, Vicia fava and V. sativa plants, used as rhizobial traps. This experiment allowed us to test the host preference hypothesis: different plant hosts select specific sub-populations of rhizobia from the available population present in a given soil. We have observed that plant-selected sub-populations are different at the single nucleotide polymorphism (SNP) level. We have selected individual isolates from each sub-population (9 fava-bean isolates, 14 pea isolates 9 vetch isolates and 9 lentil isolates) and sequenced their genomes at draft level (ca. 30x, 90 contigs). Genomic analyses have been carried out using J-species and CMG-Biotools. All the isolates had similar genome size (7.5 Mb) and number of genes (7,300). The resulting Average Nucleotide Identity (ANIm) tree showed that Rhizobium leguminosarum bv viciae is a highly diverse group. Each plant-selected subpopulation showed a closed pangenome and core genomes of similar size (11,500 and 4,800 genes, respectively). The addition of all four sub-population results in a larger, closed pangenome of 19,040 genes and a core genome of similar size (4,392 genes). Each sub-population contains a characteristic set of genes but no universal, plant-specific genes were found. The core genome obtained from all four sub-populations is probably a representative core genome for Rhizobium leguminosarum, given that the reference genome (Rhizobium leguminosarum bv. viciae strain 3841) contains most of the core genome. We have also analyzed the symbiotic cluster (nod), and different nod cluster genotypes were found in each sub-population. Supported by MINECO (Consolider-Ingenio 2010, MICROGEN Project, CSD2009-00006).

Estudo de associação entre disfunção neurocognitiva, estresse oxidativa e polimorfismos em pacientes jovens com Transtornos Bipolar tipo I

O Transtorno Bipolar (TB) tipo I é uma doença caracterizada por episódios de mania e depressão recorrentes com importante prejuízo do funcionamento global e comprometimento das funções cognitivas. Além disso, sabe-se que o número de episódios de humor patológico ao longo da vida pode também influenciar o funcionamento cognitivo destes sujeitos. Neste cenário, ocorreu a necessidade de se investigar marcadores genéticos para disfunção cognitiva no TB com o objetivo de estudar este fenômeno. Dentre os potenciais genes responsáveis por influenciar a cognição destacam-se os polimorfismos funcionais do fator neurotrófico derivado do cérebro (BDNF), da catecol-O-metiltransferase (COMT), da apolipoproteína-E (APOE) e do canal de cálcio de baixa voltagem subunidade 1-C (CACNA1C). Sabe-se, também, que no TB os marcadores de estresse oxidativo estão aumentados durante todas as fases da doença, entretanto, não é claro qual impacto destes na disfunção cognitiva de indivíduos com TB. O objetivo dessa tese foi avaliar o desempenho cognitivo de pacientes jovens com bipolaridade tipo I e sua associação com o genótipo de BDNF, COMT, APOE e CACNA1C e também com os níveis plasmáticos de oxidação da guanosina (8-OHdG) e citosina (5-Mec) durante os episódios de humor, eutimia e em controles. Para investigar essa associação foram incluídos 116 pacientes (79 em episódio de humor patológico e 37 eutímicos) com diagnóstico de TB tipo I (DSMIV-TR); 97 controles saudáveis foram submetidos à avaliação neuropsicológica e coleta de sangue para extração de DNA visando genotipagem para BDNF (rs6265), COMT (rs4680; rs165599), APOE (rs429358 e rs7412), CACNA1C (rs1006737), 8-OhdG e 5-Mec. A análise dos dados obtidos revelou que pacientes portadores do genótipo Met/Met rs4680/rs165599 do COMT apresentam comprometimento cognitivo mais grave (função executiva, fluência verbal, memória e inteligência) comparado ao genótipo Val/Met ou Val/Val durante episódios maníacos ou mistos. Na mesma direção destes resultados, verificou-se que pacientes portadores do alelo Met rs4680 do COMT apresentam comprometimento do reconhecimento de emoções faciais em episódios de mania e depressão. Nenhum efeito do COMT foi observado em controles. O alelo de risco Met do CACNA1C se associou a um pior comprometimento executivo independente dos sintomas maníacos ou depressivos no TB, porém nenhum efeito se observou nos controles. O alelo Met do BDNF rs6265 ou a presença do alelo 4 da APOE não representa um fator que identifique um grupo com desempenho cognitivo diferenciado durante as fases do TB ou em controles. Sujeitos com TB apresentaram níveis mais elevados de 8-OHdG e tais níveis eram diretamente proporcionais ao número de episódios maníacos ao longo da vida, sugerindo um papel dos episódios hiperdopaminérgicos na oxidação das bases de DNA. Concluiu-se que a genotipagem para COMT e CACNA1C em pacientes com TB pode identificar um grupo de pacientes associados a pior disfunção cognitiva durante as fases maníacas e mistas do TB. Tal dado pode ser um indicador do envolvimento do sistema dopaminérgico e dos canais de cálcio de baixa voltagem na fisiopatologia da disfunção cognitiva no TB e deve ser explorado em outros estudos

Pulse Pressure magnifies the effect of COMT Val158Met on 15 Year Episodic Memory Trajectories

We investigated whether a physiological marker of cardiovascular health, pulse pressure (PP), and age magnified the effect of the functional COMT Val158Met (rs4680) polymorphism on 15-years cognitive trajectories [episodic memory (EM), visuospatial ability, and semantic memory] using data from 1585 non-demented adults from the Betula study. A multiple-group latent growth curve model was specified to gauge individual differences in change, and average trends therein. The allelic variants showed negligible differences across the cognitive markers in average trends. The older portion of the sample selectively age-magnified the effects of Val158Met on EM changes, resulting in greater decline in Val compared to homozygote Met carriers. This effect was attenuated by statistical control for PP. Further, PP moderated the effects of COMT on 15-years EM trajectories, resulting in greater decline in Val carriers, even after accounting for the confounding effects of sex, education, cardiovascular diseases (diabetes, stroke, and hypertension), and chronological age, controlled for practice gains. The effect was still present after excluding individuals with a history of cardiovascular diseases. The effects of cognitive change were not moderated by any other covariates. This report underscores the importance of addressing synergistic effects in normal cognitive aging, as the addition thereof may place healthy individuals at greater risk for memory decline.

Analysis of RNA-protein interactions by flow cytometry

Flow cytometry, in combination with advances in bead coding technologies, is maturing as a powerful high-throughput approach for analyzing molecular interactions. Applications of this technology include antibody assays and single nucleotide polymorphism mapping. This review describes the recent development of a microbead flow cytometric approach to analyze RNA-protein interactions and discusses emerging bead coding strategies that together will allow genome-wide identification of RNA-protein complexes. The microbead flow cytometric approach is flexible and provides new opportunities for functional genomic studies and small-molecule screening.

Polymorphisms in the 5 '-untranslated region of the human serotonin receptor 1B (HTR1B) gene affect gene expression

We present evidence of complex balancing regulation of HTR1B transcription by common polymorphisms in its promoter. Computational analysis of the HTR1B gene predicted that a 50 segment, spanning common DNA sequence variations, T-261G, A-161T, and -182INS/DEL-181, contained a putative functional promoter. Using a secreted alkaline phosphatase (SEAP) reporter gene system, we found that the haplotype -261G_-182INS-181_A-161 enhanced transcriptional activity 2.3-fold compared with the haplotype T-261_-182INS-181_A-161. Conversely, -161T reversed this, and the net effect when -261G and -161T were in the same haplotype (-261G_-182INS-181_-161T) was equivalent to the major haplotype (T-261_-182INS-181_A-161). Electrophoretic mobility shift experiments showed that -261G and -161T modify the binding of transcription factors (TFs): -261G generates a new AP2 binding site, while alleles A-161 and -161T exhibit different binding characteristics to AP1. T-261G and A-161T were found to be in linkage disequilibrium (LD) with G861C in a European ancestry population. Interestingly, G861C has been reported to be associated with several psychiatric disorders. Our results indicate that HTR1B is the target of substantial transcriptional genetic regulation by common haplotypes, which are in LD with the HTR1B single-nucleotide polymorphism (SNP) most commonly used in association studies.

Evaluation of haplotypes associated with copper toxicosis in Bedlington Terriers in Australia

Objective-To evaluate the haplotype distribution associated with the copper toxicosis gene and the segregation of the mutated allele in a Bedlington Terrier population in Australia. Animals-131 Bedlington Terriers. Procedure-Samples of DNA and RNA were obtained from each dog. Genetic status of each dog was evaluated by use of the DNA markers C04107; single nucleotide polymorphism (SNP), which was adjacent to exon 2 of Murr1; and a deletion marker for exon 2. A subgroup of the study population was evaluated by use of biochemical and histologic techniques to elucidate the correlation between genotype and phenotype. Results-We identified a recombination between the C04107 marker and Murr1 and a variation in a nucleotide in the splice site of exon 2 in our Bedlington Terrier cohort. Furthermore, we identified a novel haplotype associated with copper toxicosis in this cohort. Conclusions and Clinical Relevance-Our findings indicate that the deletion of exon 2 was not the sole cause of copper toxicosis, although only exon 2 deletion of Murr1 has been responsible for copper toxicosis in Bedlington Terriers. Although we failed to find a novel mutation in our cohort, we identified an affected dog family with an intact exon 2. Furthermore, we found that an SNP in the 5' splicing site of exon 2 may or may not be associated with a novel mutation of the Murr1 gene or other genes. Loss of linkage between the C04107 marker and the Murr1 gene was also identified in a certain family of dogs.

mecA locus diversity in methicillin-resistant Staphylococcus aureus isolates in Brisbane, Australia, and the development of a novel diagnostic procedure for the Western Samoan phage pattern clone

An emerging public health phenomenon is the increasing incidence of methicillin-resistant Staphylococcus aureus (MRSA) infections that are acquired outside of health care facilities. One lineage of community-acquired MRSA (CA-MRSA) is known as the Western Samoan phage pattern (WSPP) clone. The central aim of this study was to develop an efficient genotyping procedure for the identification of WSPP isolates. The approach taken was to make use of the highly variable region downstream of mecA in combination with a single nucleotide polymorphism (SNP) defined by the S. aureus multilocus sequence typing (MLST) database. The premise was that a combinatorial genotyping method that interrogated both a highly variable region and the genomic backbone would deliver a high degree of informative power relative to the number of genetic polymorphisms-interrogated. Thirty-five MRSA isolates were used for this study, and their gene contents and order downstream of mecA were determined. The CA-MRSA isolates were found to contain a truncated mecA downstream region consisting of mecA-HVR-IS431 mec-dcs-Ins117, and a PCR-based method for identifying this structure was developed. The hospital-acquired isolates were found to contain eight different mecA downstream regions, three of which were novel. The Minimum SNPs computer software program was used to mine the S. aureus MLST database, and the arcC 2726 polymorph was identified as 82% discriminatory for ST-30. A real-time PCR assay was developed to interrogate this SNP. We found that the assay for the truncated mecA downstream region in combination with the interrogation of arcC position 272 provided an unambiguous identification of WSPP isolates.

Evaluation of IDDM8 susceptibility locus in a Russian simplex family data set

Type 1 diabetes (TID) susceptibility locus, IDDM8, has been accurately mapped to 200 kilobases at the terminal end of chromosome 6q27. This is within the region which harbours a cluster of three genes encoding proteasome subunit beta 1 (PMSB1), TATA-box binding protein (TBP) and a homologue of mouse programming cell death activator 2 (PDCD2). In this study, we evaluated whether these genes contribute to TID susceptibility using the transmission disequilibrium test of the data set from 114 affected Russian simplex families. The A allele of the G/A1180 single nucleotide polymorphism (SNP) at the PDCD2 gene, which was significant in its preferential transfer from parents to diabetic children (75 transmissions vs. 47 non-transmissionS, x(2) = 12.85, P corrected = 0.0038), was found to be associated with T1D. G/A1180 dimorphism and two other SNPs, C/T771 TBP and G/T(-271) PDCD2, were shown to share three common haplotypes, two of which (A-T-G and A-T-T) have been associated with higher development risk of TID. The third haplotype (G-T-G) was related to having a lower risk of disease. These findings suggest that the PDCD2 gene is a likely susceptibility gene for TID within IDDM8. However, it was not possible to exclude the TBP gene from being another putative susceptibility gene in this region. (c) 2005 Elsevier Ltd. All rights reserved.

Double-strand break repair gene polymorphisms and risk of breast or ovarian cancer

Deficiencies in DNA repair have been hypothesized to increase cancer risk and excess cancer incidence is a feature of inherited diseases caused by defects in DNA damage recognition and repair. We investigated, using a case-control design, whether the double-strand break repair gene polymorphisms RAD51 5' untranslated region -135 G > C, XRCC2 R188H G > A, and XRCC3 T241M C > T were associated with risk of breast or ovarian cancer in Australian women. Sample sets included 1,456 breast cancer cases and 793 age-matched controls ages under 60 years of age, 549 incident ovarian cancer cases, and 335 controls of similar age distribution. For the total sample and the subsample of Caucasian women, there were no significant differences in genotype distribution between breast cancer cases and controls or between ovarian cancer cases and combined control groups. The crude odds ratios (OR) and 95% confidence intervals (95% CI) associated with the RAD51 GC/CC genotype frequency was OR, 1.10; 95% CI, 0.80-1.41 for breast cancer and OR, 1.22; 95% CI, 0.92-1.62 for ovarian cancer. Similarly, there were no increased risks associated with the XRCC2 GA/AA genotype (OR, 0.98; 95% CI, 0.76-1.26 for breast cancer and OR, 0.93; 95% CI, 0.69-1.25 for ovarian cancer) or the XRCC3 CT/TT genotype (OR, 0.92; 95% Cl, 0.77-1.10 for breast cancer and OR, 0.87; 95% CI, 0.71-1.08 for ovarian cancer). Results were little changed after adjustment for age and other measured risk factors. Although there was little statistical power to detect modest increases in risk for the homozygote variant genotypes, particularly for the rare RAD51 and XRCC2 variants, the data suggest that none of these variants play a major role in the etiology of breast or ovarian cancer.