基于叶绿体基因组探讨桃金娘目及其近缘类群的系统发育关系

该研究基于叶绿体基因组数据,对桃金娘目(6科44属97种)及其近缘类群(牻牛儿苗目2科5属25种)的系统发育关系进行了分析.结果表明:(1)桃金娘目基因组大小为152～171 kb,包括的蛋白质编码基因数目为74～90个;牻牛儿苗目基因组大小为116～242 kb,包括的蛋白质编码基因数目为75～132个.(2)对比叶...     

A chromosome-level genome assembly of Ostrea denselamellosa provides initial insights into its evolution

The oyster Ostrea denselamellosa is a live-bearing species with a sharp decline in the natural population. Despite recent breakthroughs in long-read sequencing, high quality genomic data are very limited in O. denselamellosa. Here, we carried out the first whole genome sequencing at the chromosome-level in O. denselamellosa. Our studies yielded a 636 Mb assembly with scaffold N50 around 71.80 Mb. 608.3 Mb (95.6% of the assembly) were anchored to 10 chromosomes. A total of 26,412 protein-coding genes were predicted, of which 22,636 (85.7%) were functionally annotated. By comparative genomics, we found that long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) made up a larger proportion in O. denselamellosa genome than in other oysters'. Moreover, gene family analysis showed some initial insight into its evolution. This high-quality genome of O. denselamellosa provides a valuable genomic resource for studies of evolution, adaption and conservation in oysters.     

The chromosome-scale assembly of endive (Cichorium endivia) genome provides insights into the sesquiterpenoid biosynthesis

Endive (Cichorium endivia L.) is a leafy vegetable in the Asteraceae family. Sesquiterpene lactones (STLs) in endive leaves bring a bitter taste that varies between varieties. Despite their importance in breeding varieties with unique flavours, sesquiterpenoid biosynthesis pathways in endive are poorly understood. We assembled a chromosome-scale endive genome of 641 Mb with a contig N50 of 5.16 Mb and annotated 46,711 protein-coding genes. Several gene families, especially terpene synthases (TPS) genes, expanded significantly in the C. endivia genome. STLs biosynthesis-related genes and TPS genes in more bitter varieties have shown a higher level of expression, which could be attributed to genomic variations. Our results penetrate the origin and diversity of bitter taste and facilitate the molecular breeding of endive varieties with unique bitter tastes. The high-quality endive assembly would provide a reference genome for studying the evolution and diversity of Asteraceae.     

Chromosome‐level genome assembly of a cyprinid fish Onychostoma macrolepis by integration of nanopore sequencing,Bionano and Hi‐C technology

Onychostoma macrolepis is an emerging commercial cyprinid fish species. It is a model system for studies of sexual dimorphism and genome evolution. Here, we report the chromosome‐level assembly of the O.macrolepis genome obtained from the integration of nanopore long‐read sequencing with physical maps produced using Bionano and Hi‐C technology. A total of 87.9 Gb of nanopore sequence provided approximately 100‐fold coverage of the genome. The preliminary genome assembly was 883.2 Mb in size with a contig N50 size of 11.2 Mb. The 969 corrected contigs obtained from Bionano optical mapping were assembled into 853 scaffolds and produced an assembly of 886.5 Mb with a scaffold N50 of 16.5 Mb. Finally, using the Hi‐C data, 881.3 Mb (99.4% of genome) in 526 scaffolds were anchored and oriented in 25 chromosomes ranging in size from 25.27 to 56.49 Mb. In total, 24,770 protein‐coding genes were predicted in the genome, and ~96.85% of the genes were functionally annotated. The annotated assembly contains 93.3% complete genes from the BUSCO reference set. In addition, we identified 409 Mb (46.23% of the genome) of repetitive sequence, and 11,213 non‐coding RNAs, in the genome. Evolutionary analysis revealed that O. macrolepis diverged from common carp approximately 24.25 million years ago. The chromosomes of O. macrolepis showed an unambiguous correspondence to the chromosomes of zebrafish. The high‐quality genome assembled in this work provides a valuable genomic resource for further biological and evolutionary studies of O. macrolepis.     

Chromosome-level genome assembly of Scapharca kagoshimensis reveals the expanded molecular basis of heme biosynthesis in ark shells

Ark shells are commercially important clam species that inhabit in muddy sediments of shallow coasts in East Asia. For a long time, the lack of genome resources has hindered scientific research of ark shells. Here, we report a high-quality chromosome-level genome assembly of Scapharca kagoshimensis, with an aim to unravel the molecular basis of heme biosynthesis, and develop genomic resources for genetic breeding and population genetics in ark shells. Nineteen scaffolds corresponding to 19 chromosomes were constructed from 938 contigs (contig N50 = 2.01 Mb) to produce a final high-quality assembly with a total length of 1.11 Gb and scaffold N50 around 60.64 Mb. The genome assembly represents 93.4% completeness via matching 303 eukaryota core conserved genes. A total of 24,908 protein-coding genes were predicted and 24,551 genes (98.56%) of which were functionally annotated. The enrichment analyses suggested that genes in heme biosynthesis pathways were expanded and positive selection of the haemoglobin genes was also found in the genome of S. kagoshimensis, which gives important insights into the molecular mechanisms and evolution of the heme biosynthesis in mollusca. The valuable genome assembly of S. kagoshimensis would provide a solid foundation for investigating the molecular mechanisms that underlie the diverse biological functions and evolutionary adaptations of S. kagoshimensis.     

Genomic architecture of three newly isolated unclassified Butyrivibrio species elucidate their potential role in the rumen ecosystem

One cellulose-degrading strain CB08 and two xylan-degrading strains XB500–5 and X503 were isolated from buffalo rumen. All the strains were designated as putative novel species of Butyrivibrio based on phylogeny, phylogenomy, digital DNA-DNA hybridization, and average nucleotide identity with their closest type strains. The draft genome length of CB08 was ~3.54 Mb, while X503 and XB500-5 genome sizes were ~3.24 Mb and ~3.27 Mb, respectively. Only 68.28% of total orthologous clusters were shared among three genomes, and 40–44% of genes were identified as hypothetical proteins. The presence of genes encoding diverse carbohydrate-active enzymes (CAZymes) exhibited the lignocellulolytic potential of these strains. Further, the genome annotations revealed the metabolic pathways for monosaccharide fermentation to acetate, butyrate, lactate, ethanol, and hydrogen. The presence of genes for chemotaxis, antibiotic resistance, antimicrobial activity, synthesis of vitamins, and essential fatty acid suggested the versatile metabolic nature of these Butyrivibrio strains in the rumen environment.     

Genome assembly of Scorias spongiosa and comparative genomics provide insights into ecological adaptation of honeydew-dependent sooty mould fungi

Sooty moulds are fungi of economic importance and with unique lifestyle mainly growing on insect honeydew. However, the lack of genomic data hinders investigation of genetic mechanisms underlying their ecological adaptation. With long-read sequencing technology, we generated the genome of Scorias spongiosa, an extraordinary sooty mould fungus associated with honeydew of colony aphids and producing large fruiting bodies. A 24.21 Mb high-quality genome assembly with a N50 length of 3.37 Mb was obtained. The genome contained 7758 protein coding genes, 97.13% of which were homologous to known genes, and approximately 0.29 Mb repeat sequences. Comparative genomics showed S. spongiosa lost relatively more gene families and contained fewer species-specific genes and gene families, with many CAZyme families and sugar transporters reduced or absent. This study not only promotes understanding of the ecological adaptation of sooty moulds, but also provides valuable genomic data resource for future comparative genomic and genetic studies.     

The Chimonanthus salicifolius genome provides insight into magnoliid evolution and flavonoid biosynthesis

Chimonanthus salicifolius, a member of the Calycanthaceae of magnoliids, is one of the most famous medicinal plants in Eastern China. Here, we report a chromosome‐level genome assembly of C. salicifolius, comprising 820.1 Mb of genomic sequence with a contig N50 of 2.3 Mb and containing 36 651 annotated protein‐coding genes. Phylogenetic analyses revealed that magnoliids were sister to the eudicots. Two rounds of ancient whole‐genome duplication were inferred in the C. salicifolious genome. One is shared by Calycanthaceae after its divergence with Lauraceae, and the other is in the ancestry of Magnoliales and Laurales. Notably, long genes with > 20 kb in length were much more prevalent in the magnoliid genomes compared with other angiosperms, which could be caused by the length expansion of introns inserted by transposon elements. Homologous genes within the flavonoid pathway for C. salicifolius were identified, and correlation of the gene expression and the contents of flavonoid metabolites revealed potential critical genes involved in flavonoids biosynthesis. This study not only provides an additional whole‐genome sequence from the magnoliids, but also opens the door to functional genomic research and molecular breeding of C. salicifolius.     

The genome of the marine medaka Oryzias melastigma

Marine medaka (Oryzias melastigma) is considered to be a useful fish model for marine and estuarine ecotoxicology studies and has good potential for field‐based population genomics because of its geographical distribution in Asian estuarine and coastal areas. In this study, we present the first whole‐genome draft of O. melastigma. The genome assembly consists of 8,602 scaffolds (N50 = 23.737 Mb) and a total genome length of 779.4 Mb. A total of 23,528 genes were predicted, and 12,670 gene families shared with three teleost species (Japanese medaka, mangrove killifish and zebrafish) were identified. Genome analyses revealed that the O. melastigma genome is highly heterozygous and contains a large number of repeat sequences. This assembly represents a useful genomic resource for fish scientists.     

A high‐quality chromosome‐level genome assembly of a generalist herbivore,Trichoplusia ni

The cabbage looper, Trichoplusia ni, is a globally distributed highly polyphagous herbivore and an important agricultural pest. T. ni has evolved resistance to various chemical insecticides, and is one of the only two insect species that have evolved resistance to the biopesticide Bacillus thuringiensis (Bt) in agricultural systems and has been selected for resistance to baculovirus infections. We report a 333‐Mb high‐quality T. ni genome assembly, which has N50 lengths of scaffolds and contigs of 4.6 Mb and 140 Kb, respectively, and contains 14,384 protein‐coding genes. High‐density genetic maps were constructed to anchor 305 Mb (91.7%) of the assembly to 31 chromosomes. Comparative genomic analysis of T. ni with Bombyx mori showed enrichment of tandemly duplicated genes in T. ni in families involved in detoxification and digestion, consistent with the broad host range of T. ni. High levels of genome synteny were found between T. ni and other sequenced lepidopterans. However, genome synteny analysis of T. ni and the T. ni derived cell line High Five (Hi5) indicated extensive genome rearrangements in the cell line. These results provided the first genomic evidence revealing the high instability of chromosomes in lepidopteran cell lines known from karyotypic observations. The high‐quality T. ni genome sequence provides a valuable resource for research in a broad range of areas including fundamental insect biology, insect‐plant interactions and co‐evolution, mechanisms and evolution of insect resistance to chemical and biological pesticides, and technology development for insect pest management.     

The genetic adaptations of Toxoptera aurantii facilitated its rapid multiple plant hosts dispersal and invasion

Toxoptera aurantii Boyer de Fonscolombe (Hemiptera: Aphididae) can attack many plant hosts, including tea (Camellia sinensis L.), citrus (Citrus spp.), lychee (Litchi chinensis Sonn.), banana (Musa spp.), and pineapple (Ananas comasus L.) among others. It is a widely distributed hexapod and one of the most destructive pests in tea plantations, causing enormous economic losses in tea production each year. A high-quality reference genome is important to study the phylogenetics and evolution of T. aurantii because its genome is highly heterozygous and repetitive. We obtained a de novo genome assembly of T. aurantii at the chromosome level using a combination of long Nanopore reads from sequencing with high-throughput chromosome conformation capture technology. When finally assembled, the genome was 318.95 Mb on four chromosomes with a 15.19 Mb scaffold N50. A total of 12,162 genes encoded proteins, while there were 22.01% repetitive sequences that totaled 67.73 Mb. Phylogenetic analyses revealed that T. aurantii and Aphis gossypii parted ways approximately 7.6 million years ago (Mya). We used a combination of long-read single-molecule sequencing with Hi-C–based chromatin interaction maps that resulted in a reference chromosomal level reference genome of T. aurantii that was high quality. Our results will enable the exploration of the genetics behind the special biological features of T. aurantii and also provide a source of data that should be useful to compare the compare genome among the Hemiptera.     

The chicken (Gallus gallus) Z chromosome contains at least three nonlinear evolutionary strata

Birds have female heterogamety with Z and W sex chromosomes. These evolved from different autosomal precursor chromosomes than the mammalian X and Y. However, previous work has suggested that the pattern and process of sex chromosome evolution show many similarities across distantly related organisms. Here we show that stepwise restriction of recombination between the protosex chromosomes of birds has resulted in regions of the chicken Z chromosome showing discrete levels of divergence from W homologs (gametologs). The 12 genes analyzed fall into three levels of estimated divergence values, with the most recent divergence (d_S = 0.18–0.21) displayed by 6 genes in a region on the Z chromosome corresponding to the interval 1–11 Mb of the assembled genome sequence. Another 4 genes show intermediate divergence (d_S = 0.27–0.38) and are located in the interval 16–53 Mb. Two genes (at positions 42 and 50 Mb) with higher d_S values are located proximal to the most distal of the 4 genes with intermediate divergence, suggesting an inversion event. The distribution of genes and their divergence indicate at least three evolutionary strata, with estimated times for cessation of recombination between Z and W of 132–150 (stratum 1), 71–99 (stratum 2), and 47–57 (stratum 3) million years ago. An inversion event, or some other form of intrachromosomal rearrangement, subsequent to the formation of strata 1 and 2 has scrambled the gene order to give rise to the nonlinear arrangement of evolutionary strata currently seen on the chicken Z chromosome. These observations suggest that the progressive restriction of recombination is an integral feature of sex chromosome evolution and occurs also in systems of female heterogamety.     

First reported chloroplast genome sequence of <Emphasis Type="Italic">Punica granatum</Emphasis> (cultivar Helow) from Jabal Al-Akhdar,Oman: phylogenetic comparative assortment with <Emphasis Type="Italic">Lagerstroemia</Emphasis>

Pomegranate (Punica granatum L.) is one of the oldest known edible fruits. It has grown in popularity and is a profitable fruit crop due to its attractive features including a bright red appearance and its biological activities. Scientific exploration of the genetics and evolution of these beneficial traits has been hampered by limited genomic information. In this study, we sequenced the complete chloroplast (cp) genome of the native P. granatum (cultivar Helow) cultivated in the mountains of Jabal Al-Akhdar, Oman. The results revealed a P. granatum cp genome length of 158,630 bp, characterized by a relatively conserved structure containing 2 inverted repeat regions of 25,466 bp, an 18,686 bp small single copy regions, and an 89,015 bp large single copy region. The 86 protein-coding genes included 37 transfer RNA genes and 8 ribosomal RNA genes. Comparison of the P. granatum whole cp genome with seven Lagerstroemia species revealed an overall high degree of sequence similarity with divergence among intergenic spacers. The location, distribution, and divergence of repeat sequences and shared genes of the Punica and Lagerstroemia species were highly similar. Analyses of nucleotide substitution, insertion/deletions, and highly variable regions in these cp genomes identified potential plastid markers for taxonomic and phylogenetic studies in Myrtales. A phylogenetic study of the cp genomes and 76 shared coding regions generated similar cladograms. The complete cp genome of P. granatum will aid in taxonomical studies of the family Lythraceae.     

Genome analysis of phage JS98 defines a fourth major subgroup of T4-like phages in Escherichia coli

Numerous T4-like Escherichia coli phages were isolated from human stool and environmental wastewater samples in Bangladesh and Switzerland. The sequences of the major head gene (g23) revealed that these coliphages could be placed into four subgroups, represented by the phages T4, RB69, RB49, and JS98. Thus, JS98 defines a new major subgroup of E. coli T4-like phages. We conducted an analysis of the 169-kb JS98 genome sequence. Overall, 198 of the 266 JS98 open reading frames (ORFs) shared amino acid sequence identity with the reference T4 phage, 41 shared identity with other T4-like phages, and 27 ORFs lacked any database matches. Genes on the plus strand encoded virion proteins, which showed moderate to high sequence identity with T4 proteins. The right genome half of JS98 showed a higher degree of sequence conservation with T4 and RB69, even for the nonstructural genes, than did the left genome half, containing exclusively nonstructural genes. Most of the JS98-specific genes were found in the left genome half. Two came as a hypervariability cluster, but most represented isolated genes, suggesting that they were acquired separately in multiple acquisition events. No evidence for DNA exchange between JS98 phage and the E. coli host genome or coliphages other than T4 was observed. No undesired genes which could compromise its medical use were detected in the JS98 genome sequence.