PacBio But Not Illumina Technology Can Achieve Fast, Accurate and Complete Closure of the High GC, Complex Burkholderia pseudomallei Two-Chromosome Genome

Distribution and conservation of simple sequence repeats in plant pathogenic species of Zymoseptoria and development of genomic resources for its orphaned species

To better understand the structure and evolution of the genomes of four plant pathogenic species of Zymoseptoria, we analyzed the occurrence, relative abundance (RA), and density (RD) of simple sequence repeats (SSRs) in their whole genome and transcriptome sequences. In this study, SSRs are defined as repeats of more than 12 bases in length. The genome and transcriptome sequences of Zymoseptoria ardabiliae show the highest RA (201.1 and 129.9) and RD (3229.4 and 1928.2) of SSRs, while those of Zymoseptoria pseudotritici show the lowest RA (167.2 and 118.5) and RD (2482.2 and 1687.0). The majority of SSRs in the genomic and transcriptome sequences of species were trinucleotide SSRs, while dinucleotide SSRs were the least common. The most common trinucleotide motifs in the transcriptomic sequences across all species were those that encoded the amino acid arginine. As per our motif conservation study, Zymoseptoria tritici (12.4%) possessed the most unique motifs, while Z. pseudotritici (3.9%) had the fewest. Overall, only 38.1% of the motifs were found to be conserved among the species. Gene enrichment studies reveal that three of the species, Z. ardabiliae, Zymoseptoria brevis, and Z. pseudotritici, have SSRs in their genes related to cellular metabolism, while the remaining Z. tritici harbors SSRs in genes related to DNA synthesis and gene expression. In an effort to improve the genetic resources for the orphan species of pathogenic Zymoseptoria, a total of 73,134 primers were created. The genomic resources developed in this study could help with analyses of genetic relatedness within the population and the development of species-specific markers.

Comparative genome analysis reveals high-level drug resistance markers in a clinical isolate of Mycobacterium fortuitum subsp. fortuitum MF GZ001

Introduction

Infections caused by non-tuberculosis mycobacteria are significantly worsening across the globe. M. fortuitum complex is a rapidly growing pathogenic species that is of clinical relevance to both humans and animals. This pathogen has the potential to create adverse effects on human healthcare.

Methods

The MF GZ001 clinical strain was collected from the sputum of a 45-year-old male patient with a pulmonary infection. The morphological studies, comparative genomic analysis, and drug resistance profiles along with variants detection were performed in this study. In addition, comparative analysis of virulence genes led us to understand the pathogenicity of this organism.

Results

Bacterial growth kinetics and morphology confirmed that MF GZ001 is a rapidly growing species with a rough morphotype. The MF GZ001 contains 6413573 bp genome size with 66.18 % high G+C content. MF GZ001 possesses a larger genome than other related mycobacteria and included 6156 protein-coding genes. Molecular phylogenetic tree, collinearity, and comparative genomic analysis suggested that MF GZ001 is a novel member of the M. fortuitum complex. We carried out the drug resistance profile analysis and found single nucleotide polymorphism (SNP) mutations in key drug resistance genes such as rpoB, katG, AAC(2')-Ib, gyrA, gyrB, embB, pncA, blaF, thyA, embC, embR, and iniA. In addition, the MF GZ001strain contains mutations in iniA, iniC, pncA, and ribD which conferred resistance to isoniazid, ethambutol, pyrazinamide, and para-aminosalicylic acid respectively, which are not frequently observed in rapidly growing mycobacteria. A wide variety of predicted putative potential virulence genes were found in MF GZ001, most of which are shared with well-recognized mycobacterial species with high pathogenic profiles such as M. tuberculosis and M. abscessus.

Discussion

Our identified novel features of a pathogenic member of the M. fortuitum complex will provide the foundation for further investigation of mycobacterial pathogenicity and effective treatment.

The Development of Technology to Prevent, Diagnose, and Manage Antimicrobial Resistance in Healthcare-Associated Infections

There is a growing risk of antimicrobial resistance (AMR) having an adverse effect on the healthcare system, which results in higher healthcare costs, failed treatments and a higher death rate. A quick diagnostic test that can spot infections resistant to antibiotics is essential for antimicrobial stewardship so physicians and other healthcare professionals can begin treatment as soon as possible. Since the development of antibiotics in the last two decades, traditional, standard antimicrobial treatments have failed to treat healthcare-associated infections (HAIs). These results have led to the development of a variety of cutting-edge alternative methods to combat multidrug-resistant pathogens in healthcare settings. Here, we provide an overview of AMR as well as the technologies being developed to prevent, diagnose, and control healthcare-associated infections (HAIs). As a result of better cleaning and hygiene practices, resistance to bacteria can be reduced, and new, quick, and accurate instruments for diagnosing HAIs must be developed. In addition, we need to explore new therapeutic approaches to combat diseases caused by resistant bacteria. In conclusion, current infection control technologies will be crucial to managing multidrug-resistant infections effectively. As a result of vaccination, antibiotic usage will decrease and new resistance mechanisms will not develop.

Whole-Genome Sequence and Comparative Analysis of Trichoderma asperellum ND-1 Reveal Its Unique Enzymatic System for Efficient Biomass Degradation

The lignocellulosic enzymes of Trichoderma asperellum have been intensely investigated toward efficient conversion of biomass into high-value chemicals/industrial products. However, lack of genome data is a remarkable hurdle for hydrolase systems studies. The secretory enzymes of newly isolated T. asperellum ND-1 during lignocellulose degradation are currently poorly known. Herein, a high-quality genomic sequence of ND-1, obtained by both Illumina HiSeq 2000 sequencing platforms and PacBio single-molecule real-time, has an assembly size of 35.75 Mb comprising 10,541 predicted genes. Secretome analysis showed that 895 proteins were detected, with 211 proteins associated with carbohydrate-active enzymes (CAZymes) responsible for biomass hydrolysis. Additionally, T. asperellum ND-1, T. atroviride IMI 206040, and T. virens Gv-298 shared 801 orthologues that were not identified in T. reesei QM6a, indicating that ND-1 may play critical roles in biological-control. In-depth analysis suggested that, compared with QM6a, the genome of ND-1 encoded a unique enzymatic system, especially hemicellulases and chitinases. Moreover, after comparative analysis of lignocellulase activities of ND-1 and other fungi, we found that ND-1 displayed higher hemicellulases (particularly xylanases) and comparable cellulases activities. Our analysis, combined with the whole-genome sequence information, offers a platform for designing advanced T. asperellum ND-1 strains for industrial utilizations, such as bioenergy production.

Long-Reads-Based Metagenomics in Clinical Diagnosis With a Special Focus on Fungal Infections

Identification of the causative infectious agent is essential in the management of infectious diseases, with the ideal diagnostic method being rapid, accurate, and informative, while remaining cost-effective. Traditional diagnostic techniques rely on culturing and cell propagation to isolate and identify the causative pathogen. These techniques are limited by the ability and the time required to grow or propagate an agent in vitro and the facts that identification based on morphological traits are non-specific, insensitive, and reliant on technical expertise. The evolution of next-generation sequencing has revolutionized genomic studies to generate more data at a cheaper cost. These are divided into short- and long-read sequencing technologies, depending on the length of reads generated during sequencing runs. Long-read sequencing also called third-generation sequencing emerged commercially through the instruments released by Pacific Biosciences and Oxford Nanopore Technologies, although relying on different sequencing chemistries, with the first one being more accurate both platforms can generate ultra-long sequence reads. Long-read sequencing is capable of entirely spanning previously established genomic identification regions or potentially small whole genomes, drastically improving the accuracy of the identification of pathogens directly from clinical samples. Long-read sequencing may also provide additional important clinical information, such as antimicrobial resistance profiles and epidemiological data from a single sequencing run. While initial applications of long-read sequencing in clinical diagnosis showed that it could be a promising diagnostic technique, it also has highlighted the need for further optimization. In this review, we show the potential long-read sequencing has in clinical diagnosis of fungal infections and discuss the pros and cons of its implementation.

Discovery of alternatively spliced isoforms and long non-coding RNA in full length brain transcriptomes of anadromous Hilsa shad, Tenualosa ilisha (Hamilton, 1822)

BACKGROUND

Full length transcriptomes, achieved through long-read sequencing, along with the isoform analysis can reveal complexities in the gene expression profiles, as well as annotate the transcriptomes of non-model organisms.

METHODS AND RESULT

Full length transcripts of brain transcriptome of Tenualosa ilisha, Hilsa shad, were generated through PacBio single molecule real-time sequencing and were characterized. A total of 8.30 Gb clean reads were generated, with PacBio RSII, which resulted in 57,651 high quality consensus transcripts. After removing redundant reads, a total of 19,220 high-quality non-redundant transcripts and 17,341 full length ORF transcripts were classified to 7522 putative ortholog groups. Genes involved in various neural pathways were identified. In addition, isoform clusters and lncRNAs were discovered, along with Hilsa specific transcripts with coding frames and 29,147 SSRs in 944 transcripts (1141 annotated).

CONCLUSION

The present study provided, for the first time, a comprehensive view of the alternative isoforms of genes and transcriptome complexity in Hilsa shad brain and forms a rich resource for functional studies in brain of this anadromous fish.

Straglr: discovering and genotyping tandem repeat expansions using whole genome long-read sequences

Tandem repeat (TR) expansion is the underlying cause of over 40 neurological disorders. Long-read sequencing offers an exciting avenue over conventional technologies for detecting TR expansions. Here, we present Straglr, a robust software tool for both targeted genotyping and novel expansion detection from long-read alignments. We benchmark Straglr using various simulations, targeted genotyping data of cell lines carrying expansions of known diseases, and whole genome sequencing data with chromosome-scale assembly. Our results suggest that Straglr may be useful for investigating disease-associated TR expansions using long-read sequencing.

Knockout of a highly GC-rich gene in Burkholderia pyrrocinia by recombineering with freeze-thawing transformation

Genetic transformation is a valuable and essential method that provides powerful insights into the gene function of microorganisms and contributes to the construction of engineered bacteria. Here, we developed a novel genetic transformation system to easily knock out a highly GC-rich gene (74.71% GC) from Burkholderia pyrrocinia JK-SH007, a biocontrol strain of poplar canker disease. This system revealed a reliable selectable marker (trimethoprim resistance gene, Tmp) and a simplified, efficient transformation method (6,363.64 CFU/μg, pHKT2) that was developed via freeze-thawing. The knockout recombineering of B. pyrrocinia JK-SH007 was achieved through a suicide plasmid with a three-fragment mutagenesis construct. The three-fragment cassette for mutagenesis was generated by overlap extension and touchdown PCRs and composed of Tmp flanked by GC-rich upstream and downstream fragments from B. pyrrocinia JK-SH007. The mutant strain (ΔBpEG), which was verified by PCR, lost 93.3% of its ability to degrade carboxymethyl cellulose over 40 days. Overall, this system may contribute to future research on B. pyrrocinia traits.

Helicobacter pylori Antimicrobial Resistance and Gene Variants in High- and Low-Gastric-Cancer-Risk Populations

Colombia, South America has one of the world's highest burdens of Helicobacter pylori infection and gastric cancer. While multidrug antibiotic regimens can effectively eradicate H. pylori, treatment efficacy is being jeopardized by the emergence of antibiotic-resistant H. pylori strains. Moreover, the spectrum of and genetic mechanisms for antibiotic resistance in Colombia is underreported. In this study, 28 H. pylori strains isolated from gastric biopsy specimens from a high-gastric-cancer-risk (HGCR) population living in the Andes Mountains in Túquerres, Colombia and 31 strains from a low-gastric-cancer-risk (LGCR) population residing on the Pacific coast in Tumaco, Colombia were subjected to antibiotic susceptibility testing for amoxicillin, clarithromycin, levofloxacin, metronidazole, rifampin, and tetracycline. Resistance-associated genes were amplified by PCR for all isolates, and 29 isolates were whole-genome sequenced (WGS). No strains were resistant to amoxicillin, clarithromycin, or rifampin. One strain was resistant to tetracycline and had an A926G mutation in its 16S rRNA gene. Levofloxacin resistance was observed in 12/59 isolates and was significantly associated with N87I/K and/or D91G/Y mutations in gyrA Most isolates were resistant to metronidazole; this resistance was significantly higher in the LGCR (31/31) group compared to the HGCR (24/28) group. Truncations in rdxA and frxA were present in nearly all metronidazole-resistant strains. There was no association between phylogenetic relationship and resistance profiles based on WGS analysis. Our results indicate H. pylori isolates from Colombians exhibit multidrug antibiotic resistance. Continued surveillance of H. pylori antibiotic resistance in Colombia is warranted in order to establish appropriate eradication treatment regimens for this population.

Rapid Methods for Antimicrobial Resistance Diagnostics

Antimicrobial resistance (AMR) is one of the most challenging threats in public health; thus, there is a growing demand for methods and technologies that enable rapid antimicrobial susceptibility testing (AST). The conventional methods and technologies addressing AMR diagnostics and AST employed in clinical microbiology are tedious, with high turnaround times (TAT), and are usually expensive. As a result, empirical antimicrobial therapies are prescribed leading to AMR spread, which in turn causes higher mortality rates and increased healthcare costs. This review describes the developments in current cutting-edge methods and technologies, organized by key enabling research domains, towards fighting the looming AMR menace by employing recent advances in AMR diagnostic tools. First, we summarize the conventional methods addressing AMR detection, surveillance, and AST. Thereafter, we examine more recent non-conventional methods and the advancements in each field, including whole genome sequencing (WGS), matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) spectrometry, Fourier transform infrared (FTIR) spectroscopy, and microfluidics technology. Following, we provide examples of commercially available diagnostic platforms for AST. Finally, perspectives on the implementation of emerging concepts towards developing paradigm-changing technologies and methodologies for AMR diagnostics are discussed.

High Prevalence and Mechanism Associated With Extended Spectrum Beta-Lactamase-Positive Phenotype in Laribacter hongkongensis

In this study, we reported the prevalence and mechanism associated with the extended-spectrum beta-lactamase (ESBL)-positive phenotype in Laribacter hongkongensis isolated from patients and fish. Using the inhibition zone enhancement test, 20 (95.2%) of the 21 patient strains and 8 (57.1%) of the 14 fish strains were tested ESBL-positive. However, ESBL genes, including SHV, TEM, CTX-M, GES, and PER, were not detected in all of these 28 L. hongkongensis isolates. No ESBL gene could be detected in either the complete genome of L. hongkongensis HLHK9 or the draft genome of PW3643. PCR and DNA sequencing revealed that all the 35 L. hongkongensis isolates (showing both ESBL-positive and ESBL-negative phenotypes) were positive for the ampC gene. When the AmpC deletion mutant, HLHK9ΔampC, was subject to the zone enhancement test, the difference of zone size between ceftazidime/clavulanate and ceftazidime was less than 5 mm. When boronic acid was added to the antibiotic disks, none of the 28 “ESBL-positive” isolates showed a ≥ 5 mm enhancement of inhibition zone size diameter between ceftazidime/clavulanate and ceftazidime and between cefotaxime/clavulanate and cefotaxime. A high prevalence (80%) of ESBL-positive phenotype is present in L. hongkongensis. Overall, our results suggested that the ESBL-positive phenotype in L. hongkongensis results from the expression of the intrinsic AmpC beta-lactamase. Confirmatory tests should be performed before issuing laboratory reports for L. hongkongensis isolates that are tested ESBL-positive by disk diffusion clavulanate inhibition test.

Genome Insights into the Novel Species Jejubacter calystegiae, a Plant Growth-Promoting Bacterium in Saline Conditions

Jejubacter calystegiae KSNA2T, a moderately halophilic, endophytic bacterium isolated from beach morning glory (Calystegia soldanella), was determined to be a novel species in a new genus in the family Enterobacteriaceae. To gain insights into the genetic basis of the salinity stress response of strain KSNA2T, we sequenced its genome using two complementary sequencing platforms (Illumina HiSeq and PacBio RSII). The genome contains a repertoire of metabolic pathways, such as those for nitrogen, phosphorus, and some amino acid metabolism pathways. Functional annotation of the KSNA2T genome revealed several genes involved in salt tolerance pathways, such as those encoding sodium transporters, potassium transporters, and osmoprotectant enzymes. Plant growth-promoting bacteria-based experiments indicated that strain KSNA2T promotes the germination of vegetable seeds in saline conditions. Overall, the genetic and biological analyses of strain KSNA2T provide valuable insights into bacteria-mediated salt tolerance in agriculture.

Combining phage display with SMRTbell next-generation sequencing for the rapid discovery of functional scFv fragments

Phage display technology in combination with next-generation sequencing (NGS) currently is a state-of-the-art method for the enrichment and isolation of monoclonal antibodies from diverse libraries. However, the current NGS methods employed for sequencing phage display libraries are limited by the short contiguous read lengths associated with second-generation sequencing platforms. Consequently, the identification of antibody sequences has conventionally been restricted to individual antibody domains or to the analysis of single domain binding moieties such as camelid VHH or cartilaginous fish IgNAR antibodies. In this study, we report the application of third-generation sequencing to address this limitation. We used single molecule real time (SMRT) sequencing coupled with hairpin adaptor loop ligation to facilitate the accurate interrogation of full-length single-chain Fv (scFv) libraries. Our method facilitated the rapid isolation and testing of scFv antibodies enriched from phage display libraries within days following panning. Two libraries against CD160 and CD123 were panned and monitored by NGS. Analysis of NGS antibody data sets led to the isolation of several functional scFv antibodies that were not identified by conventional panning and screening strategies. Our approach, which combines phage display selection of immune libraries with the full-length interrogation of scFv fragments, is an easy method to discover functional antibodies, with a range of affinities and biophysical characteristics.

Complete genome sequence and epigenetic profile of Bacillus velezensis UCMB5140 used for plant and crop protection in comparison with other plant-associated Bacillus strains

The application of biocontrol biopesticides based on plant growth-promoting rhizobacteria (PGPR), particularly members of the genus Bacillus, is considered a promising perspective to make agricultural practices sustainable and ecologically safe. Recent advances in genome sequencing by third-generation sequencing technologies, e.g., Pacific Biosciences' Single Molecule Real-Time (PacBio SMRT) platform, have allowed researchers to gain deeper insights into the molecular and genetic mechanisms of PGPR activities, and to compare whole genome sequences and global patterns of epigenetic modifications. In the current work, this approach was used to sequence and compare four Bacillus strains that exhibited various PGPR activities including the strain UCMB5140, which is used in the commercial biopesticide Phytosubtil. Whole genome comparison and phylogenomic inference assigned the strain UCMB5140 to the species Bacillus velezensis. Strong biocontrol activities of this strain were confirmed in several bioassays. Several factors that affect the evolution of active PGPR B. velezensis strains were identified: (1) horizontal acquisition of novel non-ribosomal peptide synthetases (NRPS) and adhesion genes; (2) rearrangements of functional modules of NRPS genes leading to strain specific combinations of their encoded products; (3) gain and loss of methyltransferases that can cause global alterations in DNA methylation patterns, which eventually may affect gene expression and regulate transcription. Notably, we identified a horizontally transferred NRPS operon encoding an uncharacterized polypeptide antibiotic in B. velezensis UCMB5140. Other horizontally acquired genes comprised a possible adhesin and a methyltransferase, which may explain the strain-specific methylation pattern of the chromosomal DNA of UCMB5140. KEY POINTS: • Whole genome sequence of the active PGPR Bacillus velezensis UCMB5140. • Identification of genetic determinants responsible for PGPR activities. • Role of methyltransferases and epigenetic mechanisms in evolution of bacteria.

Tsukamurella asaccharolytica sp. nov., Tsukamurella conjunctivitidis sp. nov. and Tsukamurella sputi sp. nov., isolated from patients with bacteraemia, conjunctivitis and respiratory infection in Hong Kong

Three bacterial strains, HKU70^T, HKU71^T and HKU72^T, were isolated from the conjunctival swab, blood and sputum samples of three patients with conjunctivitis, bacteraemia and respiratory infection, respectively, in Hong Kong. The three strains were aerobic, Gram-stain positive, catalase-positive, non-sporulating and non-motile bacilli and exhibited unique biochemical profiles distinguishable from currently recognized Tsukamurella species. 16S rRNA, secA, rpoB and groEL gene sequence analyses revealed that the three strains shared 99.6-99.9, 94.5-96.8, 95.7-97.8 and 97.7-98.9 % nucleotide identities with their corresponding closest Tsukamurella species respectively. DNA-DNA hybridization confirmed that they were distinct from other known species of the genus Tsukamurella (26.2±2.4 to 36.8±1.2 % DNA-DNA relatedness), in line with results of in silico genome-to-genome comparison (32.2-40.9 % Genome-to-Genome Distance Calculator and 86.3-88.9 % average nucleotide identity values]. Fatty acids, mycolic acids, cell-wall sugars and peptidoglycan analyses showed that they were typical of members of Tsukamurella. The G+C content determined based on the genome sequence of strains HKU70^T, HKU71^T and HKU72^T were 69.9, 70.2 and 70.5 mol%, respectively. Taken together, our results supported the proposition and description of three new species, i.e. Tsukamurella sputi HKU70^T (=JCM 33387^T=DSM 109106^T) sp. nov., Tsukamurella asaccharolytica HKU71^T (=JCM 33388^T=DSM 109107^T) sp. nov. and Tsukamurella conjunctivitidis HKU72^T (=JCM 33389^T=DSM 109108^T) sp. nov.

Comparative genomic analysis of Chinese human leptospirosis vaccine strain and circulating isolate

LEPTOSPIRA INTERROGANS

serogroup Canicola is one of the most important pathogens causing leptospirosis and is used as a vaccine strain of the current Chinese human leptospirosis vaccine. To characterize leptospiral pathogens, L. interrogans serogroup Canicola vaccine strain 611 and circulating isolate LJ178 from different hosts at different periods were sequenced using a combined strategy of Illumina X10 and PacBio technologies, and a comprehensive comparative analysis with other published Leptospira strains was conducted in this study. High levels of genomic similarities were observed between vaccine strain 611 and circulating isolate LJ178; both had two circular chromosomes and two circular extrachromosomal replicons. Compared with the strain 611 genome, 132 single nucleotide polymorphisms and 92 indels were found in strain LJ178. The larger lipopolysaccharide biosynthesis locus of serogroup Canicola was identified in both genomes. The phylogenetic analysis based on whole-genome sequences revealed that serogroup Canicola was not restricted to a specific host or geographic location, suggesting adaptive evolution associated with the ecologic diversity. In summary, our findings provide insights into a better molecular understanding of the component strains of human leptospirosis vaccine in China. Furthermore, these data detail the genetic composition and evolutionary relatedness of Leptospira strains that pose a health risk to humans.

Functional Genomic Identification of Cadmium Resistance Genes from a High GC Clone Library by Coupling the Sanger and PacBio Sequencing Strategies

Functional (meta) genomics allows the high-throughput identification of functional genes in a premise-free way. However, it is still difficult to perform Sanger sequencing for high GC DNA templates, which hinders the functional genomic exploration of a high GC genomic library. Here, we developed a procedure to resolve this problem by coupling the Sanger and PacBio sequencing strategies. Identification of cadmium (Cd) resistance genes from a small-insert high GC genomic library was performed to test the procedure. The library was generated from a high GC (75.35%) bacterial genome. Nineteen clones that conferred Cd resistance to Escherichia coli subject to Sanger sequencing directly. The positive clones were in parallel subject to in vivo amplification in host cells, from which recombinant plasmids were extracted and linearized by selected restriction endonucleases. PacBio sequencing was performed to obtain the full-length sequences. As the identities, partial sequences from Sanger sequencing were aligned to the full-length sequences from PacBio sequencing, which led to the identification of seven unique full-length sequences. The unique sequences were further aligned to the full genome sequence of the source strain. Functional screening showed that the identified positive clones were all able to improve Cd resistance of the host cells. The functional genomic procedure developed here couples the Sanger and PacBio sequencing methods and overcomes the difficulties in PCR approaches for high GC DNA. The procedure can be a promising option for the high-throughput sequencing of functional genomic libraries, and realize a cost-effective and time-efficient identification of the positive clones, particularly for high GC genetic materials.

TRANSNAP: a web database providing comprehensive information on Japanese pear transcriptome

Japanese pear (Pyrus pyrifolia) is a major fruit tree in the family Rosaceae and is bred for fruit production. To promote the development of breeding strategies and molecular research for Japanese pear, we sequenced the transcripts of Japanese pear variety 'Hosui'. To exhaustively collect information of total gene expression, RNA samples from various organs and stages of Japanese pear were sequenced by three technologies, single-molecule real-time (SMRT) sequencing, 454 pyrosequencing, and Sanger sequencing. Using all those reads, we determined comprehensive reference sequences of Japanese pear. Then, their protein sequences were predicted, and biological functional annotations were assigned. Finally, we developed a web database, TRANSNAP (http://plantomics.mind.meiji.ac.jp/nashi), which is the first web resource of Japanese pear omics information. This database provides highly reliable information via a user-friendly web interface: the reference sequences, gene functional annotations, and gene expression profiles from microarray experiments. In addition, based on sequence comparisons among Japanese, Chinese and European pears, similar protein sequences among the pears and species-specific proteins in Japanese pear can be quickly and efficiently identified. TRANSNAP will aid molecular research and breeding in Japanese pear, and its information is available for comparative analysis among other pear species and families.

Sequence properties of certain GC rich avian genes, their origins and absence from genome assemblies: case studies

Background

More and more eukaryotic genomes are sequenced and assembled, most of them presented as a complete model in which missing chromosomal regions are filled by Ns and where a few chromosomes may be lacking. Avian genomes often contain sequences with high GC content, which has been hypothesized to be at the origin of many missing sequences in these genomes. We investigated features of these missing sequences to discover why some may not have been integrated into genomic libraries and/or sequenced.

Results

The sequences of five red jungle fowl cDNA models with high GC content were used as queries to search publicly available datasets of Illumina and Pacbio sequencing reads. These were used to reconstruct the leptin, TNFα, MRPL52, PCP2 and PET100 genes, all of which are absent from the red jungle fowl genome model. These gene sequences displayed elevated GC contents, had intron sizes that were sometimes larger than non-avian orthologues, and had non-coding regions that contained numerous tandem and inverted repeat sequences with motifs able to assemble into stable G-quadruplexes and intrastrand dyadic structures. Our results suggest that Illumina technology was unable to sequence the non-coding regions of these genes. On the other hand, PacBio technology was able to sequence these regions, but with dramatically lower efficiency than would typically be expected.

Conclusions

High GC content was not the principal reason why numerous GC-rich regions of avian genomes are missing from genome assembly models. Instead, it is the presence of tandem repeats containing motifs capable of assembling into very stable secondary structures that is likely responsible.

Integrative Meta-Assembly Pipeline (IMAP): Chromosome-level genome assembler combining multiple de novo assemblies

Background

Genomic data have become major resources to understand complex mechanisms at fine-scale temporal and spatial resolution in functional and evolutionary genetic studies, including human diseases, such as cancers. Recently, a large number of whole genomes of evolving populations of yeast (Saccharomyces cerevisiae W303 strain) were sequenced in a time-dependent manner to identify temporal evolutionary patterns. For this type of study, a chromosome-level sequence assembly of the strain or population at time zero is required to compare with the genomes derived later. However, there is no fully automated computational approach in experimental evolution studies to establish the chromosome-level genome assembly using unique features of sequencing data.

Methods and results

In this study, we developed a new software pipeline, the integrative meta-assembly pipeline (IMAP), to build chromosome-level genome sequence assemblies by generating and combining multiple initial assemblies using three de novo assemblers from short-read sequencing data. We significantly improved the continuity and accuracy of the genome assembly using a large collection of sequencing data and hybrid assembly approaches. We validated our pipeline by generating chromosome-level assemblies of yeast strains W303 and SK1, and compared our results with assemblies built using long-read sequencing and various assembly evaluation metrics. We also constructed chromosome-level sequence assemblies of S. cerevisiae strain Sigma1278b, and three commonly used fungal strains: Aspergillus nidulans A713, Neurospora crassa 73, and Thielavia terrestris CBS 492.74, for which long-read sequencing data are not yet available. Finally, we examined the effect of IMAP parameters, such as reference and resolution, on the quality of the final assembly of the yeast strains W303 and SK1.

Conclusions

We developed a cost-effective pipeline to generate chromosome-level sequence assemblies using only short-read sequencing data. Our pipeline combines the strengths of reference-guided and meta-assembly approaches. Our pipeline is available online at http://github.com/jkimlab/IMAP including a Docker image, as well as a Perl script, to help users install the IMAP package, including several prerequisite programs. Users can use IMAP to easily build the chromosome-level assembly for the genome of their interest.

A new full-length circular DNA sequencing method for viral-sized genomes reveals that RNAi transgenic plants provoke a shift in geminivirus populations in the field

We present a new method, CIDER-Seq (Circular DNA Enrichment sequencing) for the unbiased enrichment and long-read sequencing of viral-sized circular DNA molecules. We used CIDER-Seq to produce single-read full-length virus genomes for the first time. CIDER-Seq combines PCR-free virus enrichment with Single Molecule Real Time sequencing and a new sequence de-concatenation algorithm. We apply our technique to produce >1200 full-length, highly accurate geminivirus genomes from RNAi-transgenic and control plants in a field trial in Kenya. Using CIDER-Seq we can demonstrate for the first time that the expression of antiviral double-stranded RNA (dsRNA) in transgenic plants causes a consistent shift in virus populations towards species sharing low homology to the transgene derived dsRNA. Our method and its application in an economically important crop plant opens new possibilities in periodic virus sequence surveillance and accurate profiling of diverse circular DNA elements.

RNA Sequencing Data: Hitchhiker's Guide to Expression Analysis

Gene expression is the fundamental level at which the results of various genetic and regulatory programs are observable. The measurement of transcriptome-wide gene expression has convincingly switched from microarrays to sequencing in a matter of years. RNA sequencing (RNA-seq) provides a quantitative and open system for profiling transcriptional outcomes on a large scale and therefore facilitates a large diversity of applications, including basic science studies, but also agricultural or clinical situations. In the past 10 years or so, much has been learned about the characteristics of the RNA-seq data sets, as well as the performance of the myriad of methods developed. In this review, we give an overview of the developments in RNA-seq data analysis, including experimental design, with an explicit focus on the quantification of gene expression and statistical approachesfor differential expression. We also highlight emerging data types, such as single-cell RNA-seq and gene expression profiling using long-read technologies.

Complete Genome Sequences of Colwellia sp. Arc7-635, a Denitrifying Bacterium Isolated from Arctic Seawater

Colwellia sp. Arc7-635, a psychrophilic denitrifying bacterium isolated from Arctic seawater, uses NO₃^- or NH₄⁺ as the sole nitrogen source to grow at low temperatures. In this article, we describe the complete genome of Colwellia sp. Arc7-635. The genome has one circular chromosome of 4,741,350 bp (38.41 mol% G+C content), consisting of 3841 coding genes, 91 tRNA genes, as well as seven rRNA operons of 16S-23S-5S rRNA, and one operon of 16S-23S-5S-5S rRNA. According to the genomic annotation results, strain Colwellia sp. Arc7-635 encodes a complete denitrifying pathway consisting of genes affiliated with nitrate reductase, nitrite reductase, nitric oxide reductase, and nitrous oxide reductase. Genes affiliated with nitrate reduction to ammonia including nitrate reductases (NapA and NapB) and nitrite reductases (NirA, NirB, and NirD) were also identified. The whole genome sequences of Arc7-635 provide information that is useful for further clarifying nitrogen metabolisms and facilitate its potential applications in the bioremediation of nitrogen pollutions.

100 Years of evolving gene-disease complexities and scientific debutants

It's been over 100 years since the word `gene' is around and progressively evolving in several scientific directions. Time-to-time technological advancements have heavily revolutionized the field of genomics, especially when it's about, e.g. triple code development, gene number proposition, genetic mapping, data banks, gene-disease maps, catalogs of human genes and genetic disorders, CRISPR/Cas9, big data and next generation sequencing, etc. In this manuscript, we present the progress of genomics from pea plant genetics to the human genome project and highlight the molecular, technical and computational developments. Studying genome and epigenome led to the fundamentals of development and progression of human diseases, which includes chromosomal, monogenic, multifactorial and mitochondrial diseases. World Health Organization has classified, standardized and maintained all human diseases, when many academic and commercial online systems are sharing information about genes and linking to associated diseases. To efficiently fathom the wealth of this biological data, there is a crucial need to generate appropriate gene annotation repositories and resources. Our focus has been how many gene-disease databases are available worldwide and which sources are authentic, timely updated and recommended for research and clinical purposes. In this manuscript, we have discussed and compared 43 such databases and bioinformatics applications, which enable users to connect, explore and, if possible, download gene-disease data.

Genomic repeats, misassembly and reannotation: a case study with long-read resequencing of Porphyromonas gingivalis reference strains

BACKGROUND

Without knowledge of their genomic sequences, it is impossible to make functional models of the bacteria that make up human and animal microbiota. Unfortunately, the vast majority of publicly available genomes are only working drafts, an incompleteness that causes numerous problems and constitutes a major obstacle to genotypic and phenotypic interpretation. In this work, we began with an example from the class Bacteroidia in the phylum Bacteroidetes, which is preponderant among human orodigestive microbiota. We successfully identify the genetic loci responsible for assembly breaks and misassemblies and demonstrate the importance and usefulness of long-read sequencing and curated reannotation.

RESULTS

We showed that the fragmentation in Bacteroidia draft genomes assembled from massively parallel sequencing linearly correlates with genomic repeats of the same or greater size than the reads. We also demonstrated that some of these repeats, especially the long ones, correspond to misassembled loci in three reference Porphyromonas gingivalis genomes marked as circularized (thus complete or finished). We prove that even at modest coverage (30X), long-read resequencing together with PCR contiguity verification (rrn operons and an integrative and conjugative element or ICE) can be used to identify and correct the wrongly combined or assembled regions. Finally, although time-consuming and labor-intensive, consistent manual biocuration of three P. gingivalis strains allowed us to compare and correct the existing genomic annotations, resulting in a more accurate interpretation of the genomic differences among these strains.

CONCLUSIONS

In this study, we demonstrate the usefulness and importance of long-read sequencing in verifying published genomes (even when complete) and generating assemblies for new bacterial strains/species with high genomic plasticity. We also show that when combined with biological validation processes and diligent biocurated annotation, this strategy helps reduce the propagation of errors in shared databases, thus limiting false conclusions based on incomplete or misleading information.

The Pelagic Bacterium Paraphotobacterium marinum Has the Smallest Complete Genome Within the Family Vibrionaceae

Members of the family Vibrionaceae are metabolically versatile and ubiquitous in natural environments, with extraordinary genome feature of two chromosomes. Here we reported the complete genome of Paraphotobacterium marinum NSCS20N07D^T, a recently described novel genus-level species in the family Vibrionaceae. It contained two circular chromosomes with a size of 2,593,992 bp with G+C content of 31.2 mol%, and a plasmid with a size of 5,539 bp. The larger chromosome (Chr. I) had a genome size of 1,426,504 bp with G+C content of 31.6 mol%, and the smaller one (Chr. II) had a genome size of 1,161,949 bp with G+C content of 30.8 mol%. The two chromosomes have strikingly similar G+C contents with difference of <1% and similar percentages of coding regions. Interestingly, by comparison to 134 species affiliated with seven genera within the family Vibrionaceae, P. marinum NSCS20N07D^T possessed the smallest genome size and lowest G+C content. Clusters of orthologous groups of proteins functional categories revealed that the two chromosomes had different distributions of functional classes, indicating they take different cellular functions. Surprisingly, Chr. II had a large proportion of unknown genes than Chr. I. Metabolic characteristics predicted that Chr. I performed the essential metabolism, which can be complemented by the Chr. II, such as amino acids biosynthesis. Microbial community analysis of in situ surface seawater revealed that P. marinum accounted for one to four sequences among more than 20,000 of 16S ribosomal RNA gene V4 contigs, representing it apparently appeared as a rare species. What’s more, P. marinum was anticipated to be specific to the pelagic ocean. This study will provide new insight into more understanding the genomic and metabolic features of multiple chromosomes in prokaryote and emphasize the ecological distribution of the members in the family Vibrionaceae as a rare species.