Mosaic Genome Research Articles

Virus-induced gene silencing for in planta validation of gene function in cucurbits.

Virus-induced gene silencing (VIGS) is a powerful tool for high-throughput analysis of gene function. Here, we developed the VIGS vector pCF93, from which expression of the cucumber fruit mottle mosaic virus genome is driven by the cauliflower mosaic virus 35S promoter to produce viral transcripts in inoculated plants. To test the utility of the pCF93 vector, we identified candidate genes related to male sterility (MS) in watermelon (Citrullus lanatus), which is recalcitrant to genetic transformation. Specifically, we exploited previously reported reference-based and de novo transcriptome data to define 38 differentially expressed genes between a male-sterile line and its fertile near-isogenic line in the watermelon cultivar DAH. We amplified 200- to 300-bp fragments of these genes, cloned them into pCF93, and inoculated DAH with the resulting VIGS clones. The small watermelon cultivar DAH enabled high-throughput screening using a small cultivation area. We simultaneously characterized the phenotypes associated with each of the 38 candidate genes in plants grown in a greenhouse. Silencing of 8 of the 38 candidate genes produced male-sterile flowers with abnormal stamens and no pollen. We confirmed the extent of gene silencing in inoculated flowers using reverse transcription-qPCR. Histological analysis of stamens from male-fertile and male-sterile floral buds and mature flowers revealed developmental defects and shrunken pollen sacs. Based on these findings, we propose that the pCF93 vector and our VIGS system will facilitate high-throughput analysis for the study of gene function in watermelons.

Sympatric Recombination in Zoonotic Cryptosporidium Leads to Emergence of Populations with Modified Host Preference.

Genetic recombination plays a critical role in the emergence of pathogens with phenotypes such as drug resistance, virulence, and host adaptation. Here, we tested the hypothesis that recombination between sympatric ancestral populations leads to the emergence of divergent variants of the zoonotic parasite Cryptosporidium parvum with modified host ranges. Comparative genomic analyses of 101 isolates have identified seven subpopulations isolated by distance. They appear to be descendants of two ancestral populations, IIa in northwestern Europe and IId from southwestern Asia. Sympatric recombination in areas with both ancestral subtypes and subsequent selective sweeps have led to the emergence of new subpopulations with mosaic genomes and modified host preference. Subtelomeric genes could be involved in the adaptive selection of subpopulations, while copy number variations of genes encoding invasion-associated proteins are potentially associated with modified host ranges. These observations reveal ancestral origins of zoonotic C. parvum and suggest that pathogen import through modern animal farming might promote the emergence of divergent subpopulations of C. parvum with modified host preference.

Characterization of Mungbean yellow mosaic India virus genome with a recombinant DNA-B in Southern Peninsular India.

Mungbean yellow mosaic India virus (MYMIV) is a representative of the genus begomovirus/Begomoviridae, which is prevalent in the northern part of Indian subcontinent causing yellow mosaic disease (YMD). This virus is rapidly evolving and breaking the resistance in the advanced lines causing huge economic losses in the pulse production. In this context, the present investigation on characterization of the causal organism of YMD was undertaken METHODS AND RESULTS: A novel recombinant isolate (YMV-BG-BPT) causing YMD was identified from blackgram in Andhra Pradesh, southern peninsular region of India. The association of a bipartite begomovirus with the disease was done by sequence analyses of the cloned full-length genome. The full length genome sequences were submitted in NCBI GenBank with accession numbers MZ235792 (DNA-A) and MZ356197 (DNA-B). The sequence analysis of DNA-A of YMV-BG-BPT showed maximum of 99.12% similarity at nucleotide level with Mungbean yellow mosaic India virus (MYMIV) isolate reported from Tamil Nadu (KC911719), India which is also confirmed by clustering pattern in phylogenic analysis and DNA-B showed 95.79% with Mungbean yellow mosaic virus (MYMV) isolate reported from Tamil Nadu (KP319016) and 95.05% with MYMIV isolate reported from Karnataka (MT027037). The huge variation in DNA-B lead us to suspect a recombination in DNA-B, where a recombination event in the CR, region coding for nuclear shuttle protein and movement protein of DNA B was detected in which MYMV-BG-AP-IND (KF928962) and MYMIV-GG-CH-IND (MN020536) have been identified as major and minor parents, respectively. Overall, the present study revealed occurrence of MYMIV with recombinant DNA B component in southern peneinsular India.

The mosaic oat genome gives insights into a uniquely healthy cereal crop

Cultivated oat (Avena sativa L.) is an allohexaploid (AACCDD, 2n = 6x = 42) thought to have been domesticated more than 3,000 years ago while growing as a weed in wheat, emmer and barley fields in Anatolia1,2. Oat has a low carbon footprint, substantial health benefits and the potential to replace animal-based food products. However, the lack of a fully annotated reference genome has hampered efforts to deconvolute its complex evolutionary history and functional gene dynamics. Here we present a high-quality reference genome of A. sativa and close relatives of its diploid (Avena longiglumis, AA, 2n = 14) and tetraploid (Avena insularis, CCDD, 2n = 4x = 28) progenitors. We reveal the mosaic structure of the oat genome, trace large-scale genomic reorganizations in the polyploidization history of oat and illustrate a breeding barrier associated with the genome architecture of oat. We showcase detailed analyses of gene families implicated in human health and nutrition, which adds to the evidence supporting oat safety in gluten-free diets, and we perform mapping-by-sequencing of an agronomic trait related to water-use efficiency. This resource for the Avena genus will help to leverage knowledge from other cereal genomes, improve understanding of basic oat biology and accelerate genomics-assisted breeding and reanalysis of quantitative trait studies.

Evolution and origin of bread wheat.

Bread wheat (Triticum aestivum, genome BBAADD) is a young hexaploid species formed only 8,500–9,000 years ago through hybridization between a domesticated free-threshing tetraploid progenitor, genome BBAA, and Aegilops tauschii, the diploid donor of the D subgenome. Very soon after its formation, it spread globally from its cradle in the fertile crescent into new habitats and climates, to become a staple food of humanity. This extraordinary global expansion was probably enabled by allopolyploidy that accelerated genetic novelty through the acquisition of new traits, new intergenomic interactions, and buffering of mutations, and by the attractiveness of bread wheat’s large, tasty, and nutritious grain with high baking quality. New genome sequences suggest that the elusive donor of the B subgenome is a distinct (unknown or extinct) species rather than a mosaic genome. We discuss the origin of the diploid and tetraploid progenitors of bread wheat and the conflicting genetic and archaeological evidence on where it was formed and which species was its free-threshing tetraploid progenitor. Wheat experienced many environmental changes throughout its evolution, therefore, while it might adapt to current climatic changes, efforts are needed to better use and conserve the vast gene pool of wheat biodiversity on which our food security depends.

Genomic Sequencing and Analysis of a Novel Human Cowpox Virus With Mosaic Sequences From North America and Old World Orthopoxvirus.

Orthopoxviruses (OPXVs) not only infect their natural hosts, but some OPXVs can also cause disease in humans. Previously, we partially characterized an OPXV isolated from an 18-year-old male living in Northern Norway. Restriction enzyme analysis and partial genome sequencing characterized this virus as an atypical cowpox virus (CPXV), which we named CPXV-No-H2. In this study, we determined the complete genome sequence of CPXV-No-H2 using Illumina and Nanopore sequencing. Our results showed that the whole CPXV-No-H2 genome is 220,276 base pairs (bp) in length, with inverted terminal repeat regions of approximately 7 kbp, containing 217 predicted genes. Seventeen predicted CPXV-No-H2 proteins were most similar to OPXV proteins from the Old World, including Ectromelia virus (ECTV) and Vaccinia virus, and North America, Alaskapox virus (AKPV). CPXV-No-H2 has a mosaic genome with genes most similar to other OPXV genes, and seven potential recombination events were identified. The phylogenetic analysis showed that CPXV-No-H2 formed a separate clade with the German CPXV isolates CPXV_GerMygEK938_17 and CPXV_Ger2010_MKY, sharing 96.4 and 96.3% nucleotide identity, respectively, and this clade clustered closely with the ECTV-OPXV Abatino clade. CPXV-No-H2 is a mosaic virus that may have arisen out of several recombination events between OPXVs, and its phylogenetic clustering suggests that ECTV-Abatino-like cowpox viruses form a distinct, new clade of cowpox viruses.

Genomic Bootstrap Barcodes and Their Application to Study the Evolution of Sarbecoviruses.

Recombination creates mosaic genomes containing regions with mixed ancestry, and the accumulation of such events over time can complicate greatly many aspects of evolutionary inference. Here, we developed a sliding window bootstrap (SWB) method to generate genomic bootstrap (GB) barcodes to highlight the regions supporting phylogenetic relationships. The method was applied to an alignment of 56 sarbecoviruses, including SARS-CoV and SARS-CoV-2, responsible for the SARS epidemic and COVID-19 pandemic, respectively. The SWB analyses were also used to construct a consensus tree showing the most reliable relationships and better interpret hidden phylogenetic signals. Our results revealed that most relationships were supported by just a few genomic regions and confirmed that three divergent lineages could be found in bats from Yunnan: SCoVrC, which groups SARS-CoV related coronaviruses from China; SCoV2rC, which includes SARS-CoV-2 related coronaviruses from Southeast Asia and Yunnan; and YunSar, which contains a few highly divergent viruses recently described in Yunnan. The GB barcodes showed evidence for ancient recombination between SCoV2rC and YunSar genomes, as well as more recent recombination events between SCoVrC and SCoV2rC genomes. The recombination and phylogeographic patterns suggest a strong host-dependent selection of the viral RNA-dependent RNA polymerase. In addition, SARS-CoV-2 appears as a mosaic genome composed of regions sharing recent ancestry with three bat SCoV2rCs from Yunnan (RmYN02, RpYN06, and RaTG13) or related to more ancient ancestors in bats from Yunnan and Southeast Asia. Finally, our results suggest that viral circular RNAs may be key molecules for the mechanism of recombination.

A myxozoan genome reveals mosaic evolution in a parasitic cnidarian

BackgroundParasite evolution has been conceptualized as a process of genetic loss and simplification. Contrary to this model, there is evidence of expansion and conservation of gene families related to essential functions of parasitism in some parasite genomes, reminiscent of widespread mosaic evolution—where subregions of a genome have different rates of evolutionary change. We found evidence of mosaic genome evolution in the cnidarian Myxobolus honghuensis, a myxozoan parasite of fish, with extremely simple morphology.ResultsWe compared M. honghuensis with other myxozoans and free-living cnidarians, and determined that it has a relatively larger myxozoan genome (206 Mb), which is less reduced and less compact due to gene retention, large introns, transposon insertion, but not polyploidy. Relative to other metazoans, the M. honghuensis genome is depleted of neural genes and has only the simplest animal immune components. Conversely, it has relatively more genes involved in stress resistance, tissue invasion, energy metabolism, and cellular processes compared to other myxozoans and free-living cnidarians. We postulate that the expansion of these gene families is the result of evolutionary adaptations to endoparasitism. M. honghuensis retains genes found in free-living Cnidaria, including a reduced nervous system, myogenic components, ANTP class Homeobox genes, and components of the Wnt and Hedgehog pathways.ConclusionsOur analyses suggest that the M. honghuensis genome evolved as a mosaic of conservative, divergent, depleted, and enhanced genes and pathways. These findings illustrate that myxozoans are not as genetically simple as previously regarded, and the evolution of some myxozoans is driven by both genomic streamlining and expansion.

Bat coronaviruses related to SARS-CoV-2 and infectious for human cells.

The animal reservoir of SARS-CoV-2 is unknown despite reports of SARS-CoV-2-related viruses in Asian Rhinolophus bats1-4, including the closest virus from R. affinis, RaTG13 (refs. 5,6), and pangolins7-9. SARS-CoV-2 has a mosaic genome, to which different progenitors contribute. The spike sequence determines the binding affinity and accessibility of its receptor-binding domain to the cellular angiotensin-converting enzyme 2 (ACE2) receptor and is responsible for host range10-12. SARS-CoV-2 progenitor bat viruses genetically close to SARS-CoV-2 and able to enter human cells through a human ACE2 (hACE2) pathway have not yet been identified, although they would be key in understanding the origin of the epidemic. Here we show that such viruses circulate in cave bats living in the limestone karstic terrain in northern Laos, in the Indochinese peninsula. We found that the receptor-binding domains of these viruses differ from that of SARS-CoV-2 by only one or two residues at the interface with ACE2, bind more efficiently to the hACE2 protein than that of the SARS-CoV-2 strain isolated in Wuhan from early human cases, and mediate hACE2-dependent entry and replication in human cells, which is inhibited by antibodies that neutralize SARS-CoV-2. None of these bat viruses contains a furin cleavage site in the spike protein. Our findings therefore indicate that bat-borne SARS-CoV-2-like viruses that arepotentially infectious for humans circulate in Rhinolophus spp. in the Indochinese peninsula.

Genomic interrogation of a MAGIC population highlights genetic factors controlling fiber quality traits in cotton

Cotton (Gossypium hirsutum L.) fiber is the most important resource of natural and renewable fiber for the textile industry. However, the understanding of genetic components and their genome-wide interactions controlling fiber quality remains fragmentary. Here, we sequenced a multiple-parent advanced-generation inter-cross (MAGIC) population, consisting of 550 individuals created by inter-crossing 11 founders, and established a mosaic genome map through tracing the origin of haplotypes that share identity-by-descent (IBD). We performed two complementary GWAS methods—SNP-based GWAS (sGWAS) and IBD-based haplotype GWAS (hGWAS). A total of 25 sQTLs and 14 hQTLs related to cotton fiber quality were identified, of which 26 were novel QTLs. Two major QTLs detected by both GWAS methods were responsible for fiber strength and length. The gene Ghir_D11G020400 (GhZF14) encoding the MATE efflux family protein was identified as a novel candidate gene for fiber length. Beyond the additive QTLs, we detected prevalent epistatic interactions that contributed to the genetics of fiber quality, pinpointing another layer for trait variance. This study provides new targets for future molecular design breeding of superior fiber quality.

Comparative genomic analysis of dwarf Vibrio myoviruses defines a conserved gene cluster for successful phage infection.

Tailed bacteriophages have been at the center of attention, not only for their ability to infect and kill pathogenic bacteria but also due to their peculiar and intriguing complex contractile tail structure. Tailed bacteriophages with contractile tails are known to have a Myoviridae morphotype and are members of the order Caudovirales. Large bacteriophages with a genome larger than 150 kbp have been studied for their ability to use multiple infection and lysis strategies to replicate more efficiently. On the other hand, smaller bacteriophages with fewer genes are represented in the GenBank database in greater numbers, and have several genes with unknown function. Isolation and molecular characterization of a newly reported bacteriophage named Athena1 revealed that it is a strongly lytic bacteriophage with a genome size of 39,826 bp. This prompted us to perform a comparative genomic analysis of Vibrio myoviruses with a genome size of no more than 50 kbp. The results revealed a pattern of genomic organization that includes sets of genes responsible for virion morphogenesis, replication/recombination of DNA, and lysis/lysogeny switching. By studying phylogenetic gene markers, we were able to draw conclusions about evolutionary events that shaped the genomic mosaicism of these phages, pinpointing the importance of a conserved organization of the genomic region encoding the baseplate protein for successful infection of Gram-negative bacteria. In addition, we propose the creation of new genera for dwarf Vibrio myoviruses. Comparative genomics of phages infecting aquatic bacteria could provide information that is useful for combating fish pathogens in aquaculture, using novel strategies.

Genomic patterns and the evolutionary origin of an invasive fungal pathogen (Hymenoscyphus fraxineus) in Europe

In the 1990s, a highly pathogenic and invasive lineage of the Asian ascomycete Hymenoscyphus fraxineus appeared in Europe. Spreading across the continent, the fungal lineage devastates common ash (Fraxinus excelsior) populations and is a threat to ash species worldwide.In contrast to highly diverse reference populations in Japan and East Russia, the expanding European population is characterized by the persistence of a mosaic of two ancestral haplotypes, and presumed constraints to recombination. Four hypotheses were proposed for the evolutionary origin of the European lineage. The first hypothesis H1 states that two, and only two, pathogenic founder individuals were transferred from Asia to Europe and successfully established. Alternatively, in addition to the fungus arriving in Europe, a change in function might have occurred and only then the lineage became invasive. The cause of such a change might have been H2: one to two mutations during meiosis; H3: the lateral transfer of a European cytoplasmatic element modifying the behavior of H. fraxineus; or H4: a hybridization event with a closely related European ascomycete.Only the hybridization hypothesis seems to explain all of the genome-wide characteristics observed for the nuclear genome. Existing reconstructions of phylogenetic relationships at different taxonomic levels provide the evolutionary context of the lineages potentially involved. Vital information was gained from three configurations of nucleotide diversity and population differentiation associated with the genomic mosaic of monomorphic tracts and diversity islands.This synthesis points to the importance of identifying the taxonomic and geographic origins of the conserved two genomic ancestries within the European lineage of H. fraxineus and of reconstructing the mosaic of ancestral parental haplotypes (ancestry blocks) across the genome. Only with this information at hand will it be possible, to reliably interpret population-genomic data in the search for the functional elements underlying this unbalanced host-pathogen system.

RNA In Situ Hybridization of Detecting Cucumber Mosaic Virus in Shoots of Nicotiana benthamiana Plants.

RNA in situ hybridization, a histological technique derived from Southern blotting and northern blotting, has been an important approach in biology studies for many years. In the studies of virus-plant interactions, RNA in situ hybridization provides a direct visualization of viral RNA in host plants. Here, we provide a detailed protocol for viral RNA in situ hybridization that has been successfully used to detect Cucumber mosaic virus genome (CMV) RNAs in shoots of N. benthamiana plants.

Nitrate-inducible MdBT2 acts as a restriction factor to limit apple necrotic mosaic virus genome replication in Malus domestica.

Apple necrotic mosaic virus (ApNMV) is highly associated with the occurrence of apple mosaic disease in China. However, ApNMV–host interactions and defence mechanisms of host plants against this virus are poorly studied. Here, we report that nitrate treatment restrains ApNMV genomic RNA accumulation by destabilizing viral replication protein 1a through the MdBT2‐mediated ubiquitin‐proteasome pathway. MdBT2, a nitrate‐responsive BTB/TAZ domain‐containing protein, was identified in a yeast two‐hybrid screen of an apple cDNA library using viral protein 1a as bait, and 1a was further confirmed to interact with MdBT2 both in vivo and in vitro. It was further verified that MdBT2 promoted the ubiquitination and degradation of viral protein 1a through the ubiquitin‐proteasome pathway in an MdCUL3A‐independent manner. Viral genomic RNA accumulation was reduced in MdBT2‐overexpressing transgenic apple leaves but enhanced in MdBT2‐antisense leaves compared to the wild type. Moreover, MdBT2 was found to interfere with the interaction between viral replication proteins 1a and 2apol by competitively interacting with 1a. Taken together, our results demonstrate that nitrate‐inducible MdBT2 functions as a limiting factor in ApNMV viral RNA accumulation by promoting the ubiquitination and degradation of viral protein 1a and interfering with interactions between viral replication proteins.

Colistin Resistance Onset Strategies and Genomic Mosaicism in Clinical Acinetobacter baumannii Lineages.

The treatment of multidrug-resistant Gram-negative infections is based on colistin. As result, COL-resistance (COL-R) can develop and spread. In Acinetobacter baumannii, a crucial step is to understand COL-R onset and stability, still far to be elucidated. COL-R phenotypic stability, onset modalities, and phylogenomics were investigated in a clinical A. baumannii sample showing a COL resistant (COLR) phenotype at first isolation. COL-R was confirmed by Minimum-Inhibitory-Concentrations as well as investigated by Resistance-Induction assays and Population-Analysis-Profiles (PAPs) to determine: (i) stability; (ii) inducibility; (iii) heteroresistance. Genomics was performed by Mi-Seq Whole-Genome-Sequencing, Phylogenesis, and Genomic Epidemiology by bioinformatics. COLR A. baumannii were subdivided as follows: (i) 3 A. baumannii with stable and high COL MICs defining the “homogeneous-resistant” onset phenotype; (ii) 6 A. baumannii with variable and lower COL MICs displaying a “COL-inducible” onset phenotype responsible for adaptive-resistance or a “subpopulation” onset phenotype responsible for COL-heteroresistance. COL-R stability and onset strategies were not uniquely linked to the amount of LPS and cell envelope charge. Phylogenomics categorized 3 lineages clustering stable and/or unstable COL-R phenotypes with increasing genomic complexity. Likewise, different nsSNP profiling in genes already associated with COL-R marked the stable and/or unstable COL-R phenotypes. Our investigation finds out that A. baumannii can range through unstable or stable COLR phenotypes emerging via different “onset strategies” within phylogenetic lineages displaying increasing genomic mosaicism.

MoG+: a database of genomic variations across three mouse subspecies for biomedical research

Laboratory mouse strains have mosaic genomes derived from at least three major subspecies that are distributed in Eurasia. Here, we describe genomic variations in ten inbred strains: Mus musculus musculus-derived BLG2/Ms, NJL/Ms, CHD/Ms, SWN/Ms, and KJR/Ms; M. m. domesticus-derived PGN2/Ms and BFM/Ms; M. m. castaneus-derived HMI/Ms; and JF1/Ms and MSM/Ms, which were derived from a hybrid between M. m. musculus and M. m. castaneus. These strains were established by Prof. Moriwaki in the 1980s and are collectively named the “Mishima Battery”. These strains show large phenotypic variations in body size and in many physiological traits. We resequenced the genomes of the Mishima Battery strains and performed a comparative genomic analysis with dbSNP data. More than 81 million nucleotide coordinates were identified as variant sites due to the large genetic distances among the mouse subspecies; 8,062,070 new SNP sites were detected in this study, and these may underlie the large phenotypic diversity observed in the Mishima Battery. The new information was collected in a reconstructed genome database, termed MoG+ that includes new application software and viewers. MoG+ intuitively visualizes nucleotide variants in genes and intergenic regions, and amino acid substitutions across the three mouse subspecies. We report statistical data from the resequencing and comparative genomic analyses and newly collected phenotype data of the Mishima Battery, and provide a brief description of the functions of MoG+, which provides a searchable and unique data resource of the numerous genomic variations across the three mouse subspecies. The data in MoG+ will be invaluable for research into phenotype-genotype links in diverse mouse strains.

Genome Mosaicism in Field Strains of Mycoplasma bovis as Footprints of In-Host Horizontal Chromosomal Transfer.

Horizontal gene transfer was long thought to be marginal in Mollicutes, but the capacity of some of these wall-less bacteria to exchange large chromosomal regions has been recently documented. Mycoplasma chromosomal transfer (MCT) is an unconventional mechanism that relies on the presence of a functional integrative conjugative element (ICE) in at least one partner and involves the horizontal acquisition of small and large chromosomal fragments from any part of the donor genome, which results in progenies composed of an infinite variety of mosaic genomes. The present study focuses on Mycoplasma bovis, an important pathogen of cattle responsible for major economic losses worldwide. By combining phylogenetic tree reconstructions and detailed comparative genome analyses of 36 isolates collected in Spain (2016 to 2018), we confirmed the mosaic nature of 16 field isolates and mapped chromosomal transfers exchanged between their hypothetical ancestors. This study provides evidence that MCT can take place in the field, most likely during coinfections by multiple strains. Because mobile genetic elements (MGEs) are classical contributors of genome plasticity, the presence of phages, insertion sequences (ISs), and ICEs was also investigated. Data revealed that these elements are widespread within the M. bovis species and evidenced classical horizontal transfer of phages and ICEs in addition to MCT. These events contribute to wide-genome diversity and reorganization within this species and may have a tremendous impact on diagnostic and disease control. IMPORTANCE Mycoplasma bovis is a major pathogen of cattle that has significant detrimental effects on economics and animal welfare in cattle rearing worldwide. Understanding the evolution and the adaptative potential of pathogenic mycoplasma species in the natural host is essential to combating them. In this study, we documented the occurrence of mycoplasma chromosomal transfer, an atypical mechanism of horizontal gene transfer, in field isolates of M. bovis that provide new insights into the evolution of this pathogenic species in their natural host. Although these events are expected to occur at low frequency, their impact is accountable for genome-wide variety and reorganization within M. bovis species, which may compromise both diagnostic and disease control.

Genome and Methylome analysis of a phylogenetic novel Campylobacter coli cluster with C. jejuni introgression.

The intriguing recent discovery of Campylobacter coli strains, especially of clade 1, that (i) possess mosaic C. coli / C. jejuni alleles, (ii) demonstrate mixed multilocus sequence types (MLSTs) and (iii) have undergone genome-wide introgression has led to the speculation that these two species may be involved in an accelerated rate of horizontal gene transfer that is progressively leading to the merging of both species in a process coined ‘despeciation’. In an MLST-based neighbour-joining tree of a number of C. coli and C. jejuni isolates of different clades, three prominent Campylobacter isolates formed a seemingly separate cluster besides the previously described C. coli and C. jejuni clades. In the light of the suspected, ongoing genetic introgression between the C. coli and C. jejuni species, this cluster of Campylobacter isolates is proposed to present one of the hybrid clonal complexes in the despeciation process of the genus. Specific DNA methylation as well as restriction modification systems are known to be involved in selective uptake of external DNA and their role in such genetic introgression remains to be further investigated. In this study, the phylogeny and DNA methylation of these putative C. coli / C. jejuni hybrid strains were explored, their genomic mosaic structure caused by C. jejuni introgression was demonstrated and basic phenotypic assays were used to characterize these isolates. The genomes of the three hybrid Campylobacter strains were sequenced using PacBio SMRT sequencing, followed by methylome analysis by Restriction-Modification Finder and genome analysis by Parsnp, Smash++ and blast. Additionally, the strains were phenotypically characterized with respect to growth behaviour, motility, eukaryotic cell invasion and adhesion, autoagglutination, biofilm formation, and water survival ability. Our analyses show that the three hybrid Campylobacter strains are clade 1 C . coli strains, which have acquired between 8.1 and 9.1 % of their genome from C. jejuni . The C. jejuni genomic segments acquired are distributed over the entire genome and do not form a coherent cluster. Most of the genes originating from C. jejuni are involved in chemotaxis and motility, membrane transport, cell signalling, or the resistance to toxic compounds such as bile acids. Interspecies gene transfer from C. jejuni has contributed 8.1–9.1% to the genome of three C. coli isolates and initiated the despeciation between C. jejuni and C. coli . Based on their functional annotation, the genes originating from C. jejuni enable the adaptation of the three strains to an intra-intestinal habitat. The transfer of a fused type II restriction-modification system that recognizes the CAYNNNNNCTC/GAGNNNNNRTG motif seems to be the key for the recombination of the C. jejuni genetic material with C. coli genomes.

Taxonomy-aware, sequence similarity ranking reliably predicts phage\u2013host relationships

BackgroundCharacterizing phage–host interactions is critical to understanding the ecological role of both partners and effective isolation of phage therapeuticals. Unfortunately, experimental methods for studying these interactions are markedly slow, low-throughput, and unsuitable for phages or hosts difficult to maintain in laboratory conditions. Therefore, a number of in silico methods emerged to predict prokaryotic hosts based on viral sequences. One of the leading approaches is the application of the BLAST tool that searches for local similarities between viral and microbial genomes. However, this prediction method has three major limitations: (i) top-scoring sequences do not always point to the actual host; (ii) mosaic virus genomes may match to many, typically related, bacteria; and (iii) viral and host sequences may diverge beyond the point where their relationship can be detected by a BLAST alignment.ResultsWe created an extension to BLAST, named Phirbo, that improves host prediction quality beyond what is obtainable from standard BLAST searches. The tool harnesses information concerning sequence similarity and bacteria relatedness to predict phage–host interactions. Phirbo was evaluated on three benchmark sets of known virus–host pairs, and it improved precision and recall by 11–40 percentage points over currently available, state-of-the-art, alignment-based, alignment-free, and machine-learning host prediction tools. Moreover, the discriminatory power of Phirbo for the recognition of virus–host relationships surpassed the results of other tools by at least 10 percentage points (area under the curve = 0.95), yielding a mean host prediction accuracy of 57% and 68% at the genus and family levels, respectively, and drops by 12 percentage points when using only a fraction of viral genome sequences (3 kb). Finally, we provide insights into a repertoire of protein and ncRNA genes that are shared between phages and hosts and may be prone to horizontal transfer during infection.ConclusionsOur results suggest that Phirbo is a simple and effective tool for predicting phage–host relationships.

Contribution of uniparental disomy in a clinical trio exome cohort of 2675 patients.

BackgroundUniparental disomy (UPD) is the inheritance of two homologous chromosomes from the same parent. UPD may result in clinical phenotypes when occurring on chromosomes with specific imprinting pattern, when leading to homozygosity of a deleterious recessive allele inherited from one carrier parent, or when associated with a mosaic aneuploidy. Due to the importance of UPD in genetic disease etiology, UPD analysis has started to be implemented in the context of exome sequencing (ES) or genome sequencing.MethodsWe developed an in‐house algorithm TRIPS (Trio Parentage/UPD Studies) to identify UPD events in trio ES cases. This method identifies regions with uniparental inheritance by utilizing the trio genotyping data obtained from the concurrent SNP array to delineate the parental origin of the SNPs in the proband.ResultsWe identified 16 UPD events from 2675 ES trios. Among those, four events led to imprinting disorders, seven unmasked a pathogenic/likely pathogenic variant in a recessive disease gene, and two were consistent with a mosaic genome wide paternal UPD pattern. Twelve of these UPD events directly contributed to the molecular diagnosis of the patients.ConclusionOur study demonstrated the contribution of UPD to the molecular diagnosis in one clinical ES cohort, thus UPD analysis should be incorporated into routine clinical ES interpretation.