Whole Genome Sequencing Reveals Local Transmission Patterns of Mycobacterium bovis in Sympatric Cattle and Badger Populations

Whole genome sequencing (WGS) technology holds great promise as a tool for the forensic epidemiology of bacterial pathogens. It is likely to be particularly useful for studying the transmission dynamics of an observed epidemic involving a largely unsampled 'reservoir' host, as for bovine tuberculosis (bTB) in British and Irish cattle and badgers. BTB is caused by Mycobacterium bovis, a member of the M. tuberculosis complex that also includes the aetiological agent for human TB. In this study, we identified a spatio-temporally linked group of 26 cattle and 4 badgers infected with the same Variable Number Tandem Repeat (VNTR) type of M. bovis. Single-nucleotide polymorphisms (SNPs) between sequences identified differences that were consistent with bacterial lineages being persistent on or near farms for several years, despite multiple clear whole herd tests in the interim. Comparing WGS data to mathematical models showed good correlations between genetic divergence and spatial distance, but poor correspondence to the network of cattle movements or within-herd contacts. Badger isolates showed between zero and four SNP differences from the nearest cattle isolate, providing evidence for recent transmissions between the two hosts. This is the first direct genetic evidence of M. bovis persistence on farms over multiple outbreaks with a continued, ongoing interaction with local badgers. However, despite unprecedented resolution, directionality of transmission cannot be inferred at this stage. Despite the often notoriously long timescales between time of infection and time of sampling for TB, our results suggest that WGS data alone can provide insights into TB epidemiology even where detailed contact data are not available, and that more extensive sampling and analysis will allow for quantification of the extent and direction of transmission between cattle and badgers. © 2012 Biek et al.

Human-directed gazing behaviour in puppies and adult dogs, Canis lupus familiaris

Recent evidence indicates that dogs' sociocognitive abilities and behaviour in a test situation are shaped by both genetic factors and life experiences. We used the 'unsolvable task' paradigm to investigate the effect of breed and age/experience on the use of human-directed gazing behaviour. Following a genetic classification based on recent genome analyses, dogs were allocated to three breed groups, namely Primitive, Hunting/Herding and Molossoid. Furthermore, we tested dogs at 2 months, 4.5. months and as adults. The test consisted of three solvable trials in which dogs could obtain food by manipulating a plastic container followed by an unsolvable trial in which obtaining the food became impossible. The dogs' behaviour towards the apparatus and the people present was analysed. At 2 months no breed group differences emerged and although human-directed gazing behaviour was observed in approximately half of the pups, it occurred for brief periods, suggesting that the aptitude to use human-directed gazing as a request for obtaining help probably develops at a later date when dogs have had more experience with human communication. Breed group differences, however, did emerge strongly in adult dogs and, although less pronounced, also in 4.5-month-old subjects, with dogs in the Hunting/Herding group showing significantly more human-directed gazing behaviour than dogs in the other two breed groups. These results suggest that, although the domestication process may have shaped the dog's human-directed communicative abilities, the later selection for specific types of work might also have had a significant impact on their emergence. © 2011 The Association for the Study of Animal Behaviour.

Cellular responses to genome mistranslation

Low level protein synthesis errors can have profound effects on normal cell physiology and disease development, namely neurodegeneration, cancer and aging. The biology of errors introduced into proteins during mRNA translation, herein referred as mistranslation, is not yet fully understood. In order to shed new light into this biological phenomenon, we have engineered constitutive codon misreading in S. cerevisiae, using a mutant tRNA that misreads leucine CUG codons as serine, representing a 240 fold increase in mRNA translational error relative to typical physiological error (0.0001%). Our studies show that mistranslation induces autophagic activity, increases accumulation of insoluble proteins, production of reactive oxygen species, and morphological disruption of the mitochondrial network. Mistranslation also up-regulates the expression of the longevity gene PNC1, which is a regulator of Sir2p deacetylase activity. We show here that both PNC1 and SIR2 are involved in the regulation of autophagy induced by mistranslation, but not by starvation-induced autophagy. Mistranslation leads to P-body but not stress-granule assembly, down-regulates the expression of ribosomal protein genes and increases slightly the selective degradation of ribosomes (ribophagy). The study also indicates that yeast cells are much more resistant to mistranslation than expected and highlights the importance of autophagy in the cellular response to mistranslation. Morpho-functional alterations of the mitochondrial network are the most visible phenotype of mistranslation. Since most of the basic cellular processes are conserved between yeast and humans, this study reinforces the importance of yeast as a model system to study mistranslation and suggests that oxidative stress and accumulation of misfolded proteins arising from aberrant protein synthesis are important causes of the cellular degeneration observed in human diseases associated to mRNA mistranslation.

Phenotypic consequences of genome mistranslation in Candida albicans

Várias espécies do género Candida traduzem o codão CUG de leucine como serina. Em C. albicans este codão é traduzido pelo tRNACAG Ser de serina que é reconhecido por leucil- e seril-tRNA sintetases (LeuRS e SerRS), permitindo a incorporação de leucina ou serina em posições com CUG. Em condições padrão de crescimento os codões CUG é incorporam 3% de leucina e 97% de serina, no entanto estes valores são flexíveis uma vez que a incorporação de serina pode variar entre 0.6% e 5% em resposta a condições de stress. Estudos anteriores realizados in vivo em Escherichia coli sugeriram que a ambiguidade em codões CUG é regulada pela SerRS. De facto, o gene da SerRS de C. albicans tem um codão CUG na posição 197 (Ser197) cuja descodificação ambígua resulta na produção de duas isoformas de SerRS. A isoforma SerRS_Leu197 é mais ativa, apesar de menos estável, que a isoforma SerRS_Ser197, suportando a ideia da existência de um feedback loop negativo, envolvendo estas duas isoformas de SerRS, a enzima LeuRS e o tRNACAG Ser, que mantem os níveis de incorporação de leucina no codões CUG baixos. Nesta tese demonstramos que tal mecanismo não é operacional nas células de C. albicans. De facto, os níveis de incorporação de leucina em codões CUG flutuam drasticamente em resposta a alterações ambientais. Por exemplo, a incorporação de leucina pode chegar a níveis de 49.33% na presença de macrófagos e anfotericina B, mostrando a notória tolerância de C. albicans à ambiguidade. Para compreender a relevância biológica da ambiguidade do código genético em C. albicans construímos estirpes que incorporam serina em vários codões. Apesar da taxa crescimento ter sido negativamente afetada em condições padrão de crescimento, as estirpes construídas crescem favoravelmente em várias condições de stresse, sugerindo que a ambiguidade desempenha um papel importante na adaptação a novos nichos ecológicos. O transcriptoma das estirpes construídas de C. albicans e Saccharomyces. cerevisiae mostram que as leveduras respondem à ambiguidade dos codões de modo distinto. A ambiguidade induziu uma desregulação moderada da expressão génica de C. albicans, mas ativou uma resposta comum ao stresse em S. cerevisiae. O único processo celular que foi induzido na maioria das estirpes foi a oxidação redução. De salientar, que enriquecimento em elementos cis de fatores de transcrição que regulam a resposta à ambiguidade em ambas as leveduras foi distinta, sugerindo que ambas respondem ao stresse de modo diferente. Na globalidade, o nosso estudo aprofunda o conhecimento da elevada tolerância à ambiguidade de codões em C. albicans. Os resultados sugerem que este fungo usa a ambiguidade do codão CUG durante infeção, possivelmente para modular a sua interação com o hospedeiro e a resposta a drogas antifúngicas.

Quantifying diet-induced metabolic changes of the human gut microbiome

The human gut microbiome is known to be associated with various human disorders, but a major challenge is to go beyond association studies and elucidate causalities. Mathematical modeling of the human gut microbiome at a genome scale is a useful tool to decipher microbe-microbe, diet-microbe and microbe-host interactions. Here, we describe the CASINO (Community And Systems-level INteractive Optimization) toolbox, a comprehensive computational platform for analysis of microbial communities through metabolic modeling. We first validated the toolbox by simulating and testing the performance of single bacteria and whole communities in vitro. Focusing on metabolic interactions between the diet, gut microbiota, and host metabolism, we demonstrated the predictive power of the toolbox in a diet-intervention study of 45 obese and overweight individuals and validated our predictions by fecal and blood metabolomics data. Thus, modeling could quantitatively describe altered fecal and serum amino acid levels in response to diet intervention.

Evolutionary dynamics of Chlamydia trachomatis genome and identification of molecular patterns of hypothetical protein coding genes

The obligate intracellular bacterium Chlamydia trachomatis is a human pathogen of major public health significance. Strains can be classified into 15 main serovars (A to L3) that preferentially cause ocular infections (A-C), genital infections (D-K) or lymphogranuloma venereum (LGV) (L1-L3), but the molecular basis behind their distinct tropism, ecological success and pathogenicity is not welldefined. Most chlamydial research demands culture in eukaryotic cell lines, but it is not known if stains become laboratory adapted. By essentially using genomics and transcriptomics, we aimed to investigate the evolutionary patterns underlying the adaptation of C. trachomatis to the different human tissues, given emphasis to the identification of molecular patterns of genes encoding hypothetical proteins, and to understand the adaptive process behind the C. trachomatis in vivo to in vitro transition. Our results highlight a positive selection-driven evolution of C. trachomatis towards nichespecific adaptation, essentially targeting host-interacting proteins, namely effectors and inclusion membrane proteins, where some of them also displayed niche-specific expression patterns. We also identified potential "ocular-specific" pseudogenes, and pointed out the major gene targets of adaptive mutations associated with LGV infections. We further observed that the in vivo-derived genetic makeup of C. trachomatis is not significantly compromised by its long-term laboratory propagation. In opposition, its introduction in vitro has the potential to affect the phenotype, likely yielding virulence attenuation. In fact, we observed a "genital-specific" rampant inactivation of the virulence gene CT135, which may impact the interpretation of data derived from studies requiring culture. Globally, the findings presented in this Ph.D. thesis contribute for the understanding of C.trachomatis adaptive evolution and provides new insights into the biological role of C. trachomatishypothetical proteins. They also launch research questions for future functional studies aiming toclarify the determinants of tissue tropism, virulence or pathogenic dissimilarities among C. trachomatisstrains.

Side effects of genome structural changes.

The first extensive catalog of structural human variation was recently released. It showed that large stretches of genomic DNA that vary considerably in copy number were extremely abundant. Thus it is conceivable that they play a major role in functional variation. Consistently, genomic insertions and deletions were shown to contribute to phenotypic differences by modifying not only the expression levels of genes within the aneuploid segments but also of normal copy-number neighboring genes. In this report, we review the possible mechanisms behind this latter effect.

Genetic variation in IL28B is associated with chronic hepatitis C and treatment failure: a genome-wide association study.

BACKGROUND & AIMS: Hepatitis C virus (HCV) induces chronic infection in 50% to 80% of infected persons; approximately 50% of these do not respond to therapy. We performed a genome-wide association study to screen for host genetic determinants of HCV persistence and response to therapy. METHODS: The analysis included 1362 individuals: 1015 with chronic hepatitis C and 347 who spontaneously cleared the virus (448 were coinfected with human immunodeficiency virus [HIV]). Responses to pegylated interferon alfa and ribavirin were assessed in 465 individuals. Associations between more than 500,000 single nucleotide polymorphisms (SNPs) and outcomes were assessed by multivariate logistic regression. RESULTS: Chronic hepatitis C was associated with SNPs in the IL28B locus, which encodes the antiviral cytokine interferon lambda. The rs8099917 minor allele was associated with progression to chronic HCV infection (odds ratio [OR], 2.31; 95% confidence interval [CI], 1.74-3.06; P = 6.07 x 10(-9)). The association was observed in HCV mono-infected (OR, 2.49; 95% CI, 1.64-3.79; P = 1.96 x 10(-5)) and HCV/HIV coinfected individuals (OR, 2.16; 95% CI, 1.47-3.18; P = 8.24 x 10(-5)). rs8099917 was also associated with failure to respond to therapy (OR, 5.19; 95% CI, 2.90-9.30; P = 3.11 x 10(-8)), with the strongest effects in patients with HCV genotype 1 or 4. This risk allele was identified in 24% of individuals with spontaneous HCV clearance, 32% of chronically infected patients who responded to therapy, and 58% who did not respond (P = 3.2 x 10(-10)). Resequencing of IL28B identified distinct haplotypes that were associated with the clinical phenotype. CONCLUSIONS: The association of the IL28B locus with natural and treatment-associated control of HCV indicates the importance of innate immunity and interferon lambda in the pathogenesis of HCV infection.

A glycine-arginine domain in control of the human MRE11 DNA repair protein.

Human MRE11 is a key enzyme in DNA double-strand break repair and genome stability. Human MRE11 bears a glycine-arginine-rich (GAR) motif that is conserved among multicellular eukaryotic species. We investigated how this motif influences MRE11 function. Human MRE11 alone or a complex of MRE11, RAD50, and NBS1 (MRN) was methylated in insect cells, suggesting that this modification is conserved during evolution. We demonstrate that PRMT1 interacts with MRE11 but not with the MRN complex, suggesting that MRE11 arginine methylation occurs prior to the binding of NBS1 and RAD50. Moreover, the first six methylated arginines are essential for the regulation of MRE11 DNA binding and nuclease activity. The inhibition of arginine methylation leads to a reduction in MRE11 and RAD51 focus formation on a unique double-strand break in vivo. Furthermore, the MRE11-methylated GAR domain is sufficient for its targeting to DNA damage foci and colocalization with gamma-H2AX. These studies highlight an important role for the GAR domain in regulating MRE11 function at the biochemical and cellular levels during DNA double-strand break repair.

An atlas of human gene expression from massively parallel signature sequencing (MPSS).

We have used massively parallel signature sequencing (MPSS) to sample the transcriptomes of 32 normal human tissues to an unprecedented depth, thus documenting the patterns of expression of almost 20,000 genes with high sensitivity and specificity. The data confirm the widely held belief that differences in gene expression between cell and tissue types are largely determined by transcripts derived from a limited number of tissue-specific genes, rather than by combinations of more promiscuously expressed genes. Expression of a little more than half of all known human genes seems to account for both the common requirements and the specific functions of the tissues sampled. A classification of tissues based on patterns of gene expression largely reproduces classifications based on anatomical and biochemical properties. The unbiased sampling of the human transcriptome achieved by MPSS supports the idea that most human genes have been mapped, if not functionally characterized. This data set should prove useful for the identification of tissue-specific genes, for the study of global changes induced by pathological conditions, and for the definition of a minimal set of genes necessary for basic cell maintenance. The data are available on the Web at http://mpss.licr.org and http://sgb.lynxgen.com.

Genome-wide association analysis identifies 20 loci that influence adult height.

Adult height is a model polygenic trait, but there has been limited success in identifying the genes underlying its normal variation. To identify genetic variants influencing adult human height, we used genome-wide association data from 13,665 individuals and genotyped 39 variants in an additional 16,482 samples. We identified 20 variants associated with adult height (P < 5 x 10(-7), with 10 reaching P < 1 x 10(-10)). Combined, the 20 SNPs explain approximately 3% of height variation, with a approximately 5 cm difference between the 6.2% of people with 17 or fewer 'tall' alleles compared to the 5.5% with 27 or more 'tall' alleles. The loci we identified implicate genes in Hedgehog signaling (IHH, HHIP, PTCH1), extracellular matrix (EFEMP1, ADAMTSL3, ACAN) and cancer (CDK6, HMGA2, DLEU7) pathways, and provide new insights into human growth and developmental processes. Finally, our results provide insights into the genetic architecture of a classic quantitative trait.

Functional and structural basis for a bacteriophage homolog of human RAD52.

In eukaryotes, homologous recombination proteins such as RAD51 and RAD52 play crucial roles in DNA repair and genome stability. Human RAD52 is a member of a large single-strand annealing protein (SSAP) family [1] and stimulates Rad51-dependent recombination [2, 3]. In prokaryotes and phages, it has been difficult to establish the presence of RAD52 homologs with conserved sequences. Putative SSAPs were recently found in several phages that infect strains of Lactococcus lactis[4]. One of these SSAPs was identified as Sak and was found in the virulent L. lactis phage ul36, which belongs to the Siphoviridae family [4, 5]. In this study, we show that Sak is homologous to the N terminus of human RAD52. Purified Sak binds single-stranded DNA (ssDNA) preferentially over double-stranded DNA (dsDNA) and promotes the renaturation of long complementary ssDNAs. Sak also binds RecA and stimulates homologous recombination reactions. Mutations shown to modulate RAD52 DNA binding [6] affect Sak similarly. Remarkably, electron-microscopic reconstruction of Sak reveals an undecameric (11) subunit ring, similar to the crystal structure of the N-terminal fragment of human RAD52 [7, 8]. For the first time, we propose a viral homolog of RAD52 at the amino acid, phylogenic, functional, and structural levels.

The European dimension for the mouse genome mutagenesis program.

The European Mouse Mutagenesis Consortium is the European initiative contributing to the international effort on functional annotation of the mouse genome. Its objectives are to establish and integrate mutagenesis platforms, gene expression resources, phenotyping units, storage and distribution centers and bioinformatics resources. The combined efforts will accelerate our understanding of gene function and of human health and disease.

A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance.

Recent genome-wide association studies have described many loci implicated in type 2 diabetes (T2D) pathophysiology and β-cell dysfunction but have contributed little to the understanding of the genetic basis of insulin resistance. We hypothesized that genes implicated in insulin resistance pathways might be uncovered by accounting for differences in body mass index (BMI) and potential interactions between BMI and genetic variants. We applied a joint meta-analysis approach to test associations with fasting insulin and glucose on a genome-wide scale. We present six previously unknown loci associated with fasting insulin at P < 5 × 10(-8) in combined discovery and follow-up analyses of 52 studies comprising up to 96,496 non-diabetic individuals. Risk variants were associated with higher triglyceride and lower high-density lipoprotein (HDL) cholesterol levels, suggesting a role for these loci in insulin resistance pathways. The discovery of these loci will aid further characterization of the role of insulin resistance in T2D pathophysiology.