Menu
September 22, 2019

Whole-genome resequencing and pan-transcriptome reconstruction highlight the impact of genomic structural Variation on secondary metabolite gene clusters in the grapevine Esca pathogen Phaeoacremonium minimum.

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.


September 22, 2019

Phenotypic and genomic properties of Brachybacterium vulturis sp. nov. and Brachybacterium avium sp. nov.

Two strains, VM2412T and VR2415T, were isolated from the feces of an Andean condor (Vultur gryphus) living in Seoul Grand Park, Gyeonggi-do, South Korea. Cells of both strains were observed to be Gram-stain positive, non-motile, aerobic, catalase positive and oxidase negative. Growth was found to occur at 10-30°C, showing optimum growth at 30°C. The strains could tolerate up to 15% (w/v) NaCl concentration and grow at pH 6-9. The strains shared 99.3% 16S rRNA gene sequence similarity to each other but were identified as two distinct species based on 89.0-89.2% ANIb, 90.3% ANIm, 89.7% OrthoANI and 38.0% dDDH values calculated using whole genome sequences. Among species with validly published names, Brachybacterium ginsengisoli DCY80T shared high 16S rRNA gene sequence similarities with strains VM2412T (98.7%) and VR2415T (98.4%) and close genetic relatedness with strains VM2412T (83.3-83.5% ANIb, 87.0% ANIm, 84.3% OrthoANI and 27.8% dDDH) and VR2415T (82.8-83.2% ANIb, 86.7% ANIm, 83.9% OrthoANI and 27.2% dDDH). The major fatty acid of the two strains was identified as anteiso-C15:0 and the polar lipids consisted of phosphatidylglycerol, diphosphatidylglycerol, presumptively phosphatidylethanolamine and three unidentified glycolipids. Strain VR2415T also produced an unidentified phospholipid. The cell walls of the two strains contained meso-diaminopimelic acid as diagnostic diamino acid and the whole cell sugars were ribose, glucose, and galactose. The strains contained MK-7 as their predominant menaquinone. The genomes of strains VM2412T, VR2415T, and B. ginsengisoli DCY80T were sequenced in this study. The genomic G+C contents of strains VM2412T and VR2415T were determined to be 70.8 and 70.4 mol%, respectively. A genome-based phylogenetic tree constructed using an up-to-date bacterial core gene set (UBCG) showed that the strains formed a clade with members of the genus Brachybacterium, supporting their taxonomic classification into the genus Brachybacterium. Based on phenotypic and genotypic analyses in this study, strains VM2412T and VR2415T are considered to represent two novel species of the genus Brachybacterium and the names Brachybacterium vulturis sp. nov. and Brachybacterium avium sp. nov. are proposed for strains VM2412T (=KCTC 39996T = JCM 32142T) and VR2415T (=KCTC 39997T = JCM 32143T), respectively.


September 22, 2019

Complete genome of streamlined marine actinobacterium Pontimonas salivibrio strain CL-TW6T adapted to coastal planktonic lifestyle.

Pontimonas salivibrio strain CL-TW6T (=KCCM 90105?=?JCM18206) was characterized as the type strain of a new genus within the Actinobacterial family Microbacteriaceae. It was isolated from a coastal marine environment in which members of Microbactericeae have not been previously characterized.The genome of P. salivibrio CL-TW6T was a single chromosome of 1,760,810 bp. Genomes of this small size are typically found in bacteria growing slowly in oligotrophic zones and said to be streamlined. Phylogenetic analysis showed it to represent a lineage originating in the Microbacteriaceae radiation occurring before the snowball Earth glaciations, and to have a closer relationship with some streamlined bacteria known through metagenomic data. Several genomic characteristics typical of streamlined bacteria are found: %G?+?C is lower than non-streamlined members of the phylum; there are a minimal number of rRNA and tRNA genes, fewer paralogs in most gene families, and only two sigma factors; there is a noticeable absence of some nonessential metabolic pathways, including polyketide synthesis and catabolism of some amino acids. There was no indication of any phage genes or plasmids, however, a system of active insertion elements was present. P. salivibrio appears to be unusual in having polyrhamnose-based cell wall oligosaccharides instead of mycolic acid or teichoic acid-based oligosaccharides. Oddly, it conducts sulfate assimilation apparently for sulfating cell wall components, but not for synthesizing amino acids. One gene family it has more of, rather than fewer of, are toxin/antitoxin systems, which are thought to down-regulate growth during nutrient deprivation or other stressful conditions.Because of the relatively small number of paralogs and its relationship to the heavily characterized Mycobacterium tuberculosis, we were able to heavily annotate the genome of P. salivibrio CL-TW6T. Its streamlined status and relationship to streamlined metagenomic constructs makes it an important reference genome for study of the streamlining concept. The final evolutionary trajectory of CL-TW6 T was to adapt to growth in a non-oligotrophic coastal zone. To understand that adaptive process, we give a thorough accounting of gene content, contrasting with both oligotrophic streamlined bacteria and large genome bacteria, and distinguishing between genes derived by vertical and horizontal descent.


September 22, 2019

Temperature responses of mutation rate and mutational spectrum in an Escherichia coli strain and the correlation with metabolic rate.

Temperature is a major determinant of spontaneous mutation, but the precise mode, and the underlying mechanisms, of the temperature influences remain less clear. Here we used a mutation accumulation approach combined with whole-genome sequencing to investigate the temperature dependence of spontaneous mutation in an Escherichia coli strain. Experiments were performed under aerobic conditions at 25, 28 and 37 °C, three temperatures that were non-stressful for the bacterium but caused significantly different bacterial growth rates.Mutation rate did not differ between 25 and 28 °C, but was higher at 37 °C. Detailed analyses of the molecular spectrum of mutations were performed; and a particularly interesting finding is that higher temperature led to a bias of mutation to coding, relative to noncoding, DNA. Furthermore, the temperature response of mutation rate was extremely similar to that of metabolic rate, consistent with an idea that metabolic rate predicts mutation rate.Temperature affects mutation rate and the types of mutation supply, both being crucial for the opportunity of natural selection. Our results help understand how temperature drives evolutionary speed of organisms and thus the global patterns of biodiversity. This study also lend support to the metabolic theory of ecology for linking metabolic rate and molecular evolution rate.


September 22, 2019

Reprogramming of the antimycin NRPS-PKS assembly lines inspired by gene evolution.

Reprogramming of the NRPS/PKS assembly line is an attractive method for the production of new bioactive molecules. However, it is usually hampered by the loss of intimate domain/module interactions required for the precise control of chain transfer and elongation reactions. In this study, we first establish heterologous expression systems of the unique antimycin-type cyclic depsipeptides: JBIR-06 (tri-lactone) and neoantimycin (tetra-lactone), and engineer their biosyntheses by taking advantage of bioinformatic analyses and evolutionary insights. As a result, we successfully accomplish three manipulations: (i) ring contraction of neoantimycin (from tetra-lactone to tri-lactone), (ii) ring expansion of JBIR-06 (from tri-lactone to tetra-lactone), and (iii) alkyl chain diversification of JBIR-06 by the incorporation of various alkylmalonyl-CoA extender units, to generate a set of unnatural derivatives in practical yields. This study presents a useful strategy for engineering NRPS-PKS module enzymes, based on nature’s diversification of the domain and module organizations.


September 22, 2019

Orphan legumes growing in dry environments: Marama bean as a case study.

Plants have developed morphological, physiological, biochemical, cellular, and molecular mechanisms to survive in drought-stricken environments with little or no water caused by below-average precipitation. In this mini-review, we highlight the characteristics that allows marama bean [Tylosema esculentum (Burchell) Schreiber], an example of an orphan legume native to arid regions of southwestern Southern Africa, to flourish under an inhospitable climate and dry soil conditions where no other agricultural crop competes in this agro-ecological zone. Orphan legumes are often better suited to withstand such harsh growth environments due to development of survival strategies using a combination of different traits and responses. Recent findings on questions on marama bean speciation, hybridization, population dynamics, and the evolutionary history of the bean and mechanisms by which the bean is able to extract and conserve water and nutrients from its environment as well as aspects of morphological and physiological adaptation will be reviewed. The importance of the soil microbiome and the genetic diversity in this species, and their interplay, as a reservoir for improvement will also be considered. In particular, the application of the newly established marama bean genome sequence will facilitate both the identification of important genes involved in the interaction with the soil microbiome and the identification of the diversity within the wild germplasm for genes involved drought tolerance. Since predicted future changes in climatic conditions, with less water availability for plant growth, will severely affect agricultural productivity, an understanding of the mechanisms of unique adaptations in marama bean to such conditions may also provide insights as to how to improve the performance of the major crops.


September 22, 2019

Discovery of new genes involved in curli production by a uropathogenic Escherichia coli strain from the highly virulent O45:K1:H7 lineage.

Curli are bacterial surface-associated amyloid fibers that bind to the dye Congo red (CR) and facilitate uropathogenic Escherichia coli (UPEC) biofilm formation and protection against host innate defenses. Here we sequenced the genome of the curli-producing UPEC pyelonephritis strain MS7163 and showed it belongs to the highly virulent O45:K1:H7 neonatal meningitis-associated clone. MS7163 produced curli at human physiological temperature, and this correlated with biofilm growth, resistance of sessile cells to the human cationic peptide cathelicidin, and enhanced colonization of the mouse bladder. We devised a forward genetic screen using CR staining as a proxy for curli production and identified 41 genes that were required for optimal CR binding, of which 19 genes were essential for curli synthesis. Ten of these genes were novel or poorly characterized with respect to curli synthesis and included genes involved in purine de novo biosynthesis, a regulator that controls the Rcs phosphorelay system, and a novel repressor of curli production (referred to as rcpA). The involvement of these genes in curli production was confirmed by the construction of defined mutants and their complementation. The mutants did not express the curli major subunit CsgA and failed to produce curli based on CR binding. Mutation of purF (the first gene in the purine biosynthesis pathway) and rcpA also led to attenuated colonization of the mouse bladder. Overall, this work has provided new insight into the regulation of curli and the role of these amyloid fibers in UPEC biofilm formation and pathogenesis.IMPORTANCE Uropathogenic Escherichia coli (UPEC) strains are the most common cause of urinary tract infection, a disease increasingly associated with escalating antibiotic resistance. UPEC strains possess multiple surface-associated factors that enable their colonization of the urinary tract, including fimbriae, curli, and autotransporters. Curli are extracellular amyloid fibers that enhance UPEC virulence and promote biofilm formation. Here we examined the function and regulation of curli in a UPEC pyelonephritis strain belonging to the highly virulent O45:K1:H7 neonatal meningitis-associated clone. Curli expression at human physiological temperature led to increased biofilm formation, resistance of sessile cells to the human cationic peptide LL-37, and enhanced bladder colonization. Using a comprehensive genetic screen, we identified multiple genes involved in curli production, including several that were novel or poorly characterized with respect to curli synthesis. In total, this study demonstrates an important role for curli as a UPEC virulence factor that promotes biofilm formation, resistance, and pathogenesis. Copyright © 2018 Nhu et al.


September 22, 2019

Comprehensive evaluation of the host responses to infection with differentially virulent classical swine fever virus strains in pigs.

Classical swine fever virus (CSFV) infection causes most variable clinical syndromes from chronic or latent infection to acute death, and it is generally acknowledged that the course of disease is affected by both virus and host factors. To compare host immune responses to differentially virulent CSFV strains in pigs, fifteen 8-week-old specific-pathogen-free pigs were randomly divided into four groups and inoculated with the CSFV Shimen strain (a highly virulent strain), the HLJZZ2014 strain (a moderately virulent strains), C-strain (an avirulent strain), and DMEM (mock control), respectively. Infection with the Shimen or HLJZZ2014 strain resulted in fever, clinical signs and histopathological lesions, which were not observed in the C-strain-inoculated pigs, though low viral genome copies were detected in the peripheral blood and tissue samples. The data showed that the virulence of the strains affected the outcome of duration and intensity of the disease rather than the tissue tropism of the virus. Furthermore, leukopenia, lymphocytopenia, differentiation of T-cells, and the secretion of cytokines associated with inflammation or apoptosis such as interferon alpha (IFN-a), tumor necrosis factor alpha (TNF-a), interleukin 2 (IL-2), IL-4, IL-6, and IL-10 were induced by the virulent CSFV infection, the differences reflected in onset and extent of the regulation. Taken together, our results revealed that the major differences among the three strains resided in the kinetics of host response to the infection: severe and immediate with the highly virulent strain, while progressive and delayed with the moderately virulent one. This comparative study will help to dissect the pathogenesis of CSFV. Copyright © 2018 Elsevier B.V. All rights reserved.


September 22, 2019

Opposite polarity monospore genome de novo sequencing and comparative analysis reveal the possible heterothallic life cycle of Morchella importuna.

Morchella is a popular edible fungus worldwide due to its rich nutrition and unique flavor. Many research efforts were made on the domestication and cultivation of Morchella all over the world. In recent years, the cultivation of Morchella was successfully commercialized in China. However, the biology is not well understood, which restricts the further development of the morel fungus cultivation industry. In this paper, we performed de novo sequencing and assembly of the genomes of two monospores with a different mating type (M04M24 and M04M26) isolated from the commercially cultivated strain M04. Gene annotation and comparative genome analysis were performed to study differences in CAZyme (Carbohydrate-active enzyme) enzyme content, transcription factors, duplicated sequences, structure of mating type sites, and differences at the gene and functional levels between the two monospore strains of M. importuna. Results showed that the de novo assembled haploid M04M24 and M04M26 genomes were 48.98 and 51.07 Mb, respectively. A complete fine physical map of M. importuna was obtained from genome coverage and gene completeness evaluation. A total of 10,852 and 10,902 common genes and 667 and 868 endemic genes were identified from the two monospore strains, respectively. The Gene Ontology (GO) and KAAS (KEGG Automatic Annotation Serve) enrichment analyses showed that the endemic genes performed different functions. The two monospore strains had 99.22% collinearity with each other, accompanied with certain position and rearrangement events. Analysis of complete mating-type loci revealed that the two monospore M. importuna strains contained an independent mating-type structure and remained conserved in sequence and location. The phylogenetic and divergence time of M. importuna was analyzed at the whole-genome level for the first time. The bifurcation time of morel and tuber was estimated to be 201.14 million years ago (Mya); the two monospore strains with a different mating type represented the evolution of different nuclei, and the single copy homologous genes between them were also different due to a genetic differentiation distance about 0.65 Mya. Compared with truffles, M. importuna had an extension of 28 clusters of orthologous genes (COGs) and a contraction of two COGs. The two different polar nuclei with different degrees of contraction and expansion suggested that they might have undergone different evolutionary processes. The different mating-type structures, together with the functional clustering and enrichment analysis results of the endemic genes of the two different polar nuclei, imply that M. importuna might be a heterothallic fungus and the interaction between the endemic genes may be necessary for its complete life history. Studies on the genome of M. importuna facilitate a better understanding of morel biology and evolution.


September 22, 2019

Molecular basis for the final oxidative rearrangement steps in chartreusin biosynthesis.

Oxidative rearrangements play key roles in introducing structural complexity and biological activities of natural products biosynthesized by type II polyketide synthases (PKSs). Chartreusin (1) is a potent antitumor polyketide that contains a unique rearranged pentacyclic aromatic bilactone aglycone derived from a type II PKS. Herein, we report an unprecedented dioxygenase, ChaP, that catalyzes the final a-pyrone ring formation in 1 biosynthesis using flavin-activated oxygen as an oxidant. The X-ray crystal structures of ChaP and two homologues, docking studies, and site-directed mutagenesis provided insights into the molecular basis of the oxidative rearrangement that involves two successive C-C bond cleavage steps followed by lactonization. ChaP is the first example of a dioxygenase that requires a flavin-activated oxygen as a substrate despite lacking flavin binding sites, and represents a new class in the vicinal oxygen chelate enzyme superfamily.


September 22, 2019

Characterization of Haemophilus parasuis serovar 2 CL120103, a moderately virulent strain in China

Haemophilus parasuis is an important bacterium affecting pigs, causing Glässer’s disease. To further characterize this species, we determined the complete genomic sequence of H. parasuis CL120103, which was isolated from diseased pigs. The strain H. parasuis CL120103 was identified as serovar 2. The size of the largest scaffold is 2,326,318 bp and contains 145 large contigs, with the N50 contig being 20,573 bp in length. The complete genome of H. parasuis CL120103 is 2,305,354 bp in length with 39.97% GC content and contains 2227 protein-coding genes, 19 ribosomal rRNA operons and 60 tRNA genes. Sequence similarity of the genome of H. parasuis CL120103 to the previously sequenced genome of H. parasuis was up to 96% and query cover to 86%. Annotation of the genome of H. parasuis CL120103 identified a number of genes encoding potential virulence factors. These virulence factors are involved in metabolism, adhesion, secretion and LPS biosynthesis. These related genes pave the way to better understand mechanisms underlying metabolic capabilities. The comprehensive genetic and phylogenetic analysis shows that H. parasuis is closely related to Actinobacillus pleuropneumoniae and provides a foundation for future experimental confirmation of the virulence and pathogen-host interactions in H. parasuis.


September 22, 2019

Comparative genomics of Salmonella enterica serovar Montevideo reveals lineage-specific gene differences that may influence ecological niche association.

Salmonella enterica serovar Montevideo has been linked to recent foodborne illness outbreaks resulting from contamination of products such as fruits, vegetables, seeds and spices. Studies have shown that Montevideo also is frequently associated with healthy cattle and can be isolated from ground beef, yet human salmonellosis outbreaks of Montevideo associated with ground beef contamination are rare. This disparity fuelled our interest in characterizing the genomic differences between Montevideo strains isolated from healthy cattle and beef products, and those isolated from human patients and outbreak sources. To that end, we sequenced 13 Montevideo strains to completion, producing high-quality genome assemblies of isolates from human patients (n=8) or from healthy cattle at slaughter (n=5). Comparative analysis of sequence data from this study and publicly available sequences (n=72) shows that Montevideo falls into four previously established clades, differentially occupied by cattle and human strains. The results of these analyses reveal differences in metabolic islands, environmental adhesion determinants and virulence factors within each clade, and suggest explanations for the infrequent association between bovine isolates and human illnesses.


September 22, 2019

Distinct genomic features characterize two clades of Corynebacterium diphtheriae: Proposal of Corynebacterium diphtheriae subsp. diphtheriae subsp. nov. and Corynebacterium diphtheriae subsp. lausannense subsp. nov.

Corynebacterium diphtheriae is the etiological agent of diphtheria, a disease caused by the presence of the diphtheria toxin. However, an increasing number of records report non-toxigenic C. diphtheriae infections. Here, a C. diphtheriae strain was recovered from a patient with a past history of bronchiectasis who developed a severe tracheo-bronchitis with multiple whitish lesions of the distal trachea and the mainstem bronchi. Whole-genome sequencing (WGS), performed in parallel with PCR targeting the toxin gene and the Elek test, provided clinically relevant results in a short turnaround time, showing that the isolate was non-toxigenic. A comparative genomic analysis of the new strain (CHUV2995) with 56 other publicly available genomes of C. diphtheriae revealed that the strains CHUV2995, CCUG 5865 and CMCNS703 share a lower average nucleotide identity (ANI) (95.24 to 95.39%) with the C. diphtheriae NCTC 11397T reference genome than all other C. diphtheriae genomes (>98.15%). Core genome phylogeny confirmed the presence of two monophyletic clades. Based on these findings, we propose here two new C. diphtheriae subspecies to replace the lineage denomination used in previous multilocus sequence typing studies: C. diphtheriae subsp. lausannense subsp. nov. (instead of lineage-2), regrouping strains CHUV2995, CCUG 5865, and CMCNS703, and C. diphtheriae subsp. diphtheriae subsp. nov, regrouping all other C. diphtheriae in the dataset (instead of lineage-1). Interestingly, members of subspecies lausannense displayed a larger genome size than subspecies diphtheriae and were enriched in COG categories related to transport and metabolism of lipids (I) and inorganic ion (P). Conversely, they lacked all genes involved in the synthesis of pili (SpaA-type, SpaD-type and SpaH-type), molybdenum cofactor and of the nitrate reductase. Finally, the CHUV2995 genome is particularly enriched in mobility genes and harbors several prophages. The genome encodes a type II-C CRISPR-Cas locus with 2 spacers that lacks csn2 or cas4, which could hamper the acquisition of new spacers and render strain CHUV2995 more susceptible to bacteriophage infections and gene acquisition through various mechanisms of horizontal gene transfer.


September 22, 2019

The energy-coupling factor transporter module EcfAA’T, a novel candidate for the genetic basis of fatty acid-auxotrophic small-colony variants of Staphylococcus aureus.

Staphylococcal small-colony variants (SCVs) are invasive and persistent due to their ability to thrive intracellularly and to evade the host immune response. Thus, the course of infections due to this phenotype is often chronic, relapsing, and therapy-refractory. In order to improve treatment of patients suffering from SCV-associated infections, it is of major interest to understand triggers for the development of this phenotype, in particular for strains naturally occurring in clinical settings. Within this study, we comprehensively characterized two different Staphylococcus aureus triplets each consisting of isogenic strains comprising (i) clinically derived SCV phenotypes with auxotrophy for unsaturated fatty acids, (ii) the corresponding wild-types (WTs), and (iii) spontaneous in vitro revertants displaying the normal phenotype (REVs). Comparison of whole genomes revealed that clinical SCV isolates were closely related to their corresponding WTs and REVs showing only seven to eight alterations per genome triplet. However, both SCVs carried a mutation within the energy-coupling factor (ECF) transporter-encoding ecf module (EcfAA’T) resulting in truncated genes. In both cases, these mutations were shown to be naturally restored in the respective REVs. Since ECF transporters are supposed to be essential for optimal bacterial growth, their dysfunction might constitute another mechanism for the formation of naturally occurring SCVs. Another three triplets analyzed revealed neither mutations in the EcfAA’T nor in other FASII-related genes underlining the high diversity of mechanisms leading to the fatty acid-dependent phenotype. This is the first report on the ECF transporter as genetic basis of fatty acid-auxotrophic staphylococcal SCVs.


September 22, 2019

Assembly and comparative analysis of the complete mitochondrial genome sequence of Sophora japonica ‘JinhuaiJ2’.

Sophora japonica L. (Faboideae, Leguminosae) is an important traditional Chinese herb with a long history of cultivation. Its flower buds and fruits contain abundant flavonoids, and therefore, the plants are cultivated for the industrial extraction of rutin. Here, we determined the complete nucleotide sequence of the mitochondrial genome of S. japonica ‘JinhuaiJ2’, the most widely planted variety in Guangxi region of China. The total length of the mtDNA sequence is 484,916 bp, with a GC content of 45.4%. Sophora japonica mtDNA harbors 32 known protein-coding genes, 17 tRNA genes, and three rRNA genes with 17 cis-spliced and five trans-spliced introns disrupting eight protein-coding genes. The gene coding and intron regions, and intergenic spacers account for 7.5%, 5.8% and 86.7% of the genome, respectively. The gene profile of S. japonica mitogenome differs from that of the other Faboideae species by only one or two gene gains or losses. Four of the 17 cis-spliced introns showed distinct length variations in the Faboideae, which could be attributed to the homologous recombination of the short repeats measuring a few bases located precisely at the edges of the putative deletions. This reflects the importance of small repeats in the sequence evolution in Faboideae mitogenomes. Repeated sequences of S. japonica mitogenome are mainly composed of small repeats, with only 20 medium-sized repeats, and one large repeat, adding up to 4% of its mitogenome length. Among the 25 pseudogene fragments detected in the intergenic spacer regions, the two largest ones and their corresponding functional gene copies located in two different sets of medium-sized repeats, point to their origins from homologous recombinations. As we further observed the recombined reads associated with the longest repeats of 2,160 bp with the PacBio long read data set of just 15 × in depth, repeat mediated homologous recombinations may play important role in the mitogenomic evolution of S. japonica. Our study provides insightful knowledge to the genetic background of this important herb species and the mitogenomic evolution in the Faboideae species.


Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.