Is reliance on an inaccurate genome sequence sabotaging your experiments?

Klumpy: A tool to evaluate the integrity of long-read genome assemblies and illusive sequence motifs

The improvement and decreasing costs of third-generation sequencing technologies has widened the scope of biological questions researchers can address with de novo genome assemblies. With the increasing number of reference genomes, validating their integrity with minimal overhead is vital for establishing confident results in their applications. Here, we present Klumpy, a tool for detecting and visualizing both misassembled regions in a genome assembly and genetic elements (e.g. genes) of interest in a set of sequences. By leveraging the initial raw reads in combination with their respective genome assembly, we illustrate Klumpy's utility by investigating antifreeze glycoprotein (afgp) loci across two icefishes, by searching for a reported absent gene in the northern snakehead fish, and by scanning the reference genomes of a mudskipper and bumblebee for misassembled regions. In the two former cases, we were able to provide support for the noncanonical placement of an afgp locus in the icefishes and locate the missing snakehead gene. Furthermore, our genome scans were able identify an unmappable locus in the mudskipper reference genome and identify a putative repetitive element shared among several species of bees.

Multicopy subtelomeric genes underlie animal infectivity of divergent Cryptosporidium hominis subtypes

The anthroponotic Cryptosporidium hominis differs from the zoonotic C. parvum in its lack of infectivity to animals, but several divergent subtypes have recently been found in nonhuman primates and equines. Here, we sequence 17 animal C. hominis isolates and generate a new IbA12G3 genome at the chromosome level. Comparative analysis with 222 human isolates shows significant genetic divergence of the animal isolates, with genetic recombination among them. They have additional subtelomeric insulinase and MEDLE genes. In interferon-γ knockout mice, three monkey isolates show differences in infectivity and induce higher and longer oocyst shedding than a reference C. parvum isolate. Deletion of the MEDLE genes significantly reduces the growth and pathogenicity of a virulent strain in mice. Co-infection of two fluorescence-tagged C. hominis subtypes produces bicolored oocysts, supporting the conclusion that mixed subtype infections can lead to genetic recombination. These data provide insight into potential determinants of host infectivity in Cryptosporidium, and a convenient animal model for biological studies of C. hominis.

A practical guide and Galaxy workflow to avoid inter-plasmidic repeat collapse and false gene loss in Unicycler's hybrid assemblies

Generating complete, high-quality genome assemblies is key for any downstream analysis, such as comparative genomics. For bacterial genome assembly, various algorithms and fully automated pipelines exist, which are free-of-charge and easily accessible. However, these assembly tools often cannot unambiguously resolve a bacterial genome, for example due to the presence of sequence repeat structures on the chromosome or on plasmids. Then, a more sophisticated approach and/or manual curation is needed. Such modifications can be challenging, especially for non-bioinformaticians, because they are generally not considered as a straightforward process. In this study, we propose a standardized approach for manual genome completion focusing on the popular hybrid assembly pipeline Unicycler. The provided Galaxy workflow addresses two weaknesses in Unicycler's hybrid assemblies: (i) collapse of inter-plasmidic repeats and (ii) false loss of single-copy sequences. To demonstrate and validate how to detect and resolve these assembly errors, we use two genomes from the Bacillus cereus group. By applying the proposed pipeline following an automated assembly, the genome sequence quality can be significantly improved.

From contigs towards chromosomes: automatic improvement of long read assemblies (ILRA)

Recent advances in long read technologies not only enable large consortia to aim to sequence all eukaryotes on Earth, but they also allow individual laboratories to sequence their species of interest with relatively low investment. Long read technologies embody the promise of overcoming scaffolding problems associated with repeats and low complexity sequences, but the number of contigs often far exceeds the number of chromosomes and they may contain many insertion and deletion errors around homopolymer tracts. To overcome these issues, we have implemented the ILRA pipeline to correct long read-based assemblies. Contigs are first reordered, renamed, merged, circularized, or filtered if erroneous or contaminated. Illumina short reads are used subsequently to correct homopolymer errors. We successfully tested our approach by improving the genome sequences of Homo sapiens, Trypanosoma brucei, and Leptosphaeria spp., and by generating four novel Plasmodium falciparum assemblies from field samples. We found that correcting homopolymer tracts reduced the number of genes incorrectly annotated as pseudogenes, but an iterative approach seems to be required to correct more sequencing errors. In summary, we describe and benchmark the performance of our new tool, which improved the quality of novel long read assemblies up to 1 Gbp. The pipeline is available at GitHub: https://github.com/ThomasDOtto/ILRA.

New T2T assembly of Cryptosporidium parvum IOWA annotated with reference genome gene identifiers

Cryptosporidium parvum is a significant pathogen causing gastrointestinal infections in humans and animals, that is spread through the ingestion of contaminated food and water. Despite its global impact on public health, generating a C. parvum genome sequence has always been challenging due to a lack of in vitro cultivation systems and challenging sub-telomeric gene families. A gapless telomere to telomere genome assembly has been created for Cryptosporidium parvum IOWA obtained from Bunch Grass Farms, named here as CpBGF. There are 8 chromosomes that total 9,259,183 bp. The new hybrid assembly which was generated with Illumina and Oxford Nanopore resolves complex sub-telomeric regions of chromosomes 1, 7 and 8. To facilitate ease of use and consistency with the literature, whenever possible, chromosomes have been oriented and genes in this annotation have been given the same gene IDs used in the current reference genome sequence generated in 2004. The annotation of this assembly utilized considerable RNA expression evidence, thus, untranslated regions, long noncoding RNAs and antisense RNAs are annotated. The CpBGF genome assembly serves as a valuable resource for understanding the biology, pathogenesis, and transmission of C. parvum, and it facilitates the development of diagnostics, drugs, and vaccines against cryptosporidiosis.

Noise reduction strategies in metagenomic chromosome confirmation capture to link antibiotic resistance genes to microbial hosts

The gut microbiota is a reservoir for antimicrobial resistance genes (ARGs). With current sequencing methods, it is difficult to assign ARGs to their microbial hosts, particularly if these ARGs are located on plasmids. Metagenomic chromosome conformation capture approaches (meta3C and Hi-C) have recently been developed to link bacterial genes to phylogenetic markers, thus potentially allowing the assignment of ARGs to their hosts on a microbiome-wide scale. Here, we generated a meta3C dataset of a human stool sample and used previously published meta3C and Hi-C datasets to investigate bacterial hosts of ARGs in the human gut microbiome. Sequence reads mapping to repetitive elements were found to cause problematic noise in, and may importantly skew interpretation of, meta3C and Hi-C data. We provide a strategy to improve the signal-to-noise ratio by discarding reads that map to insertion sequence elements and to the end of contigs. We also show the importance of using spike-in controls to quantify whether the cross-linking step in meta3C and Hi-C protocols has been successful. After filtering to remove artefactual links, 87 ARGs were assigned to their bacterial hosts across all datasets, including 27 ARGs in the meta3C dataset we generated. We show that commensal gut bacteria are an important reservoir for ARGs, with genes coding for aminoglycoside and tetracycline resistance being widespread in anaerobic commensals of the human gut.

Cloning, Expression Analysis and SNP Screening of the kiss1 Gene in Male Schizothorax biddulphi

Schizothorax biddulphi is an endemic fish distributed only in southern Xinjiang, China. Due to overfishing, water conservancy facilities, and other factors, as well as inherent biological limitations, resource recovery is quite difficult. For endangered fish with slow growth, late sexual maturity, and insufficient natural population supplementation, large-scale artificial reproduction and breeding are important for restoring resources. Therefore, it is urgent to optimize the reproductive regulation methods of the fish. The kiss1 gene is a key regulator of the reproductive regulation cascade, and identifying and analyzing the role of kiss1 are important for further elucidating the reproductive mechanism of S. biddulphi. To understand the characteristics of the kiss1 of S. biddulphi, the full-length cDNA sequence of kiss1 was obtained in this study, and its tissue expression specificity and association with phenotypic traits were analyzed in male fish. The full-length cDNA sequence of kiss1 in S. biddulphi was 658 bp, with an ORF of 327 bp, and encoded a 108-amino acid, unstable protein. Homology results indicated that kiss1 was highly conserved. qPCR showed kiss1 expression in different tissues in male S. biddulphi, with the highest expression in the gonads, followed by muscle, and significantly lower expression in the swim bladder, pituitary gland, heart, hypothalamus, gill, fin, liver, eye, and mid-kidney. qPCR revealed three SNP loci in the exonic region of kiss1. The c.3G>T locus was significantly correlated (p < 0.05) with gonad mass and the maturation coefficient in S. biddulphi. These results will help uncover the reproductive endocrinology network of S. biddulphi, improve artificial breeding technology for fish, and unveil new directions for breeding excellent strains of S. biddulphi and molecular marker-assisted breeding.

A Nile rat transcriptomic landscape across 22 organs by ultra-deep sequencing and comparative RNA-seq pipeline (CRSP)

RNA sequencing (RNA-seq) has been a widely used high-throughput method to characterize transcriptomic dynamics spatiotemporally. However, RNA-seq data analysis pipelines typically depend on either a sequenced genome and/or corresponding reference transcripts. This limitation is a challenge for species lacking sequenced genomes and corresponding reference transcripts. The Nile rat (Arvicanthis niloticus) has two key features - it is daytime active, and it is prone to diet-induced diabetes, which makes it more similar to humans than regular laboratory rodents. However, at the time of this study, neither a Nile rat genome nor a reference transcript set were available, making it technically challenging to perform large-scale RNA-seq based transcriptomic studies. This genome-independent work progressed concurrently with our generation of a Nile rat genome. A well-annotated genome requires several iterations of manually reviewing curated transcripts and takes years to achieve. Here, we developed a Comparative RNA-Seq Pipeline (CRSP), integrating a comparative species strategy independent of a specific sequenced genome or species-matched reference transcripts. We performed benchmarking to validate that our CRSP tool can accurately quantify gene expression levels. In this study, we generated the first ultra-deep (2.3 billion × 2 paired-end) Nile rat RNA-seq data from 59 biopsy samples representing 22 major organs, providing a unique resource and spatial gene expression reference for Nile rat researchers. Importantly, CRSP is not limited to the Nile rat species and can be applied to any species without prior genomic knowledge. To facilitate a general use of CRSP, we also characterized the number of RNA-seq reads required for accurate estimation via simulation studies. CRSP and documents are available at: https://github.com/pjiang1105/CRSP.

Accessing the Variability of Multicopy Genes in Complex Genomes using Unassembled Next-Generation Sequencing Reads: The Case of Trypanosoma cruzi Multigene Families

Repetitive elements cause assembly fragmentation in complex eukaryotic genomes, limiting the study of their variability. The genome of Trypanosoma cruzi, the parasite that causes Chagas disease, has a high repetitive content, including multigene families. Although many T. cruzi multigene families encode surface proteins that play pivotal roles in host-parasite interactions, their variability is currently underestimated, as their high repetitive content results in collapsed gene variants. To estimate sequence variability and copy number variation of multigene families, we developed a read-based approach that is independent of gene-specific read mapping and de novo assembly. This methodology was used to estimate the copy number and variability of MASP, TcMUC, and Trans-Sialidase (TS), the three largest T. cruzi multigene families, in 36 strains, including members of all six parasite discrete typing units (DTUs). We found that these three families present a specific pattern of variability and copy number among the distinct parasite DTUs. Inter-DTU hybrid strains presented a higher variability of these families, suggesting that maintaining a larger content of their members could be advantageous. In addition, in a chronic murine model and chronic Chagasic human patients, the immune response was focused on TS antigens, suggesting that targeting TS conserved sequences could be a potential avenue to improve diagnosis and vaccine design against Chagas disease. Finally, the proposed approach can be applied to study multicopy genes in any organism, opening new avenues to access sequence variability in complex genomes. IMPORTANCE Sequences that have several copies in a genome, such as multicopy-gene families, mobile elements, and microsatellites, are among the most challenging genomic segments to study. They are frequently underestimated in genome assemblies, hampering the correct assessment of these important players in genome evolution and adaptation. Here, we developed a new methodology to estimate variability and copy numbers of repetitive genomic regions and employed it to characterize the T. cruzi multigene families MASP, TcMUC, and transsialidase (TS), which are important virulence factors in this parasite. We showed that multigene families vary in sequence and content among the parasite's lineages, whereas hybrid strains have a higher sequence variability that could be advantageous to the parasite's survivability. By identifying conserved sequences within multigene families, we showed that the mammalian host immune response toward these multigene families is usually focused on the TS multigene family. These TS conserved and immunogenic peptides can be explored in future works as diagnostic targets or vaccine candidates for Chagas disease. Finally, this methodology can be easily applied to any organism of interest, which will aid in our understanding of complex genomic regions.

Nematode chromosomes

The nematode Caenorhabditis elegans has shed light on many aspects of eukaryotic biology, including genetics, development, cell biology, and genomics. A major factor in the success of C. elegans as a model organism has been the availability, since the late 1990s, of an essentially gap-free and well-annotated nuclear genome sequence, divided among 6 chromosomes. In this review, we discuss the structure, function, and biology of C. elegans chromosomes and then provide a general perspective on chromosome biology in other diverse nematode species. We highlight malleable chromosome features including centromeres, telomeres, and repetitive elements, as well as the remarkable process of programmed DNA elimination (historically described as chromatin diminution) that induces loss of portions of the genome in somatic cells of a handful of nematode species. An exciting future prospect is that nematode species may enable experimental approaches to study chromosome features and to test models of chromosome evolution. In the long term, fundamental insights regarding how speciation is integrated with chromosome biology may be revealed.

Long-read assembly and comparative evidence-based reanalysis of Cryptosporidium genome sequences reveal expanded transporter repertoire and duplication of entire chromosome ends including subtelomeric regions

Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community has only had access to a good, but incomplete, Cryptosporidium parvum IOWA reference genome sequence. Incomplete reference sequences hamper annotation, experimental design, and interpretation. We have generated a new C. parvum IOWA genome assembly supported by Pacific Biosciences (PacBio) and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species: C. parvum, Cryptosporidium hominis, and Cryptosporidium tyzzeri We made 1926 C. parvum annotation updates based on experimental evidence. They include new transporters, ncRNAs, introns, and altered gene structures. The new assembly and annotation revealed a complete Dnmt2 methylase ortholog. Comparative annotation between C. parvum, C. hominis, and C. tyzzeri revealed that most "missing" orthologs are found, suggesting that the biological differences between the species must result from gene copy number variation, differences in gene regulation, and single-nucleotide variants (SNVs). Using the new assembly and annotation as reference, 190 genes are identified as evolving under positive selection, including many not detected previously. The new C. parvum IOWA reference genome assembly is larger, gap free, and lacks ambiguous bases. This chromosomal assembly recovers all 16 chromosome ends, 13 of which are contiguously assembled. The three remaining chromosome ends are provisionally placed. These ends represent duplication of entire chromosome ends including subtelomeric regions revealing a new level of genome plasticity that will both inform and impact future research.

Overview and New Insights Into the Diversity, Evolution, Role, and Regulation of Kisspeptins and Their Receptors in Teleost Fish

In the last two decades, kisspeptin (Kiss) has been identified as an important player in the regulation of reproduction and other physiological functions in vertebrates, including several fish species. To date, two ligands (Kiss1, Kiss2) and three kisspeptin receptors (Kissr1, Kissr2, Kissr3) have been identified in teleosts, likely due to whole-genome duplication and loss of genes that occurred early in teleost evolution. Recent results in zebrafish and medaka mutants have challenged the notion that the kisspeptin system is essential for reproduction in fish, in marked contrast to the situation in mammals. In this context, this review focuses on the role of kisspeptins at three levels of the reproductive, brain-pituitary-gonadal (BPG) axis in fish. In addition, this review compiled information on factors controlling the Kiss/Kissr system, such as photoperiod, temperature, nutritional status, sex steroids, neuropeptides, and others. In this article, we summarize the available information on the molecular diversity and evolution, tissue expression and neuroanatomical distribution, functional significance, signaling pathways, and gene regulation of Kiss and Kissr in teleost fishes. Of particular note are recent advances in understanding flatfish kisspeptin systems, which require further study to reveal their structural and functional diversity.

Comparative Genomics across Three Ensifer Species Using a New Complete Genome Sequence of the Medicago Symbiont Sinorhizobium (Ensifer) meliloti WSM1022

Here, we report an improved and complete genome sequence of Sinorhizobium (Ensifer) meliloti strain WSM1022, a microsymbiont of Medicago species, revealing its tripartite structure. This improved genome sequence was generated combining Illumina and Oxford nanopore sequencing technologies to better understand the symbiotic properties of the bacterium. The 6.75 Mb WSM1022 genome consists of three scaffolds, corresponding to a chromosome (3.70 Mb) and the pSymA (1.38 Mb) and pSymB (1.66 Mb) megaplasmids. The assembly has an average GC content of 62.2% and a mean coverage of 77X. Genome annotation of WSM1022 predicted 6058 protein coding sequences (CDSs), 202 pseudogenes, 9 rRNAs (3 each of 5S, 16S, and 23S), 55 tRNAs, and 4 ncRNAs. We compared the genome of WSM1022 to two other rhizobial strains, closely related Sinorhizobium (Ensifer) meliloti Sm1021 and Sinorhizobium (Ensifer) medicae WSM419. Both WSM1022 and WSM419 species are high-efficiency rhizobial strains when in symbiosis with Medicago truncatula, whereas Sm1021 is ineffective. Our findings report significant genomic differences across the three strains with some similarities between the meliloti strains and some others between the high efficiency strains WSM1022 and WSM419. The addition of this high-quality rhizobial genome sequence in conjunction with comparative analyses will help to unravel the features that make a rhizobial symbiont highly efficient for nitrogen fixation.

High-quality reannotation of the king scallop genome reveals no 'gene-rich' feature and evolution of toxin resistance

The king scallop, Pecten maximus is a well-known, commercially important scallop species and is featured with remarkable tolerance to potent phytotoxins such as domoic acid. A high-quality genome can shed light on its biology and innovative evolution of toxin resistance. A reference genome has recently been published for P. maximus, however, it is suspicious that over 67,700 genes are annotated in this genome, which is unexpectedly larger than its close relatives of pectinids. Herein, we provide an improved high-quality chromosome-level reference genome assembly and annotation for the king scallop P. maximus. A final set of 26,995 genes is annotated after carefully checking and curation of the predicted gene models, which significantly improves the accuracy of gene structure information. The large number of gene duplicates in the previous genome is mainly distorted by the fragmented annotation. Through integrated genomic, evolutionary and transcriptomic analyses, we reveal that the Phi subfamily of ionotropic glutamate receptors (iGluRs) are well preserved in molluscs, and P. maximus experienced the rapid expansion of the Phi class of iGluR (GluF) gene family. The GluF genes exhibit ubiquitously high expression and altered sequence characteristics for ligand selectivity, which may contribute to the remarkable tolerance to neurotoxins in P. maximus. Taken together, our study disapproves the previous claim of the 'gene-rich' genome of this species and provides a high-quality genome assembly for further understanding of its biology and evolution of toxin resistance.

Comparative genomics of Leishmania isolates from Brazil confirms the presence of Leishmania major in the Americas

Leishmania (Leishmania) major is an important agent of cutaneous leishmaniasis, having as a vector sandflies belonging to the genus Phlebotomus. Although this species has been described as restricted to the Old World, parasites similar to L. major have been isolated from South American patients who have never travelled abroad. These parasites were named "L. major-like", and several studies have been carried out to characterise them biochemically, molecularly, and biologically. However, the phylogenetic origin of these isolates is still unknown. In the present study we characterised three L. major-like isolates, named BH49, BH121 and BH129, using comparative genomics approaches. We evaluated the presence of gene and segmental duplications/deletions and the presence of aneuploidies that could explain the differences in infectivity observed in the BH49 and BH121 isolates. All isolates presented a pattern of mosaic aneuploidy and gene copy number variation, which are common in the genus Leishmania. Virulence factors such as phosphatases and peptidases were found to have increased gene copy numbers in the infective isolate, which could explain the difference in infectivity previously observed between BH121 and BH49. Phylogenetic analyses revealed that BH49, BH121 and BH129 L. major-like grouped with L. major isolates, and suggest they were imported from the Old World in at least two independent events. We suggest that new epidemiological inquiries should also evaluate L. major infections in South America, to assess the epidemiological importance of this species in the New World.

Cultivation and Genome Sequencing of Bacteria Isolated From the Coffee Berry Borer (Hypothenemus hampei), With Emphasis on the Role of Caffeine Degradation

The coffee berry borer, the most economically important insect pest of coffee worldwide, is the only insect capable of feeding and reproducing solely on the coffee seed, a food source containing the purine alkaloid caffeine. Twenty-one bacterial species associated with coffee berry borers from Hawai’i, Mexico, or a laboratory colony in Maryland (Acinetobacter sp. S40, S54, S55, Bacillus aryabhattai, Delftia lacustris, Erwinia sp. S38, S43, S63, Klebsiella oxytoca, Ochrobactrum sp. S45, S46, Pantoea sp. S61, Pseudomonas aeruginosa, P. parafulva, and Pseudomonas sp. S30, S31, S32, S37, S44, S60, S75) were found to have at least one of five caffeine N-demethylation genes (ndmA, ndmB, ndmC, ndmD, ndmE), with Pseudomonas spp. S31, S32, S37, S60 and P. parafulva having the full complement of these genes. Some of the bacteria carrying the ndm genes were detected in eggs, suggesting possible vertical transmission, while presence of caffeine-degrading bacteria in frass, e.g., P. parafulva (ndmABCDE) and Bacillus aryabhattai (ndmA) could result in horizontal transmission to all insect life stages. Thirty-five bacterial species associated with the insect (Acinetobacter sp. S40, S54, S55, B. aryabhattai, B. cereus group, Bacillus sp. S29, S70, S71, S72, S73, D. lacustris, Erwinia sp. S38, S43, S59, S63, K. oxytoca, Kosakonia cowanii, Ochrobactrum sp. S45, S46, Paenibacillus sp. S28, Pantoea sp. S61, S62, P. aeruginosa, P. parafulva, Pseudomonas sp. S30, S31, S32, S37, S44, S60, S75, Stenotrophomonas sp. S39, S41, S48, S49) might contribute to caffeine breakdown using the C-8 oxidation pathway, based on presence of genes required for this pathway. It is possible that caffeine-degrading bacteria associated with the coffee berry borer originated as epiphytes and endophytes in the coffee plant microbiota.

Continuous Microevolution Accelerates Disease Progression during Sequential Episodes of Infection

Bacteria adapt to dynamic changes in the host during chronic and recurrent infections. Bacterial microevolution is one type of adaptation that imparts a selective advantage. We hypothesize that recurrent episodes of disease promote microevolution through genetic mutations that modulate disease severity. We use a pre-clinical model of otitis media (OM) to determine the potential role for microevolution of nontypeable Haemophilus influenzae (NTHI) during sequential episodes of disease. Whole genome sequencing reveals microevolution of hemoglobin binding and lipooligosaccharide (LOS) biosynthesis genes, suggesting that adaptation of these systems is critical for infection. These OM-adapted strains promote increased biofilm formation, inflammation, stromal fibrosis, and an increased propensity to form intracellular bacterial communities (IBCs). Remarkably, IBCs remain for at least one month following clinical resolution of infection, suggesting an intracellular reservoir as a nidus for recurrent OM. Additional approaches for therapeutic design tailored to combat this burdensome disease will arise from these studies.

Harrison et al. develop a sequential model of otitis media (OM) to investigate microevolution through genetic mutations that modulate disease severity. OM-adapted strains promote increased biofilm, inflammation, stromal fibrosis, and intracellular bacterial community (IBC) development. IBCs remain one month following clinical resolution of infection, suggesting a nidus for recurrent OM.