Consensus generation and variant detection by Celera Assembler

Advancements in long-read genome sequencing technologies and algorithms

The recent advent of long read sequencing technologies, such as Pacific Biosciences (PacBio) and Oxford Nanopore technology (ONT), have led to substantial improvements in accuracy and computational cost in sequencing genomes. However, de novo whole-genome assembly still presents significant challenges related to the quality of the results. Pursuing de novo whole-genome assembly remains a formidable challenge, underscored by intricate considerations surrounding computational demands and result quality. As sequencing accuracy and throughput steadily advance, a continuous stream of innovative assembly tools floods the field. Navigating this dynamic landscape necessitates a reasonable choice of sequencing platform, depth, and assembly tools to orchestrate high-quality genome reconstructions. This comprehensive review delves into the intricate interplay between cutting-edge long read sequencing technologies, assembly methodologies, and the ever-evolving field of genomics. With a focus on addressing the pivotal challenges and harnessing the opportunities presented by these advancements, we provide an in-depth exploration of the crucial factors influencing the selection of optimal strategies for achieving robust and insightful genome assemblies.

Demographic dynamics and molecular evolution of the rare and endangered subsect. Gerardianae of Pinus: insights from chloroplast genomes and mitochondrial DNA markers

MAIN CONCLUSION

The divergence of subsect. Gerardianae was likely triggered by the uplift of the Qinghai-Tibetan Plateau and adjacent mountains. Pinus bungeana might have probably experienced expansion since Last Interglacial period. Historical geological and climatic oscillations have profoundly affected patterns of nucleotide variability, evolutionary history, and species divergence in numerous plants of the Northern Hemisphere. However, how long-lived conifers responded to geological and climatic fluctuations in East Asia remain poorly understood. Here, based on paternally inherited chloroplast genomes and maternally inherited mitochondrial DNA markers, we investigated the population demographic history and molecular evolution of subsect. Gerardianae (only including three species, Pinus bungeana, P. gerardiana, and P. squamata) of Pinus. A low level of nucleotide diversity was found in P. bungeana (π was 0.00016 in chloroplast DNA sequences, and 0.00304 in mitochondrial DNAs). The haplotype-based phylogenetic topology and unimodal distributions of demographic analysis suggested that P. bungeana probably originated in the southern Qinling Mountains and experienced rapid population expansion since Last Interglacial period. Phylogenetic analysis revealed that P. gerardiana and P. squamata had closer genetic relationship. The species divergence of subsect. Gerardianae occurred about 27.18 million years ago (Mya) during the middle to late Oligocene, which was significantly associated with the uplift of the Qinghai-Tibetan Plateau and adjacent mountains from the Eocene to the mid-Pliocene. The molecular evolutionary analysis showed that two chloroplast genes (psaI and ycf1) were under positive selection, the genetic lineages of P. bungeana exhibited higher transition and nonsynonymous mutations, which were involved with the strongly environmental adaptation. These findings shed light on the population evolutionary history of white pine species and provide striking insights for comprehension of their species divergence and molecular evolution.

Unraveling the phylogenomic diversity of Methanomassiliicoccales and implications for mitigating ruminant methane emissions

BACKGROUND

Methanomassiliicoccales are a recently identified order of methanogens that are diverse across global environments particularly the gastrointestinal tracts of animals; however, their metabolic capacities are defined via a limited number of cultured strains.

RESULTS

Here, we profile and analyze 243 Methanomassiliicoccales genomes assembled from cultured representatives and uncultured metagenomes recovered from various biomes, including the gastrointestinal tracts of different animal species. Our analyses reveal the presence of numerous undefined genera and genetic variability in metabolic capabilities within Methanomassiliicoccales lineages, which is essential for adaptation to their ecological niches. In particular, gastrointestinal tract Methanomassiliicoccales demonstrate the presence of co-diversified members with their hosts over evolutionary timescales and likely originated in the natural environment. We highlight the presence of diverse clades of vitamin transporter BtuC proteins that distinguish Methanomassiliicoccales from other archaeal orders and likely provide a competitive advantage in efficiently handling B12. Furthermore, genome-centric metatranscriptomic analysis of ruminants with varying methane yields reveal elevated expression of select Methanomassiliicoccales genera in low methane animals and suggest that B12 exchanges could enable them to occupy ecological niches that possibly alter the direction of H2 utilization.

CONCLUSIONS

We provide a comprehensive and updated account of divergent Methanomassiliicoccales lineages, drawing from numerous uncultured genomes obtained from various habitats. We also highlight their unique metabolic capabilities involving B12, which could serve as promising targets for mitigating ruminant methane emissions by altering H2 flow.

Thermococcus thermotolerans sp. nov., a hyperthermophilic archaeon isolated from a chimney in the Southwest Indian Ocean

An anaerobic hyperthermophilic archaeon was isolated from a black smoker chimney with a snail attachment at a water depth of 2 739 m in the Southwest Indian Ocean. The sample was taken from the chimney exterior wall. The enrichment was conducted under a continuous culture with temperature fluctuation of 80-130 °C over 24 h for 42 days at 28 MPa. The isolation was performed at 90 °C at 0.1 MPa. Cells of the isolated strain 813A4T were irregular cocci. Strain 813A4T grew at 60-94 °C (optimal growth at 85 °C) at 0.1 MPa, and growth was detected at up to 99 °C at 28 MPa. At 85 °C, the strain was able to grow at pressures ranging from 0.1 to 110 MPa (optimal pressure, 0.1-40 MPa). At 85 °C, the cells of 813A4T grew at pH 5.5-9 (optimal, pH 7.0) and a NaCl concentration of 1.0-4.0 % (w/v; optimum concentration, 2.5 % NaCl). Strain 813A4T utilized yeast extract, tryptone and peptone as single carbon sources for growth. Elemental sulphur stimulated its growth. The G+C content of the complete genome was 53.48 mol%. Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain 813A4T belonged to the genus Thermococcus, with the highest sequence similarity to Thermococcus barossii SHCK-94T (99.73 %). The average nucleotide identity between strains 813A4T and SHCK-94T was 82.56 %. All these data indicated that strain 813A4T should be classified as representing a novel species of the genus Thermococcus, for which Thermococcus thermotolerans sp. nov. is proposed. The type strain is 813A4T (=JCM 39367T=MCCC M28628T).

Convergent reductive evolution and host adaptation in Mycoavidus bacterial endosymbionts of Mortierellaceae fungi

Intimate associations between fungi and intracellular bacterial endosymbionts are becoming increasingly well understood. Phylogenetic analyses demonstrate that bacterial endosymbionts of Mucoromycota fungi are related either to free-living Burkholderia or Mollicutes species. The so-called Burkholderia-related endosymbionts or BRE comprise Mycoavidus, Mycetohabitans and Candidatus Glomeribacter gigasporarum. These endosymbionts are marked by genome contraction thought to be associated with intracellular selection. However, the conclusions drawn thus far are based on a very small subset of endosymbiont genomes, and the mechanisms leading to genome streamlining are not well understood. The purpose of this study was to better understand how intracellular existence shapes Mycoavidus and BRE functionally at the genome level. To this end we generated and analyzed 14 novel draft genomes for Mycoavidus living within the hyphae of Mortierellomycotina fungi. We found that our novel Mycoavidus genomes were significantly reduced compared to free-living Burkholderiales relatives. Using a genome-scale phylogenetic approach including the novel and available existing genomes of Mycoavidus, we show that the genus is an assemblage composed of two independently derived lineages including three well supported clades of Mycoavidus. Using a comparative genomic approach, we shed light on the functional implications of genome reduction, documenting shared and unique gene loss patterns between the three Mycoavidus clades. We found that many endosymbiont isolates demonstrate patterns of vertical transmission and host-specificity, but others are present in phylogenetically disparate hosts. We discuss how reductive evolution and host specificity reflect convergent adaptation to the intrahyphal selective landscape, and commonalities of eukaryotic endosymbiont genome evolution.

Systematic analyses with genomic and metabolomic insights reveal a new species, Ophiocordyceps indica sp. nov. from treeline area of Indian Western Himalayan region

Ophiocordyceps is a species-rich genus in the order Hypocreales (Sordariomycetes, Ascomycota) depicting a fascinating relationship between microbes and insects. In the present study, a new species, Ophiocordyceps indica sp. nov., is discovered infecting lepidopteran larvae from tree line locations (2,202-2,653 m AMSL) of the Kullu District, Himachal Pradesh, Indian Western Himalayan region, using combinations of morphological and molecular phylogenetic analyses. A phylogeny for Ophiocordyceps based on a combined multigene (nrSSU, nrLSU, tef-1α, and RPB1) dataset is provided, and its taxonomic status within Ophiocordycipitaceae is briefly discussed. Its genome size (~59 Mb) revealed 94% genetic similarity with O. sinensis; however, it differs from other extant Ophiocordyceps species based on morphological characteristics, molecular phylogenetic relationships, and genetic distance. O. indica is identified as the second homothallic species in the family Ophiocordycipitaceae, after O. sinensis. The presence of targeted marker components, viz. nucleosides (2,303.25 μg/g), amino acids (6.15%), mannitol (10.13%), and biological activity data, suggests it to be a new potential source of nutraceutical importance. Data generated around this economically important species will expand our understanding regarding the diversity of Ophiocordyceps-like taxa from new locations, thus providing new research avenues.

Identification of a deep-branching thermophilic clade sheds light on early bacterial evolution

It has been proposed that early bacteria, or even the last universal common ancestor of all cells, were thermophilic. However, research on the origin and evolution of thermophily is hampered by the difficulties associated with the isolation of deep-branching thermophilic microorganisms in pure culture. Here, we isolate a deep-branching thermophilic bacterium from a deep-sea hydrothermal vent, using a two-step cultivation strategy ("Subtraction-Suboptimal", StS) designed to isolate rare organisms. The bacterium, which we name Zhurongbacter thermophilus 3DAC, is a sulfur-reducing heterotroph that is phylogenetically related to Coprothermobacterota and other thermophilic bacterial groups, forming a clade that seems to represent a major, early-diverging bacterial lineage. The ancestor of this clade might be a thermophilic, strictly anaerobic, motile, hydrogen-dependent, and mixotrophic bacterium. Thus, our study provides insights into the early evolution of thermophilic bacteria.

Description of Brachybacterium sillae sp. nov., a thermophilic bacterium isolated from a hot spring

The taxonomic position of strain EF45031^T, isolated from the Neungam Carbonate hot spring, was examined using the polyphasic taxonomic approach. Strain EF45031^T shared the highest percentage of 16S rRNA gene sequence with Brachybacterium nesterenkovii CIP 104813 ^T (97.7%). The average nucleotide identity (ANI), average amino acid identity (AAI), and digital DNA-DNA hybridization (dDDH) values between strain EF45031^T and the type strains B. nesterenkovii CIP 104813 ^T and B. phenoliresistens Phenol-A^T were 77.0%, 69.15%, 21.9% and 75.73%, 68.81%, 20.5%, respectively. Phylogenomic analysis using an up-to-date bacterial core gene (UBCG) set revealed that strain EF45031^T belonged to the genus Brachybacterium. Growth occurred between 25 and 50 ℃ at pH 6.0-9.0 and could tolerate salinity up to 5% (w/v). Strain had anteiso-C_15:0 and anteiso-C_17:0 as major fatty acids. Menaquinone-7 (MK-7) was the predominant respiratory menaquinone. The polar lipids were diphosphatidylglycerol, phosphatidylglycerol, three aminolipids, and two unidentified glycolipids. The cell-wall peptidoglycan contained meso-diaminopimelic acid as a diagnostic diamino acid. The genome comprised 2,663,796 bp, with a G + C content of 70.9%. Stress-responsive periplasmic chaperone/protease coding genes were identified in the genome of EF45031^T and were not detected in other Brachybacterium species. The polyphasic taxonomic properties indicate that the strain represents a novel species within the genus Brachybacterium, for which the name Brachybacterium sillae sp. nov. is proposed. The type strain is EF45031^T (= KCTC 49702 ^T = NBRC 115869 ^T).

Mass Spectrometry-Guided Discovery of Multi-N-Methylated Cyclodecapeptides Auyuittuqamides E-H from Sesquicillium sp. QL0466

Mass spectrometry-based dereplication and prioritization led to the discovery of four multi-N-methylated cyclodecapeptides, auyuittuqamides E-H (1-4), from a soil-derived Sesquicillium sp. The planar structures of these compounds were elucidated based on analysis of HRESIMS and NMR data. Absolute configurations of the chiral amino acid residues were assigned by a combination of the advanced Marfey's method, chiral-phase LC-MS analysis, and J-based configuration analysis, revealing that 1-4 contain both d- and l-isomers of N-methylleucine (MeLeu). Differentiation of d- and l-MeLeu in the sequence was achieved by advanced Marfey's analysis of the diagnostic peptide fragments generated from partial hydrolysis of 1. Bioinformatic analysis identified a putative biosynthetic gene cluster (auy) for auyuittuqamides E-H, and a plausible biosynthetic pathway was proposed. These newly identified fungal cyclodecapeptides (1-4) displayed in vitro growth inhibitory activity against vancomycin-resistant Enterococcus faecium with MIC values of 8 μg/mL.

Comparative genome analysis of Sesbania cannabina-nodulating Rhizobium spp. revealing the symbiotic and transferrable characteristics of symbiosis plasmids

Symbiotic nitrogen fixation between legumes and rhizobia makes a great contribution to the terrestrial ecosystem. The successful symbiosis between the partners mainly depends on the nod and nif genes in rhizobia, while the specific symbiosis is mainly determined by the structure of Nod factors and the corresponding secretion systems (type III secretion system; T3SS), etc. These symbiosis genes are usually located on symbiotic plasmids or a chromosomal symbiotic island, both could be transferred interspecies. In our previous studies, Sesbania cannabina-nodulating rhizobia across the world were classified into 16 species of four genera and all the strains, especially those of Rhizobium spp., harboured extraordinarily highly conserved symbiosis genes, suggesting that horizontal transfer of symbiosis genes might have happened among them. In order to learn the genomic basis of diversification of rhizobia under the selection of host specificity, we performed this study to compare the complete genome sequences of four Rhizobium strains associated with S. cannabina, YTUBH007, YTUZZ027, YTUHZ044 and YTUHZ045. Their complete genomes were sequenced and assembled at the replicon level. Each strain represents a different species according to the average nucleotide identity (ANI) values calculated using the whole-genome sequences; furthermore, except for YTUBH007, which was classified as Rhizobium binae, the remaining three strains were identified as new candidate species. A single symbiotic plasmid sized 345-402 kb containing complete nod, nif, fix, T3SS and conjugal transfer genes was detected in each strain. The high ANI and amino acid identity (AAI) values, as well as the close phylogenetic relationships among the entire symbiotic plasmid sequences, indicate that they have the same origin and the entire plasmid has been transferred among different Rhizobium species. These results indicate that S. cannabina stringently selects a certain symbiosis gene background of the rhizobia for nodulation, which might have forced the symbiosis genes to transfer from some introduced rhizobia to the related native or local-condition-adapted bacteria. The existence of almost complete conjugal transfer related elements, but not the gene virD, indicated that the self-transfer of the symbiotic plasmid in these rhizobial strains may be realized via a virD-independent pathway or through another unidentified gene. This study provides insight for the better understanding of high-frequency symbiotic plasmid transfer, host-specific nodulation and the host shift for rhizobia.

Identification of a 17β-estradiol-degrading Microbacterium hominis SJTG1 with high adaptability and characterization of the genes for estrogen degradation

Environmental estrogen contamination poses severe threat to wildlife and human. Biodegradation is an efficient strategy to remove the wide-spread natural estrogen, while strains suitable for hostile environments and fit for practical application are rare. In this work, Microbacterium hominis SJTG1 was isolated and identified with high degrading efficiency for 17β-estradiol (E2) and great environment fitness. It could degrade nearly 100% of 10 mg/L E2 in minimal medium in 6 days, and remove 93% of 1 mg/L E2 and 74% of 10 mg/L E2 in the simulated E2-polluted solid soil in 10 days. It maintained stable E2-degrading efficiency in various harsh conditions like non-neutral pH, high salinity, stress of heavy metals and surfactants. Genome mining and comparative genome analysis revealed that there are multiple genes potentially associated with steroid degradation in strain SJTG1. One 3β/17β-hydroxysteroid dehydrogenase HSD-G129 induced by E2 catalyzed the 3β/17β-dehydrogenation of E2 and other steroids efficiently. The transcription of hsd-G129 gene was negatively regulated by the adjacent LysR-type transcriptional regulator LysR-G128, through specific binding to the conserved site. E2 can release this binding and initiate the degradation process. This work provides an efficient and adaptive E2-degrading strain and promotes the biodegrading mechanism study and actual remediation application.

Engineering Vibrio alginolyticus as a novel chassis for PHB production from starch

Vibrio alginolyticus LHF01 was engineered to efficiently produce poly-3-hydroxybutyrate (PHB) from starch in this study. Firstly, the ability of Vibrio alginolyticus LHF01 to directly accumulate PHB using soluble starch as the carbon source was explored, and the highest PHB titer of 2.06 g/L was obtained in 18 h shake flask cultivation. Then, with the analysis of genomic information of V. alginolyticus LHF01, the PHB synthesis operon and amylase genes were identified. Subsequently, the effects of overexpressing PHB synthesis operon and amylase on PHB production were studied. Especially, with the co-expression of PHB synthesis operon and amylase, the starch consumption rate was improved and the PHB titer was more than doubled. The addition of 20 g/L insoluble corn starch could be exhausted in 6-7 h cultivation, and the PHB titer was 4.32 g/L. To the best of our knowledge, V. alginolyticus was firstly engineered to produce PHB with the direct utilization of starch, and this stain can be considered as a novel host to produce PHB using starch as the raw material.

Era of gapless plant genomes: innovations in sequencing and mapping technologies revolutionize genomics and breeding

Whole-genome sequencing and assembly have revolutionized plant genetics and molecular biology over the last two decades. However, significant shortcomings in first- and second-generation technology resulted in imperfect reference genomes: numerous and large gaps of low quality or undeterminable sequence in areas of highly repetitive DNA along with limited chromosomal phasing restricted the ability of researchers to characterize regulatory noncoding elements and genic regions that underwent recent duplication events. Recently, advances in long-read sequencing have resulted in the first gapless, telomere-to-telomere (T2T) assemblies of plant genomes. This leap forward has the potential to increase the speed and confidence of genomics and molecular experimentation while reducing costs for the research community.

Description of Microbacterium neungamense sp. nov. isolated from a hot spring

The Gram-positive, nonmotile, rod-shaped bacterium EF45044^T was isolated from a hot spring in Chungju, South Korea. The strain was able to grow at concentrations of 0‒5% (w/v) NaCl, at pH 6.0‒10.0 and in the temperature range of 18‒50 °C. Strain EF45044^T showed the highest 16S rRNA gene sequence similarity (98.2%) with Microbacterium ketosireducens DSM 12510^T, and the digital DNA‒DNA hybridization (dDDH), average amino acid identity (AAI), and average nucleotide identity (ANI) values were all lower than the accepted species threshold. Strain EF45044^T contained MK‒12 and MK‒13 as the predominant respiratory quinones and anteiso‒C_17:0, anteiso‒C_15:0, and iso‒C_16:0 as the major fatty acids. Diphosphatidylglycerol, phosphatidylglycerol, and glycolipid were detected as the major polar lipids. The cell-wall peptidoglycan contained ornithine. The DNA G + C content was 71.4 mol%. Based on the polyphasic data, strain EF45044^T (= KCTC 49703^T) presents a novel species of the genus Microbacterium, for which the name Microbacterium neungamense sp. nov. is proposed.

Complete chloroplast genome of Ilex dabieshanensis: Genome structure, comparative analyses with three traditional Ilex tea species, and its phylogenetic relationships within the family Aquifoliaceae

Ilex dabieshanensis K. Yao & M. B. Deng is not only a highly valued tree species for landscaping, it is also a good material for making kuding tea due to its anti-inflammatory and lipid-lowering medicinal properties. Utilizing next-generation and long-read sequencing technologies, we assembled the whole chloroplast genome of I. dabieshanensis. The genome was 157,218 bp in length, exhibiting a typical quadripartite structure with a large single copy (LSC: 86,607 bp), a small single copy (SSC: 18,427 bp) and a pair of inverted repeat regions (IRA and IRB: each of 26,092 bp). A total of 121 predicted genes were encoded, including 113 distinctive (79 protein-coding genes, 30 tRNAs, and 4 rRNAs) and 8 duplicated (8 protein-coding genes) located in the IR regions. Overall, 132 SSRs and 43 long repeats were detected and could be used as potential molecular markers. Comparative analyses of four traditional Ilex tea species (I. dabieshanensis, I. paraguariensis, I. latifolia and I. cornuta) revealed seven divergent regions: matK-rps16, trnS-psbZ, trnT-trnL, atpB-rbcL, petB-petD, rpl14-rpl16, and rpl32-trnL. These variations might be applicable for distinguishing different species within the genus Ilex. Phylogenetic reconstruction strongly suggested that I. dabieshanensis formed a sister clade to I. cornuta and also showed a close relationship to I. latifolia. The generated chloroplast genome information in our study is significant for Ilex tea germplasm identification, phylogeny and genetic improvement.

Parallel decay of vision genes in subterranean water beetles

In the framework of neutral theory of molecular evolution, genes specific to the development and function of eyes in subterranean animals living in permanent darkness are expected to evolve by relaxed selection, ultimately becoming pseudogenes. However, definitive empirical evidence for the role of neutral processes in the loss of vision over evolutionary time remains controversial. In previous studies, we characterized an assemblage of independently-evolved water beetle (Dytiscidae) species from a subterranean archipelago in Western Australia, where parallel vision and eye loss have occurred. Using a combination of transcriptomics and exon capture, we present evidence of parallel coding sequence decay, resulting from the accumulation of frameshift mutations and premature stop codons, in eight phototransduction genes (arrestins, opsins, ninaC and transient receptor potential channel genes) in 32 subterranean species in contrast to surface species, where these genes have open reading frames. Our results provide strong evidence to support neutral evolutionary processes as a major contributing factor to the loss of phototransduction genes in subterranean animals, with the ultimate fate being the irreversible loss of a light detection system.

Natural selenium stress influences the changes of antibiotic resistome in seleniferous forest soils

BACKGROUND

Metal(loid)s can promote the spread and enrichment of antibiotic resistance genes (ARGs) in the environment through a co-selection effect. However, it remains unclear whether exposure of microorganisms to varying concentrations of selenium (Se), an essential but potentially deleterious metal(loid) to living organisms, can influence the migration and distribution of ARGs in forest soils.

RESULTS

Precisely 235 ARGs conferring resistance to seven classes of antibiotics were detected along a Se gradient (0.06-20.65 mg kg^-1) across 24 forest soils. (flor)/(chlor)/(am)phenicol resistance genes were the most abundant in all samples. The total abundance of ARGs first increased and then decreased with an elevated available Se content threshold of 0.034 mg kg^-1 (P = 2E-05). A structural equation model revealed that the dominant mechanism through which Se indirectly influences the vertical migration of ARGs is by regulating the abundance of the bacterial community. In addition, the methylation of Se (mediated by tehB) and the repairing of DNA damages (mediated by ruvB and recG) were the dominant mechanisms involved in Se resistance in the forest soils. The co-occurrence network analysis revealed a significant correlated cluster between Se-resistance genes, MGEs and ARGs, suggesting the co-transfer potential. Lelliottia amnigena YTB01 isolated from the soil was able to tolerate 50 μg mL^-1 ampicillin and 1000 mg kg^-1 sodium selenite, and harbored both Se resistant genes and ARGs in the genome.

CONCLUSIONS

Our study demonstrated that the spread and enrichment of ARGs are enhanced under moderate Se pressure but inhibited under severe Se pressure in the forest soil (threshold at 0.034 mg kg^-1 available Se content). The data generated in this pilot study points to the potential health risk associated with Se contamination and its associated influence on ARGs distribution in soil.

SLDMS: A Tool for Calculating the Overlapping Regions of Sequences

In the field of genome assembly, contig assembly is one of the most important parts. Contig assembly requires the processing of overlapping regions of a large number of DNA sequences and this calculation usually takes a lot of time. The time consumption of contig assembly algorithms is an important indicator to evaluate the degree of algorithm superiority. Existing methods for processing overlapping regions of sequences consume too much in terms of running time. Therefore, we propose a method SLDMS for processing sequence overlapping regions based on suffix array and monotonic stack, which can effectively improve the efficiency of sequence overlapping regions processing. The running time of the SLDMS is much less than that of Canu and Flye in dealing with the sequence overlap interval and in some data with most sequencing errors occur at both the ends of the sequencing data, the running time of the SLDMS is only about one-tenth of the other two methods.

Complete Chloroplast Genome Analysis of Two Important Medicinal Alpinia Species: Alpinia galanga and Alpinia kwangsiensis

Most Alpinia species are valued as foods, ornamental plants, or plants with medicinal properties. However, morphological characteristics and commonly used DNA barcode fragments are not sufficient for accurately identifying Alpinia species. Difficulties in species identification have led to confusion in the sale and use of Alpinia for medicinal use. To mine resources and improve the molecular methods for distinguishing among Alpinia species, we report the complete chloroplast (CP) genomes of Alpinia galanga and Alpinia kwangsiensis species, obtained via high-throughput Illumina sequencing. The CP genomes of A. galanga and A. kwangsiensis exhibited a typical circular tetramerous structure, including a large single-copy region (87,565 and 87,732 bp, respectively), a small single-copy region (17,909 and 15,181 bp, respectively), and a pair of inverted repeats (27,313 and 29,705 bp, respectively). The guanine-cytosine content of the CP genomes is 36.26 and 36.15%, respectively. Furthermore, each CP genome contained 133 genes, including 87 protein-coding genes, 38 distinct tRNA genes, and 8 distinct rRNA genes. We identified 110 and 125 simple sequence repeats in the CP genomes of A. galanga and A. kwangsiensis, respectively. We then combined these data with publicly available CP genome data from four other Alpinia species (A. hainanensis, A. oxyphylla, A. pumila, and A. zerumbet) and analyzed their sequence characteristics. Nucleotide diversity was analyzed based on the alignment of the complete CP genome sequences, and five candidate highly variable site markers (trnS-trnG, trnC-petN, rpl32-trnL, psaC-ndhE, and ndhC-trnV) were found. Twenty-eight complete CP genome sequences belonging to Alpinieae species were used to construct phylogenetic trees. The results fully demonstrated the phylogenetic relationship among the genera of the Alpinieae, and further proved that Alpinia is a non-monophyletic group. The complete CP genomes of the two medicinal Alpinia species provides lays the foundation for the use of CP genomes in species identification and phylogenetic analyses of Alpinia species.

Molecular evolution of chloroplast genomes in subfamily Zingiberoideae (Zingiberaceae)

BACKGROUND

Zingiberoideae is a large and diverse subfamily of the family Zingiberaceae. Four genera in subfamily Zingiberoideae each possess 50 or more species, including Globba (100), Hedychium (> 80), Kaempferia (50) and Zingiber (150). Despite the agricultural, medicinal and horticultural importance of these species, genomic resources and suitable molecular markers for them are currently sparse.

RESULTS

Here, we have sequenced, assembled and analyzed ten complete chloroplast genomes from nine species of subfamily Zingiberoideae: Globba lancangensis, Globba marantina, Globba multiflora, Globba schomburgkii, Globba schomburgkii var. angustata, Hedychium coccineum, Hedychium neocarneum, Kaempferia rotunda 'Red Leaf', Kaempferia rotunda 'Silver Diamonds' and Zingiber recurvatum. These ten chloroplast genomes (size range 162,630-163,968 bp) possess typical quadripartite structures that consist of a large single copy (LSC, 87,172-88,632 bp), a small single copy (SSC, 15,393-15,917 bp) and a pair of inverted repeats (IRs, 29,673-29,833 bp). The genomes contain 111-113 different genes, including 79 protein coding genes, 28-30 tRNAs and 4 rRNA genes. The dynamics of the genome structures, gene contents, amino acid frequencies, codon usage patterns, RNA editing sites, simple sequence repeats and long repeats exhibit similarities, with slight differences observed among the ten genomes. Further comparative analysis of seventeen related Zingiberoideae species, 12 divergent hotspots are identified. Positive selection is observed in 14 protein coding genes, including accD, ccsA, ndhA, ndhB, psbJ, rbcL, rpl20, rpoC1, rpoC2, rps12, rps18, ycf1, ycf2 and ycf4. Phylogenetic analyses, based on the complete chloroplast-derived single-nucleotide polymorphism data, strongly support that Globba, Hedychium, and Curcuma I + "the Kaempferia clade" consisting of Curcuma II, Kaempferia and Zingiber, form a nested evolutionary relationship in subfamily Zingiberoideae.

CONCLUSIONS

Our study provides detailed information on ten complete Zingiberoideae chloroplast genomes, representing a valuable resource for future studies that seek to understand the molecular evolutionary dynamics in family Zingiberaceae. The identified divergent hotspots can be used for development of molecular markers for phylogenetic inference and species identification among closely related species within four genera of Globba, Hedychium, Kaempferia and Zingiber in subfamily Zingiberoideae.

Genome-resolved metagenomics using environmental and clinical samples

Recent advances in high-throughput sequencing technologies and computational methods have added a new dimension to metagenomic data analysis i.e. genome-resolved metagenomics. In general terms, it refers to the recovery of draft or high-quality microbial genomes and their taxonomic classification and functional annotation. In recent years, several studies have utilized the genome-resolved metagenome analysis approach and identified previously unknown microbial species from human and environmental metagenomes. In this review, we describe genome-resolved metagenome analysis as a series of four necessary steps: (i) preprocessing of the sequencing reads, (ii) de novo metagenome assembly, (iii) genome binning and (iv) taxonomic and functional analysis of the recovered genomes. For each of these four steps, we discuss the most commonly used tools and the currently available pipelines to guide the scientific community in the recovery and subsequent analyses of genomes from any metagenome sample. Furthermore, we also discuss the tools required for validation of assembly quality as well as for improving quality of the recovered genomes. We also highlight the currently available pipelines that can be used to automate the whole analysis without having advanced bioinformatics knowledge. Finally, we will highlight the most widely adapted and actively maintained tools and pipelines that can be helpful to the scientific community in decision making before they commence the analysis.

Genomic insights into the origin, domestication and diversification of Brassica juncea

Despite early domestication around 3000 BC, the evolutionary history of the ancient allotetraploid species Brassica juncea (L.) Czern & Coss remains uncertain. Here, we report a chromosome-scale de novo assembly of a yellow-seeded B. juncea genome by integrating long-read and short-read sequencing, optical mapping and Hi-C technologies. Nuclear and organelle phylogenies of 480 accessions worldwide supported that B. juncea is most likely a single origin in West Asia, 8,000-14,000 years ago, via natural interspecific hybridization. Subsequently, new crop types evolved through spontaneous gene mutations and introgressions along three independent routes of eastward expansion. Selective sweeps, genome-wide trait associations and tissue-specific RNA-sequencing analysis shed light on the domestication history of flowering time and seed weight, and on human selection for morphological diversification in this versatile species. Our data provide a comprehensive insight into the origin and domestication and a foundation for genomics-based breeding of B. juncea.

Complete genome sequencing and comparative genome analysis of the extremely halophilic archaea, Haloterrigena daqingensis

In the present study, we report the complete genome sequencing of Haloterrigena daqingensis species. The genome of H. daqingensis JX313^T consisted of a circular chromosome with three plasmids. The genome size and G+C content were estimated to be 3835796 bp and 61.7%, respectively. A total of 4158 genes were predicted with six rRNAs and 45 tRNAs. Metabolic pathway analysis suggests that H. daqingensis JX313^T codes for all the necessary genes responsible to sustain its life at saline environment. The pan-genome analysis suggests that the number of singleton-gene between H. daqingensis and other Haloterrigena species varied. The study not only helps us understand H. daqingensis strategy for dealing with high stress, but it also provides an overview of its genomic makeup.

Genome and transcriptome sequencing of a newly isolated 2,4-dinitrophenol-degrading strain Rhodococcus imtechensis XM24D

BACKGROUND

2,4-Dinitrophenol (2,4-DNP) is an important organic environmental pollutant that is highly toxic to all forms of living organisms. A gram-positive strain (designated XM24D) was isolated from 2,4-DNP-contaminated soil by an enrichment technique.

OBJECTIVE

The study was designed to analyze the ability of XM24D to degrade 2,4-DNP and its analogs and to reveal the degradation pathways of these aromatic compounds.

METHODS

The degradation ability of XM24D was tested by a growth experiment. 2,4-DNP and its analog degradation pathways were predicted by genome and comparative transcriptome sequencing.

RESULTS

Growth profiles showed that XM24D was able to utilize 2,4-DNP as the sole source of carbon, nitrogen and energy. Analogs of 2,4-DNP, including 4-nitrophenol (PNP) and 2-chloro-4-nitrophenol (2C4NP), can also be degraded by XM24D. Genome analysis showed that the XM24D genome contains two chromosomes with a combined size of 9.08 Mb and an average GC content of 67.07 %. Average nucleotide identity analysis indicated that Rhodococcus imtechensis RKJ300 is the most closely related strain to XM24D. Comparative transcriptome analysis revealed that the 2,4-DNP/PNP/2C4NP degradation pathway in XM24D is highly similar in sequence and organization to the 2,4-DNP degradation pathway in Rhodococcus opacus HL PM-1, the PNP degradation pathway in Rhodococcus opacus SAO101 and the 2C4NP degradation pathway in Rhodococcus imtechensis RKJ300. These results suggested that 2,4-DNP/PNP/2C4NP was degraded via the 2,4-dinitrocyclohexanone/4-nitrocatechol/hydroxyquinol pathway in XM24D.

CONCLUSIONS

Genomic and transcriptomic information on XM24D provides a valuable reference for further investigating the evolutionary characteristics of nitrophenol degradation pathways in microorganisms.

Complete Genome Sequence of Curtobacterium sp. Strain YC1, Isolated from the Surface of Nostoc flagelliforme Colonies in Yinchuan, Ningxia, China

The genus Curtobacterium belongs to the family Microbacteriaceae, within the phylum Actinobacteria. This genus includes a wide range of Gram-positive species associated with plants and soils. Here, we report the genome sequence of a new strain, Curtobacterium sp. strain YC1, which was isolated from the surface of Nostoc flagelliforme colonies. The genome of this strain contains one chromosome and one plasmid, and its size is 3.4 Mb.

Complete and Draft Genome Sequences of Aerobic Methanotrophs Isolated from a Riparian Wetland

Wetlands are important sources of methane emissions, and the impacts of these emissions can be mitigated by methanotrophic bacteria. The genomes of methanotrophs Methylomonas sp. strain LL1 and Methylosinus sp. strain H3A, as well as Methylocystis sp. strains H4A, H15, H62, and L43, were sequenced and are reported here.

Wetlands are important sources of methane emissions, and the impacts of these emissions can be mitigated by methanotrophic bacteria. The genomes of methanotrophs Methylomonas sp. strain LL1 and Methylosinus sp. strain H3A, as well as Methylocystis sp. strains H4A, H15, H62, and L43, were sequenced and are reported here.

Genome sequences of Tropheus moorii and Petrochromis trewavasae, two eco-morphologically divergent cichlid fishes endemic to Lake Tanganyika

With more than 1000 species, East African cichlid fishes represent the fastest and most species-rich vertebrate radiation known, providing an ideal model to tackle molecular mechanisms underlying recurrent adaptive diversification. We add high-quality genome reconstructions for two phylogenetic key species of a lineage that diverged about ~ 3-9 million years ago (mya), representing the earliest split of the so-called modern haplochromines that seeded additional radiations such as those in Lake Malawi and Victoria. Along with the annotated genomes we analysed discriminating genomic features of the study species, each representing an extreme trophic morphology, one being an algae browser and the other an algae grazer. The genomes of Tropheus moorii (TM) and Petrochromis trewavasae (PT) comprise 911 and 918 Mbp with 40,300 and 39,600 predicted genes, respectively. Our DNA sequence data are based on 5 and 6 individuals of TM and PT, and the transcriptomic sequences of one individual per species and sex, respectively. Concerning variation, on average we observed 1 variant per 220 bp (interspecific), and 1 variant per 2540 bp (PT vs PT)/1561 bp (TM vs TM) (intraspecific). GO enrichment analysis of gene regions affected by variants revealed several candidates which may influence phenotype modifications related to facial and jaw morphology, such as genes belonging to the Hedgehog pathway (SHH, SMO, WNT9A) and the BMP and GLI families.

Comparative genomics of Klebsiella michiganensis BD177 and related members of Klebsiella sp. reveal the symbiotic relationship with Bactrocera dorsalis

Background

Bactrocera dorsalis is a destructive polyphagous and highly invasive insect pest of tropical and subtropical species of fruit and vegetable crops. The sterile insect technique (SIT) has been used for decades to control insect pests of agricultural, veterinary, and human health importance. Irradiation of pupae in SIT can reduce the ecological fitness of the sterile insects. Our previous study has shown that a gut bacterial strain BD177 that could restore ecological fitness by promoting host food intake and metabolic activities.

Results

Using long-read sequence technologies, we assembled the complete genome of K. michiganensis BD177 strain. The complete genome of K. michiganensis BD177 comprises one circular chromosome and four plasmids with a GC content of 55.03%. The pan-genome analysis was performed on 119 genomes (strain BD177 genome and 118 out of 128 published Klebsiella sp. genomes since ten were discarded). The pan-genome includes a total of 49305 gene clusters, a small number of 858 core genes, and a high number of accessory (10566) genes. Pan-genome and average nucleotide identity (ANI) analysis showed that BD177 is more similar to the type strain K. michiganensis DSM2544, while away from the type strain K. oxytoca ATCC13182. Comparative genome analysis with 21 K. oxytoca and 12 K. michiganensis strains, identified 213 unique genes, several of them related to amino acid metabolism, metabolism of cofactors and vitamins, and xenobiotics biodegradation and metabolism in BD177 genome.

Conclusions

Phylogenomics analysis reclassified strain BD177 as a member of the species K. michiganensis. Comparative genome analysis suggested that K. michiganensis BD177 has the strain-specific ability to provide three essential amino acids (phenylalanine, tryptophan and methionine) and two vitamins B (folate and riboflavin) to B. dorsalis. The clear classification status of BD177 strain and identification of unique genetic characteristics may contribute to expanding our understanding of the symbiotic relationship of gut microbiota and B. dorsalis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-020-00945-0.

Complete chloroplast genomes of Zingiber montanum and Zingiber zerumbet: Genome structure, comparative and phylogenetic analyses

Zingiber montanum (Z. montanum) and Zingiber zerumbet (Z. zerumbet) are important medicinal and ornamental herbs in the genus Zingiber and family Zingiberaceae. Chloroplast-derived markers are useful for species identification and phylogenetic studies, but further development is warranted for these two Zingiber species. In this study, we report the complete chloroplast genomes of Z. montanum and Z. zerumbet, which had lengths of 164,464 bp and 163,589 bp, respectively. These genomes had typical quadripartite structures with a large single copy (LSC, 87,856-89,161 bp), a small single copy (SSC, 15,803-15,642 bp), and a pair of inverted repeats (IRa and IRb, 29,393-30,449 bp). We identified 111 unique genes in each chloroplast genome, including 79 protein-coding genes, 28 tRNAs and 4 rRNA genes. We analyzed the molecular structures, gene information, amino acid frequencies, codon usage patterns, RNA editing sites, simple sequence repeats (SSRs) and long repeats from the two chloroplast genomes. A comparison of the Z. montanum and Z. zerumbet chloroplast genomes detected 489 single-nucleotide polymorphisms (SNPs) and 172 insertions/deletions (indels). Thirteen highly divergent regions, including ycf1, rps19, rps18-rpl20, accD-psaI, psaC-ndhE, psbA-trnK-UUU, trnfM-CAU-rps14, trnE-UUC-trnT-UGU, ccsA-ndhD, psbC-trnS-UGA, start-psbA, petA-psbJ, and rbcL-accD, were identified and might be useful for future species identification and phylogeny in the genus Zingiber. Positive selection was observed for ATP synthase (atpA and atpB), RNA polymerase (rpoA), small subunit ribosomal protein (rps3) and other protein-coding genes (accD, clpP, ycf1, and ycf2) based on the Ka/Ks ratios. Additionally, chloroplast SNP-based phylogeny analyses found that Zingiber was a monophyletic sister branch to Kaempferia and that chloroplast SNPs could be used to identify Zingiber species. The genome resources in our study provide valuable information for the identification and phylogenetic analysis of the genus Zingiber and family Zingiberaceae.

Whole-genome, transcriptome, and methylome analyses provide insights into the evolution of platycoside biosynthesis in Platycodon grandiflorus, a medicinal plant

Triterpenoid saponins (TSs) are common plant defense phytochemicals with potential pharmaceutical properties. Platycodon grandiflorus (Campanulaceae) has been traditionally used to treat bronchitis and asthma in East Asia. The oleanane-type TSs, platycosides, are a major component of the P. grandiflorus root extract. Recent studies show that platycosides exhibit anti-inflammatory, antiobesity, anticancer, antiviral, and antiallergy properties. However, the evolutionary history of platycoside biosynthesis genes remains unknown. In this study, we sequenced the genome of P. grandiflorus and investigated the genes involved in platycoside biosynthesis. The draft genome of P. grandiflorus is 680.1 Mb long and contains 40,017 protein-coding genes. Genomic analysis revealed that the CYP716 family genes play a major role in platycoside oxidation. The CYP716 gene family of P. grandiflorus was much larger than that of other Asterid species. Orthologous gene annotation also revealed the expansion of β-amyrin synthases (bASs) in P. grandiflorus, which was confirmed by tissue-specific gene expression. In these expanded gene families, we identified key genes showing preferential expression in roots and association with platycoside biosynthesis. In addition, whole-genome bisulfite sequencing showed that CYP716 and bAS genes are hypomethylated in P. grandiflorus, suggesting that epigenetic modification of these two gene families affects platycoside biosynthesis. Thus whole-genome, transcriptome, and methylome data of P. grandiflorus provide novel insights into the regulation of platycoside biosynthesis by CYP716 and bAS gene families.

Decolorization and biodegradation of textile di-azo dye Acid Blue 113 by Pseudomonas stutzeri AK6

Textile industry is one of the anthropogenic activities that consume a large amount of water and pollute water bodies. It uses a massive amount of dyes, which is one of the main constituents of polluting textile effluent. In the present research, biodegradation of Acid Blue 113 dye, a commonly used textile di-azo dye, has been studied exploiting Pseudomonas stutzeri, strain AK6. The dye (300 ppm) was decolorized up to 86.2% within 96 h. The metabolites of Acid Blue 113 obtained after biodegradation were identified by various analytical techniques viz. HPLC (high-performance liquid chromatography) and GC-MS (gas chromatography-mass spectrometry). Genome analysis of isolate AK6 using IMG/M (Integrated Microbial Genomes and Microbiomes) system supported the role of azoreductase and laccase for the decolorization and degradation of azo dye. The ability of P. stutzeri AK6 to tolerate high amount of dye makes it a potential candidate for bioremediation and pre-processing to remove dyes from textile effluents.

Characterization of an efficient estrogen-degrading bacterium Stenotrophomonas maltophilia SJTH1 in saline-, alkaline-, heavy metal-contained environments or solid soil and identification of four 17β-estradiol-oxidizing dehydrogenases

The efficient bioremediation of estrogen contamination in complex environments is of great concern. Here the strain Stenotrophomonas maltophilia SJTH1 was found with great and stable estrogen-degradation efficiency even under stress environments. The strain could utilize 17β-estradiol (E2) as a carbon source and degrade 90% of 10 mg/L E2 in a week; estrone (E1) was the first degrading intermediate of E2. Notably, diverse pH conditions (3.0-11.0) and supplements of 4% salinity, 6.25 mg/L of heavy metal (Cd²⁺ or Cu²⁺), or 1 CMC of surfactant (Tween 80/ Triton X-100) had little effect on its cell growth and estrogen degradation. The addition of low concentrations of copper and Tween 80 even promoted its E2 degradation. Bioaugmentation of strain SJTH1 into solid clay soil achieved over 80% removal of E2 contamination (10 mg/kg) within two weeks. Further, the whole genome sequence of S. maltophilia SJTH1 was obtained, and a series of potential genes participating in stress-tolerance and estrogen-degradation were predicted. Four dehydrogenases similar to 17β-hydroxysteroid dehydrogenases (17β-HSDs) were found to be induced by E2, and the four heterogenous-expressed enzymes could oxidize E2 into E1 efficiently. This work could promote bioremediation appliance potential with microorganisms and biodegradation mechanism study of estrogens in complex real environments.

Complete Chloroplast Genomes of Three Medicinal Alpinia Species: Genome Organization, Comparative Analyses and Phylogenetic Relationships in Family Zingiberaceae

Alpinia katsumadai (A. katsumadai), Alpinia oxyphylla (A. oxyphylla) and Alpinia pumila (A. pumila), which belong to the family Zingiberaceae, exhibit multiple medicinal properties. The chloroplast genome of a non-model plant provides valuable information for species identification and phylogenetic analysis. Here, we sequenced three complete chloroplast genomes of A. katsumadai, A. oxyphylla sampled from Guangdong and A. pumila, and analyzed the published chloroplast genomes of Alpinia zerumbet (A. zerumbet) and A. oxyphylla sampled from Hainan to retrieve useful chloroplast molecular resources for Alpinia. The five Alpinia chloroplast genomes possessed typical quadripartite structures comprising of a large single copy (LSC, 87,248-87,667 bp), a small single copy (SSC, 15,306-18,295 bp) and a pair of inverted repeats (IR, 26,917-29,707 bp). They had similar gene contents, gene orders and GC contents, but were slightly different in the numbers of small sequence repeats (SSRs) and long repeats. Interestingly, fifteen highly divergent regions (rpl36, ycf1, rps15, rpl22, infA, psbT-psbN, accD-psaI, petD-rpoA, psaC-ndhE, ccsA-ndhD, ndhF-rpl32, rps11-rpl36, infA-rps8, psbC-psbZ, and rpl32-ccsA), which could be suitable for species identification and phylogenetic studies, were detected in the Alpinia chloroplast genomes. Comparative analyses among the five chloroplast genomes indicated that 1891 mutational events, including 304 single nucleotide polymorphisms (SNPs) and 118 insertion/deletions (indels) between A. pumila and A. katsumadai, 367 SNPs and 122 indels between A. pumila and A. oxyphylla sampled from Guangdong, 331 SNPs and 115 indels between A. pumila and A. zerumbet, 371 SNPs and 120 indels between A. pumila and A. oxyphylla sampled from Hainan, and 20 SNPs and 23 indels between the two accessions of A. oxyphylla, were accurately located. Additionally, phylogenetic relationships based on SNP matrix among 28 whole chloroplast genomes showed that Alpinia was a sister branch to Amomum in the family Zingiberaceae, and that the five Alpinia accessions were divided into three groups, one including A. pumila, another including A. zerumbet and A. katsumadai, and the other including two accessions of A. oxyphylla. In conclusion, the complete chloroplast genomes of the three medicinal Alpinia species in this study provided valuable genomic resources for further phylogeny and species identification in the family Zingiberaceae.

Complete Genome Sequencing of Mycoplasma synoviae Strain HN01, Isolated from Chicken in Henan Province, China

Here, we report the complete genome sequence of Mycoplasma synoviae HN01, a virulent epidemic strain isolated from a sick chicken with synovitis in Henan Province, China. HN01 is the Asian source of an M. synoviae strain that is completely sequenced, genome annotated, and published with relevant data.

Here, we report the complete genome sequence of Mycoplasma synoviae HN01, a virulent epidemic strain isolated from a sick chicken with synovitis in Henan Province, China. HN01 is the Asian source of an M. synoviae strain that is completely sequenced, genome annotated, and published with relevant data.

Characterization of an 17β-estradiol-degrading bacterium Stenotrophomonas maltophilia SJTL3 tolerant to adverse environmental factors

Bioremediation of environmental estrogens requires microorganisms with stable degradation efficiency and great stress tolerance in complex environments. In this work, Stenotrophomonas maltophilia SJTL3 isolated from wastewater was found to be able to degrade over 90% of 10 μg/mL 17β-estradiol (E2) in a week and the degradation dynamic was fitted by the first-order kinetic equations. Estrone was the first and major intermediate of E2 biodegradation. Strain SJTL3 exhibited strong tolerance to several adverse conditions like extreme pH (3.0-11.0), high osmolality (2%), co-existing heavy metals (6.25 μg/mL of Cu²⁺) and surfactants (5 CMC of Tween 80), and retained normal cell vitality and stable E2-degradaing efficiency. In solid soil, strain SJTL3 could remove nearly 100% of 1 μg/mL of E2 after the bacteria inoculation and 8-day culture. As to the contamination of 10 μg/mL E2 in soil, the biodegradation efficiency was about 90%. The further obtainment of the whole genome of strain SJTL3 and genome analysis revealed that this strain contained not only the potential genes responsible for estrogen degradation, but also the genes encoding proteins involved in stress tolerance. This work could promote the estrogen-biodegrading mechanism study and provide insights into the bioremediation application.

Comparative genomic analysis of six Glossina genomes, vectors of African trypanosomes

BACKGROUND

Tsetse flies (Glossina sp.) are the vectors of human and animal trypanosomiasis throughout sub-Saharan Africa. Tsetse flies are distinguished from other Diptera by unique adaptations, including lactation and the birthing of live young (obligate viviparity), a vertebrate blood-specific diet by both sexes, and obligate bacterial symbiosis. This work describes the comparative analysis of six Glossina genomes representing three sub-genera: Morsitans (G. morsitans morsitans, G. pallidipes, G. austeni), Palpalis (G. palpalis, G. fuscipes), and Fusca (G. brevipalpis) which represent different habitats, host preferences, and vectorial capacity.

RESULTS

Genomic analyses validate established evolutionary relationships and sub-genera. Syntenic analysis of Glossina relative to Drosophila melanogaster shows reduced structural conservation across the sex-linked X chromosome. Sex-linked scaffolds show increased rates of female-specific gene expression and lower evolutionary rates relative to autosome associated genes. Tsetse-specific genes are enriched in protease, odorant-binding, and helicase activities. Lactation-associated genes are conserved across all Glossina species while male seminal proteins are rapidly evolving. Olfactory and gustatory genes are reduced across the genus relative to other insects. Vision-associated Rhodopsin genes show conservation of motion detection/tracking functions and variance in the Rhodopsin detecting colors in the blue wavelength ranges.

CONCLUSIONS

Expanded genomic discoveries reveal the genetics underlying Glossina biology and provide a rich body of knowledge for basic science and disease control. They also provide insight into the evolutionary biology underlying novel adaptations and are relevant to applied aspects of vector control such as trap design and discovery of novel pest and disease control strategies.

Comparative genomic analysis of six Glossina genomes, vectors of African trypanosomes

BACKGROUND

Tsetse flies (Glossina sp.) are the vectors of human and animal trypanosomiasis throughout sub-Saharan Africa. Tsetse flies are distinguished from other Diptera by unique adaptations, including lactation and the birthing of live young (obligate viviparity), a vertebrate blood-specific diet by both sexes, and obligate bacterial symbiosis. This work describes the comparative analysis of six Glossina genomes representing three sub-genera: Morsitans (G. morsitans morsitans, G. pallidipes, G. austeni), Palpalis (G. palpalis, G. fuscipes), and Fusca (G. brevipalpis) which represent different habitats, host preferences, and vectorial capacity.

RESULTS

Genomic analyses validate established evolutionary relationships and sub-genera. Syntenic analysis of Glossina relative to Drosophila melanogaster shows reduced structural conservation across the sex-linked X chromosome. Sex-linked scaffolds show increased rates of female-specific gene expression and lower evolutionary rates relative to autosome associated genes. Tsetse-specific genes are enriched in protease, odorant-binding, and helicase activities. Lactation-associated genes are conserved across all Glossina species while male seminal proteins are rapidly evolving. Olfactory and gustatory genes are reduced across the genus relative to other insects. Vision-associated Rhodopsin genes show conservation of motion detection/tracking functions and variance in the Rhodopsin detecting colors in the blue wavelength ranges.

CONCLUSIONS

Expanded genomic discoveries reveal the genetics underlying Glossina biology and provide a rich body of knowledge for basic science and disease control. They also provide insight into the evolutionary biology underlying novel adaptations and are relevant to applied aspects of vector control such as trap design and discovery of novel pest and disease control strategies.

Metabolic engineering of Corynebacterium glutamicum S9114 based on whole-genome sequencing for efficient N-acetylglucosamine synthesis

Glucosamine (GlcN) and its acetylated derivative N-acetylglucosamine (GlcNAc) are widely used in the pharmaceutical industries. Here, we attempted to achieve efficient production of GlcNAc via genomic engineering of Corynebacterium glutamicum. Specifically, we ligated the GNA1 gene, which converts GlcN-6-phosphate to GlcNAc-6-phosphate by transferring the acetyl group in Acetyl-CoA to the amino group of GlcN-6-phosphate, into the plasmid pJYW4 and then transformed this recombinant vector into the C. glutamicum ATCC 13032, ATCC 13869, ATCC 14067, and S9114 strains, and we assessed the GlcNAc titers at 0.5 g/L, 1.2 g/L, 0.8 g/L, and 3.1 g/L from each strain, respectively. This suggested that there were likely to be significant differences among the key genes in the glutamate and GlcNAc synthesis pathways of these C. glutamicum strains. Therefore, we performed whole genome sequencing of the S9114 strain, which has not been previously published, and found that there are many differences among the genes in the glutamate and GlcNAc synthesis pathways among the four strains tested. Next, nagA (encoding GlcNAc-6-phosphate deacetylase) and gamA (encoding GlcN-6-phosphate deaminase) were deleted in C. glutamicum S9114 to block the catabolism of intracellular GlcNAc, leading to a 54.8% increase in GlcNAc production (from 3.1 to 4.8 g/L) when grown in a shaker flask. In addition, lactate synthesis was blocked by knockout of ldh (encoding lactate dehydrogenase); thus, further increasing the GlcNAc titer to 5.4 g/L. Finally, we added a key gene of the GlcN synthetic pathway, glmS, from different sources into the expression vector pJYW-4-ceN, and the resulting recombinant strain CGGN2-GNA1-CgglmS produced the GlcNAc titer of 6.9 g/L. This is the first report concerning the metabolic engineering of C. glutamicum, and the results of this study provide a good starting point for further metabolic engineering to achieve industrial-scale production of GlcNAc.

The putative functions of lysogeny in mediating the survivorship of Escherichia coli in seawater and marine sediment

Escherichia coli colonizes various body parts of animal hosts as a commensal and a pathogen. It can also persist in the external environment in the absence of fecal pollution. It remains unclear how this species has evolved to adapt to such contrasting habitats. Lysogeny plays pivotal roles in the diversification of the phenotypic and ecologic characters of E. coli as a symbiont. We hypothesized that lysogeny could also confer fitness to survival in the external environment. To test this hypothesis, we used the induced phages of an E. coli strain originating from marine sediment to infect a fecal E. coli strain to obtain an isogenic lysogen of the latter. The three strains were tested for survivorship in microcosms of seawater, marine sediment and sediment interstitial water as well as swimming motility, glycogen accumulation, biofilm formation, substrate utilization and stress resistance. The results indicate that lysogenic infection led to tractable changes in many of the ecophysiological attributes tested. Particularly, the lysogen had better survivorship in the microcosms and had a substrate utilization profile resembling the sediment strain more than the wild type fecal strain. Our findings provide new insights into the understanding of how E. coli survives in the natural environment.

Complete Chloroplast Genome Sequences of Kaempferia Galanga and Kaempferia Elegans: Molecular Structures and Comparative Analysis

Kaempferia galanga and Kaempferia elegans, which belong to the genus Kaempferia family Zingiberaceae, are used as valuable herbal medicine and ornamental plants, respectively. The chloroplast genomes have been used for molecular markers, species identification and phylogenetic studies. In this study, the complete chloroplast genome sequences of K. galanga and K. elegans are reported. Results show that the complete chloroplast genome of K. galanga is 163,811 bp long, having a quadripartite structure with large single copy (LSC) of 88,405 bp and a small single copy (SSC) of 15,812 bp separated by inverted repeats (IRs) of 29,797 bp. Similarly, the complete chloroplast genome of K. elegans is 163,555 bp long, having a quadripartite structure in which IRs of 29,773 bp length separates 88,020 bp of LSC and 15,989 bp of SSC. A total of 111 genes in K. galanga and 113 genes in K. elegans comprised 79 protein-coding genes and 4 ribosomal RNA (rRNA) genes, as well as 28 and 30 transfer RNA (tRNA) genes in K. galanga and K. elegans, respectively. The gene order, GC content and orientation of the two Kaempferia chloroplast genomes exhibited high similarity. The location and distribution of simple sequence repeats (SSRs) and long repeat sequences were determined. Eight highly variable regions between the two Kaempferia species were identified and 643 mutation events, including 536 single-nucleotide polymorphisms (SNPs) and 107 insertion/deletions (indels), were accurately located. Sequence divergences of the whole chloroplast genomes were calculated among related Zingiberaceae species. The phylogenetic analysis based on SNPs among eleven species strongly supported that K. galanga and K. elegans formed a cluster within Zingiberaceae. This study identified the unique characteristics of the entire K. galanga and K. elegans chloroplast genomes that contribute to our understanding of the chloroplast DNA evolution within Zingiberaceae species. It provides valuable information for phylogenetic analysis and species identification within genus Kaempferia.

Complete Genome Sequence of Saccharospirillum mangrovi HK-33T Sheds Light on the Ecological Role of a Bacterium in Mangrove Sediment Environment

We present the genome sequence of Saccharospirillum mangrovi HK-33^T, isolated from a mangrove sediment sample in Haikou, China. The complete genome of S. mangrovi HK-33^T consisted of a single-circular chromosome with the size of 3,686,911 bp as well as an average G + C content of 57.37%, and contained 3,383 protein-coding genes, 4 operons of 16S-23S-5S rRNA genes, and 52 tRNA genes. Genomic annotation indicated that the genome of S. mangrovi HK-33^T had many genes related to oligosaccharide and polysaccharide degradation and utilization of polyhydroxyalkanoate. For nitrogen cycle, genes encoding nitrate and nitrite reductase, glutamate dehydrogenase, glutamate synthase, and glutamine synthetase could be found. For phosphorus cycle, genes related to polyphosphate kinases (ppk1 and ppk2), the high-affinity phosphate-specific transport (Pst) system, and the low-affinity inorganic phosphate transporter (pitA) were predicted. For sulfur cycle, cysteine synthase and type III acyl coenzyme A transferase (dddD) coding genes were searched out. This study provides evidence about carbon, nitrogen, phosphorus, and sulfur metabolic patterns of S. mangrovi HK-33^T and broadens our understandings about ecological roles of this bacterium in the mangrove sediment environment.

Hierarchical Clustering of DNA k-mer Counts in RNAseq Fastq Files Identifies Sample Heterogeneities

We apply hierarchical clustering (HC) of DNA k-mer counts on multiple Fastq files. The tree structures produced by HC may reflect experimental groups and thereby indicate experimental effects, but clustering of preparation groups indicates the presence of batch effects. Hence, HC of DNA k-mer counts may serve as a diagnostic device. In order to provide a simple applicable tool we implemented sequential analysis of Fastq reads with low memory usage in an R package (seqTools) available on Bioconductor. The approach is validated by analysis of Fastq file batches containing RNAseq data. Analysis of three Fastq batches downloaded from ArrayExpress indicated experimental effects. Analysis of RNAseq data from two cell types (dermal fibroblasts and Jurkat cells) sequenced in our facility indicate presence of batch effects. The observed batch effects were also present in reads mapped to the human genome and also in reads filtered for high quality (Phred > 30). We propose, that hierarchical clustering of DNA k-mer counts provides an unspecific diagnostic tool for RNAseq experiments. Further exploration is required once samples are identified as outliers in HC derived trees.

The pathogenic mechanisms of Tilletia horrida as revealed by comparative and functional genomics

Tilletia horrida is a soil-borne, mononucleate basidiomycete fungus with a biotrophic lifestyle that causes rice kernel smut, a disease that is distributed throughout hybrid rice growing areas worldwide. Here we report on the high-quality genome sequence of T. horrida; it is composed of 23.2 Mb that encode 7,729 predicted genes and 6,973 genes supported by RNA-seq. The genome contains few repetitive elements that account for 8.45% of the total. Evolutionarily, T. horrida lies close to the Ustilago fungi, suggesting grass species as potential hosts, but co-linearity was not observed between T. horrida and the barley smut Ustilago hordei. Genes and functions relevant to pathogenicity were presumed. T. horrida possesses a smaller set of carbohydrate-active enzymes and secondary metabolites, which probably reflect the specific characteristics of its infection and biotrophic lifestyle. Genes that encode secreted proteins and enzymes of secondary metabolism, and genes that are represented in the pathogen-host interaction gene database genes, are highly expressed during early infection; this is consistent with their potential roles in pathogenicity. Furthermore, among the 131 candidate pathogen effectors identified according to their expression patterns and functionality, we validated two that trigger leaf cell death in Nicotiana benthamiana. In summary, we have revealed new molecular mechanisms involved in the evolution, biotrophy, and pathogenesis of T. horrida.

Identification and genome analysis of Deinococcus actinosclerus SJTR1, a novel 17β-estradiol degradation bacterium

Biodegradation with microorganisms is considered as an efficient strategy to remove the environmental pollutants. In this work, Deinococcus actinosclerus SJTR1 isolated from the wastewater was confirmed with great degradation capability to 17β-estradiol, one typical estrogen chemical. It could degrade nearly 90% of 17β-estradiol (10 mg/L) in 5 days and transform it into estrone; its degradation kinetics fitted for the first-order kinetic equation. The whole genome sequence of D. actinosclerus SJTR1 was obtained and annotated, containing one chromosome (3,315,586 bp) and four plasmids (ranging from 17,267 bp to 460,244 bp). A total of 3913 CDSs and 73 RNA genes (including 12 rRNA genes, 50 tRNA genes, and 11 ncRNA genes) were identified in its whole genome sequence. On this basis, a series of potential genes involved in steroid metabolism and stress responses of D. actinosclerus SJTR1 were predicted. It is the first report of Deinococcus strain with the degradation capability to estrogens. This work could enrich the genome sources of the estrogen-degrading strains and promote the degradation mechanism study of 17β-estradiol in bacteria.

Genome Sequencing and Assembly by Long Reads in Plants

Plant genomes generated by Sanger and Next Generation Sequencing (NGS) have provided insight into species diversity and evolution. However, Sanger sequencing is limited in its applications due to high cost, labor intensity, and low throughput, while NGS reads are too short to resolve abundant repeats and polyploidy, leading to incomplete or ambiguous assemblies. The advent and improvement of long-read sequencing by Third Generation Sequencing (TGS) methods such as PacBio and Nanopore have shown promise in producing high-quality assemblies for complex genomes. Here, we review the development of sequencing, introducing the application as well as considerations of experimental design in TGS of plant genomes. We also introduce recent revolutionary scaffolding technologies including BioNano, Hi-C, and 10× Genomics. We expect that the informative guidance for genome sequencing and assembly by long reads will benefit the initiation of scientists' projects.