Cartilaginous fishes offer unique insights into the evolution of the nuclear receptor gene repertoire in gnathostomes

Nuclear receptors (NRs) are key transcription factors that originated in the common ancestor of metazoans. The vast majority of NRs are triggered by binding to either endogenous (e.g. retinoic acid) or exogenous (e.g. xenobiotics) ligands, and their evolution and expansion is tightly linked to the function of endocrine systems. Importantly, they represent classic targets of physiological exploitation by endocrine disrupting chemicals. The NR gene repertoire in different lineages has been shaped by gene loss, duplication and mutation, denoting a dynamic evolutionary route. As the earliest diverging class of gnathostomes (jawed vertebrates), cartilaginous fishes offer an exceptional opportunity to address the early diversification of NR gene families and the evolution of the endocrine system in jawed vertebrates. Here we provide an exhaustive analysis into the NR gene composition in five elasmobranch (sharks and rays) and two holocephalan (chimaeras) species. For this purpose, we generated also a low coverage draft genome assembly of the chimaera small-eyed rabbitfish, Hydrolagus affinis. We show that cartilaginous fish retain an archetypal NR gene repertoire, similar to that of mammals and coincident with the two rounds of whole genome duplication that occurred in the gnathostome ancestor. Furthermore, novel gene members of the non-canonical NR0B receptors were found in the genomes of this lineage. Our findings provide an essential view into the early diversification of NRs in gnathostomes, paving the way for functional studies.

Sequencing the Genome of Indian Flying Fox, Natural Reservoir of Nipah Virus, Using Hybrid Assembly and Conservative Secondary Scaffolding

Indian fruit bats, flying fox Pteropus medius was identified as an asymptomatic natural host of recently emerged Nipah virus, which is known to induce a severe infectious disease in humans. The absence of P. medius genome sequence presents an important obstacle for further studies of virus–host interactions and better understanding of mechanisms of zoonotic viral emergence. Generation of the high-quality genome sequence is often linked to a considerable effort associated to elevated costs. Although secondary scaffolding methods have reduced sequencing expenses, they imply the development of new tools for the integration of different data sources to achieve more reliable sequencing results. We initially sequenced the P. medius genome using the combination of Illumina paired-end and Nanopore sequencing, with a depth of 57.4x and 6.1x, respectively. Then, we introduced the novel scaff2link software to integrate multiple sources of information for secondary scaffolding, allowing to remove the association with discordant information among two sources. Different quality metrics were next produced to validate the benefits from secondary scaffolding. The P. medius genome, assembled by this method, has a length of 1,985 Mb and consists of 33,613 contigs and 16,113 scaffolds with an NG50 of 19 Mb. At least 22.5% of the assembled sequences is covered by interspersed repeats already described in other species and 19,823 coding genes are annotated. Phylogenetic analysis demonstrated the clustering of P. medius genome with two other Pteropus bat species, P. alecto and P. vampyrus, for which genome sequences are currently available. SARS-CoV entry receptor ACE2 sequence of P. medius was 82.7% identical with ACE2 of Rhinolophus sinicus bats, thought to be the natural host of SARS-CoV. Altogether, our results confirm that a lower depth of sequencing is enough to obtain a valuable genome sequence, using secondary scaffolding approaches and demonstrate the benefits of the scaff2link application. The genome sequence is now available to the scientific community to (i) proceed with further genomic analysis of P. medius, (ii) to characterize the underlying mechanism allowing Nipah virus maintenance and perpetuation in its bat host, and (iii) to monitor their evolutionary pathways toward a better understanding of bats’ ability to control viral infections.

Gut microorganisms act together to exacerbate inflammation in spinal cords

Accumulating evidence indicates that gut microorganisms have a pathogenic role in autoimmune diseases, including in multiple sclerosis¹. Studies of experimental autoimmune encephalomyelitis (an animal model of multiple sclerosis)^2,3, as well as human studies^4-6, have implicated gut microorganisms in the development or severity of multiple sclerosis. However, it remains unclear how gut microorganisms act on the inflammation of extra-intestinal tissues such as the spinal cord. Here we show that two distinct signals from gut microorganisms coordinately activate autoreactive T cells in the small intestine that respond specifically to myelin oligodendrocyte glycoprotein (MOG). After induction of experimental autoimmune encephalomyelitis in mice, MOG-specific CD4⁺ T cells are observed in the small intestine. Experiments using germ-free mice that were monocolonized with microorganisms from the small intestine demonstrated that a newly isolated strain in the family Erysipelotrichaceae acts similarly to an adjuvant to enhance the responses of T helper 17 cells. Shotgun sequencing of the contents of the small intestine revealed a strain of Lactobacillus reuteri that possesses peptides that potentially mimic MOG. Mice that were co-colonized with these two strains showed experimental autoimmune encephalomyelitis symptoms that were more severe than those of germ-free or monocolonized mice. These data suggest that the synergistic effects that result from the presence of these microorganisms should be considered in the pathogenicity of multiple sclerosis, and that further study of these microorganisms may lead to preventive strategies for this disease.

Comparative genomic analysis of Flavobacteriaceae: insights into carbohydrate metabolism, gliding motility and secondary metabolite biosynthesis

BACKGROUND

Members of the bacterial family Flavobacteriaceae are widely distributed in the marine environment and often found associated with algae, fish, detritus or marine invertebrates. Yet, little is known about the characteristics that drive their ubiquity in diverse ecological niches. Here, we provide an overview of functional traits common to taxonomically diverse members of the family Flavobacteriaceae from different environmental sources, with a focus on the Marine clade. We include seven newly sequenced marine sponge-derived strains that were also tested for gliding motility and antimicrobial activity.

RESULTS

Comparative genomics revealed that genome similarities appeared to be correlated to 16S rRNA gene- and genome-based phylogeny, while differences were mostly associated with nutrient acquisition, such as carbohydrate metabolism and gliding motility. The high frequency and diversity of genes encoding polymer-degrading enzymes, often arranged in polysaccharide utilization loci (PULs), support the capacity of marine Flavobacteriaceae to utilize diverse carbon sources. Homologs of gliding proteins were widespread among all studied Flavobacteriaceae in contrast to members of other phyla, highlighting the particular presence of this feature within the Bacteroidetes. Notably, not all bacteria predicted to glide formed spreading colonies. Genome mining uncovered a diverse secondary metabolite biosynthesis arsenal of Flavobacteriaceae with high prevalence of gene clusters encoding pathways for the production of antimicrobial, antioxidant and cytotoxic compounds. Antimicrobial activity tests showed, however, that the phenotype differed from the genome-derived predictions for the seven tested strains.

CONCLUSIONS

Our study elucidates the functional repertoire of marine Flavobacteriaceae and highlights the need to combine genomic and experimental data while using the appropriate stimuli to unlock their uncharted metabolic potential.

Metagenomic approaches in microbial ecology: an update on whole-genome and marker gene sequencing analyses

Metagenomics and marker gene approaches, coupled with high-throughput sequencing technologies, have revolutionized the field of microbial ecology. Metagenomics is a culture-independent method that allows the identification and characterization of organisms from all kinds of samples. Whole-genome shotgun sequencing analyses the total DNA of a chosen sample to determine the presence of micro-organisms from all domains of life and their genomic content. Importantly, the whole-genome shotgun sequencing approach reveals the genomic diversity present, but can also give insights into the functional potential of the micro-organisms identified. The marker gene approach is based on the sequencing of a specific gene region. It allows one to describe the microbial composition based on the taxonomic groups present in the sample. It is frequently used to analyse the biodiversity of microbial ecosystems. Despite its importance, the analysis of metagenomic sequencing and marker gene data is quite a challenge. Here we review the primary workflows and software used for both approaches and discuss the current challenges in the field.

Genome Sequence Resource of Magnaporthe oryzae Laboratory Strain 2539

Magnaporthe oryzae causes blast disease on more than 50 species of monocot plants, including important crops such as rice, millet, and most recently wheat. Additionally, it is an important model system for studying host-pathogen interaction. Here, we report a high-quality genome assembly and annotation of a laboratory strain 2539 of M. oryzae, which is a widely used progeny of a rice-infecting isolate and a grass-infecting isolate. The genome sequence of strain 2539 will be useful for studying the evolution, host adaption, and pathogenicity of M. oryzae, which will be beneficial for a better understanding of the mechanisms of host-pathogen interaction.

More Accurate Transcript Assembly via Parameter Advising

Computational tools used for genomic analyses are becoming more accurate but also increasingly sophisticated and complex. This introduces a new problem in that these pieces of software have a large number of tunable parameters that often have a large influence on the results that are reported. We quantify the impact of parameter choice on transcript assembly and take some first steps toward generating a truly automated genomic analysis pipeline by developing a method for automatically choosing input-specific parameter values for reference-based transcript assembly using the Scallop tool. By choosing parameter values for each input, the area under the receiver operator characteristic curve (AUC) when comparing assembled transcripts to a reference transcriptome is increased by an average of 28.9% over using only the default parameter choices on 1595 RNA-Seq samples in the Sequence Read Archive. This approach is general, and when applied to StringTie, it increases the AUC by an average of 13.1% on a set of 65 RNA-Seq experiments from ENCODE. Parameter advisors for both Scallop and StringTie are available on Github.

Pathogenic Allodiploid Hybrids of Aspergillus Fungi

Interspecific hybridization substantially alters genotypes and phenotypes and can give rise to new lineages. Hybrid isolates that differ from their parental species in infection-relevant traits have been observed in several human-pathogenic yeasts and plant-pathogenic filamentous fungi but have yet to be found in human-pathogenic filamentous fungi. We discovered 6 clinical isolates from patients with aspergillosis originally identified as Aspergillus nidulans (section Nidulantes) that are actually allodiploid hybrids formed by the fusion of Aspergillus spinulosporus with an unknown close relative of Aspergillus quadrilineatus, both in section Nidulantes. Evolutionary genomic analyses revealed that these isolates belong to Aspergillus latus, an allodiploid hybrid species. Characterization of diverse infection-relevant traits further showed that A. latus hybrid isolates are genomically and phenotypically heterogeneous but also differ from A. nidulans, A. spinulosporus, and A. quadrilineatus. These results suggest that allodiploid hybridization contributes to the genomic and phenotypic diversity of filamentous fungal pathogens of humans.

Complete chloroplast genome sequence of Chrysosplenium ramosum and Chrysosplenium alternifolium (Saxifragaceae)

The complete chloroplast genome of Chrysosplenium ramosum Maxim. and Chrysosplenium alternifolium L. were reported in this study. The chloroplast genomes were 153,460 bp for C. ramosum and 152,619 bp for C. alternifolium. LSC and SSC of 83,670 bp and 17,342 bp were separated by two IRs of 26,224 bp each in C. ramosum. While C. alternifolium contained IRs of 25,992 bp, LSC of 83,524 bp and SSC of 17,111 bp. The chloroplast genome of C. ramosum contains 112 unique genes, including 79 protein-coding genes, four ribosomal RNA genes, and 30 transfer RNA genes. And the chloroplast genome of C. alternifolium contains 112 unique genes, including 79 protein-coding genes, four ribosomal RNA genes, and 30 transfer RNA genes. In addition, the rps12 gene was recognized as a trans-spliced gene and 17 intron-containing genes were also detected.

Detection of extended-spectrum beta-lactamase (ESBL) genes and plasmid replicons in Enterobacteriaceae using PlasmidSPAdes assembly of short-read sequence data

Knowledge of the epidemiology of plasmids is essential for understanding the evolution and spread of antimicrobial resistance. PlasmidSPAdes attempts to reconstruct plasmids using short-read sequence data. Accurate detection of extended-spectrum beta-lactamase (ESBL) genes and plasmid replicon genes is a prerequisite for the use of plasmid assembly tools to investigate the role of plasmids in the spread and evolution of ESBL production in Enterobacteriaceae. This study evaluated the performance of PlasmidSPAdes plasmid assembly for Enterobacteriaceae in terms of detection of ESBL-encoding genes, plasmid replicons and chromosomal wgMLST genes, and assessed the effect of k-mer size. Short-read sequence data for 59 ESBL-producing Enterobacteriaceae were assembled with PlasmidSPAdes using different k-mer sizes (21, 33, 55, 77, 99 and 127). For every k-mer size, the presence of ESBL genes, plasmid replicons and a selection of chromosomal wgMLST genes in the plasmid assembly was determined. Out of 241 plasmid replicons and 66 ESBL genes detected by whole-genome assembly, 213 plasmid replicons [88 %; 95 % confidence interval (CI): 83.9-91.9] and 43 ESBL genes (65 %; 95 % CI: 53.1-75.6) were detected in the plasmid assemblies obtained by PlasmidSPAdes. For most ESBL genes (83.3 %) and plasmid replicons (72.0 %), detection results did not differ between the k-mer sizes used in the plasmid assembly. No optimal k-mer size could be defined for the number of ESBL genes and plasmid replicons detected. For most isolates, the number of chromosomal wgMLST genes detected in the plasmid assemblies decreased with increasing k-mer size. Based on our findings, PlasmidSPAdes is not a suitable plasmid assembly tool for short-read sequence data for ESBL-encoding plasmids of Enterobacteriaceae.

Lateral Gene Transfer Shapes Diversity of Gardnerella spp

Gardnerella spp. are pathognomonic for bacterial vaginosis, which increases the risk of preterm birth and the transmission of sexually transmitted infections. Gardnerella spp. are genetically diverse, comprising what have recently been defined as distinct species with differing functional capacities. Disease associations with Gardnerella spp. are not straightforward: patients with BV are usually infected with multiple species, and Gardnerella spp. are also found in the vaginal microbiome of healthy women. Genome comparisons of Gardnerella spp. show evidence of lateral gene transfer (LGT), but patterns of LGT have not been characterized in detail. Here we sought to define the role of LGT in shaping the genetic structure of Gardnerella spp. We analyzed whole genome sequencing data for 106 Gardnerella strains and used these data for pan genome analysis and to characterize LGT in the core and accessory genomes, over recent and remote timescales. In our diverse sample of Gardnerella strains, we found that both the core and accessory genomes are clearly differentiated in accordance with newly defined species designations. We identified putative competence and pilus assembly genes across most species; we also found them to be differentiated between species. Competence machinery has diverged in parallel with the core genome, with selection against deleterious mutations as a predominant influence on their evolution. By contrast, the virulence factor vaginolysin, which encodes a toxin, appears to be readily exchanged among species. We identified five distinct prophage clusters in Gardnerella genomes, two of which appear to be exchanged between Gardnerella species. Differences among species are apparent in their patterns of LGT, including their exchange with diverse gene pools. Despite frequent LGT and co-localization in the same niche, our results show that Gardnerella spp. are clearly genetically differentiated and yet capable of exchanging specific genetic material. This likely reflects complex interactions within bacterial communities associated with the vaginal microbiome. Our results provide insight into how such interactions evolve and are maintained, allowing these multi-species communities to colonize and invade human tissues and adapt to antibiotics and other stressors.

The Basis for Natural Multiresistance to Phage in Pseudomonas aeruginosa

Pseudomonas aeruginosa is responsible for long-term infections and is particularly resistant to treatments when hiding inside the extracellular matrix or biofilms. Phage therapy might represent an alternative to antibiotic treatment, but up to 10% of clinical strains appear to resist multiple phages. We investigated the characteristics of P. aeruginosa clinical strains naturally resistant to phages and compared them to highly susceptible strains. The phage-resistant strains were defective in lipopolysaccharide (LPS) biosynthesis, were nonmotile and displayed an important degree of autolysis, releasing phages and pyocins. Complete genome sequencing of three resistant strains showed the existence of a large accessory genome made of multiple insertion elements, genomic islands, pyocins and prophages, including two phages performing lateral transduction. Mutations were found in genes responsible for the synthesis of LPS and/or type IV pilus, the major receptors for most phages. CRISPR-Cas systems appeared to be absent or inactive in phage-resistant strains, confirming that they do not play a role in the resistance to lytic phages but control the insertion of exogenous sequences. We show that, despite their apparent weakness, the multiphage-resistant strains described in this study displayed selective advantages through the possession of various functions, including weapons to eliminate other strains of the same or closely related species.

Genome mining of the citrus pathogen Elsinoë fawcettii; prediction and prioritisation of candidate effectors, cell wall degrading enzymes and secondary metabolite gene clusters

Elsinoë fawcettii, a necrotrophic fungal pathogen, causes citrus scab on numerous citrus varieties around the world. Known pathotypes of E. fawcettii are based on host range; additionally, cryptic pathotypes have been reported and more novel pathotypes are thought to exist. E. fawcettii produces elsinochrome, a non-host selective toxin which contributes to virulence. However, the mechanisms involved in potential pathogen-host interactions occurring prior to the production of elsinochrome are unknown, yet the host-specificity observed among pathotypes suggests a reliance upon such mechanisms. In this study we have generated a whole genome sequencing project for E. fawcettii, producing an annotated draft assembly 26.01 Mb in size, with 10,080 predicted gene models and low (0.37%) coverage of transposable elements. A small proportion of the assembly showed evidence of AT-rich regions, potentially indicating genomic regions with increased plasticity. Using a variety of computational tools, we mined the E. fawcettii genome for potential virulence genes as candidates for future investigation. A total of 1,280 secreted proteins and 276 candidate effectors were predicted and compared to those of other necrotrophic (Botrytis cinerea, Parastagonospora nodorum, Pyrenophora tritici-repentis, Sclerotinia sclerotiorum and Zymoseptoria tritici), hemibiotrophic (Leptosphaeria maculans, Magnaporthe oryzae, Rhynchosporium commune and Verticillium dahliae) and biotrophic (Ustilago maydis) plant pathogens. Genomic and proteomic features of known fungal effectors were analysed and used to guide the prioritisation of 120 candidate effectors of E. fawcettii. Additionally, 378 carbohydrate-active enzymes were predicted and analysed for likely secretion and sequence similarity with known virulence genes. Furthermore, secondary metabolite prediction indicated nine additional genes potentially involved in the elsinochrome biosynthesis gene cluster than previously described. A further 21 secondary metabolite clusters were predicted, some with similarity to known toxin producing gene clusters. The candidate virulence genes predicted in this study provide a comprehensive resource for future experimental investigation into the pathogenesis of E. fawcettii.

An in silico survey of Clostridioides difficile extrachromosomal elements

The Gram-positive enteropathogen Clostridioides difficile (Clostridium difficile) is the major cause of healthcare-associated diarrhoea and is also an important cause of community-acquired infectious diarrhoea. Considering the burden of the disease, many studies have employed whole-genome sequencing of bacterial isolates to identify factors that contribute to virulence and pathogenesis. Though extrachromosomal elements (ECEs) such as plasmids are important for these processes in other bacteria, the few characterized plasmids of C. difficile have no relevant functions assigned and no systematic identification of plasmids has been carried out to date. Here, we perform an in silico analysis of publicly available sequence data to show that ~13 % of all C. difficile strains contain ECEs, with 1–6 elements per strain. Our approach identifies known plasmids (e.g. pCD6, pCD630 and cloning plasmids) and six novel putative plasmid families. Our study shows that plasmids are abundant and may encode functions that are relevant for C. difficile physiology. The newly identified plasmids may also form the basis for the construction of novel cloning plasmids for C. difficile that are compatible with existing tools.

Identification of gene fusion events in Mycobacterium tuberculosis that encode chimeric proteins

Mycobacterium tuberculosis is a facultative intracellular pathogen responsible for causing tuberculosis. The harsh environment in which M. tuberculosis survives requires this pathogen to continuously adapt in order to maintain an evolutionary advantage. However, the apparent absence of horizontal gene transfer in M. tuberculosis imposes restrictions in the ways by which evolution can occur. Large-scale changes in the genome can be introduced through genome reduction, recombination events and structural variation. Here, we identify a functional chimeric protein in the ppe38-71 locus, the absence of which is known to have an impact on protein secretion and virulence. To examine whether this approach was used more often by this pathogen, we further develop software that detects potential gene fusion events from multigene deletions using whole genome sequencing data. With this software we could identify a number of other putative gene fusion events within the genomes of M. tuberculosis isolates. We were able to demonstrate the expression of one of these gene fusions at the protein level using mass spectrometry. Therefore, gene fusions may provide an additional means of evolution for M. tuberculosis in its natural environment whereby novel chimeric proteins and functions can arise.

Identification of key genes related to seedlessness by genome-wide detection of structural variation and transcriptome analysis in 'Shijiwuhe' pear

Seedless fruits are highly marketable because they are easier to eat than fruits with seeds. 'Shijiwuhe' is a seedless pear cultivar that is a mutant derived from an F1 hybridization population ('Bartlett' x 'Yali'). Little is known about the key genes controlling seedless pear fruit. In this study, field experiments revealed that seedless 'Shijiwuhe' pear was not due to parthenocarpy, and that it was self-incompatible. Single nucleotide polymorphisms (SNPs), small insertions and deletions (InDels) and structural variations (SVs) were characterized using DNA sequencing data between 'Shijiwuhe' and parental cultivars. A total of 1498 genes were found to be affected by SV and over 50% of SVs were located in promoter regions. Transcriptome analysis was conducted at three time points (4, 8, and 12 days after cross-pollination) during early fruit development of 'Shijiwuhe', 'Bartlett', and 'Yali'. In total, 1438 differentially expressed genes (DEGs) were found between 'Shijiwuhe' and parental cultivars 'Bartlett' and 'Yali'. We found 1193 SVs that caused differential expression of genes at 4 DACP. Among them, over 100 genes were in pathways related to seed nutrition and energy storage and 41 candidate genes encoded several important transcription factors, such as MYB, WRKY, NAC, and bHLH, which might play important roles in seed development. The qRT-PCR results also confirmed that the candidate genes with SVs showed differential expression between 'Shijiwuhe' pear and 'Bartlett' or 'Yali'. This study, which combined field experiments, SV detection, and transcriptome analysis might provide an effective way to predict the candidate genes regulating the seedless trait and important gene resources for genetic improvement of pear.

Deconvolute individual genomes from metagenome sequences through short read clustering

Metagenome assembly from short next-generation sequencing data is a challenging process due to its large scale and computational complexity. Clustering short reads by species before assembly offers a unique opportunity for parallel downstream assembly of genomes with individualized optimization. However, current read clustering methods suffer either false negative (under-clustering) or false positive (over-clustering) problems. Here we extended our previous read clustering software, SpaRC, by exploiting statistics derived from multiple samples in a dataset to reduce the under-clustering problem. Using synthetic and real-world datasets we demonstrated that this method has the potential to cluster almost all of the short reads from genomes with sufficient sequencing coverage. The improved read clustering in turn leads to improved downstream genome assembly quality.

Nosocomial transmission of extensively drug resistant Acinetobacter baumannii strains in a tertiary level hospital

Acinetobacter baumannii is an opportunistic infectious agent that affects primarily immunocompromised individuals. A. baumannii is highly prevalent in hospital settings being commonly associated with nosocomial transmission and drug resistance. Here, we report the identification and genetic characterization of A. baumannii strains among patients in a tertiary level hospital in Mexico. Whole genome sequencing analysis was performed to establish their genetic relationship and drug resistance mutations profile. Ten genetically different, extensively drug resistant strains were identified circulating among seven wards. The genetic profiles showed resistance primarily against aminoglycosides and beta-lactam antibiotics. Importantly, no mutants conferring resistance to colistin were observed. The results highlight the importance of implementing robust classification schemes for advanced genetic characterization of A. baumannii clinical isolates and simultaneous detection of drug resistance markers for adequate patient's management in clinical settings.

A High-Quality Genome Assembly of the North American Song Sparrow, Melospiza melodia

The song sparrow, Melospiza melodia, is one of the most widely distributed species of songbirds found in North America. It has been used in a wide range of behavioral and ecological studies. This species' pronounced morphological and behavioral diversity across populations makes it a favorable candidate in several areas of biomedical research. We have generated a high-quality de novo genome assembly of M. melodia using Illumina short read sequences from genomic and in vitro proximity-ligation libraries. The assembled genome is 978.3 Mb, with a physical coverage of 24.9×, N50 scaffold size of 5.6 Mb and N50 contig size of 31.7 Kb. Our genome assembly is highly complete, with 87.5% full-length genes present out of a set of 4,915 universal single-copy orthologs present in most avian genomes. We annotated our genome assembly and constructed 15,086 gene models, a majority of which have high homology to related birds, Taeniopygia guttata and Junco hyemalis In total, 83% of the annotated genes are assigned with putative functions. Furthermore, only ∼7% of the genome is found to be repetitive; these regions and other non-coding functional regions are also identified. The high-quality M. melodia genome assembly and annotations we report will serve as a valuable resource for facilitating studies on genome structure and evolution that can contribute to biomedical research and serve as a reference in population genomic and comparative genomic studies of closely related species.

GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes

An important assessment prior to genome assembly and related analyses is genome profiling, where the k-mer frequencies within raw sequencing reads are analyzed to estimate major genome characteristics such as size, heterozygosity, and repetitiveness. Here we introduce GenomeScope 2.0 (https://github.com/tbenavi1/genomescope2.0), which applies combinatorial theory to establish a detailed mathematical model of how k-mer frequencies are distributed in heterozygous and polyploid genomes. We describe and evaluate a practical implementation of the polyploid-aware mixture model that quickly and accurately infers genome properties across thousands of simulated and several real datasets spanning a broad range of complexity. We also present a method called Smudgeplot (https://github.com/KamilSJaron/smudgeplot) to visualize and estimate the ploidy and genome structure of a genome by analyzing heterozygous k-mer pairs. We successfully apply the approach to systems of known variable ploidy levels in the Meloidogyne genus and the extreme case of octoploid Fragaria × ananassa.

Gene clustering and copy number variation in alkaloid metabolic pathways of opium poppy

Genes in plant secondary metabolic pathways enable biosynthesis of a range of medically and industrially important compounds, and are often clustered on chromosomes. Here, we study genomic clustering in the benzylisoquinoline alkaloid (BIA) pathway in opium poppy (Papaver somniferum), exploring relationships between gene expression, copy number variation, and metabolite production. We use Hi-C to improve the existing draft genome assembly, yielding chromosome-scale scaffolds that include 35 previously unanchored BIA genes. We find that co-expression of BIA genes increases within clusters and identify candidates with unknown function based on clustering and covariation in expression and alkaloid production. Copy number variation in critical BIA genes correlates with stark differences in alkaloid production, linking noscapine production with an 11-gene deletion, and increased thebaine/decreased morphine production with deletion of a T6ODM cluster. Our results show that the opium poppy genome is still dynamically evolving in ways that contribute to medically and industrially important phenotypes.

Whole genome sequencing and comparative genomic analysis of oleaginous red yeast Sporobolomyces pararoseus NGR identifies candidate genes for biotechnological potential and ballistospores-shooting

BACKGROUND

Sporobolomyces pararoseus is regarded as an oleaginous red yeast, which synthesizes numerous valuable compounds with wide industrial usages. This species hold biotechnological interests in biodiesel, food and cosmetics industries. Moreover, the ballistospores-shooting promotes the colonizing of S. pararoseus in most terrestrial and marine ecosystems. However, very little is known about the basic genomic features of S. pararoseus. To assess the biotechnological potential and ballistospores-shooting mechanism of S. pararoseus on genome-scale, the whole genome sequencing was performed by next-generation sequencing technology.

RESULTS

Here, we used Illumina Hiseq platform to firstly assemble S. pararoseus genome into 20.9 Mb containing 54 scaffolds and 5963 predicted genes with a N50 length of 2,038,020 bp and GC content of 47.59%. Genome completeness (BUSCO alignment: 95.4%) and RNA-seq analysis (expressed genes: 98.68%) indicated the high-quality features of the current genome. Through the annotation information of the genome, we screened many key genes involved in carotenoids, lipids, carbohydrate metabolism and signal transduction pathways. A phylogenetic assessment suggested that the evolutionary trajectory of the order Sporidiobolales species was evolved from genus Sporobolomyces to Rhodotorula through the mediator Rhodosporidiobolus. Compared to the lacking ballistospores Rhodotorula toruloides and Saccharomyces cerevisiae, we found genes enriched for spore germination and sugar metabolism. These genes might be responsible for the ballistospores-shooting in S. pararoseus NGR.

CONCLUSION

These results greatly advance our understanding of S. pararoseus NGR in biotechnological potential and ballistospores-shooting, which help further research of genetic manipulation, metabolic engineering as well as its evolutionary direction.

Evidence of Recent Genomic Evolution in Gonococcal Strains With Decreased Susceptibility to Cephalosporins or Azithromycin in the United States, 2014-2016

BACKGROUND

Given the lack of new antimicrobials or a vaccine, understanding the evolutionary dynamics of Neisseria gonorrhoeae is a significant public and global health priority. We investigated the emergence and spread of gonococcal strains with decreased susceptibility to cephalosporins and azithromycin using detailed genomic analyses of gonococcal isolates collected in the United States, 2014-2016.

METHODS

We sequenced genomes of 649 isolates collected through the Gonococcal Isolate Surveillance Project. We examined the genetic relatedness of isolates and assessed associations between clades and various genotypic and phenotypic combinations.

RESULTS

We identified a large and clonal lineage of strains (MLST ST9363) associated with elevated azithromycin minimum inhibitory concentration (AZIem), characterized by a mosaic mtr locus (C substitution in the mtrR promoter, mosaic mtrR and mtrD). Mutations in 23S rRNA were sporadically distributed among AZIem strains. Another clonal group (MLST ST1901) possessed 7 unique PBP2 patterns, and it shared common mutations in other genes associated with cephalosporin resistance.

CONCLUSIONS

Whole-genome sequencing methods can enhance monitoring of antimicrobial resistant gonococcal strains by identifying gonococcal populations containing mutations of concern. These methods could inform the development of point-of-care diagnostic tests designed to determine the specific antibiotic susceptibility profile of a gonococcal infection in a patient.

Comparative Genomic Characterization of the Multimammate Mouse Mastomys coucha

Mastomys are the most widespread African rodent and carriers of various diseases such as the plague or Lassa virus. In addition, mastomys have rapidly gained a large number of mammary glands. Here, we generated a genome, variome, and transcriptomes for Mastomys coucha. As mastomys diverged at similar times from mouse and rat, we demonstrate their utility as a comparative genomic tool for these commonly used animal models. Furthermore, we identified over 500 mastomys accelerated regions, often residing near important mammary developmental genes or within their exons leading to protein sequence changes. Functional characterization of a noncoding mastomys accelerated region, located in the HoxD locus, showed enhancer activity in mouse developing mammary glands. Combined, our results provide genomic resources for mastomys and highlight their potential both as a comparative genomic tool and for the identification of mammary gland number determining factors.

Genome Assembly and Analysis of the North American Mountain Goat (Oreamnos americanus) Reveals Species-Level Responses to Extreme Environments

The North American mountain goat (Oreamnos americanus) is an iconic alpine species that faces stressors from climate change, industrial development, and recreational activities. This species’ phylogenetic position within the Caprinae lineage has not been resolved and their phylogeographic history is dynamic and controversial. Genomic data could be used to address these questions and provide valuable insights to conservation and management initiatives. We sequenced short-read genomic libraries constructed from a DNA sample of a 2.5-year-old female mountain goat at 80X coverage. We improved the short-read assembly by generating Chicago library data and scaffolding using the HiRise approach. The final assembly was 2,506 Mbp in length with an N50 of 66.6 Mbp, which is within the length range and in the upper quartile for N50 published ungulate genome assemblies. Comparative analysis identified 84 gene families unique to the mountain goat. The species demographic history in terms of effective population size generally mirrored climatic trends over the past one hundred thousand years and showed a sharp decline during the last glacial maximum. This genome assembly will provide a reference basis for future population and comparative genomic analyses.

First Report of Isolation and Genome Sequence of L. petauri Strain from a Rainbow Trout Lactococcosis Outbreak

Lactococcosis is a disease encountered in a wide variety of fish species causing mortalities and having great economic impact on farmed fish. In this study, we report for the first time the isolation of a strain of the recently described novel species Lactococcus petauri, from rainbow trout suffering from lactococcosis. The aim of this study was to determine the complete genome sequence of L. petauri strain LG_SAV_20 and to characterize its antimicrobial resistance and virulence. The genome of L. petauri LG_SAV_20 consists of 2,078,949 base pair (bp) with a GC content of 38.05%, 1950 predicted coding sequence (CDS), and 60 RNAs (51 tRNAs, 3 ncRNAs, and 6 rRNAs). Phylogenetic analysis revealed that L. petauri LG_SAV_20 shares most of its genome with L. garvieae strains isolated from rainbow trout. Detection of genes associated with antimicrobial resistance indicated that the isolate possesses the multidrug transporter mdt(A) gene, while using comparative analysis we identified several genes that might be related to bacterial pathogenesis. This genomic information provides new insights into the role of this novel species as an etiological agent of lactococcosis.

Whole-genome sequencing of wild Siberian musk deer (Moschus moschiferus) provides insights into its genetic features

Background

Siberian musk deer, one of the seven species, is distributed in coniferous forests of Asia. Worldwide, the population size of Siberian musk deer is threatened by severe illegal poaching for commercially valuable musk and meat, habitat losses, and forest fire. At present, this species is categorized as Vulnerable on the IUCN Red List. However, the genetic information of Siberian musk deer is largely unexplored.

Results

Here, we produced 3.10 Gb draft assembly of wild Siberian musk deer with a contig N50 of 29,145 bp and a scaffold N50 of 7,955,248 bp. We annotated 19,363 protein-coding genes and estimated 44.44% of the genome to be repetitive. Our phylogenetic analysis reveals that wild Siberian musk deer is closer to Bovidae than to Cervidae. Comparative analyses showed that the genetic features of Siberian musk deer adapted in cold and high-altitude environments. We sequenced two additional genomes of Siberian musk deer constructed demographic history indicated that changes in effective population size corresponded with recent glacial epochs. Finally, we identified several candidate genes that may play a role in the musk secretion based on transcriptome analysis.

Conclusions

Here, we present a high-quality draft genome of wild Siberian musk deer, which will provide a valuable genetic resource for further investigations of this economically important musk deer.

Enrichment of Anaerobic Syngas-Converting Communities and Isolation of a Novel Carboxydotrophic Acetobacterium wieringae Strain JM

Syngas is a substrate for the anaerobic bioproduction of fuels and valuable chemicals. In this study, anaerobic sludge was used for microbial enrichments with synthetic syngas and acetate as main substrates. The objectives of this study were to identify microbial networks (in enrichment cultures) for the conversion of syngas to added-value products, and to isolate robust, non-fastidious carboxydotrophs. Enrichment cultures produced methane and propionate, this last one an unusual product from syngas fermentation. A bacterium closely related to Acetobacterium wieringae was identified as most prevalent (87% relative abundance) in the enrichments. Methanospirillum sp. and propionate-producing bacteria clustering within the genera Anaerotignum and Pelobacter were also found. Further on, strain JM, was isolated and was found to be 99% identical (16S rRNA gene) to A. wieringae DSM 1911^T. Digital DNA-DNA hybridization (dDDH) value between the genomes of strain JM and A. wieringae was 77.1%, indicating that strain JM is a new strain of A. wieringae. Strain JM can grow on carbon monoxide (100% CO, total pressure 170 kPa) without yeast extract or formate, producing mainly acetate. Remarkably, conversion of CO by strain JM showed shorter lag phase than in cultures of A. wieringae DSM 1911^T, and about four times higher amount of CO was consumed in 7 days. Genome analysis suggests that strain JM uses the Wood-Ljungdahl pathway for the conversion of one carbon compounds (CO, formate, CO₂/H₂). Genes encoding bifurcational enzyme complexes with similarity to the bifurcational formate dehydrogenase (Fdh) of Clostridium autoethanogenum are present, and possibly relate to the higher tolerance to CO of strain JM compared to other Acetobacterium species. A. wieringae DSM 1911^T grew on CO in medium containing 1 mM formate.

High quality 3C de novo assembly and annotation of a multidrug resistant ST-111 Pseudomonas aeruginosa genome: Benchmark of hybrid and non-hybrid assemblers

Genotyping methods and genome sequencing are indispensable to reveal genomic structure of bacterial species displaying high level of genome plasticity. However, reconstruction of genome or assembly is not straightforward due to data complexity, including repeats, mobile and accessory genetic elements of bacterial genomes. Moreover, since the solution to this problem is strongly influenced by sequencing technology, bioinformatics pipelines, and selection criteria to assess assemblers, there is no systematic way to select a priori the optimal assembler and parameter settings. To assembly the genome of Pseudomonas aeruginosa strain AG1 (PaeAG1), short reads (Illumina) and long reads (Oxford Nanopore) sequencing data were used in 13 different non-hybrid and hybrid approaches. PaeAG1 is a multiresistant high-risk sequence type 111 (ST-111) clone that was isolated from a Costa Rican hospital and it was the first report of an isolate of P. aeruginosa carrying both blaVIM-2 and blaIMP-18 genes encoding for metallo-β-lactamases (MBL) enzymes. To assess the assemblies, multiple metrics regard to contiguity, correctness and completeness (3C criterion, as we define here) were used for benchmarking the 13 approaches and select a definitive assembly. In addition, annotation was done to identify genes (coding and RNA regions) and to describe the genomic content of PaeAG1. Whereas long reads and hybrid approaches showed better performances in terms of contiguity, higher correctness and completeness metrics were obtained for short read only and hybrid approaches. A manually curated and polished hybrid assembly gave rise to a single circular sequence with 100% of core genes and known regions identified, >98% of reads mapped back, no gaps, and uniform coverage. The strategy followed to obtain this high-quality 3C assembly is detailed in the manuscript and we provide readers with an all-in-one script to replicate our results or to apply it to other troublesome cases. The final 3C assembly revealed that the PaeAG1 genome has 7,190,208 bp, a 65.7% GC content and 6,709 genes (6,620 coding sequences), many of which are included in multiple mobile genomic elements, such as 57 genomic islands, six prophages, and two complete integrons with blaVIM-2 and blaIMP-18 MBL genes. Up to 250 and 60 of the predicted genes are anticipated to play a role in virulence (adherence, quorum sensing and secretion) or antibiotic resistance (β-lactamases, efflux pumps, etc). Altogether, the assembly and annotation of the PaeAG1 genome provide new perspectives to continue studying the genomic diversity and gene content of this important human pathogen.

Draft Genome Sequences of Three Strains of Verticillium nonalfalfae Exhibiting Different Levels of Aggressiveness on Ailanthus altissima

Verticillium nonalfalfae, a soilborne vascular fungus, shows promise for biocontrol of highly invasive Ailanthus altissima strains. This announcement provides draft genome sequences of the aggressive isolate G1/5 (wild-type strain), the highly aggressive isolate Vert56 (improved strain), and the mildly aggressive isolate I3/2, all obtained from symptomatic A. altissima trees in Austria.

Verticillium nonalfalfae, a soilborne vascular fungus, shows promise for biocontrol of highly invasive Ailanthus altissima strains. This announcement provides draft genome sequences of the aggressive isolate G1/5 (wild-type strain), the highly aggressive isolate Vert56 (improved strain), and the mildly aggressive isolate I3/2, all obtained from symptomatic A. altissima trees in Austria.

Organohalide-respiring Desulfoluna species isolated from marine environments

The genus Desulfoluna comprises two anaerobic sulfate-reducing strains, D. spongiiphila AA1^T and D. butyratoxydans MSL71^T, of which only the former was shown to perform organohalide respiration (OHR). Here we isolated a third strain, designated D. spongiiphila strain DBB, from marine intertidal sediment using 1,4-dibromobenzene and sulfate as the electron acceptors and lactate as the electron donor. Each strain harbors three reductive dehalogenase gene clusters (rdhABC) and corrinoid biosynthesis genes in their genomes, and dehalogenated brominated but not chlorinated organohalogens. The Desulfoluna strains maintained OHR in the presence of 20 mM sulfate or 20 mM sulfide, which often negatively affect other organohalide-respiring bacteria. Strain DBB sustained OHR with 2% oxygen in the gas phase, in line with its genetic potential for reactive oxygen species detoxification. Reverse transcription-quantitative PCR revealed differential induction of rdhA genes in strain DBB in response to 1,4-dibromobenzene or 2,6-dibromophenol. Proteomic analysis confirmed expression of rdhA1 with 1,4-dibromobenzene, and revealed a partially shared electron transport chain from lactate to 1,4-dibromobenzene and sulfate, which may explain accelerated OHR during concurrent sulfate reduction. Versatility in using electron donors, de novo corrinoid biosynthesis, resistance to sulfate, sulfide and oxygen, and concurrent sulfate reduction and OHR may confer an advantage to marine Desulfoluna strains.

The complete mitochondrial genome of a Chinese endemic species, Paraprenanthes diversifolia (Vaniot) N. Kilian (Cichorieae; Asteraceae)

The first complete mitochondrial genome (mt) of Paraprenanthes diversifolia (Vaniot) N. Kilian (Cichorieae; Asteraceae) was sequenced and successfully assembled in this study. The full length of the mt genome is 360,751 bp, containing 73 genes (33 protein-coding genes, 29 tRNA genes, 6 rRNA genes, and 5 protein-coding genes containing internal stop codons). There are two pairs of long (over 1000 bp) repeat regions in the mt genome of P. diversifolia. The phylogenetic analysis indicated that P. diversifolia has a close relationship with other Lactucinae species.

A draft genome sequence of the elusive giant squid, Architeuthis dux

BACKGROUND

The giant squid (Architeuthis dux; Steenstrup, 1857) is an enigmatic giant mollusc with a circumglobal distribution in the deep ocean, except in the high Arctic and Antarctic waters. The elusiveness of the species makes it difficult to study. Thus, having a genome assembled for this deep-sea-dwelling species will allow several pending evolutionary questions to be unlocked.

FINDINGS

We present a draft genome assembly that includes 200 Gb of Illumina reads, 4 Gb of Moleculo synthetic long reads, and 108 Gb of Chicago libraries, with a final size matching the estimated genome size of 2.7 Gb, and a scaffold N50 of 4.8 Mb. We also present an alternative assembly including 27 Gb raw reads generated using the Pacific Biosciences platform. In addition, we sequenced the proteome of the same individual and RNA from 3 different tissue types from 3 other species of squid (Onychoteuthis banksii, Dosidicus gigas, and Sthenoteuthis oualaniensis) to assist genome annotation. We annotated 33,406 protein-coding genes supported by evidence, and the genome completeness estimated by BUSCO reached 92%. Repetitive regions cover 49.17% of the genome.

CONCLUSIONS

This annotated draft genome of A. dux provides a critical resource to investigate the unique traits of this species, including its gigantism and key adaptations to deep-sea environments.

The Genome of the Blind Soil-Dwelling and Ancestrally Wingless Dipluran Campodea augens: A Key Reference Hexapod for Studying the Emergence of Insect Innovations

The dipluran two-pronged bristletail Campodea augens is a blind ancestrally wingless hexapod with the remarkable capacity to regenerate lost body appendages such as its long antennae. As sister group to Insecta (sensu stricto), Diplura are key to understanding the early evolution of hexapods and the origin and evolution of insects. Here we report the 1.2-Gb draft genome of C. augens and results from comparative genomic analyses with other arthropods. In C. augens, we uncovered the largest chemosensory gene repertoire of ionotropic receptors in the animal kingdom, a massive expansion that might compensate for the loss of vision. We found a paucity of photoreceptor genes mirroring at the genomic level the secondary loss of an ancestral external photoreceptor organ. Expansions of detoxification and carbohydrate metabolism gene families might reflect adaptations for foraging behavior, and duplicated apoptotic genes might underlie its high regenerative potential. The C. augens genome represents one of the key references for studying the emergence of genomic innovations in insects, the most diverse animal group, and opens up novel opportunities to study the under-explored biology of diplurans.

Fast Approximation of Frequent k-Mers and Applications to Metagenomics

Estimating the abundances of all k-mers in a set of biological sequences is a fundamental and challenging problem with many applications in biological analysis. Although several methods have been designed for the exact or approximate solution of this problem, they all require to process the entire data set, which can be extremely expensive for high-throughput sequencing data sets. Although in some applications it is crucial to estimate all k-mers and their abundances, in other situations it may be sufficient to report only frequent k-mers, which appear with relatively high frequency in a data set. This is the case, for example, in the computation of k-mers' abundance-based distances among data sets of reads, commonly used in metagenomic analyses. In this study, we develop, analyze, and test a sampling-based approach, called Sampling Algorithm for K-mErs approxIMAtion (SAKEIMA), to approximate the frequent k-mers and their frequencies in a high-throughput sequencing data set while providing rigorous guarantees on the quality of the approximation. SAKEIMA employs an advanced sampling scheme and we show how the characterization of the Vapnik-Chervonenkis dimension, a core concept from statistical learning theory, of a properly defined set of functions leads to practical bounds on the sample size required for a rigorous approximation. Our experimental evaluation shows that SAKEIMA allows to rigorously approximate frequent k-mers by processing only a fraction of a data set and that the frequencies estimated by SAKEIMA lead to accurate estimates of k-mer-based distances between high-throughput sequencing data sets. Overall, SAKEIMA is an efficient and rigorous tool to estimate k-mers' abundances providing significant speedups in the analysis of large sequencing data sets.

The C-Fern (Ceratopteris richardii) genome: insights into plant genome evolution with the first partial homosporous fern genome assembly

Ferns are notorious for possessing large genomes and numerous chromosomes. Despite decades of speculation, the processes underlying the expansive genomes of ferns are unclear, largely due to the absence of a sequenced homosporous fern genome. The lack of this crucial resource has not only hindered investigations of evolutionary processes responsible for the unusual genome characteristics of homosporous ferns, but also impeded synthesis of genome evolution across land plants. Here, we used the model fern species Ceratopteris richardii to address the processes (e.g., polyploidy, spread of repeat elements) by which the large genomes and high chromosome numbers typical of homosporous ferns may have evolved and have been maintained. We directly compared repeat compositions in species spanning the green plant tree of life and a diversity of genome sizes, as well as both short- and long-read-based assemblies of Ceratopteris. We found evidence consistent with a single ancient polyploidy event in the evolutionary history of Ceratopteris based on both genomic and cytogenetic data, and on repeat proportions similar to those found in large flowering plant genomes. This study provides a major stepping-stone in the understanding of land plant evolutionary genomics by providing the first homosporous fern reference genome, as well as insights into the processes underlying the formation of these massive genomes.

Novel redox-active enzymes for ligninolytic applications revealed from multiomics analyses of Peniophora sp. CBMAI 1063, a laccase hyper-producer strain

The repertoire of redox-active enzymes produced by the marine fungus Peniophora sp. CBMAI 1063, a laccase hyper-producer strain, was characterized by omics analyses. The genome revealed 309 Carbohydrate-Active Enzymes (CAZymes) genes, including 48 predicted genes related to the modification and degradation of lignin, whith 303 being transcribed under cultivation in optimized saline conditions for laccase production. The secretome confirmed that the fungus can produce a versatile ligninolytic enzyme cocktail. It secretes 56 CAZymes, including 11 oxidative enzymes classified as members of auxiliary activity families (AAs), comprising two laccases, Pnh_Lac1 and Pnh_Lac2, the first is the major secretory protein of the fungi. The Pnh_Lac1-mediator system was able to promote the depolymerization of lignin fragments and polymeric lignin removal from pretreated sugarcane bagasse, confirming viability of this fungus enzymatic system for lignocellulose-based bioproducts applications.

Ecological and functional capabilities of an uncultured Kordia sp

Cultivable bacteria represent only a fraction of the diversity in microbial communities. However, the official procedures for classification and characterization of a novel prokaryotic species still rely on isolates. Nevertheless, due to single cell genomics, it is possible to retrieve genomes from environmental samples by sequencing them individually, and to assign specific genes to a specific taxon, regardless of their ability to grow in culture. In this study, a complete description was performed for uncultured Kordia sp. TARA_039_SRF, a proposed novel species within the genus Kordia, using culture-independent techniques. The type material was a high-quality draft genome (94.97% complete, 4.65% gene redundancy) co-assembled using ten nearly identical single amplified genomes (SAGs) from surface seawater in the North Indian Ocean during the Tara Oceans Expedition. The assembly process was optimized to obtain the best possible assembly metrics and a less fragmented genome. The closest relative of the species was Kordia periserrulae, which shared 97.56% similarity of the 16S rRNA gene, 75% orthologs and 89.13% average nucleotide identity. The functional potential of the proposed novel species included proteorhodopsin, the ability to incorporate nitrate, cytochrome oxidases with high affinity for oxygen, and CAZymes that were unique features within the genus. Its abundance at different depths and size fractions was also evaluated together with its functional annotation, revealing that its putative ecological niche could be particles of phytoplanktonic origin. It could putatively attach to these particles and consume them while sinking to the deeper and oxygen depleted layers of the North Indian Ocean.

The apricot (Prunus armeniaca L.) genome elucidates Rosaceae evolution and beta-carotenoid synthesis

Apricots, scientifically known as Prunus armeniaca L, are drupes that resemble and are closely related to peaches or plums. As one of the top consumed fruits, apricots are widely grown worldwide except in Antarctica. A high-quality reference genome for apricot is still unavailable, which has become a handicap that has dramatically limited the elucidation of the associations of phenotypes with the genetic background, evolutionary diversity, and population diversity in apricot. DNA from P. armeniaca was used to generate a standard, size-selected library with an average DNA fragment size of ~20 kb. The library was run on Sequel SMRT Cells, generating a total of 16.54 Gb of PacBio subreads (N50 = 13.55 kb). The high-quality P. armeniaca reference genome presented here was assembled using long-read single-molecule sequencing at approximately 70× coverage and 171× Illumina reads (40.46 Gb), combined with a genetic map for chromosome scaffolding. The assembled genome size was 221.9 Mb, with a contig NG50 size of 1.02 Mb. Scaffolds covering 92.88% of the assembled genome were anchored on eight chromosomes. Benchmarking Universal Single-Copy Orthologs analysis showed 98.0% complete genes. We predicted 30,436 protein-coding genes, and 38.28% of the genome was predicted to be repetitive. We found 981 contracted gene families, 1324 expanded gene families and 2300 apricot-specific genes. The differentially expressed gene (DEG) analysis indicated that a change in the expression of the 9-cis-epoxycarotenoid dioxygenase (NCED) gene but not lycopene beta-cyclase (LcyB) gene results in a low β-carotenoid content in the white cultivar "Dabaixing". This complete and highly contiguous P. armeniaca reference genome will be of help for future studies of resistance to plum pox virus (PPV) and the identification and characterization of important agronomic genes and breeding strategies in apricot.

Athena: Automated Tuning of k-mer based Genomic Error Correction Algorithms using Language Models

The performance of most error-correction (EC) algorithms that operate on genomics reads is dependent on the proper choice of its configuration parameters, such as the value of k in k-mer based techniques. In this work, we target the problem of finding the best values of these configuration parameters to optimize error correction and consequently improve genome assembly. We perform this in an adaptive manner, adapted to different datasets and to EC tools, due to the observation that different configuration parameters are optimal for different datasets, i.e., from different platforms and species, and vary with the EC algorithm being applied. We use language modeling techniques from the Natural Language Processing (NLP) domain in our algorithmic suite, Athena, to automatically tune the performance-sensitive configuration parameters. Through the use of N-Gram and Recurrent Neural Network (RNN) language modeling, we validate the intuition that the EC performance can be computed quantitatively and efficiently using the “perplexity” metric, repurposed from NLP. After training the language model, we show that the perplexity metric calculated from a sample of the test (or production) data has a strong negative correlation with the quality of error correction of erroneous NGS reads. Therefore, we use the perplexity metric to guide a hill climbing-based search, converging toward the best configuration parameter value. Our approach is suitable for both de novo and comparative sequencing (resequencing), eliminating the need for a reference genome to serve as the ground truth. We find that Athena can automatically find the optimal value of k with a very high accuracy for 7 real datasets and using 3 different k-mer based EC algorithms, Lighter, Blue, and Racer. The inverse relation between the perplexity metric and alignment rate exists under all our tested conditions—for real and synthetic datasets, for all kinds of sequencing errors (insertion, deletion, and substitution), and for high and low error rates. The absolute value of that correlation is at least 73%. In our experiments, the best value of k found by Athena achieves an alignment rate within 0.53% of the oracle best value of k found through brute force searching (i.e., scanning through the entire range of k values). Athena’s selected value of k lies within the top-3 best k values using N-Gram models and the top-5 best k values using RNN models With best parameter selection by Athena, the assembly quality (NG50) is improved by a Geometric Mean of 4.72X across the 7 real datasets.

Draft genome sequence data of maqui (Aristotelia chilensis) and identification of SSR markers

Maqui (Aristotelia chilensis [Molina] Stunz) is a small dioecious tree, belonging to the Elaeocarpaceae family. Maqui fruit has high levels of antioxidant activity, which are due to elevated anthocyanin and polyphenol content. Here we describe a draft genome sequence data of maqui (A. chilensis). The genomic sequence datasets were obtained using Illumina NextSeq platform. Nucleotide sequences of raw reads and the assembled draft genome are available at NCBI's Sequence Read Archive as BioProject PRJNA544858. Also, a total of 210067 microsatellite or simple sequence repeat (SSR) markers were identified.

Data on the draft genome sequence of Caryocar brasiliense Camb. (Caryocaraceae): An important genetic resource from Brazilian savannas

Caryocar brasiliense (Caryocaraceae) is a Neotropical tree species widely distributed in Brazilian savannas. This species is very popular in central Brazil mainly due to the use of its fruits in the local cuisine and their anti-inflammatory proprieties, and indeed it is one of the candidates, among Brazilian native plants, for fast track incorporation into cropping systems. Considering the importance of Caryocar brasiliense, little is known about its genetics and genomics, and determination of a reference genome sequence could improve the understanding of its evolution, as well as the development of tools for domestication. Here, we provide the first draft genome of C. brasiliense, the raw sequencing data and some multiplex sets of high quality microsatellite primers. Data on the genome project can be obtained from the BioProject at NCBI (https://www.ncbi.nlm.nih.gov/bioproject/?term=caryocar).

doepipeline: a systematic approach to optimizing multi-level and multi-step data processing workflows

BACKGROUND

Selecting the proper parameter settings for bioinformatic software tools is challenging. Not only will each parameter have an individual effect on the outcome, but there are also potential interaction effects between parameters. Both of these effects may be difficult to predict. To make the situation even more complex, multiple tools may be run in a sequential pipeline where the final output depends on the parameter configuration for each tool in the pipeline. Because of the complexity and difficulty of predicting outcomes, in practice parameters are often left at default settings or set based on personal or peer experience obtained in a trial and error fashion. To allow for the reliable and efficient selection of parameters for bioinformatic pipelines, a systematic approach is needed.

RESULTS

We present doepipeline, a novel approach to optimizing bioinformatic software parameters, based on core concepts of the Design of Experiments methodology and recent advances in subset designs. Optimal parameter settings are first approximated in a screening phase using a subset design that efficiently spans the entire search space, then optimized in the subsequent phase using response surface designs and OLS modeling. Doepipeline was used to optimize parameters in four use cases; 1) de-novo assembly, 2) scaffolding of a fragmented genome assembly, 3) k-mer taxonomic classification of Oxford Nanopore Technologies MinION reads, and 4) genetic variant calling. In all four cases, doepipeline found parameter settings that produced a better outcome with respect to the characteristic measured when compared to using default values. Our approach is implemented and available in the Python package doepipeline.

CONCLUSIONS

Our proposed methodology provides a systematic and robust framework for optimizing software parameter settings, in contrast to labor- and time-intensive manual parameter tweaking. Implementation in doepipeline makes our methodology accessible and user-friendly, and allows for automatic optimization of tools in a wide range of cases. The source code of doepipeline is available at https://github.com/clicumu/doepipeline and it can be installed through conda-forge.

Draft genome sequence of Solanum aethiopicum provides insights into disease resistance, drought tolerance, and the evolution of the genome

BACKGROUND

The African eggplant (Solanum aethiopicum) is a nutritious traditional vegetable used in many African countries, including Uganda and Nigeria. It is thought to have been domesticated in Africa from its wild relative, Solanum anguivi. S. aethiopicum has been routinely used as a source of disease resistance genes for several Solanaceae crops, including Solanum melongena. A lack of genomic resources has meant that breeding of S. aethiopicum has lagged behind other vegetable crops.

RESULTS

We assembled a 1.02-Gb draft genome of S. aethiopicum, which contained predominantly repetitive sequences (78.9%). We annotated 37,681 gene models, including 34,906 protein-coding genes. Expansion of disease resistance genes was observed via 2 rounds of amplification of long terminal repeat retrotransposons, which may have occurred ∼1.25 and 3.5 million years ago, respectively. By resequencing 65 S. aethiopicum and S. anguivi genotypes, 18,614,838 single-nucleotide polymorphisms were identified, of which 34,171 were located within disease resistance genes. Analysis of domestication and demographic history revealed active selection for genes involved in drought tolerance in both "Gilo" and "Shum" groups. A pan-genome of S. aethiopicum was assembled, containing 51,351 protein-coding genes; 7,069 of these genes were missing from the reference genome.

CONCLUSIONS

The genome sequence of S. aethiopicum enhances our understanding of its biotic and abiotic resistance. The single-nucleotide polymorphisms identified are immediately available for use by breeders. The information provided here will accelerate selection and breeding of the African eggplant, as well as other crops within the Solanaceae family.

Population Genomic Approaches for Weed Science

Genomic approaches are opening avenues for understanding all aspects of biological life, especially as they begin to be applied to multiple individuals and populations. However, these approaches typically depend on the availability of a sequenced genome for the species of interest. While the number of genomes being sequenced is exploding, one group that has lagged behind are weeds. Although the power of genomic approaches for weed science has been recognized, what is needed to implement these approaches is unfamiliar to many weed scientists. In this review we attempt to address this problem by providing a primer on genome sequencing and provide examples of how genomics can help answer key questions in weed science such as: (1) Where do agricultural weeds come from; (2) what genes underlie herbicide resistance; and, more speculatively, (3) can we alter weed populations to make them easier to control? This review is intended as an introduction to orient weed scientists who are thinking about initiating genome sequencing projects to better understand weed populations, to highlight recent publications that illustrate the potential for these methods, and to provide direction to key tools and literature that will facilitate the development and execution of weed genomic projects.

Identification of telomerase RNAs in species of the Yarrowia clade provides insights into the co-evolution of telomerase, telomeric repeats and telomere-binding proteins

Telomeric repeats in fungi of the subphylum Saccharomycotina exhibit great inter- and intra-species variability in length and sequence. Such variations challenged telomeric DNA-binding proteins that co-evolved to maintain their functions at telomeres. Here, we compare the extent of co-variations in telomeric repeats, encoded in the telomerase RNAs (TERs), and the repeat-binding proteins from 13 species belonging to the Yarrowia clade. We identified putative TER loci, analyzed their sequence and secondary structure conservation, and predicted functional elements. Moreover, in vivo complementation assays with mutant TERs showed the functional importance of four novel TER substructures. The TER-derived telomeric repeat unit of all species, except for one, is 10 bp long and can be represented as 5′-TTNNNNAGGG-3′, with repeat sequence variations occuring primarily outside the vertebrate telomeric motif 5′-TTAGGG-3′. All species possess a homologue of the Yarrowia lipolytica Tay1 protein, YlTay1p. In vitro, YlTay1p displays comparable DNA-binding affinity to all repeat variants, suggesting a conserved role among these species. Taken together, these results add significant insights into the co-evolution of TERs, telomeric repeats and telomere-binding proteins in yeasts.

When the levee breaks: a practical guide to sketching algorithms for processing the flood of genomic data

Considerable advances in genomics over the past decade have resulted in vast amounts of data being generated and deposited in global archives. The growth of these archives exceeds our ability to process their content, leading to significant analysis bottlenecks. Sketching algorithms produce small, approximate summaries of data and have shown great utility in tackling this flood of genomic data, while using minimal compute resources. This article reviews the current state of the field, focusing on how the algorithms work and how genomicists can utilize them effectively. References to interactive workbooks for explaining concepts and demonstrating workflows are included at https://github.com/will-rowe/genome-sketching .