Explore the vast range of titles available.

Leif Andersson and colleagues report the genome sequence of the ruff, a bird species with three male morphs with different reproductive strategies. Satellite and faeder morphs differ from the common independent morph by a 4.5-Mb inversion that occurred approximately 3.8 million years ago, and multiple genetic changes within this inverted region are associated with the satellite and faeder morphs. The ruff is a Palearctic wader with a spectacular lekking behavior where highly ornamented males compete for females 1 , 2 , 3 , 4 . This bird has one of the most remarkable mating systems in the animal kingdom, comprising three different male morphs (independents, satellites and faeders) that differ in behavior, plumage color and body size. Remarkably, the satellite and faeder morphs are controlled by dominant alleles 5 , 6 . Here we have used whole-genome sequencing and resolved the enigma of how such complex phenotypic differences can have a simple genetic basis. The Satellite and Faeder alleles are both associated with a 4.5-Mb inversion that occurred about 3.8 million years ago. We propose an evolutionary scenario where the Satellite chromosome arose by a rare recombination event about 500,000 years ago. The ruff mating system is the result of an evolutionary process in which multiple genetic changes contributing to phenotypic differences between morphs have accumulated within the inverted region.

The extent to which disturbances in the resident microbiota can compromise an animal's health is poorly understood. Hydra is one of the evolutionary oldest animals with naturally occurring tumors. Here, we found a causal relationship between an environmental spirochete (Turneriella spec.) and tumorigenesis in Hydra. Unexpectedly, virulence of this pathogen requires the presence of Pseudomonas spec., a member of Hydra´s beneficial microbiome indicating that dynamic interactions between a resident bacterium and a pathogen cause tumor formation. The observation points to the crucial role of commensal bacteria in maintaining tissue homeostasis and adds support to the view that microbial community interactions are essential for disease. These findings in an organism that shares deep evolutionary connections with all animals have implications for our understanding of cancer.

The evolution of asexual organisms is driven not only by the inheritance of genetic modification but also by the acquisition of foreign DNA. The contribution of vertical and horizontal processes to genome evolution depends on their rates per year and is quantified by the ratio of recombination to mutation. These rates have been estimated for bacteria; however, no estimates have been reported for phages. Here, we delineate the contribution of mutation and recombination to dsDNA phage genome evolution. We analyzed 34 isolates of the 936 group of Siphoviridae phages using a Lactococcus lactis strain from a single dairy over 29 years. We estimate a constant substitution rate of 1.9 × 10−4 substitutions per site per year due to mutation that is within the range of estimates for eukaryotic RNA and DNA viruses. The reconstruction of recombination events reveals a constant rate of five recombination events per year and 4.5 × 10−3 nucleotide alterations due to recombination per site per year. Thus, the recombination rate exceeds the substitution rate, resulting in a relative effect of recombination to mutation (r/m) of ∼24 that is homogenous over time. Especially in the early transcriptional region, we detect frequent gene loss and regain due to recombination with phages of the 936 group, demonstrating the role of the 936 group pangenome as a reservoir of genetic variation. The observed substitution rate homogeneity conforms to the neutral theory of evolution; hence, the neutral theory can be applied to phage genome evolution and also to genetic variation brought about by recombination.

Conservation management strategies for many highly threatened species include conservation breeding to prevent extinction and enhance recovery. Pairing decisions for these conservation breeding programmes can be informed by pedigree data to minimize relatedness between individuals in an effort to avoid inbreeding, maximize diversity and maintain evolutionary potential. However, conservation breeding programmes struggle to use this approach when pedigrees are shallow or incomplete. While genetic data (i.e., microsatellites) can be used to estimate relatedness to inform pairing decisions, emerging evidence indicates this approach may lack precision in genetically depauperate species, and more effective estimates will likely be obtained from genomic data (i.e., thousands of genome‐wide single nucleotide polymorphisms, or SNPs). Here, we compare relatedness estimates and subsequent pairing decisions using pedigrees, microsatellites and SNPs from whole‐genome resequencing approaches in two critically endangered birds endemic to New Zealand: kakī/black stilt (Himantopus novaezelandiae) and kākāriki karaka/orange‐fronted parakeet (Cyanoramphus malherbi). Our findings indicate that SNPs provide more precise estimates of relatedness than microsatellites when assessing empirical parent–offspring and full sibling relationships. Further, our results show that relatedness estimates and subsequent pairing recommendations using PMx are most similar between pedigree‐ and SNP‐based approaches. These combined results indicate that in lieu of robust pedigrees, SNPs are an effective tool for informing pairing decisions, which has important implications for many poorly pedigreed conservation breeding programmes worldwide.

Differentiation into infectious cysts through the process of encystation is crucial for transmission and survival of the intestinal protozoan parasite Giardia intestinalis. Hitherto the majority of studies have focused on the early events, leaving late encystation poorly defined. In order to further study encystation, focusing on the later events, we developed a new encystation protocol that generates a higher yield of mature cysts compared to standard methods. Transcriptome changes during the entire differentiation from trophozoites to cysts were thereafter studied using RNA sequencing (RNA-seq). A high level of periodicity was observed for up- and down-regulated genes, both at the level of the entire transcriptome and putative regulators. This suggests the trajectory of differentiation to be coordinated through developmentally linked gene regulatory activities. Our study identifies a core of 13 genes that are consistently up-regulated during initial encystation. Of these, two constitute previously uncharacterized proteins that we were able to localize to a new type of encystation-specific vesicles. Interestingly, the largest transcriptional changes were seen in the late phase of encystation with the majority of the highly up-regulated genes encoding hypothetical proteins. Several of these were epitope-tagged and localized to further characterize these previously unknown genetic components of encystation and possibly excystation. Finally, we also detected a switch of variant specific surface proteins (VSPs) in the late phase of encystation. This occurred at the same time as nuclear division and DNA replication, suggesting a potential link between the processes.

The RNA world hypothesis, that RNA genomes and catalysts preceded DNA genomes and genetically-encoded protein catalysts, has been central to models for the early evolution of life on Earth. A key part of such models is continuity between the earliest stages in the evolution of life and the RNA repertoires of extant lineages. Some assessments seem consistent with a diverse RNA world, yet direct continuity between modern RNAs and an RNA world has not been demonstrated for the majority of RNA families, and, anecdotally, many RNA functions appear restricted in their distribution. Despite much discussion of the possible antiquity of RNA families, no systematic analyses of RNA family distribution have been performed. To chart the broad evolutionary history of known RNA families, we performed comparative genomic analysis of over 3 million RNA annotations spanning 1446 families from the Rfam 10 database. We report that 99% of known RNA families are restricted to a single domain of life, revealing discrete repertoires for each domain. For the 1% of RNA families/clans present in more than one domain, over half show evidence of horizontal gene transfer (HGT), and the rest show a vertical trace, indicating the presence of a complex protein synthesis machinery in the Last Universal Common Ancestor (LUCA) and consistent with the evolutionary history of the most ancient protein-coding genes. However, with limited interdomain transfer and few RNA families exhibiting demonstrable antiquity as predicted under RNA world continuity, our results indicate that the majority of modern cellular RNA repertoires have primarily evolved in a domain-specific manner.

Here, we present the genome of the industrial ethanol production strain Brettanomyces bruxellensis CBS 11270. The nuclear genome was found to be diploid, containing four chromosomes with sizes of ranging from 2.2 to 4.0 Mbp. A 75 Kbp mitochondrial genome was also identified. Comparing the homologous chromosomes, we detected that 0.32% of nucleotides were polymorphic, i.e. formed single nucleotide polymorphisms (SNPs), 40.6% of them were found in coding regions (i.e. 0.13% of all nucleotides formed SNPs and were in coding regions). In addition, 8,538 indels were found. The total number of protein coding genes was 4897, of them, 4,284 were annotated on chromosomes; and the mitochondrial genome contained 18 protein coding genes. Additionally, 595 genes, which were annotated, were on contigs not associated with chromosomes. A number of genes was duplicated, most of them as tandem repeats, including a six-gene cluster located on chromosome 3. There were also examples of interchromosomal gene duplications, including a duplication of a six-gene cluster, which was found on both chromosomes 1 and 4. Gene copy number analysis suggested loss of heterozygosity for 372 genes. This may reflect adaptation to relatively harsh but constant conditions of continuous fermentation. Analysis of gene topology showed that most of these losses occurred in clusters of more than one gene, the largest cluster comprising 33 genes. Comparative analysis against the wine isolate CBS 2499 revealed 88,534 SNPs and 8,133 indels. Moreover, when the scaffolds of the CBS 2499 genome assembly were aligned against the chromosomes of CBS 11270, many of them aligned completely, some have chunks aligned to different chromosomes, and some were in fact rearranged. Our findings indicate a highly dynamic genome within the species B. bruxellensis and a tendency towards reduction of gene number in long-term continuous cultivation.

Background During infection by intracellular pathogens, a highly complex interplay occurs between the infected cell trying to degrade the invader and the pathogen which actively manipulates the host cell to enable survival and proliferation. Many intracellular pathogens pose important threats to human health and major efforts have been undertaken to better understand the host-pathogen interactions that eventually determine the outcome of the infection. Over the last decades, the unicellular eukaryote Dictyostelium discoideum has become an established infection model, serving as a surrogate macrophage that can be infected with a wide range of intracellular pathogens. In this study, we use high-throughput RNA-sequencing to analyze the transcriptional response of D. discoideum when infected with Mycobacterium marinum and Legionella pneumophila . The results were compared to available data from human macrophages . Results The majority of the transcriptional regulation triggered by the two pathogens was found to be unique for each bacterial challenge. Hallmark transcriptional signatures were identified for each infection, e.g. induction of endosomal sorting complexes required for transport (ESCRT) and autophagy genes in response to M. marinum and inhibition of genes associated with the translation machinery and energy metabolism in response to L. pneumophila . However, a common response to the pathogenic bacteria was also identified, which was not induced by non-pathogenic food bacteria. Finally, comparison with available data sets of regulation in human monocyte derived macrophages shows that the elicited response in D. discoideum is in many aspects similar to what has been observed in human immune cells in response to Mycobacterium tuberculosis and L. pneumophila. Conclusions Our study presents high-throughput characterization of D. discoideum transcriptional response to intracellular pathogens using RNA-seq. We demonstrate that the transcriptional response is in essence distinct to each pathogen and that in many cases, the corresponding regulation is recapitulated in human macrophages after infection by mycobacteria and L. pneumophila . This indicates that host-pathogen interactions are evolutionary conserved, derived from the early interactions between free-living phagocytic cells and bacteria. Taken together, our results strengthen the use of D. discoideum as a general infection model.

Massively parallel DNA sequencing of clinical samples holds great promise for the gene-based diagnosis of human inherited diseases because it allows rapid detection of putatively causative mutations at genome-wide level. Without additional evidence complementing their initial bioinformatics evaluation, however, the clinical relevance of such candidate genetic variants often remains unclear. In consequence, dedicated 'matching' services have been established in recent years that aim at the discovery of other, comparable case reports to facilitate individual diagnoses. However, legal concerns have been raised about the global sharing of genetic data, particularly in Europe where the recently enacted General Data Protection Regulation EU-2016/679 classifies genetic data as highly sensitive. Hence, unrestricted sharing of genetic data from clinical cases on platforms outside the national jurisdiction increasingly may be perceived as problematic. To allow collaborative data producers, particularly large consortia of diagnostic laboratories, to acknowledge these concerns while still practicing efficient case matching internally, novel tools are required. To this end, we developed VarWatch, an easy-to-deploy and highly scalable case matching software that provides users with comprehensive programmatic tools and a user-friendly interface to fulfil said purpose.

Evolutionary arms races between pathogens and their hosts may be manifested as selection for rapid evolutionary change of key genes, and are sometimes detectable through sequence-level analyses. In the case of protein-coding genes, such analyses frequently predict that specific codons are under positive selection. However, detecting positive selection can be non-trivial, and false positive predictions are a common concern in such analyses. It is therefore helpful to place such predictions within a structural and functional context. Here, we focus on the p19 protein from tombusviruses. P19 is a homodimer that sequesters siRNAs, thereby preventing the host RNAi machinery from shutting down viral infection. Sequence analysis of the p19 gene is complicated by the fact that it is constrained at the sequence level by overprinting of a viral movement protein gene. Using homology modeling, in silico mutation and molecular dynamics simulations, we assess how non-synonymous changes to two residues involved in forming the dimer interface-one invariant, and one predicted to be under positive selection-impact molecular function. Interestingly, we find that both observed variation and potential variation (where a non-synonymous change to p19 would be synonymous for the overprinted movement protein) does not significantly impact protein structure or RNA binding. Consequently, while several methods identify residues at the dimer interface as being under positive selection, MD results suggest they are functionally indistinguishable from a site that is free to vary. Our analyses serve as a caveat to using sequence-level analyses in isolation to detect and assess positive selection, and emphasize the importance of also accounting for how non-synonymous changes impact structure and function.