Genome sequencing of the extinct Eurasian wild aurochs, Bos primigenius, illuminates the phylogeography and evolution of cattle

Background

Domestication of the now-extinct wild aurochs, Bos primigenius, gave rise to the two major domestic extant cattle taxa, B. taurus and B. indicus. While previous genetic studies have shed some light on the evolutionary relationships between European aurochs and modern cattle, important questions remain unanswered, including the phylogenetic status of aurochs, whether gene flow from aurochs into early domestic populations occurred, and which genomic regions were subject to selection processes during and after domestication. Here, we address these questions using whole-genome sequencing data generated from an approximately 6,750-year-old British aurochs bone and genome sequence data from 81 additional cattle plus genome-wide single nucleotide polymorphism data from a diverse panel of 1,225 modern animals.

Results

Phylogenomic analyses place the aurochs as a distinct outgroup to the domestic B. taurus lineage, supporting the predominant Near Eastern origin of European cattle. Conversely, traditional British and Irish breeds share more genetic variants with this aurochs specimen than other European populations, supporting localized gene flow from aurochs into the ancestors of modern British and Irish cattle, perhaps through purposeful restocking by early herders in Britain. Finally, the functions of genes showing evidence for positive selection in B. taurus are enriched for neurobiology, growth, metabolism and immunobiology, suggesting that these biological processes have been important in the domestication of cattle.

Conclusions

This work provides important new information regarding the origins and functional evolution of modern cattle, revealing that the interface between early European domestic populations and wild aurochs was significantly more complex than previously thought.

Electronic supplementary material

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

Mycobacterial lineages causing pulmonary and extrapulmonary tuberculosis, Ethiopia

Molecular typing of 964 specimens from patients in Ethiopia with lymph node or pulmonary tuberculosis showed a similar distribution of Mycobacterium tuberculosis strains between the 2 disease manifestations and a minimal role for M. bovis. We report a novel phylogenetic lineage of M. tuberculosis strongly associated with the Horn of Africa.

Control of reproductive floral organ identity specification in Arabidopsis by the C function regulator AGAMOUS

The floral organ identity factor AGAMOUS (AG) is a key regulator of Arabidopsis thaliana flower development, where it is involved in the formation of the reproductive floral organs as well as in the control of meristem determinacy. To obtain insights into how AG specifies organ fate, we determined the genes and processes acting downstream of this C function regulator during early flower development and distinguished between direct and indirect effects. To this end, we combined genome-wide localization studies, gene perturbation experiments, and computational analyses. Our results demonstrate that AG controls flower development to a large extent by controlling the expression of other genes with regulatory functions, which are involved in mediating a plethora of different developmental processes. One aspect of this function is the suppression of the leaf development program in emerging floral primordia. Using trichome initiation as an example, we demonstrate that AG inhibits an important aspect of leaf development through the direct control of key regulatory genes. A comparison of the gene expression programs controlled by AG and the B function regulators APETALA3 and PISTILLATA, respectively, showed that while they control many developmental processes in conjunction, they also have marked antagonistic, as well as independent activities.

Whole-transcriptome, high-throughput RNA sequence analysis of the bovine macrophage response to Mycobacterium bovis infection in vitro

Background

Mycobacterium bovis, the causative agent of bovine tuberculosis, is an intracellular pathogen that can persist inside host macrophages during infection via a diverse range of mechanisms that subvert the host immune response. In the current study, we have analysed and compared the transcriptomes of M. bovis-infected monocyte-derived macrophages (MDM) purified from six Holstein-Friesian females with the transcriptomes of non-infected control MDM from the same animals over a 24 h period using strand-specific RNA sequencing (RNA-seq). In addition, we compare gene expression profiles generated using RNA-seq with those previously generated by us using the high-density Affymetrix® GeneChip® Bovine Genome Array platform from the same MDM-extracted RNA.

Results

A mean of 7.2 million reads from each MDM sample mapped uniquely and unambiguously to single Bos taurus reference genome locations. Analysis of these mapped reads showed 2,584 genes (1,392 upregulated; 1,192 downregulated) and 757 putative natural antisense transcripts (558 upregulated; 119 downregulated) that were differentially expressed based on sense and antisense strand data, respectively (adjusted P-value ≤ 0.05). Of the differentially expressed genes, 694 were common to both the sense and antisense data sets, with the direction of expression (i.e. up- or downregulation) positively correlated for 693 genes and negatively correlated for the remaining gene. Gene ontology analysis of the differentially expressed genes revealed an enrichment of immune, apoptotic and cell signalling genes. Notably, the number of differentially expressed genes identified from RNA-seq sense strand analysis was greater than the number of differentially expressed genes detected from microarray analysis (2,584 genes versus 2,015 genes). Furthermore, our data reveal a greater dynamic range in the detection and quantification of gene transcripts for RNA-seq compared to microarray technology.

Conclusions

This study highlights the value of RNA-seq in identifying novel immunomodulatory mechanisms that underlie host-mycobacterial pathogen interactions during infection, including possible complex post-transcriptional regulation of host gene expression involving antisense RNA.

Genome of Acanthamoeba castellanii highlights extensive lateral gene transfer and early evolution of tyrosine kinase signaling

BACKGROUND

The Amoebozoa constitute one of the primary divisions of eukaryotes, encompassing taxa of both biomedical and evolutionary importance, yet its genomic diversity remains largely unsampled. Here we present an analysis of a whole genome assembly of Acanthamoeba castellanii (Ac) the first representative from a solitary free-living amoebozoan.

RESULTS

Ac encodes 15,455 compact intron-rich genes, a significant number of which are predicted to have arisen through inter-kingdom lateral gene transfer (LGT). A majority of the LGT candidates have undergone a substantial degree of intronization and Ac appears to have incorporated them into established transcriptional programs. Ac manifests a complex signaling and cell communication repertoire, including a complete tyrosine kinase signaling toolkit and a comparable diversity of predicted extracellular receptors to that found in the facultatively multicellular dictyostelids. An important environmental host of a diverse range of bacteria and viruses, Ac utilizes a diverse repertoire of predicted pattern recognition receptors, many with predicted orthologous functions in the innate immune systems of higher organisms.

CONCLUSIONS

Our analysis highlights the important role of LGT in the biology of Ac and in the diversification of microbial eukaryotes. The early evolution of a key signaling facility implicated in the evolution of metazoan multicellularity strongly argues for its emergence early in the Unikont lineage. Overall, the availability of an Ac genome should aid in deciphering the biology of the Amoebozoa and facilitate functional genomic studies in this important model organism and environmental host.

The transcriptional landscape and small RNAs of Salmonella enterica serovar Typhimurium

More than 50 y of research have provided great insight into the physiology, metabolism, and molecular biology of Salmonella enterica serovar Typhimurium (S. Typhimurium), but important gaps in our knowledge remain. It is clear that a precise choreography of gene expression is required for Salmonella infection, but basic genetic information such as the global locations of transcription start sites (TSSs) has been lacking. We combined three RNA-sequencing techniques and two sequencing platforms to generate a robust picture of transcription in S. Typhimurium. Differential RNA sequencing identified 1,873 TSSs on the chromosome of S. Typhimurium SL1344 and 13% of these TSSs initiated antisense transcripts. Unique findings include the TSSs of the virulence regulators phoP, slyA, and invF. Chromatin immunoprecipitation revealed that RNA polymerase was bound to 70% of the TSSs, and two-thirds of these TSSs were associated with σ(70) (including phoP, slyA, and invF) from which we identified the -10 and -35 motifs of σ(70)-dependent S. Typhimurium gene promoters. Overall, we corrected the location of important genes and discovered 18 times more promoters than identified previously. S. Typhimurium expresses 140 small regulatory RNAs (sRNAs) at early stationary phase, including 60 newly identified sRNAs. Almost half of the experimentally verified sRNAs were found to be unique to the Salmonella genus, and <20% were found throughout the Enterobacteriaceae. This description of the transcriptional map of SL1344 advances our understanding of S. Typhimurium, arguably the most important bacterial infection model.

Methicillin resistance alters the biofilm phenotype and attenuates virulence in Staphylococcus aureus device-associated infections

Clinical isolates of Staphylococcus aureus can express biofilm phenotypes promoted by the major cell wall autolysin and the fibronectin-binding proteins or the icaADBC-encoded polysaccharide intercellular adhesin/poly-N-acetylglucosamine (PIA/PNAG). Biofilm production in methicillin-susceptible S. aureus (MSSA) strains is typically dependent on PIA/PNAG whereas methicillin-resistant isolates express an Atl/FnBP-mediated biofilm phenotype suggesting a relationship between susceptibility to β-lactam antibiotics and biofilm. By introducing the methicillin resistance gene mecA into the PNAG-producing laboratory strain 8325-4 we generated a heterogeneously resistant (HeR) strain, from which a homogeneous, high-level resistant (HoR) derivative was isolated following exposure to oxacillin. The HoR phenotype was associated with a R₆₀₂H substitution in the DHHA1 domain of GdpP, a recently identified c-di-AMP phosphodiesterase with roles in resistance/tolerance to β-lactam antibiotics and cell envelope stress. Transcription of icaADBC and PNAG production were impaired in the 8325-4 HoR derivative, which instead produced a proteinaceous biofilm that was significantly inhibited by antibodies against the mecA-encoded penicillin binding protein 2a (PBP2a). Conversely excision of the SCCmec element in the MRSA strain BH1CC resulted in oxacillin susceptibility and reduced biofilm production, both of which were complemented by mecA alone. Transcriptional activity of the accessory gene regulator locus was also repressed in the 8325-4 HoR strain, which in turn was accompanied by reduced protease production and significantly reduced virulence in a mouse model of device infection. Thus, homogeneous methicillin resistance has the potential to affect agr- and icaADBC-mediated phenotypes, including altered biofilm expression and virulence, which together are consistent with the adaptation of healthcare-associated MRSA strains to the antibiotic-rich hospital environment in which they are frequently responsible for device-related infections in immuno-compromised patients.

The acquisition of mecA, which encodes penicillin binding protein 2a (PBP2a) and methicillin resistance, by Staphylococcus aureus has added to an already impressive array of virulence mechanisms including enzyme and toxin production, biofilm forming capacity and immune evasion. And yet clinical data does not indicate that healthcare-associated methicillin resistant S. aureus (MRSA) strains are more virulent than their methicillin-susceptible counterparts. Here our findings suggest that MRSA sacrifices virulence potential for antibiotic resistance and that expression of methicillin resistance alters the biofilm phenotype but does not interfere with the colonization of implanted medical devices in vivo. High level expression of PBP2a, which was associated with a mutation in the c-di-AMP phosphodiesterase gene gdpP, resulted in these pleiotrophic effects by blocking icaADBC-dependent polysaccharide type biofilm development and promoting an alternative PBP2a-mediated biofilm, repressing the accessory gene regulator and extracellular protease production, and attenuating virulence in a mouse device-infection model. Thus the adaptation of MRSA to the hospital environment has apparently focused on the acquisition of antibiotic resistance and retention of biofilm forming capacity, which are likely to be more advantageous than metabolically-expensive enzyme and toxin production in immunocompromised patients with implanted medical devices offering a route to infection.

Sequencing illustrates the transcriptional response of Legionella pneumophila during infection and identifies seventy novel small non-coding RNAs

Second generation sequencing has prompted a number of groups to re-interrogate the transcriptomes of several bacterial and archaeal species. One of the central findings has been the identification of complex networks of small non-coding RNAs that play central roles in transcriptional regulation in all growth conditions and for the pathogen's interaction with and survival within host cells. Legionella pneumophila is a Gram-negative facultative intracellular human pathogen with a distinct biphasic lifestyle. One of its primary environmental hosts in the free-living amoeba Acanthamoeba castellanii and its infection by L. pneumophila mimics that seen in human macrophages. Here we present analysis of strand specific sequencing of the transcriptional response of L. pneumophila during exponential and post-exponential broth growth and during the replicative and transmissive phase of infection inside A. castellanii. We extend previous microarray based studies as well as uncovering evidence of a complex regulatory architecture underpinned by numerous non-coding RNAs. Over seventy new non-coding RNAs could be identified; many of them appear to be strain specific and in configurations not previously reported. We discover a family of non-coding RNAs preferentially expressed during infection conditions and identify a second copy of 6S RNA in L. pneumophila. We show that the newly discovered putative 6S RNA as well as a number of other non-coding RNAs show evidence for antisense transcription. The nature and extent of the non-coding RNAs and their expression patterns suggests that these may well play central roles in the regulation of Legionella spp. specific traits and offer clues as to how L. pneumophila adapts to its intracellular niche. The expression profiles outlined in the study have been deposited into Genbank's Gene Expression Omnibus (GEO) database under the series accession GSE27232.

Sequencing and analysis of an Irish human genome

Background

Recent studies generating complete human sequences from Asian, African and European subgroups have revealed population-specific variation and disease susceptibility loci. Here, choosing a DNA sample from a population of interest due to its relative geographical isolation and genetic impact on further populations, we extend the above studies through the generation of 11-fold coverage of the first Irish human genome sequence.

Results

Using sequence data from a branch of the European ancestral tree as yet unsequenced, we identify variants that may be specific to this population. Through comparisons with HapMap and previous genetic association studies, we identified novel disease-associated variants, including a novel nonsense variant putatively associated with inflammatory bowel disease. We describe a novel method for improving SNP calling accuracy at low genome coverage using haplotype information. This analysis has implications for future re-sequencing studies and validates the imputation of Irish haplotypes using data from the current Human Genome Diversity Cell Line Panel (HGDP-CEPH). Finally, we identify gene duplication events as constituting significant targets of recent positive selection in the human lineage.

Conclusions

Our findings show that there remains utility in generating whole genome sequences to illustrate both general principles and reveal specific instances of human biology. With increasing access to low cost sequencing we would predict that even armed with the resources of a small research group a number of similar initiatives geared towards answering specific biological questions will emerge.

A complete mitochondrial genome sequence from a mesolithic wild aurochs (Bos primigenius)

BACKGROUND

The derivation of domestic cattle from the extinct wild aurochs (Bos primigenius) has been well-documented by archaeological and genetic studies. Genetic studies point towards the Neolithic Near East as the centre of origin for Bos taurus, with some lines of evidence suggesting possible, albeit rare, genetic contributions from locally domesticated wild aurochsen across Eurasia. Inferences from these investigations have been based largely on the analysis of partial mitochondrial DNA sequences generated from modern animals, with limited sequence data from ancient aurochsen samples. Recent developments in DNA sequencing technologies, however, are affording new opportunities for the examination of genetic material retrieved from extinct species, providing new insight into their evolutionary history. Here we present DNA sequence analysis of the first complete mitochondrial genome (16,338 base pairs) from an archaeologically-verified and exceptionally-well preserved aurochs bone sample.

METHODOLOGY

DNA extracts were generated from an aurochs humerus bone sample recovered from a cave site located in Derbyshire, England and radiocarbon-dated to 6,738+/-68 calibrated years before present. These extracts were prepared for both Sanger and next generation DNA sequencing technologies (Illumina Genome Analyzer). In total, 289.9 megabases (22.48%) of the post-filtered DNA sequences generated using the Illumina Genome Analyzer from this sample mapped with confidence to the bovine genome. A consensus B. primigenius mitochondrial genome sequence was constructed and was analysed alongside all available complete bovine mitochondrial genome sequences.

CONCLUSIONS

For all nucleotide positions where both Sanger and Illumina Genome Analyzer sequencing methods gave high-confidence calls, no discrepancies were observed. Sequence analysis reveals evidence of heteroplasmy in this sample and places this mitochondrial genome sequence securely within a previously identified aurochsen haplogroup (haplogroup P), thus providing novel insights into pre-domestic patterns of variation. The high proportion of authentic, endogenous aurochs DNA preserved in this sample bodes well for future efforts to determine the complete genome sequence of a wild ancestor of domestic cattle.

Analysis of 5'- or 3'-terminal tRNA editing: mitochondrial 5' tRNA editing in Acanthamoeba castellanii as the exemplar

Editing processes that result in the structural retailoring of the aminoacyl acceptor stems of mitochondrial tRNAs are the focus of this chapter. This type of tRNA editing is the most frequently observed and widely distributed and involves nucleotide replacement within the 5' or 3' half of the aminoacyl acceptor stem in either a template-directed or a template-independent fashion. We provide a detailed protocol that allows demarcation of 5'-terminal tRNA editing events from those occurring on the 3' side of the acceptor stem. We present the mitochondrial 5' tRNA editing system in Acanthamoeba castellanii as the exemplar of terminal tRNA editing. The methodology involves RNA ligase-mediated circularization of tRNAs, cDNA synthesis primed by tRNA-specific oligonucleotides, amplification of cDNA via polymerase chain reaction, and cloning and sequencing of multiple products. This approach permits (1) simultaneous determination of 5' and 3' acceptor stem sequences, (2) delineation of 5' or 3' editing, (3) identification of mature tRNAs (characterized by having a 3'-CCA(OH) motif), (4) identification of processing/editing intermediates, and (5) mechanistic insights.

Methods for analysis of mitochondrial tRNA editing in Acanthamoeba castellanii

We present here methodology for assaying 5'-terminal editing of mitochondrial tRNAs in the amoeboid protist Acanthamoeba castellanii. This type of editing involves replacement of one or more nucleotides within the first three positions at the 5' end of a tRNA substrate. The assay procedure involves RNA ligase-mediated joining of the 5' and 3' ends of a tRNA, use of the resulting circularized tRNA as template for cDNA synthesis primed by tRNA-specific primers over a region that encompasses the ligated 5' and 3' halves of the acceptor stem, amplification of cDNA via polymerase chain reaction, and cloning and sequencing of double-stranded cDNA product. This approach has the advantage that it simultaneously reveals potential editing events on the 5' and 3' side of an acceptor stem, as well as serves to identify mature tRNAs (characterized by having a 3'-CCAOH motif) and partially processed intermediates.

tRNAHis maturation: an essential yeast protein catalyzes addition of a guanine nucleotide to the 5' end of tRNAHis

All tRNAHis molecules are unusual in having an extra 5' GMP residue (G(-1)) that, in eukaryotes, is added after transcription and RNase P cleavage. Incorporation of this G(-1) residue is a rare example of nucleotide addition occurring at an RNA 5' end in a normal phosphodiester linkage. We show here that the essential Saccharomyces cerevisiae ORF YGR024c (THG1) is responsible for this guanylyltransferase reaction. Thg1p was identified by survey of a genomic collection of yeast GST-ORF fusion proteins for addition of [alpha-32P]GTP to tRNAHis. End analysis confirms the presence of G(-1). Thg1p is required for tRNAHis guanylylation in vivo, because cells depleted of Thg1p lack G(-1) in their tRNAHis. His6-Thg1p purified from Escherichia coli catalyzes the guanylyltransferase step of G(-1) addition using a ppp-tRNAHis substrate, and appears to catalyze the activation step using p-tRNAHis and ATP. Thg1p is highlye conserved in eukaryotes, where G(-1) addition is necessary, and is not found in eubacteria, where G(-1) is genome-encoded. Thus, Thg1p is the first member of a new family of enzymes that can catalyze phosphodiester bond formation at the 5' end of RNAs, formally in a 3'-5' direction. Surprisingly, despite its varied activities, Thg1p contains no recognizable catalytic or functional domains.

A subset of conserved tRNA genes in plastid DNA of nongreen plants

The plastid genome of the nonphotosynthetic parasitic plant Epifagus virginiana contains only 17 of the 30 tRNA genes normally found in angiosperm plastid DNA. Although this is insufficient for translation, the genome is functional, so import of cytosolic tRNAs into plastids has been suggested. This raises the question of whether the tRNA genes that remain in E. virginiana plastid DNA are active or have just fortuitously escaped deletion. We report the sequences of 20 plastid tRNA loci from Orobanche minor, which shares a nonphotosynthetic ancestor with E. virginiana. The two species have 9 intact tRNA genes in common, the others being defunct in one or both species. The intron-containing trnLUAA gene is absent from E. virginiana, but it is intact, transcribed, and spliced in O. minor. The shared intact genes are better conserved than intergenic sequences, which indicates that these genes are being maintained by natural selection and, therefore, must be functional. For the most part, the tRNA species conserved in nonphotosynthetic plastids are also those that have never been found to be imported in plant mitochondria, which suggests that the same rules may govern tRNA import in the two organelles. A small photosynthesis gene, psbI, is still intact in O. minor, and computer simulations show that some small nonessential genes have an appreciable chance of escaping deletion.

Complete DNA sequence of yeast chromosome II

In the framework of the EU genome-sequencing programmes, the complete DNA sequence of the yeast Saccharomyces cerevisiae chromosome II (807 188 bp) has been determined. At present, this is the largest eukaryotic chromosome entirely sequenced. A total of 410 open reading frames (ORFs) were identified, covering 72% of the sequence. Similarity searches revealed that 124 ORFs (30%) correspond to genes of known function, 51 ORFs (12.5%) appear to be homologues of genes whose functions are known, 52 others (12.5%) have homologues the functions of which are not well defined and another 33 of the novel putative genes (8%) exhibit a degree of similarity which is insufficient to confidently assign function. Of the genes on chromosome II, 37-45% are thus of unpredicted function. Among the novel putative genes, we found several that are related to genes that perform differentiated functions in multicellular organisms of are involved in malignancy. In addition to a compact arrangement of potential protein coding sequences, the analysis of this chromosome confirmed general chromosome patterns but also revealed particular novel features of chromosomal organization. Alternating regional variations in average base composition correlate with variations in local gene density along chromosome II, as observed in chromosomes XI and III. We propose that functional ARS elements are preferably located in the AT-rich regions that have a spacing of approximately 110 kb. Similarly, the 13 tRNA genes and the three Ty elements of chromosome II are found in AT-rich regions. In chromosome II, the distribution of coding sequences between the two strands is biased, with a ratio of 1.3:1. An interesting aspect regarding the evolution of the eukaryotic genome is the finding that chromosome II has a high degree of internal genetic redundancy, amounting to 16% of the coding capacity.

Sequence around the centromere of Saccharomyces cerevisiae chromosome II: similarity of CEN2 to CEN4

We report the sequence of a 12 kilobase region spanning the centromere of Saccharomyces cerevisiae chromosome II. The sequence from the left arm includes genes for histones H2A and H2B. The sequence from the right arm includes a gene that probably encodes a novel trehalase, as well as the COQ1 gene (for an enzyme involved in coenzyme Q biosynthesis), and an open reading frame with significant similarity to bacterial genes of unknown function. The trehalase gene (YBR0106) on chromosome II is located beside the centromere and transcribed towards it. This is identical to the arrangement of the neutral trehalase gene (NTH1) beside the centromere of chromosome IV. The centromere regions of chromosomes II and IV may therefore have arisen through a duplication of the centromere region of an ancestral chromosome. The YBR0106 and NTH1 proteins are 77% identical in predicted amino acid sequence, but there is no pronounced sequence similarity between the two centromeres (CEN2 and CEN4) outside of the universally conserved CDE I and CDE III elements. The genes flanking the centromere and trehalase genes differ between the two chromosomes, so the similarity between chromosomes II and IV may be less extensive than that recently reported between chromosomes III and XIV.