The genome sequence of the biocontrol fungus Metarhizium anisopliae and comparative genomics of Metarhizium species

Background

Metarhizium anisopliae is an important fungal biocontrol agent of insect pests of agricultural crops. Genomics can aid the successful commercialization of biopesticides by identification of key genes differentiating closely related species, selection of virulent microbial isolates which are amenable to industrial scale production and formulation and through the reduction of phenotypic variability. The genome of Metarhizium isolate ARSEF23 was recently published as a model for M. anisopliae, however phylogenetic analysis has since re-classified this isolate as M. robertsii. We present a new annotated genome sequence of M. anisopliae (isolate Ma69) and whole genome comparison to M. robertsii (ARSEF23) and M. acridum (CQMa 102).

Results

Whole genome analysis of M. anisopliae indicates significant macrosynteny with M. robertsii but with some large genomic inversions. In comparison to M. acridum, the genome of M. anisopliae shares lower sequence homology. While alignments overall are co-linear, the genome of M. acridum is not contiguous enough to conclusively observe macrosynteny. Mating type gene analysis revealed both MAT1-1 and MAT1-2 genes present in M. anisopliae suggesting putative homothallism, despite having no known teleomorph, in contrast with the putatively heterothallic M. acridum isolate CQMa 102 (MAT1-2) and M. robertsii isolate ARSEF23 (altered MAT1-1). Repetitive DNA and RIP analysis revealed M. acridum to have twice the repetitive content of the other two species and M. anisopliae to be five times more RIP affected than M. robertsii. We also present an initial bioinformatic survey of candidate pathogenicity genes in M. anisopliae.

Conclusions

The annotated genome of M. anisopliae is an important resource for the identification of virulence genes specific to M. anisopliae and development of species- and strain- specific assays. New insight into the possibility of homothallism and RIP affectedness has important implications for the development of M. anisopliae as a biopesticide as it may indicate the potential for greater inherent diversity in this species than the other species. This could present opportunities to select isolates with unique combinations of pathogenicity factors, or it may point to instability in the species, a negative attribute in a biopesticide.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-660) contains supplementary material, which is available to authorized users.     

DNA Complexity of the Root-knot Nematode (Meloidogyne spp.) Genome

Cot curves derived from renaturation kinetics of sheared denatured DNA indicated that the genome of six populations representing the four most common root-knot nematode species (Meloidogyne incognita, M. arenaria, M. javanica, and M. hapla) is composed of 20% repetitive and 80% nonrepetitive sequences of DNA. Cot curves were almost identical, indicating that all populations had a haploid genome of approximately the same size. Calculations from an average Cot curve gave an estimate of 0.51 x 108 nucleotide base pairs for the haploid genome of the four Meloidogyne species. This genome is about 12-13 times larger than the genome of the E. coli strain used as a control.     

Comparative genome analysis of rice-pathogenic Burkholderia provides insight into capacity to adapt to different environments and hosts

Background

In addition to human and animal diseases, bacteria of the genus Burkholderia can cause plant diseases. The representative species of rice-pathogenic Burkholderia are Burkholderia glumae, B. gladioli, and B. plantarii, which primarily cause grain rot, sheath rot, and seedling blight, respectively, resulting in severe reductions in rice production. Though Burkholderia rice pathogens cause problems in rice-growing countries, comprehensive studies of these rice-pathogenic species aiming to control Burkholderia-mediated diseases are only in the early stages.

Results

We first sequenced the complete genome of B. plantarii ATCC 43733^T. Second, we conducted comparative analysis of the newly sequenced B. plantarii ATCC 43733^T genome with eleven complete or draft genomes of B. glumae and B. gladioli strains. Furthermore, we compared the genome of three rice Burkholderia pathogens with those of other Burkholderia species such as those found in environmental habitats and those known as animal/human pathogens. These B. glumae, B. gladioli, and B. plantarii strains have unique genes involved in toxoflavin or tropolone toxin production and the clustered regularly interspaced short palindromic repeats (CRISPR)-mediated bacterial immune system. Although the genome of B. plantarii ATCC 43733^T has many common features with those of B. glumae and B. gladioli, this B. plantarii strain has several unique features, including quorum sensing and CRISPR/CRISPR-associated protein (Cas) systems.

Conclusions

The complete genome sequence of B. plantarii ATCC 43733^T and publicly available genomes of B. glumae BGR1 and B. gladioli BSR3 enabled comprehensive comparative genome analyses among three rice-pathogenic Burkholderia species responsible for tissue rotting and seedling blight. Our results suggest that B. glumae has evolved rapidly, or has undergone rapid genome rearrangements or deletions, in response to the hosts. It also, clarifies the unique features of rice pathogenic Burkholderia species relative to other animal and human Burkholderia species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1558-5) contains supplementary material, which is available to authorized users.     

Genome analysis reveals three genomospecies in Mycobacterium abscessus

Background

Mycobacterium abscessus complex, the third most frequent mycobacterial complex responsible for community- and health care-associated infections in developed countries, comprises of M. abscessus subsp. abscessus and M. abscessus subsp. bolletii reviously referred as Mycobacterium bolletii and Mycobacterium massiliense. The diversity of this group of opportunistic pathogens is poorly described.

Results

In-depth analysis of 14 published M. abscessus complex genomes found a pan-genome of 6,153 proteins and core-genome of 3,947 (64.1%) proteins, indicating a non-conservative genome. Analysing the average percentage of amino-acid sequence identity (from 94.19% to 98.58%) discriminates three main clusters C1, C2 and C3: C1 comprises strains belonging to M. abscessus, C2 comprises strains belonging to M. massiliense and C3 comprises strains belonging to M. bolletii; and two sub-clusters in clusters C2 and C3. The phylogenomic network confirms these three clusters. The genome length (from 4.8 to 5.51-Mb) varies from 5.07-Mb in C1, 4.89-Mb in C2A, 5.01-Mb in C2B and 5.28-Mb in C3. The mean number of prophage regions (from 0 to 7) is 2 in C1; 1.33 in C2A; 3.5 in C2B and five in C3. A total of 36 genes are uniquely present in C1, 15 in C2 and 15 in C3. These genes could be used for the detection and identification of organisms in each cluster. Further, the mean number of host-interaction factors (including PE, PPE, LpqH, MCE, Yrbe and type VII secretion system ESX3 and ESX4) varies from 70 in cluster C1, 80 in cluster C2A, 74 in cluster C2B and 93 in clusters C3A and C3B. No significant differences in antibiotic resistance genes were observed between clusters, in contrast to previously reported in-vitro patterns of drug resistance. They encode both penicillin-binding proteins targeted by β-lactam antibiotics and an Ambler class A β-lactamase for which inhibitors exist.

Conclusions

Our comparative analysis indicates that M. abscessus complex comprises three genomospecies, corresponding to M. abscessus, M. bolletii, and M. massiliense. The genomics data here reported indicate differences in virulence of medical interest; and suggest targets for the refined detection and identification of M. abscessus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-359) contains supplementary material, which is available to authorized users.     

How endogenous plant pararetroviruses shed light on Musa evolution

Background and Aims Banana genomes harbour numerous copies of viral sequences derived from banana streak viruses (BSVs) – dsDNA viruses belonging to the family Caulimoviridae. These viral integrants (eBSVs) are mostly defective, probably as a result of ‘pseudogenization’ driven by host genome evolution. However, some can give rise to infection by releasing a functional viral genome following abiotic stresses. These distinct infective eBSVs correspond to the three main widespread BSV species (BSOLV, BSGFV and BSIMV), fully described within the Musa balbisiana B genomes of the seedy diploid ‘Pisang Klutuk Wulung’ (PKW).Methods We characterize eBSV distribution among a Musa sampling including seedy BB diploids and interspecific hybrids with Musa acuminata exhibiting different levels of ploidy for the B genome (ABB, AAB, AB). We used representative samples of the two areas of sympatry between M. acuminata and M. balbisiana species representing the native area of the most widely cultivated AAB cultivars (in India and in East Asia, ranging from the Philippines to New Guinea). Seventy-seven accessions were characterized using eBSV-related PCR markers and Southern hybridization approaches. We coded both sets of results to create a common dissimilarity matrix with which to interpret eBSV distribution.Key Results We propose a Musa phylogeny driven by the M. balbisiana genome based on a dendrogram resulting from a joint neighbour-joining analysis of the three BSV species, showing for the first time lineages between BB and ABB/AAB hybrids. eBSVs appear to be relevant phylogenetic markers that can illustrate the M. balbisiana phylogeography story.Conclusion The theoretical implications of this study for further elucidation of the historical and geographical process of Musa domestication are numerous. Discovery of banana plants with B genome non-infective for eBSV opens the way to the introduction of new genitors in programmes of genetic banana improvement.     

The landscape of transposable elements in the finished genome of the fungal wheat pathogen Mycosphaerella graminicola

Background

In addition to gene identification and annotation, repetitive sequence analysis has become an integral part of genome sequencing projects. Identification of repeats is important not only because it improves gene prediction, but also because of the role that repetitive sequences play in determining the structure and evolution of genes and genomes. Several methods using different repeat-finding strategies are available for whole-genome repeat sequence analysis. Four independent approaches were used to identify and characterize the repetitive fraction of the Mycosphaerella graminicola (synonym Zymoseptoria tritici) genome. This ascomycete fungus is a wheat pathogen and its finished genome comprises 21 chromosomes, eight of which can be lost with no obvious effects on fitness so are dispensable.

Results

Using a combination of four repeat-finding methods, at least 17% of the M. graminicola genome was estimated to be repetitive. Class I transposable elements, that amplify via an RNA intermediate, account for about 70% of the total repetitive content in the M. graminicola genome. The dispensable chromosomes had a higher percentage of repetitive elements as compared to the core chromosomes. Distribution of repeats across the chromosomes also varied, with at least six chromosomes showing a non-random distribution of repetitive elements. Repeat families showed transition mutations and a CpA → TpA dinucleotide bias, indicating the presence of a repeat-induced point mutation (RIP)-like mechanism in M. graminicola. One gene family and two repeat families specific to subtelomeres also were identified in the M. graminicola genome. A total of 78 putative clusters of nested elements was found in the M. graminicola genome. Several genes with putative roles in pathogenicity were found associated with these nested repeat clusters. This analysis of the transposable element content in the finished M. graminicola genome resulted in a thorough and highly curated database of repetitive sequences.

Conclusions

This comprehensive analysis will serve as a scaffold to address additional biological questions regarding the origin and fate of transposable elements in fungi. Future analyses of the distribution of repetitive sequences in M. graminicola also will be able to provide insights into the association of repeats with genes and their potential role in gene and genome evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1132) contains supplementary material, which is available to authorized users.     

Extensive Conserved Synteny of Genes between the Karyotypes of Manduca sexta and Bombyx mori Revealed by BAC-FISH Mapping

Background

Genome sequencing projects have been completed for several species representing four highly diverged holometabolous insect orders, Diptera, Hymenoptera, Coleoptera, and Lepidoptera. The striking evolutionary diversity of insects argues a need for efficient methods to apply genome information from such models to genetically uncharacterized species. Constructing conserved synteny maps plays a crucial role in this task. Here, we demonstrate the use of fluorescence in situ hybridization with bacterial artificial chromosome probes as a powerful tool for physical mapping of genes and comparative genome analysis in Lepidoptera, which have numerous and morphologically uniform holokinetic chromosomes.

Methodology/Principal Findings

We isolated 214 clones containing 159 orthologs of well conserved single-copy genes of a sequenced lepidopteran model, the silkworm, Bombyx mori, from a BAC library of a sphingid with an unexplored genome, the tobacco hornworm, Manduca sexta. We then constructed a BAC-FISH karyotype identifying all 28 chromosomes of M. sexta by mapping 124 loci using the corresponding BAC clones. BAC probes from three M. sexta chromosomes also generated clear signals on the corresponding chromosomes of the convolvulus hawk moth, Agrius convolvuli, which belongs to the same subfamily, Sphinginae, as M. sexta.

Conclusions/Significance

Comparison of the M. sexta BAC physical map with the linkage map and genome sequence of B. mori pointed to extensive conserved synteny including conserved gene order in most chromosomes. Only a few rearrangements, including three inversions, three translocations, and two fission/fusion events were estimated to have occurred after the divergence of Bombycidae and Sphingidae. These results add to accumulating evidence for the stability of lepidopteran genomes. Generating signals on A. convolvuli chromosomes using heterologous M. sexta probes demonstrated that BAC-FISH with orthologous sequences can be used for karyotyping a wide range of related and genetically uncharacterized species, significantly extending the ability to develop synteny maps for comparative and functional genomics.     

Comparative genomics between human and animal associated subspecies of the Mycobacterium avium complex: a basis for pathogenicity

Background

A human isolate of Mycobacterium avium subsp. paratuberculosis (M. paratuberculosis 43525) was sequenced and compared genomically to other mycobacterial pathogens. M. paratuberculosis 43525 was recently isolated from a patient with ulcerative colitis and belongs to the M. avium complex, a group known to infect both humans and animals. While M. paratuberculosis is a known pathogen of livestock, there are only 20 human isolates from the last 20 years, therefore we took the opportunity to perform a whole genome comparison between human and animal mycobacterial pathogens. We also compared virulence determinants such as the mycobactin cluster, PE/PPE genes and mammalian cell entry (mce) operons between MAC subspecies that infect animals and those that infect humans. M. tuberculosis was also included in these analyses given its predominant role as a human pathogen.

Results

This genome comparison showed the PE/PPE profile of M. paratuberculosis 43525 to be largely the same as other M. paratuberculosis isolates, except that it had one PPE and one PE_PGRS protein that are only present in human MAC strains and M. tuberculosis. PE/PPE proteins that were unique to M. paratuberculosis 43525, M. avium subsp. hominissuis and a caprine M. paratuberculosis isolate, were also identified. In addition, the mycobactin cluster differed between human and animal isolates and a unique mce operon flanked by two mycobactin genes, mbtA and mbtJ, was identified in all available M. paratuberculosis genomes.

Conclusions

Despite the whole genome comparison placing M. paratuberculosis 43525 as closely related to bovine M. paratuberculosis, key virulence factors were similar to human mycobacterial pathogens. This study highlights key factors of mycobacterial pathogenesis in humans and forms the basis for future functional studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1889-2) contains supplementary material, which is available to authorized users.     

Identification of the Genomic Constitution ofMusaL. Lines (Bananas, Plantains and Hybrids) Using Molecular Cytogenetics

The genomic constitutions of someMusaL. lines (bananas, plantainsand artificial hybrids) were identified using molecular cytogenetictechniques. Double targetin situDNA:DNA hybridization to chromosomespreads using as probes, total genomic DNA isolated from diploidMusalinesof known AA (labelled with biotin-11-dUTP) and BB (labelledwith digoxigenin-11-dUTP) genome constitution was carried out.The use of 60% acetic acid combined with heating over a flamegave high quality chromosome spreads free of cytoplasm forinsituhybridization. Total genomic A DNA labelled broad centromericregions of all 22 chromosomes of the diploid line, Calcutta4 (M. acuminataColla. ssp.burmanniccoides; A genome) with somechromosomes showing stronger hybridization. Labelled DNA fromthe B genome hybridized strongly to the centromeric regionsof all 22 chromosomes of Butohan 2 (M. balbisianaColla; B genome).The two satellited chromosomes of genome B labelled stronglywith genomic A DNA.In situhybridization of labelled A and Bgenomic DNA to metaphase chromosomes of triploid AAB and ABBcultivars discriminated between A and B genome chromosomes.The plantains Agbagba, Obino l'Ewai and Mbi Egome showed 22genome A and 11 genome B chromosomes while the cooking bananasBluggoe and Fougamou showed 11 genome A and 22 genome B chromosomes.Hybridization of labelled A and B genomic DNA to chromosomesof the hybrids showed that TMP2x 2829-62 has all 22 genome Achromosomes while TMPx 4698-1 has 33 genome A and 11 genomeB chromosomes.In situhybridization of labelled total genomicDNA to chromosomes has immense potential for identificationof chromosome origin and can be used to characterize cultivarsand hybrids produced inMusabreeding.Copyright 1997 Annals ofBotany Company Genomicin situhybridization; banana; plantain; hybrids; plant breeding; genome organization; biodiversity     

Polyphasic characterization and identification of the bioremediation agent Bacillus sp. SFC 500-1E

Bacillus sp. SFC 500-1E is used for the effective treatment of tannery effluents since it consistently removes hexavalent chromium from diverse contaminated matrices. The aim of the present study was to complete identification of the strain through a polyphasic characterization, which included the pattern of carbohydrate utilization, fatty acids profile, multilocus sequence analysis, multiplex PCR profile and the analysis of the complete genome sequence. Morpho-physiological and biochemical characterization results and analysis of 16S rRNA sequences were not conclusive. The strain formed a monophyletic clade with B. toyonensis BCT-7112, B. thuringiensis MC28 and B. cereus Rock 1–3. However, genomic comparisons with type strains of B. cereus and B. thuringiensis showed that the isolated belonged to a different species. Results of this study highlight the relevance of the genome sequence of this strain, identified as Bacillus toyonensis SFC 500-1E, to expand knowledge of its bioremediation potential and to explore unknown decontamination activities.     

Genome Sequences of Populus tremula Chloroplast and Mitochondrion: Implications for Holistic Poplar Breeding

Complete Populus genome sequences are available for the nucleus (P. trichocarpa; section Tacamahaca) and for chloroplasts (seven species), but not for mitochondria. Here, we provide the complete genome sequences of the chloroplast and the mitochondrion for the clones P. tremula W52 and P. tremula x P. alba 717-1B4 (section Populus). The organization of the chloroplast genomes of both Populus clones is described. A phylogenetic tree constructed from all available complete chloroplast DNA sequences of Populus was not congruent with the assignment of the related species to different Populus sections. In total, 3,024 variable nucleotide positions were identified among all compared Populus chloroplast DNA sequences. The 5-prime part of the LSC from trnH to atpA showed the highest frequency of variations. The variable positions included 163 positions with SNPs allowing for differentiating the two clones with P. tremula chloroplast genomes (W52, 717-1B4) from the other seven Populus individuals. These potential P. tremula-specific SNPs were displayed as a whole-plastome barcode on the P. tremula W52 chloroplast DNA sequence. Three of these SNPs and one InDel in the trnH-psbA linker were successfully validated by Sanger sequencing in an extended set of Populus individuals. The complete mitochondrial genome sequence of P. tremula is the first in the family of Salicaceae. The mitochondrial genomes of the two clones are 783,442 bp (W52) and 783,513 bp (717-1B4) in size, structurally very similar and organized as single circles. DNA sequence regions with high similarity to the W52 chloroplast sequence account for about 2% of the W52 mitochondrial genome. The mean SNP frequency was found to be nearly six fold higher in the chloroplast than in the mitochondrial genome when comparing 717-1B4 with W52. The availability of the genomic information of all three DNA-containing cell organelles will allow a holistic approach in poplar molecular breeding in the future.     

Genome-wide comparative analysis of the Brassica rapa gene space reveals genome shrinkage and differential loss of duplicated genes after whole genome triplication

Background

Brassica rapa is one of the most economically important vegetable crops worldwide. Owing to its agronomic importance and phylogenetic position, B. rapa provides a crucial reference to understand polyploidy-related crop genome evolution. The high degree of sequence identity and remarkably conserved genome structure between Arabidopsis and Brassica genomes enables comparative tiling sequencing using Arabidopsis sequences as references to select the counterpart regions in B. rapa, which is a strong challenge of structural and comparative crop genomics.

Results

We assembled 65.8 megabase-pairs of non-redundant euchromatic sequence of B. rapa and compared this sequence to the Arabidopsis genome to investigate chromosomal relationships, macrosynteny blocks, and microsynteny within blocks. The triplicated B. rapa genome contains only approximately twice the number of genes as in Arabidopsis because of genome shrinkage. Genome comparisons suggest that B. rapa has a distinct organization of ancestral genome blocks as a result of recent whole genome triplication followed by a unique diploidization process. A lack of the most recent whole genome duplication (3R) event in the B. rapa genome, atypical of other Brassica genomes, may account for the emergence of B. rapa from the Brassica progenitor around 8 million years ago.

Conclusions

This work demonstrates the potential of using comparative tiling sequencing for genome analysis of crop species. Based on a comparative analysis of the B. rapa sequences and the Arabidopsis genome, it appears that polyploidy and chromosomal diploidization are ongoing processes that collectively stabilize the B. rapa genome and facilitate its evolution.     

Identification of a novel cytoplasmic male sterile line M2BS induced by partial-length <Emphasis Type="Italic">HcPDIL5-2a</Emphasis> transformation in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

A cytoplasmic male sterile line (designated as M2BS) was obtained from an indica rice maintainer M2B induced by partial-length HcPDIL5-2a (Hibiscus cannabinus protein disulfide isomerase-like) transformation. The anther of M2BS was short, slender, hygrophanous, and fissured. I2-KI staining method showed that there was typical and spherical abortion in pollen grains. M2BS was found abortive at middle and late stage of monocyte by the modified carbol fuchsin stained observation and paraffin section observation. The tapetum was observed pre-degenerated in M2BS. Hereditary analysis indicated that the male sterility of M2BS was a maternally inherited inability after six backcross generations with M2B and the combinations of M2BS hybridized with other two male fertile materials. The M2BS could be affirmed a cytoplasmic male sterile (CMS) type. Moreover, it was a transgenic plant confirmed by PCR, Southern blot and RT-PCR detection. M2BS could be distinguished from M2B and its CMS line M2A by RFLP analysis. The overall mitochondrial genome sequencing results showed, that in M2BS, the main differences of mitochondrial gene sequence were located in nad4, nad5, nad7, orf194 and intergenic region, relatively to those of M2A. The obtained results indicate that M2BS is a novel cytoplasmic male sterile line.     

Comparative Analysis between Homoeologous Genome Segments of Brassica napus and Its Progenitor Species Reveals Extensive Sequence-Level Divergence

Homoeologous regions of Brassica genomes were analyzed at the sequence level. These represent segments of the Brassica A genome as found in Brassica rapa and Brassica napus and the corresponding segments of the Brassica C genome as found in Brassica oleracea and B. napus. Analysis of synonymous base substitution rates within modeled genes revealed a relatively broad range of times (0.12 to 1.37 million years ago) since the divergence of orthologous genome segments as represented in B. napus and the diploid species. Similar, and consistent, ranges were also identified for single nucleotide polymorphism and insertion-deletion variation. Genes conserved across the Brassica genomes and the homoeologous segments of the genome of Arabidopsis thaliana showed almost perfect collinearity. Numerous examples of apparent transduplication of gene fragments, as previously reported in B. oleracea, were observed in B. rapa and B. napus, indicating that this phenomenon is widespread in Brassica species. In the majority of the regions studied, the C genome segments were expanded in size relative to their A genome counterparts. The considerable variation that we observed, even between the different versions of the same Brassica genome, for gene fragments and annotated putative genes suggest that the concept of the pan-genome might be particularly appropriate when considering Brassica genomes.     

Syntenic Relationships between the U and M Genomes of Aegilops,Wheat and the Model Species Brachypodium and Rice as Revealed by COS Markers

Diploid Aegilops umbellulata and Ae. comosa and their natural allotetraploid hybrids Ae. biuncialis and Ae. geniculata are important wild gene sources for wheat. With the aim of assisting in alien gene transfer, this study provides gene-based conserved orthologous set (COS) markers for the U and M genome chromosomes. Out of the 140 markers tested on a series of wheat-Aegilops chromosome introgression lines and flow-sorted subgenomic chromosome fractions, 100 were assigned to Aegilops chromosomes and six and seven duplications were identified in the U and M genomes, respectively. The marker-specific EST sequences were BLAST-ed to Brachypodium and rice genomic sequences to investigate macrosyntenic relationships between the U and M genomes of Aegilops, wheat and the model species. Five syntenic regions of Brachypodium identified genome rearrangements differentiating the U genome from the M genome and from the D genome of wheat. All of them seem to have evolved at the diploid level and to have been modified differentially in the polyploid species Ae. biuncialis and Ae. geniculata. A certain level of wheat–Aegilops homology was detected for group 1, 2, 3 and 5 chromosomes, while a clearly rearranged structure was showed for the group 4, 6 and 7 Aegilops chromosomes relative to wheat. The conserved orthologous set markers assigned to Aegilops chromosomes promise to accelerate gene introgression by facilitating the identification of alien chromatin. The syntenic relationships between the Aegilops species, wheat and model species will facilitate the targeted development of new markers specific for U and M genomic regions and will contribute to the understanding of molecular processes related to allopolyploidization.     

Molecular cytogenetic and genomic analyses reveal new insights into the origin of the wheat B genome

Key message

This work pinpointed the goatgrass chromosomal segment in the wheat B genome using modern cytogenetic and genomic technologies, and provided novel insights into the origin of the wheat B genome.

Abstract

Wheat is a typical allopolyploid with three homoeologous subgenomes (A, B, and D). The donors of the subgenomes A and D had been identified, but not for the subgenome B. The goatgrass Aegilops speltoides (genome SS) has been controversially considered a possible candidate for the donor of the wheat B genome. However, the relationship of the Ae. speltoides S genome with the wheat B genome remains largely obscure. The present study assessed the homology of the B and S genomes using an integrative cytogenetic and genomic approach, and revealed the contribution of Ae. speltoides to the origin of the wheat B genome. We discovered noticeable homology between wheat chromosome 1B and Ae. speltoides chromosome 1S, but not between other chromosomes in the B and S genomes. An Ae. speltoides-originated segment spanning a genomic region of approximately 10.46 Mb was detected on the long arm of wheat chromosome 1B (1BL). The Ae. speltoides-originated segment on 1BL was found to co-evolve with the rest of the B genome. Evidently, Ae. speltoides had been involved in the origin of the wheat B genome, but should not be considered an exclusive donor of this genome. The wheat B genome might have a polyphyletic origin with multiple ancestors involved, including Ae. speltoides. These novel findings will facilitate genome studies in wheat and other polyploids.

     

Complete chloroplast genome sequence of Magnolia grandiflora and comparative analysis with related species

Magnolia grandiflora is an important medicinal,ornamental and horticultural plant species.The chloroplast(cp) genome of M.grandiflora was sequenced using a 454 sequencing platform and the genome structure was compared with other related species.The complete cp genome of M.grandiflora was 159623 bp in length and contained a pair of inverted repeats(IR) of 26563 bp separated by large and small single copy(LSC,SSC) regions of 87757 and 18740 bp,respectively.A total of 129 genes were successfully annotated,18 of which included introns.The identity,number and GC content of M.grandiflora cp genes were similar to those of other Magnoliaceae species genomes.Analysis revealed 218 simple sequence repeat(SSR) loci,most composed of A or T,contributing to a bias in base composition.The types and abundances of repeat units in Magnoliaceae species were relatively conserved and these loci will be useful for developing M.grandiflora cp genome vectors.In addition,results indicated that the cp genome size in Magnoliaceae species and the position of the IR border were closely related to the length of the ycf1 gene.Phylogenetic analyses based on 66 shared genes from 30 species using maximum parsimony(MP) and maximum likelihood(ML) methods provided strong support for the phylogenetic position of Magnolia.The availability of the complete cp genome sequence of M.grandiflora provides valuable information for breeding of desirable varieties,cp genetic engineering,developing useful molecular markers and phylogenetic analyses in Magnoliaceae.     

The mitochondrial genome contribution to the phylogeny and identification of Metarhizium species and strains

The genus Metarhizium is composed of entomopathogenic fungal biological control agents (BCAs) used for invertebrate pest control. The phylogenetic relationships of species within this genus are still under scrutiny as several cryptic species can be found. In this work, the mitochondrial (mt) genome of Metarhizium brunneum ARSEF 4556 was fully sequenced and a comparative genome analysis was conducted with 7 other available mt genomes, belonging to 5 Metarhizium species: M. anisopliae, M. brunneum, M. robertsii, M. guizhouense and M. majus. Results showed that Metarhizium demonstrates greater conserved stability than other fungal mt genomes. Furthermore, this analysis located 7 diverse regions in both intergenic domains and gene fragments which were ideal for species/strain discrimination. The sequencing of these regions revealed several SNPs among 38 strains tested, 11 of which were uncharacterized. Single gene phylogenies presented variable results which may be used further for intra-species discrimination. Phylogenetic trees based on the concatenation of mt domains and the nuclear ITS1-5.8S-ITS2 region showed discrimination of the species studied and allowed the identification of uncharacterized strains. These were mostly placed within species M. anisopliae and M. brunneum. Five strains clustered together in a clade related to M. brunneum, suggesting that they comprise a cryptic species.     

Full genome SNP-based phylogenetic analysis reveals the origin and global spread of Brucella melitensis

Background

Brucellosis is an important zoonotic disease that affects both humans and animals. We sequenced the full genome and characterised the genetic diversity of two Brucella melitensis isolates from Malaysia and the Philippines. In addition, we performed a comparative whole-genome single nucleotide polymorphism (SNP) analysis of B. melitensis strains collected from around the world, to investigate the potential origin and the history of the global spread of B. melitensis.

Results

Single sequencing runs of each genome resulted in draft genome sequences of MY1483/09 and Phil1136/12, which covered 99.85% and 99.92% of the complete genome sequences, respectively. The B. melitensis genome sequences, and two B. abortus strains used as the outgroup strains, yielded a total of 13,728 SNP sites. Phylogenetic analysis using whole-genome SNPs and geographical distribution of the isolates revealed spatial clustering of the B. melitensis isolates into five genotypes, I, II, III, IV and V. The Mediterranean strains, identified as genotype I, occupied the basal node of the phylogenetic tree, suggesting that B. melitensis may have originated from the Mediterranean regions. All of the Asian B. melitensis strains clustered into genotype II with the SEA strains, including the two isolates sequenced in this study, forming a distinct clade denoted here as genotype IId. Genotypes III, IV and V of B. melitensis demonstrated a restricted geographical distribution, with genotype III representing the African lineage, genotype IV representing the European lineage and genotype V representing the American lineage.

Conclusion

We showed that SNPs retrieved from the B. melitensis draft full genomes were sufficient to resolve the interspecies relationships between B. melitensis strains and to discriminate between the vaccine and endemic strains. Phylogeographic reconstruction of the history of B. melitensis global spread at a finer scale by using whole-genome SNP analyses supported the origin of all B. melitensis strains from the Mediterranean region. The possible global distribution of B. melitensis following the ancient trade routes was also consistent with whole-genome SNP phylogeny. The whole genome SNP phylogenetics analysis, hence is a powerful tool for intraspecies discrimination of closely related species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1294-x) contains supplementary material, which is available to authorized users.     

Identification and characterization of improved nitrogen efficiency in interspecific hybridized new-type Brassica napus

Background and Aims

Oilseed rape (Brassica napus) is an important oil crop worldwide. The aim of this study was to identify the variation in nitrogen (N) efficiency of new-type B. napus (genome A^rA^rC^cC^c) genotypes, and to characterize some critical physiological and molecular mechanisms in response to N limitation.

Methods

Two genotypes with contrasting N efficiency (D4-15 and D1-1) were identified from 150 new-type B. napus lines, and hydroponic and pot experiments were conducted. Root morphology, plant biomass, N uptake parameters and seed yield of D4-15 and D1-1 were investigated. Two traditional B. napus (genome AⁿAⁿCⁿCⁿ) genotypes, QY10 and NY7, were also cultivated. Introgression of exotic genomic components in D4-15 and D1-1 was evaluated with molecular markers.

Key Results

Large genetic variation existed among traits contributing to the N efficiency of new-type B. napus. Under low N levels at the seedling stage, the N-efficient new-type D4-15 showed higher values than the N-inefficient D1-1 line and the traditional B. napus QY10 and NY7 genotypes with respect to several traits, including root and shoot biomass, root morphology, N accumulation, N utilization efficiency (NutE), N uptake efficiency (NupE), activities of nitrate reductase (NR) and glutamine synthetase (GS), and expression levels of N transporter genes and genes that are involved in N assimilation. Higher yield was produced by the N-efficient D4-15 line compared with the N-inefficient D1-1 at maturity. More exotic genome components were introgressed into the genome of D4-15 (64·97 %) compared with D1-1 (32·23 %).

Conclusions

The N-efficient new-type B. napus identified in this research had higher N efficiency (and tolerance to low-N stress) than traditional B. napus cultivars, and thus could have important potential for use in breeding N-efficient B. napus cultivars in the field.