Effective population size of an indigenous Swiss cattle breed estimated from linkage disequilibrium

Effective population size is an important parameter for the assessment of genetic diversity within a livestock population and its development over time. If pedigree information is not available, linkage disequilibrium (LD) analysis might offer an alternative perspective for the estimation of effective population size. In this study, 128 individuals of the Swiss Eringer breed were genotyped using the Illumina BovineSNP50 beadchip. We set bin size at 50 kb for LD analysis, assuming that LD for proximal single nucleotide polymorphism (SNP)-pairs reflects distant breeding history while LD from distal SNP-pairs would reflect near history. Recombination rates varied among different regions of the genome. The use of physical distances as an approximation of genetic distances (e.g. setting 1 Mb = 0.01 Morgan) led to an upward bias in LD-based estimates of effective population size for generations beyond 50, while estimates for recent history were unaffected. Correction for restricted sample size did not substantially affect these results. LD-based actual effective population size was estimated in the range of 87-149, whereas pedigree-based effective population size resulted in 321 individuals. For conservation purposes, requiring knowledge of recent history (<50 generations), approximation assuming constant recombination rate seemed adequate.

Age-modulated association between prefrontal NAA and the BDNF gene

Brain-derived neurotrophic factor (BDNF) has been implicated in the pathophysiology of psychiatric and neurological disorders and in the mechanisms of antidepressant pharmacotherapy. Psychiatric and neurological conditions have also been associated with reduced brain levels of N-acetyl-aspartate (NAA), which has been used as a putative marker of neural integrity. However, few studies have explored the relationship between BDNF polymorphisms and NAA levels directly. Here, we present data from a single-voxel proton magnetic resonance spectroscopy study of 64 individuals and explore the relationship between BDNF polymorphisms and prefrontal NAA level. Our results indicate an association between a single nucleotide polymorphism (SNP) within BDNF, known as rs1519480, and reduced NAA level (p = 0.023). NAA levels were further predicted by age and Asian ancestry. There was a significant rs1519480 × age interaction on NAA level (p = 0.031). Specifically, the effect of rs1519480 on NAA level became significant at age ⩾34.17 yr. NAA level decreased with advancing age for genotype TT (p = 0.001) but not for genotype CT (p = 0.82) or CC (p = 0.34). Additional in silico analysis of 142 post-mortem brain samples revealed an association between the same SNP and reduced BDNF mRNA expression in the prefrontal cortex. The rs1519480 SNP influences BDNF mRNA expression and has an impact on prefrontal NAA level over time. This genetic mechanism may contribute to inter-individual variation in cognitive performance seen during normal ageing, as well as contributing to the risk for developing psychiatric and neurological conditions.

Beta-D-glucoside utilization by Mycoplasma mycoides subsp. mycoides SC: possible involvement in the control of cytotoxicity towards bovine lung cells

BACKGROUND: Contagious bovine pleuropneumonia (CBPP) caused by Mycoplasma mycoides subsp. mycoides small-colony type (SC) is among the most serious threats for livestock producers in Africa. Glycerol metabolism-associated H2O2 production seems to play a crucial role in virulence of this mycoplasma. A wide number of attenuated strains of M. mycoides subsp. mycoides SC are currently used in Africa as live vaccines. Glycerol metabolism is not affected in these vaccine strains and therefore it does not seem to be the determinant of their attenuation. A non-synonymous single nucleotide polymorphism (SNP) in the bgl gene coding for the 6-phospho-beta-glucosidase (Bgl) has been described recently. The SNP differentiates virulent African strains isolated from outbreaks with severe CBPP, which express the Bgl isoform Val204, from strains to be considered less virulent isolated from CBPP outbreaks with low mortality and vaccine strains, which express the Bgl isoform Ala204. RESULTS: Strains of M. mycoides subsp. mycoides SC considered virulent and possessing the Bgl isoform Val204, but not strains with the Bgl isoform Ala204, do trigger elevated levels of damage to embryonic bovine lung (EBL) cells upon incubation with the disaccharides (i.e., beta-D-glucosides) sucrose and lactose. However, strains expressing the Bgl isoform Val204 show a lower hydrolysing activity on the chromogenic substrate p-nitrophenyl-beta-D-glucopyranoside (pNPbG) when compared to strains that possess the Bgl isoform Ala204. Defective activity of Bgl in M. mycoides subsp. mycoides SC does not lead to H2O2 production. Rather, the viability during addition of beta-D-glucosides in medium-free buffers is higher for strains harbouring the Bgl isoform Val204 than for those with the isoform Ala204. CONCLUSION: Our results indicate that the studied SNP in the bgl gene is one possible cause of the difference in bacterial virulence among strains of M. mycoides subsp. mycoides SC. Bgl does not act as a direct virulence factor, but strains possessing the Bgl isoform Val204 with low hydrolysing activity are more prone to survive in environments that contain high levels of beta-D-glucosides, thus contributing in some extent to mycoplasmaemia.

Heterosis for biomass-related traits in Arabidopsis investigated by quantitative trait loci analysis of the triple testcross design with recombinant inbred lines

Arabidopsis thaliana has emerged as a leading model species in plant genetics and functional genomics including research on the genetic causes of heterosis. We applied a triple testcross (TTC) design and a novel biometrical approach to identify and characterize quantitative trait loci (QTL) for heterosis of five biomass-related traits by (i) estimating the number, genomic positions, and genetic effects of heterotic QTL, (ii) characterizing their mode of gene action, and (iii) testing for presence of epistatic effects by a genomewide scan and marker x marker interactions. In total, 234 recombinant inbred lines (RILs) of Arabidopsis hybrid C24 x Col-0 were crossed to both parental lines and their F1 and analyzed with 110 single-nucleotide polymorphism (SNP) markers. QTL analyses were conducted using linear transformations Z1, Z2, and Z3 calculated from the adjusted entry means of TTC progenies. With Z1, we detected 12 QTL displaying augmented additive effects. With Z2, we mapped six QTL for augmented dominance effects. A one-dimensional genome scan with Z3 revealed two genomic regions with significantly negative dominance x additive epistatic effects. Two-way analyses of variance between marker pairs revealed nine digenic epistatic interactions: six reflecting dominance x dominance effects with variable sign and three reflecting additive x additive effects with positive sign. We conclude that heterosis for biomass-related traits in Arabidopsis has a polygenic basis with overdominance and/or epistasis being presumably the main types of gene action.

A MIF haplotype is associated with the outcome of patients with severe sepsis: a case control study

BACKGROUND: Macrophage migration inhibitory factor (MIF) plays an important regulatory role in sepsis. In the promoter region a C/G single nucleotide polymorphism (SNP) at position -173 (rs755622) and a CATT5-8 microsatellite at position -794 are related to modified promoter activity. The purpose of the study was to analyze their association with the incidence and outcome of severe sepsis. METHODS: Genotype distributions and allele frequencies in 169 patients with severe sepsis, 94 healthy blood donors and 183 postoperative patients without signs of infection or inflammation were analyzed by real time PCR and Sequence analysis. All included individuals were Caucasians. RESULTS: Genotype distribution and allele frequencies of severe sepsis patients were comparable to both control groups. However, the genotype and allele frequencies of both polymorphisms were associated significantly with the outcome of severe sepsis. The highest risk of dying from severe sepsis was detectable in patients carrying a haplotype with the alleles -173 C and CATT7 (p = 0.0005, fisher exact test, RR = 1,806, CI: 1.337 to 2.439). CONCLUSION: The haplotype with the combination of the -173 C allele and the -794 CATT7 allele may not serve as a marker for susceptibility to sepsis, but may help identify septic patients at risk of dying.

Quantitative sequence analysis of FBN1 premature termination codons provides evidence for incomplete NMD in leukocytes

We improved, evaluated, and used Sanger sequencing for quantification of single nucleotide polymorphism (SNP) variants in transcripts and gDNA samples. This improved assay resulted in highly reproducible relative allele frequencies (e.g., for a heterozygous gDNA 50.0+/-1.4%, and for a missense mutation-bearing transcript 46.9+/-3.7%) with a lower detection limit of 3-9%. It provided excellent accuracy and linear correlation between expected and observed relative allele frequencies. This sequencing assay, which can also be used for the quantification of copy number variations (CNVs), methylations, mosaicisms, and DNA pools, enabled us to analyze transcripts of the FBN1 gene in fibroblasts and blood samples of patients with suspected Marfan syndrome not only qualitatively but also quantitatively. We report a total of 18 novel and 19 known FBN1 sequence variants leading to a premature termination codon (PTC), 26 of which we analyzed by quantitative sequencing both at gDNA and cDNA levels. The relative amounts of PTC-containing FBN1 transcripts in fresh and PAXgene-stabilized blood samples were significantly higher (33.0+/-3.9% to 80.0+/-7.2%) than those detected in affected fibroblasts with inhibition of nonsense-mediated mRNA decay (NMD) (11.0+/-2.1% to 25.0+/-1.8%), whereas in fibroblasts without NMD inhibition no mutant alleles could be detected. These results provide evidence for incomplete NMD in leukocytes and have particular importance for RNA-based analyses not only in FBN1 but also in other genes.

Comparative genetic analyses point to HCP5 as susceptibility locus for HCV-associated hepatocellular carcinoma

BACKGROUND & AIMS: Recently, genetic variations in MICA (lead single nucleotide polymorphism [SNP] rs2596542) were identified by a genome-wide association study (GWAS) to be associated with hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC) in Japanese patients. In the present study, we sought to determine whether this SNP is predictive of HCC development in the Caucasian population as well. METHODS: An extended region around rs2596542 was genotyped in 1924 HCV-infected patients from the Swiss Hepatitis C Cohort Study (SCCS). Pair-wise correlation between key SNPs was calculated both in the Japanese and European populations (HapMap3: CEU and JPT). RESULTS: To our surprise, the minor allele A of rs2596542 in proximity of MICA appeared to have a protective impact on HCC development in Caucasians, which represents an inverse association as compared to the one observed in the Japanese population. Detailed fine-mapping analyses revealed a new SNP in HCP5 (rs2244546) upstream of MICA as strong predictor of HCV-related HCC in the SCCS (univariable p=0.027; multivariable p=0.0002, odds ratio=3.96, 95% confidence interval=1.90-8.27). This newly identified SNP had a similarly directed effect on HCC in both Caucasian and Japanese populations, suggesting that rs2244546 may better tag a putative true variant than the originally identified SNPs. CONCLUSIONS: Our data confirms the MICA/HCP5 region as susceptibility locus for HCV-related HCC and identifies rs2244546 in HCP5 as a novel tagging SNP. In addition, our data exemplify the need for conducting meta-analyses of cohorts of different ethnicities in order to fine map GWAS signals.

Multiple Hypersensitivities Including Recurrent Airway Obstruction, Insect Bite Hypersensitivity, and Urticaria in 2 Warmblood Horse Populations.

BACKGROUND Multiple hypersensitivities (MHS) have been described in humans, cats, and dogs, but not horses. HYPOTHESES Horses suffering from recurrent airway obstruction (RAO), insect bite hypersensitivity (IBH), or urticaria (URT) will have an increased risk of also being affected by another one of these hypersensitivities. This predisposition for MHS also will be associated with decreased shedding of strongylid eggs in feces and with a single nucleotide polymorphism (SNP BIEC2-224511), previously shown to be associated with RAO. ANIMALS The first population (P1) included 119 randomly sampled horses representative of the Swiss sporthorse population; the replication population (P2) included 210 RAO-affected Warmblood horses and 264 RAO-unaffected controls. All horses were Warmbloods, 14 years or older. METHODS Associations between disease phenotypes (RAO, IBH, URT, MHS) fecal egg counts, the SNP BIEC2-224511 as well as management and environmental factors were investigated. RESULTS In P1, RAO-affected horses had a 13.1 times higher odds ratio (OR) of also suffering from IBH (P = .004). In P2, the respective OR was 7.4 (P = .002) and IBH-affected horses also showed a 7.1 times increased OR of concomitantly suffering from URT (P < .001). IBH, URT, and MHS phenotypes were significantly associated with the absence of nematode eggs in the feces. CONCLUSIONS AND CLINICAL IMPORTANCE This is the first report of MHS in horses. Specifically, an increased risk for IBH should be expected in RAO-affected horses.

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).

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.