Electrophoretic karyotyping of the lignin-degrading basidiomycete Phanerochaete chrysosporium

Electrophoretic karyotyping of the two most widely studied strains of Phanerochaete chrysosporium, BKMF-1767 and ME-446, has been determined using transverse alternating field etectrophoresis. The genomic DNA of BKMF-1767 was resolved into 10 chromosomes ranging in size from 1.8–5.0 Mb, amounting to a total genome size of about 29 Mb. The genomic DNA of strain ME-446, on the other hand, was resolved into 11 chromosomes, amounting to a total genome size of about 32Mb. Lignin peroxidase genes have been localized to five chromosomes in strain BKMF-1767 and to four chromosomes in strain ME-446.     

The predatory mite Metaseiulus occidentalis: mitey small and mitey large genomes

Metaseiulus occidentalis is a representative of an important family of mites (Arthropoda: Chelicerata: Acari: Phytoseiidae) that are effective predators of pest mites in agricultural crops around the world. Like many arthropods, this mite contains multiple genomes, including the genomes of several microbial symbionts as well as its own mitochondrial and nuclear genomes. The mitochondrial genome is “mitey” large at 25 kb, due to duplication and triplication of genes. By contrast, the nuclear genome is “mitey” small at 88 Mb. This mite has a parahaploid genetic system, tolerates inbreeding, and has a haploid chromosome number of 3. This predator was genetically improved for use in agriculture by developing strains that lacked the ability to overwinter in diapause or were resistant to multiple pesticides, and can be genetically modified using recombinant DNA methods. Sequencing the nuclear genome would provide useful insights that could enhance genetic improvement programs that would result in improved pest management, could provide genes needed to resolve the evolutionary relationships of this family, and could serve as a model for understanding the evolution and genetics of chelicerate arthropod predators.     

Chromosome‐level genome assembly of Triplophysa tibetana,a fish adapted to the harsh high‐altitude environment of the Tibetan Plateau

Triplophysa is an endemic fish genus of the Tibetan Plateau in China. Triplophysa tibetana, which lives at a recorded altitude of ~4,000 m and plays an important role in the highland aquatic ecosystem, serves as an excellent model for investigating high‐altitude environmental adaptation. However, evolutionary and conservation studies of T. tibetana have been limited by scarce genomic resources for the genus Triplophysa. In the present study, we applied PacBio sequencing and the Hi‐C technique to assemble the T. tibetana genome. A 652‐Mb genome with 1,325 contigs with an N50 length of 3.1 Mb was obtained. The 1,137 contigs were further assembled into 25 chromosomes, representing 98.7% and 80.47% of all contigs at the base and sequence number level, respectively. Approximately 260 Mb of sequence, accounting for ~39.8% of the genome, was identified as repetitive elements. DNA transposons (16.3%), long interspersed nuclear elements (12.4%) and long terminal repeats (11.0%) were the most repetitive types. In total, 24,372 protein‐coding genes were predicted in the genome, and ~95% of the genes were functionally annotated via a search in public databases. Using whole genome sequence information, we found that T. tibetana diverged from its common ancestor with Danio rerio ~121.4 million years ago. The high‐quality genome assembled in this work not only provides a valuable genomic resource for future population and conservation studies of T. tibetana, but it also lays a solid foundation for further investigation into the mechanisms of environmental adaptation of endemic fishes in the Tibetan Plateau.     

Genetic characterization of the yeast Pichia anomala (strain K), an antagonist of postharvest diseases of apple

AIMS: To obtain information about the genomic organization of Pichia anomala (strain K) and about its genomic diversity at species and intraspecies level. METHODS AND RESULTS: The PFGE karyotype of strain K was composed of four bands ranging in size from 1.1 to 3.2 Mb. The number of chromosomes was estimated at six since bands 2 and 3 seemed to result from the comigration of two chromosomes with similar size. A comparison of strain K and Hansenulawingeii migration profiles led to the estimate of K strain genome size at 11.7 Mb. Comparison with isogenic strains, resulting from the sporulation of strain K, highlighted some major karyotypic differences. Two segregants (KH6 and KH7) showed supernumerary chromosomes and one (KH9) displayed chromosomal length polymorphism. This genomic instability was confirmed by molecular hybridization with four probes, consisting of URA3, LEU2, PAEXG1 and PAEXG2 genes of P. anomala. URA3 and LEU2 probes showed second hybridization signals on supernumerary chromosomes of strain KH7 and on chromosome 6 of strain K for LEU2 only. Karyotypic comparison of seven non-isogenic P. anomala strains revealed chromosomal length polymorphism, a sign of intraspecies variation. CONCLUSIONS: This work has supplied information about genome size and chromosome number of strain K of P. anomala. The strain seems to be aneuploid because of the presence of supernumerary chromosomes and additional hybridization signals for URA3 and LEU2 probes in the chromosomal profile of some segregants. The work also highlighted genomic diversity within the P. anomala species. SIGNIFICANCE AND IMPACT OF THE STUDY: Results obtained here increase information about the aneuploidy of P. anomala (strain K). Information about the genomic diversity of the segregants will be of great interest for further studies on strain K mode of action. The genome size and chromosomal profile of P. anomala presented here are different from the results obtained elsewhere for Hansenula anomala, while Hansenula is included as a synonym of Pichia. This warrants further studies to investigate this taxonomic relationship.     

The chromosome-scale reference genome of Rubus chingii Hu provides insight into the biosynthetic pathway of hydrolyzable tannins

Rubus chingii Hu (Fu-Pen-Zi), a perennial woody plant in the Rosaceae family, is a characteristic traditional Chinese medicinal plant because of its unique pharmacological effects. There are abundant hydrolyzable tannin (HT) components in R. chingii that provide health benefits. Here, an R. chingii chromosome-scale genome and related functional analysis provide insights into the biosynthetic pathway of HTs. In total, sequence data of 231.21 Mb (155 scaffolds with an N50 of 8.2 Mb) were assembled into seven chromosomes with an average length of 31.4 Mb, and 33 130 protein-coding genes were predicted, 89.28% of which were functionally annotated. Evolutionary analysis showed that R. chingii was most closely related to Rubus occidentalis, from which it was predicted to have diverged 22.46 million years ago (Table S8). Comparative genomic analysis showed that there was a tandem gene cluster of UGT, carboxylesterase (CXE) and SCPL genes on chromosome 02 of R. chingii, including 11 CXE, eight UGT, and six SCPL genes, which may be critical for the synthesis of HTs. In vitro enzyme assays indicated that the proteins encoded by the CXE (LG02.4273) and UGT (LG02.4102) genes have tannin hydrolase and gallic acid glycosyltransferase functions, respectively. The genomic sequence of R. chingii will be a valuable resource for comparative genomic analysis within the Rosaceae family and will be useful for understanding the biosynthesis of HTs.     

Common ancestry and novel genetic traits of Francisella novicida-like isolates from North America and Australia as revealed by comparative genomic analyses

Francisella novicida is a close relative of Francisella tularensis, the causative agent of tularemia. The genomes of F. novicida-like clinical isolates 3523 (Australian strain) and Fx1 (Texas strain) were sequenced and compared to F. novicida strain U112 and F. tularensis strain Schu S4. The strain 3523 chromosome is 1,945,310 bp and contains 1,854 protein-coding genes. The strain Fx1 chromosome is 1,913,619 bp and contains 1,819 protein-coding genes. NUCmer analyses revealed that the genomes of strains Fx1 and U112 are mostly colinear, whereas the genome of strain 3523 has gaps, translocations, and/or inversions compared to genomes of strains Fx1 and U112. Using the genome sequence data and comparative analyses with other members of the genus Francisella, several strain-specific genes that encode putative proteins involved in RTX toxin production, polysaccharide biosynthesis/modification, thiamine biosynthesis, glucuronate utilization, and polyamine biosynthesis were identified. The RTX toxin synthesis and secretion operon of strain 3523 contains four open reading frames (ORFs) and was named rtxCABD. Based on the alignment of conserved sequences upstream of operons involved in thiamine biosynthesis from various bacteria, a putative THI box was identified in strain 3523. The glucuronate catabolism loci of strains 3523 and Fx1 contain a cluster of nine ORFs oriented in the same direction that appear to constitute an operon. Strains U112 and Schu S4 appeared to have lost the loci for RTX toxin production, thiamine biosynthesis, and glucuronate utilization as a consequence of host adaptation and reductive evolution. In conclusion, comparative analyses provided insights into the common ancestry and novel genetic traits of these strains.     

Chromosome profile ofLeishmania donovani: Interstrain and interspecific variations

The genome ofLeishmania donovani AG83, a virulent strain causing kala-azar, was resolved into 29 chromosomal bands by pulsed field gel electrophoresis (pFGE) under standardized conditions. Comparison of the karyotype with those of other strains and species revealed variations. By Southern hybridization, specific genes were localized to individual chromosomes. Twenty-two copies ofβ-tubulin genes are located on band 27 (1.63 Mb); minor copies are present in band 16 (850 kb) and band 9 (650 kb). Aβ-tubulin related nontranscribed locus was isolated from a genomic library and shown to contain repetitive sequences hybridizing throughout the genome. Single chromosomes contain multicopy clusters of gp63 and rnini-exon-derived RNA genes, but interspecific variations were observed in each case. The results emphasize the importance of using a standard reference strain ofLeishmania donovani for coordinated genome mapping of this clinically important organism.     

Chromosome-level genome assembly and characterization of Sophora Japonica

Sophora japonica is a medium-size deciduous tree belonging to Leguminosae family and famous for its high ecological, economic and medicinal value. Here, we reveal a draft genome of S. japonica, which was ∼511.49 Mb long (contig N50 size of 17.34 Mb) based on Illumina, Nanopore and Hi-C data. We reliably assembled 110 contigs into 14 chromosomes, representing 91.62% of the total genome, with an improved N50 size of 31.32 Mb based on Hi-C data. Further investigation identified 271.76 Mb (53.13%) of repetitive sequences and 31,000 protein-coding genes, of which 30,721 (99.1%) were functionally annotated. Phylogenetic analysis indicates that S. japonica separated from Arabidopsis thaliana and Glycine max ∼107.53 and 61.24 million years ago, respectively. We detected evidence of species-specific and common-legume whole-genome duplication events in S. japonica. We further found that multiple TF families (e.g. BBX and PAL) have expanded in S. japonica, which might have led to its enhanced tolerance to abiotic stress. In addition, S. japonica harbours more genes involved in the lignin and cellulose biosynthesis pathways than the other two species. Finally, population genomic analyses revealed no obvious differentiation among geographical groups and the effective population size continuously declined since 2 Ma. Our genomic data provide a powerful comparative framework to study the adaptation, evolution and active ingredients biosynthesis in S. japonica. More importantly, our high-quality S. japonica genome is important for elucidating the biosynthesis of its main bioactive components, and improving its production and/or processing.     

Organization and evolution of mitochondrial gene clusters in human

Currently, the spatial patterns of mitochondrial genes and how the genomic localization of (pseudo)genes originated from mitochondrial DNA remain largely unexplained. The aim of this study was to elucidate the organization of mitochondrial (pseudo)genes given their evolutionary origin. We used a keyword finding method and a bootstrapping method to estimate parameter values that represent the distribution pattern of mitochondrial genes in the nuclear genome. Almost half of mitochondrial genes showing physical clusters were located in the pericentromeric and subtelomeric regions of the chromosome. Most interestingly, the size of these clusters ranged from 0.085 to 3.2 Mb (average ± SD 1.3 ± 0.73 Mb), which coincides with the size of the evolutionary pocket, or the average size of evolutionary breakpoint regions. Our findings imply that the localization of mitochondrial genes in the human genome is determined independent of adaptation.     

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.     

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.     

一株太平洋表层海水来源阿拉伯海假交替单胞菌N1230-9的全基因组测序及比较基因组学分析

【目的】阐述阿拉伯海假交替单胞菌(Pseudoalteromonas arabiensis)的分子进化与生态适应策略。【方法】借助Illumina HiSeq X Ten和Oxford Nanopore PromethION测序平台对一株分离自太平洋表层海水的菌株Pseudoalteromonas arabiensis N1230-9进行全基因组测序,利用相关生物信息学分析软件对原始数据进行组装和基因组注释,并与一株分离自深海沉积物环境的模式菌株Pseudoalteromonas arabiensis JCM 17292进行比较基因组分析。【结果】菌株N1230-9基因组由2条染色体组成,基因组大小为4 627 470 bp,G+C含量为40.85%,共编码4 202个蛋白。功能基因注释结果显示,菌株N1230-9具有适应海洋环境的多种类型功能基因,包括重金属抗性基因、多种铁离子摄取系统编码基因、多种噬菌体防御系统编码基因、种类丰富的水解酶编码基因、多种碳水化合物代谢相关基因,以及数量众多的二元信号传导系统编码基因。通过比较基因组分析发现,菌株N1230-9和菌株JCM 17292拥有适应不同生态位的特有基因,这些基因主要涉及血红素摄取、重金属抗性、噬菌体防御、二元信号系统传导和横向基因水平转移。【结论】来自表层海水的菌株P. arabiensis N1230-9已演化出了适应其生态位的特有基因。     

Mitogenomics reveals two cryptic species in Ciona intestinalis

Individual mitochondrial genes or genomic features are commonly used as phylogenetic markers at many taxonomic levels. We used a mitogenomics approach to demonstrate the existence of two cryptic species in the ascidian Ciona intestinalis, a model chordate whose status as a single species has recently been questioned. Comprehensive comparative analysis of the mitochondrial genome of the two cryptic species revealed significant differences in gene order, size and number of noncoding regions, compositional features and divergence of protein-coding genes.     

六妹羊肚菌PacBio基因组测序和DNA 6mA甲基化分析

羊肚菌属于子囊门真菌,是世界范围内最受欢迎的食用菌之一。本研究通过PacBio单分子实时测序技术对我国四川省成功栽培的六妹羊肚菌Morchella sextelata SCLS菌株进行全基因组测序,获得其高质量核基因组组装,大小为53.57Mb,重复序列含量为17.03%,包含42条重叠群（contigs）,重叠群N50高达1.82Mb,其中13条重叠群两端均含有端粒重复序列,为完整的染色体。通过链特异性RNA-seq测序和转录本拼接,并结合多种基因预测策略,预测到13 182个蛋白编码基因,包括267个碳水化合物活性酶,11个次生代谢产物合成基因簇。通过与内蒙古地区栽培的六妹羊肚菌菌株NZTD180501373基因组比较发现,六妹羊肚菌进化过程中可能发生过染色体重组事件,二者具有33 055个SNP和48 726个InDel位点差异,并且各自拥有超过6Mb的特有序列。此外,相比于梯棱羊肚菌M. importuna,六妹羊肚菌SCLS菌株中与逆转录转座酶有关的orthogroup发生了扩张,并且拥有5个成员数超过100的特有逆转录转座酶基因家族。对SCLS菌株中的DNA N⁶腺嘌呤甲基化（6mA）修饰进行了鉴定,发现SCLS菌株基因组中0.42%的腺嘌呤被6mA甲基化,其含量显著高于已报道的6种双核亚界（Dikarya）真菌。6mA甲基化位点在逆转录转座子上显著富集,表明六妹羊肚菌中6mA甲基化可能调控逆转录转座子的活性,这也是首次在真菌中报道6mA甲基化与转座子相关。