To Be or Not to Be Expressed: The First Evidence of a Nucleolar Dominance Tissue-Specificity in Brachypodium hybridum

Nucleolar dominance (ND) is an epigenetic, developmentally regulated phenomenon that describes the selective inactivation of 35S rDNA loci derived from one progenitor of a hybrid or allopolyploid. The presence of ND was documented in an allotetraploid grass, Brachypodium hybridum (genome composition DDSS), which is a polyphyletic species that arose from crosses between two putative ancestors that resembled the modern B. distachyon (DD) and B. stacei (SS). In this work, we investigated the developmental stability of ND in B. hybridum genotype 3-7-2 and compared it with the reference genotype ABR113. We addressed the question of whether the ND is established in generative tissues such as pollen mother cells (PMC). We examined condensation of rDNA chromatin by fluorescence in situ hybridization employing state-of-art confocal microscopy. The transcription of rDNA homeologs was determined by reverse-transcription cleaved amplified polymorphic sequence analysis. In ABR113, the ND was stable in all tissues analyzed (primary and adventitious root, leaf, and spikes). In contrast, the 3-7-2 individuals showed a strong upregulation of the S-genome units in adventitious roots but not in other tissues. Microscopic analysis of the 3-7-2 PMCs revealed extensive decondensation of the D-genome loci and their association with the nucleolus in meiosis. As opposed, the S-genome loci were always highly condensed and localized outside the nucleolus. These results indicate that genotype-specific loss of ND in B. hybridum occurs probably after fertilization during developmental processes. This finding supports our view that B. hybridum is an attractive model to study ND in grasses.

Properties of Two Broad Host Range Phages of Yersinia enterocolitica Isolated from Wild Animals

Yersinia (Y.) enterocolitica and Y. pseudotuberculosis are important zoonotic agents which can infect both humans and animals. To combat these pathogens, the application of strictly lytic phages may be a promising tool. Since only few Yersinia phages have been described yet, some of which demonstrated a high specificity for certain serotypes, we isolated two phages from game animals and characterized them in terms of their morphology, host specificity, lytic activity on two bio-/serotypes and genome composition. The T7-related podovirus vB_YenP_Rambo and the myovirus vB_YenM_P281, which is very similar to a previously described phage PY100, showed a broad host range. Together, they lysed all the 62 tested pathogenic Y. enterocolitica strains belonging to the most important bio-/serotypes in Europe. A cocktail containing these two phages strongly reduced cultures of a bio-/serotype B4/O:3 and a B2/O:9 strain, even at very low MOIs (multiplicity of infection) and different temperatures, though, lysis of bio-/serotype B2/O:9 by vB_YenM_P281 and also by the related phage PY100 only occurred at 37 °C. Both phages were additionally able to lyse various Y. pseudotuberculosis strains at 28 °C and 37 °C, but only when the growth medium was supplemented with calcium and magnesium cations.

A replication-linked mutational gradient drives somatic mutation accumulation and influences germline polymorphisms and genome composition in mitochondrial DNA

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational ‘hot-spots’ or ‘cold-spots’. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.

Partners in space: Discordant population structure between legume hosts and rhizobium symbionts in their native range

Coevolution is predicted to depend on how the genetic diversity of interacting species is geographically structured. Plant-microbe symbioses such as the legume-rhizobium mutualism are ecologically and economically important, but distinct life history and dispersal mechanisms for these host and microbial partners, plus dynamic genome composition in bacteria, present challenges for understanding spatial genetic processes in these systems. Here we study the model rhizobium Ensifer meliloti using a hierarchically-structured sample of 191 strains from 21 sites in the native range and compare its population structure to that of its host plant Medicago truncatula. We find high local genomic variation and minimal isolation by distance across the rhizobium genome, particularly at the two symbiosis elements pSymA and pSymB, which have evolutionary histories and population structures that are similar to each other but distinct from both the chromosome and the host. While the chromosome displays weak isolation by distance, it is uncorrelated with hosts. Patterns of discordant population structure among elements with the bacterial genome has implications for bacterial adaptation to life in the soil versus symbiosis, while discordant population genetic structure of hosts and microbes might restrict local adaptation of species to each other and give rise to phenotypic mismatches in coevolutionary traits.

Anaerobic Fungal Mevalonate Pathway Genomic Biases Lead to Heterologous Toxicity Underpredicted by Codon Adaptation Indices

Anaerobic fungi are emerging biotechnology platforms with genomes rich in biosynthetic potential. Yet, the heterologous expression of their biosynthetic pathways has had limited success in model hosts like E. coli. We find one reason for this is that the genome composition of anaerobic fungi like P. indianae are extremely AT-biased with a particular preference for rare and semi-rare AT-rich tRNAs in E coli, which are not explicitly predicted by standard codon adaptation indices (CAI). Native P. indianae genes with these extreme biases create drastic growth defects in E. coli (up to 69% reduction in growth), which is not seen in genes from other organisms with similar CAIs. However, codon optimization rescues growth, allowing for gene evaluation. In this manner, we demonstrate that anaerobic fungal homologs such as PI.atoB are more active than S. cerevisiae homologs in a hybrid pathway, increasing the production of mevalonate up to 2.5 g/L (more than two-fold) and reducing waste carbon to acetate by ~90% under the conditions tested. This work demonstrates the bioproduction potential of anaerobic fungal enzyme homologs and how the analysis of codon utilization enables the study of otherwise difficult to express genes that have applications in biocatalysis and natural product discovery.

Unraveling phylogenetic relationships, reticulate evolution, and genome composition of polyploid plant complexes by RAD-Seq and Hyb-Seq

Complex genome evolution of young polyploid complexes is poorly understood. Besides challenges caused by hybridization, polyploidization, and incomplete lineage sorting, bioinformatic analyses are often exacerbated by missing information on progenitors, ploidy, and reproduction modes. By using a comprehensive, self-developed bioinformatic pipeline covering tree, structure, network, and SNP-origin analyses, we for the first time unraveled polyploid phylogenetic relationships and genome evolution within the large Eurasian Ranunculus auricomus species complex comprising more than 840 taxa. Our results rely on 97,312 genomic RADseq loci, target enrichment of 576 nuclear genes (48 phased), and 71 plastid regions (Hybseq; OMICS-data) derived from the 75 most widespread polyploid apomictic taxa and four di- and one tetraploid potential sexual progenitor species. Phylogenetic tree and structure analyses consistently showed 3-5 supported polyploid groups, each containing sexual progenitor species. In total, analyses revealed four diploid sexual progenitors and a one unknown, probably extinct progenitor, contributing to the genome composition of R. auricomus polyploids. Phylogenetic network, structure, and SNP-origin analyses based on RADseq loci and phased nuclear genes completed by plastid data demonstrated predominantly allopolyploid origins, each involving 2-3 different diploid sexual subgenomes. Allotetraploid genomes were characterized by subgenome dominance and large proportions of interspecific, non-hybrid SNPs, indicating an enormous degree of post-origin evolution (i.e., Mendelian segregation of the diploid hybrid generations, back-crossings, and gene flow due to facultative sexuality of apomicts), but only low proportions of lineage-specific SNPs. The R. auricomus model system is the first large European polyploid species complex studied with reduced representation OMICS data. Our bioinformatic pipeline underlines the importance of combining different approaches and datasets to successfully unveil how reticulate evolution and post-origin processes shape the diversity of polyploid plant complexes.

Gradual evolution of allopolyploidy in Arabidopsis suecica

Most diploid organisms have polyploid ancestors. The evolutionary process of polyploidization is poorly understood but has frequently been conjectured to involve some form of ‘genome shock’, such as genome reorganization and subgenome expression dominance. Here we study polyploidization in Arabidopsis suecica, a post-glacial allopolyploid species formed via hybridization of Arabidopsis thaliana and Arabidopsis arenosa. We generated a chromosome-level genome assembly of A. suecica and complemented it with polymorphism and transcriptome data from all species. Despite a divergence around 6 million years ago (Ma) between the ancestral species and differences in their genome composition, we see no evidence of a genome shock: the A. suecica genome is colinear with the ancestral genomes; there is no subgenome dominance in expression; and transposon dynamics appear stable. However, we find changes suggesting gradual adaptation to polyploidy. In particular, the A. thaliana subgenome shows upregulation of meiosis-related genes, possibly to prevent aneuploidy and undesirable homeologous exchanges that are observed in synthetic A. suecica, and the A. arenosa subgenome shows upregulation of cyto-nuclear processes, possibly in response to the new cytoplasmic environment of A. suecica, with plastids maternally inherited from A. thaliana. These changes are not seen in synthetic hybrids, and thus are likely to represent subsequent evolution.

The five Urochola spp. used in development of tropical forage cultivars originate from defined subpopulations with differentiated gene pools

Background and Aims: Urochola (syn. Brachiaria, and including some Panicum and Megathyrus) is a genus of tropical and subtropical grasses widely sown as forage to feed ruminants in the tropics. A better understanding of the diversity among Urochola spp. allow us to leverage its varying ploidy levels and genome composition to accelerate its improvement, following the example from other crop genera. Methods: We explored the genetic make-up and population structure in 111 accessions, which comprise the five Urochola species used for the development of commercial cultivars. These accessions are conserved from wild materials from collection sites at their centre of origin in Africa. We used RNA-seq, averaging 40M reads per accession, to generate 1,167,542 stringently selected SNP markers that tentatively encompassed the complete Urochola gene pool used in breeding. Key Results: We identified ten subpopulations, which had no relation with geographical origin and represented ten independent gene pools, and two groups of admixed accessions. Our results support a division in U. decumbens by ploidy, with a diploid subpopulation closely related to U. ruziziensis, and a tetraploid subpopulation closely related to U. brizantha. We observed highly differentiated gene pools in U. brizantha, which were not related with origin or ploidy. Particularly, one U. brizantha subpopulation clustered distant from the other U. brizantha and U. decumbens subpopulations, so likely containing unexplored alleles. We also identified a well-supported subpopulation containing both polyploid U. decumbens and U. brizantha accessions; this was the only group containing more than one species and tentatively constitutes an independent "mixed" gene pool for both species. We observed two gene pools in U. humidicola. One subpopulation, "humidicola-2", was much less common but likely includes the only known sexual accession in the species. Conclusions: Our results offered a definitive picture of the available diversity in Urochola to inform breeding and resolve questions raised by previous studies. It also allowed us identifying prospective founders to enrich the breeding gene pool and to develop genotyping and genotype-phenotype association mapping experiments.

Whole-Genome Sequences of Mycobacterium abscessus subsp. massiliense Isolates from Brazil

The Mycobacterium abscessus complex comprises multidrug-resistant, opportunistic, and rapidly growing pathogens responsible for severe infections. Here, we report the genome composition of four Mycobacterium abscessus subsp. massiliense isolates from three sources: two from the lung of a cystic fibrosis patient, one from a mammary cyst, and one from a gutter system.