The Draft Assembly of the Radically Organized Stylonychia lemnae Macronuclear Genome.

Stylonychia lemnae is a classical model single-celled eukaryote, and a quintessential ciliate typified by dimorphic nuclei: A small, germline micronucleus and a massive, vegetative macronucleus. The genome within Stylonychia's macronucleus has a very unusual architecture, comprised variably and highly amplified "nanochromosomes," each usually encoding a single gene with a minimal amount of surrounding noncoding DNA. As only a tiny fraction of the Stylonychia genes has been sequenced, and to promote research using this organism, we sequenced its macronuclear genome. We report the analysis of the 50.2-Mb draft S. lemnae macronuclear genome assembly, containing in excess of 16,000 complete nanochromosomes, assembled as less than 20,000 contigs. We found considerable conservation of fundamental genomic properties between S. lemnae and its close relative, Oxytricha trifallax, including nanochromosomal gene synteny, alternative fragmentation, and copy number. Protein domain searches in Stylonychia revealed two new telomere-binding protein homologs and the presence of linker histones. Among the diverse histone variants of S. lemnae and O. trifallax, we found divergent, coexpressed variants corresponding to four of the five core nucleosomal proteins (H1.2, H2A.6, H2B.4, and H3.7) suggesting that these ciliates may possess specialized nucleosomes involved in genome processing during nuclear differentiation. The assembly of the S. lemnae macronuclear genome demonstrates that largely complete, well-assembled highly fragmented genomes of similar size and complexity may be produced from one library and lane of Illumina HiSeq 2000 shotgun sequencing. The provision of the S. lemnae macronuclear genome sets the stage for future detailed experimental studies of chromatin-mediated, RNA-guided developmental genome rearrangements.

Genome-wide analysis of genetic and epigenetic control of programmed DNA deletion

During the development of the somatic genome from the Paramecium germline genome the bulk of the copies of ∼45 000 unique, internal eliminated sequences (IESs) are deleted. IES targeting is facilitated by two small RNA (sRNA) classes: scnRNAs, which relay epigenetic information from the parental nucleus to the developing nucleus, and iesRNAs, which are produced and used in the developing nucleus. Why only certain IESs require sRNAs for their removal has been enigmatic. By analyzing the silencing effects of three genes: PGM (responsible for DNA excision), DCL2/3 (scnRNA production) and DCL5 (iesRNA production), we identify key properties required for IES elimination. Based on these results, we propose that, depending on the exact combination of their lengths and end bases, some IESs are less efficiently recognized or excised and have a greater requirement for targeting by scnRNAs and iesRNAs. We suggest that the variation in IES retention following silencing of DCL2/3 is not primarily due to scnRNA density, which is comparatively uniform relative to IES retention, but rather the genetic properties of IESs. Taken together, our analyses demonstrate that in Paramecium the underlying genetic properties of developmentally deleted DNA sequences are essential in determining the sensitivity of these sequences to epigenetic control.

Relationship between genome and epigenome - challenges and requirements for future research

Understanding the links between genetic, epigenetic and non-genetic factors throughout the lifespan and across generations and their role in disease susceptibility and disease progression offer entirely new avenues and solutions to major problems in our society. To overcome the numerous challenges, we have come up with nine major conclusions to set the vision for future policies and research agendas at the European level.

Whole genome scan identifies several chromosomal regions linked to equine sarcoids.

Despite the evidence for a genetic predisposition to develop equine sarcoids (ES), no whole genome scan for ES has been performed to date. The objective of this explorative study was to identify chromosome regions associated with ES. The studied population was comprised of two half-sibling sire families, involving a total of 222 horses. Twenty-six of these horses were affected with ES. All horses had been previously genotyped with 315 microsatellite markers. Quantitative trait locus (QTL) signals were suggested where the F statistic exceeded chromosome-wide significance at P < 0.05. The QTL analyses revealed significant signals reaching P < 0.05 on equine chromosome (ECA) 20, 23 and 25, suggesting a polygenic character for this trait. The candidate regions identified on ECA 20, 23 and 25 include genes regulating virus replication and host immune response. Further investigation of the chromosome regions associated with ES and of genes potentially responsible for the development of ES could form the basis for early identification of susceptible animals, breeding selection or the development of new therapeutic targets.

Impaired genome encapsidation restricts the in vitro propagation of human parvovirus B19.

The lack of a permissive cell culture system hampers the study of human parvovirus B19 (B19V). UT7/Epo is one of the few established cell lines that can be infected with B19V but generates none or few infectious progeny. Recently, hypoxic conditions or the use of primary CD36+ erythroid progenitor cells (CD36+ EPCs) have been shown to improve the infection. These novel approaches were evaluated in infection and transfection experiments. Hypoxic conditions or the use of CD36+ EPCs resulted in a significant acceleration of the infection/transfection and a modest increase in the yield of capsid progeny. However, under all tested conditions, genome encapsidation was impaired seriously. Further analysis of the cell culture virus progeny revealed that differently to the wild-type virus, the VP1 unique region (VP1u) was exposed partially and was unable to become further externalized upon heat treatment. The fivefold axes pore, which is used for VP1u externalization and genome encapsidation, might be constricted by the atypical VP1u conformation explaining the packaging failure. Although CD36+ EPCs and hypoxia facilitate B19V infection, large quantities of infectious progeny cannot be generated due to a failure in genome encapsidation, which arises as a major limiting factor for the in vitro propagation of B19V.

Tracking a tuberculosis outbreak over 21 years: strain-specific single nucleotide polymorphism-typing combined with targeted whole genome sequencing.

BACKGROUND  Whole genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single nucleotide polymorphism (SNP)-typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland. METHODS  Based on genome sequences of three historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1,642 patient isolates, and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected. RESULTS  We identified 68 patients associated with the outbreak strain. Most were diagnosed in 1991-1995, but cases were observed until 2011. Two thirds belonged to the homeless and substance abuser milieu. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into three sub-clusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs. CONCLUSIONS  Strain-specific SNP-genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real-time and at high resolution.

Fine-scale population structure analysis of seven local Swiss sheep breeds using genome-wide SNP data

As part of the global sheep Hapmap project, 24 individuals from each of seven indigenous Swiss sheep breeds (Bundner Oberländer sheep (BOS), Engadine Red sheep (ERS), Swiss Black-Brown Mountain sheep (SBS), Swiss Mirror sheep (SMS), Swiss White Alpine (SWA) sheep, Valais Blacknose sheep (VBS) and Valais Red sheep (VRS)), were genotyped using Illumina’s Ovine SNP50 BeadChip. In total, 167 animals were subjected to a detailed analysis for genetic diversity using 45 193 informative single nucleotide polymorphisms. The results of the phylogenetic analyses supported the known proximity between populations such as VBS and VRS or SMS and SWA. Average genomic relatedness within a breed was found to be 12 percent (BOS), 5 percent (ERS), 9 percent (SBS), 10 percent (SMS), 9 percent (SWA), 12 percent (VBS) and 20 percent (VRS). Furthermore, genomic relationships between breeds were found for single individuals from SWA and SMS, VRS and VBS as well as VRS and BOS. In addition, seven out of 40 indicated parent–offspring pairs could not be confirmed. These results were further supported by results from the genome-wide population cluster analysis. This study provides a better understanding of fine-scale population structures within and between Swiss sheep breeds. This relevant information will help to increase the conservation activities of the local Swiss sheep breeds.

Complete Genome Sequences of Virulent Mycoplasma capricolum subsp. capripneumoniae Strains F38 and ILRI181

Contagious caprine pleuropneumonia (CCPP) caused by Mycoplasma capricolum subsp. capripneumoniae is a severe epidemic affecting mainly domestic Caprinae species but also affects wild Caprinae species. M. capricolum subsp. capripneumoniae belongs to the "Mycoplasma mycoides cluster." The disease features prominently in East Africa, in particular Kenya, Tanzania, and Ethiopia. CCPP also endangers wildlife and thus affects not only basic nutritional resources of large populations but also expensively built-up game resorts in affected countries. Here, we report the complete sequences of two M. capricolum subsp. capripneumoniae strains: the type strain F38 and strain ILRI181 isolated druing a recent outbreak in Kenya. Both genomes have a G+C content of 24% with sizes of 1,016,760 bp and 1,017,183 bp for strains F38 and ILRI181, respectively.

Genome-wide identification of pathogenicity factors of the free-living amoeba Naegleria fowleri.

BACKGROUND The free-living amoeba Naegleria fowleri is the causative agent of the rapidly progressing and typically fatal primary amoebic meningoencephalitis (PAM) in humans. Despite the devastating nature of this disease, which results in > 97% mortality, knowledge of the pathogenic mechanisms of the amoeba is incomplete. This work presents a comparative proteomic approach based on an experimental model in which the pathogenic potential of N. fowleri trophozoites is influenced by the compositions of different media. RESULTS As a scaffold for proteomic analysis, we sequenced the genome and transcriptome of N. fowleri. Since the sequence similarity of the recently published genome of Naegleria gruberi was far lower than the close taxonomic relationship of these species would suggest, a de novo sequencing approach was chosen. After excluding cell regulatory mechanisms originating from different media compositions, we identified 22 proteins with a potential role in the pathogenesis of PAM. Functional annotation of these proteins revealed, that the membrane is the major location where the amoeba exerts its pathogenic potential, possibly involving actin-dependent processes such as intracellular trafficking via vesicles. CONCLUSION This study describes for the first time the 30 Mb-genome and the transcriptome sequence of N. fowleri and provides the basis for the further definition of effective intervention strategies against the rare but highly fatal form of amoebic meningoencephalitis.

Whole genome RNAi screens reveal a critical role of REV3 in coping with replication stress

REV3, the catalytic subunit of translesion polymerase zeta (polζ), is commonly associated with DNA damage bypass and repair. Despite sharing accessory subunits with replicative polymerase δ, very little is known about the role of polζ in DNA replication. We previously demonstrated that inhibition of REV3 expression induces persistent DNA damage and growth arrest in cancer cells. To reveal determinants of this sensitivity and obtain insights into the cellular function of REV3, we performed whole human genome RNAi library screens aimed at identification of synthetic lethal interactions with REV3 in A549 lung cancer cells. The top confirmed hit was RRM1, the large subunit of ribonucleotide reductase (RNR), a critical enzyme of de novo nucleotide synthesis. Treatment with the RNR-inhibitor hydroxyurea (HU) synergistically increased the fraction of REV3-deficient cells containing single stranded DNA (ssDNA) as indicated by an increase in replication protein A (RPA). However, this increase was not accompanied by accumulation of the DNA damage marker γH2AX suggesting a role of REV3 in counteracting HU-induced replication stress (RS). Consistent with a role of REV3 in DNA replication, increased RPA staining was confined to HU-treated S-phase cells. Additionally, we found genes related to RS to be significantly enriched among the top hits of the synthetic sickness/lethality (SSL) screen further corroborating the importance of REV3 for DNA replication under conditions of RS.

Evaluation of genome-wide loci of iron metabolism in hereditary hemochromatosis identifies PCSK7 as a host risk factor of liver cirrhosis

Genome-wide association studies (GWAS) have revealed genetic determinants of iron metabolism, but correlation of these with clinical phenotypes is pending. Homozygosity for HFE C282Y is the predominant genetic risk factor for hereditary hemochromatosis (HH) and may cause liver cirrhosis. However, this genotype has a low penetrance. Thus, detection of yet unknown genetic markers that identify patients at risk of developing severe liver disease is necessary for better prevention. Genetic loci associated with iron metabolism (TF, TMPRSS6, PCSK7, TFR2 and Chr2p14) in recent GWAS and liver fibrosis (PNPLA3) in recent meta-analysis were analyzed for association with either liver cirrhosis or advanced fibrosis in 148 German HFE C282Y homozygotes. Replication of associations was sought in additional 499 Austrian/Swiss and 112 HFE C282Y homozygotes from Sweden. Only variant rs236918 in the PCSK7 gene (proprotein convertase subtilisin/kexin type 7) was associated with cirrhosis or advanced fibrosis (P = 1.02 × 10(-5)) in the German cohort with genotypic odds ratios of 3.56 (95% CI 1.29-9.77) for CG heterozygotes and 5.38 (95% CI 2.39-12.10) for C allele carriers. Association between rs236918 and cirrhosis was confirmed in Austrian/Swiss HFE C282Y homozygotes (P = 0.014; ORallelic = 1.82 (95% CI 1.12-2.95) but not in Swedish patients. Post hoc combined analyses of German/Swiss/Austrian patients with available liver histology (N = 244, P = 0.00014, ORallelic = 2.84) and of males only (N = 431, P = 2.17 × 10(-5), ORallelic = 2.54) were consistent with the premier finding. Association between rs236918 and cirrhosis was not confirmed in alcoholic cirrhotics, suggesting specificity of this genetic risk factor for HH. PCSK7 variant rs236918 is a risk factor for cirrhosis in HH patients homozygous for the HFE C282Y mutation.

Complete Genome Sequence of Bovine Pestivirus Strain PG-2, a Second Member of the Tentative Pestivirus Species Giraffe

We report the complete genome sequence of bovine pestivirus strain PG-2. The sequence data from this virus showed that PG-2 is closely related to the giraffe pestivirus strain H138. PG-2 and H138 belong to one pestivirus species that should be considered an approved member of the genus Pestivirus.

A Genome-Wide Association Study Reveals Loci Influencing Height and Other Conformation Traits in Horses

The molecular analysis of genes influencing human height has been notoriously difficult. Genome-wide association studies (GWAS) for height in humans based on tens of thousands to hundreds of thousands of samples so far revealed ∼200 loci for human height explaining only 20% of the heritability. In domestic animals isolated populations with a greatly reduced genetic heterogeneity facilitate a more efficient analysis of complex traits. We performed a genome-wide association study on 1,077 Franches-Montagnes (FM) horses using ∼40,000 SNPs. Our study revealed two QTL for height at withers on chromosomes 3 and 9. The association signal on chromosome 3 is close to the LCORL/NCAPG genes. The association signal on chromosome 9 is close to the ZFAT gene. Both loci have already been shown to influence height in humans. Interestingly, there are very large intergenic regions at the association signals. The two detected QTL together explain ∼18.2% of the heritable variation of height in horses. However, another large fraction of the variance for height in horses results from ECA 1 (11.0%), although the association analysis did not reveal significantly associated SNPs on this chromosome. The QTL region on ECA 3 associated with height at withers was also significantly associated with wither height, conformation of legs, ventral border of mandible, correctness of gaits, and expression of the head. The region on ECA 9 associated with height at withers was also associated with wither height, length of croup and length of back. In addition to these two QTL regions on ECA 3 and ECA 9 we detected another QTL on ECA 6 for correctness of gaits. Our study highlights the value of domestic animal populations for the genetic analysis of complex traits.