Unraveling the nuclear and chloroplast genomes of an agar producing red macroalga, Gracilaria changii (Rhodophyta,Gracilariales)

Agar and agarose have wide applications in food and pharmaceutical industries. Knowledge on the genome of red seaweeds that produce them is still lacking. To fill the gap in genome analyses of these red algae, we have sequenced the nuclear and organellar genomes of an agarophyte, Gracilaria changii. The partial nuclear genome sequence of G. changii has a total length of 35.8 Mb with 10,912 predicted protein coding sequences. Only 39.4% predicted proteins were found to have significant matches to protein sequences in SwissProt. The chloroplast genome of G. changii is 183,855 bp with a total of 201 open reading frames (ORFs), 29 tRNAs and 3 rRNAs predicted. Five genes: ssrA, leuC and leuD CP76_p173 (orf139) and pbsA were absent in the chloroplast genome of G. changii. The genome information is valuable in accelerating functional studies of individual genes and resolving evolutionary relationship of red seaweeds.     

EST analysis of Ostreococcus lucimarinus, the most compact eukaryotic genome, shows an excess of introns in highly expressed genes

Background

The genome of the pico-eukaryotic (bacterial-sized) prasinophyte green alga Ostreococcus lucimarinus has one of the highest gene densities known in eukaryotes, yet it contains many introns. Phylogenetic studies suggest this unusually compact genome (13.2 Mb) is an evolutionarily derived state among prasinophytes. The presence of introns in the highly reduced O. lucimarinus genome appears to be in opposition to simple explanations of genome evolution based on unidirectional tendencies, either neutral or selective. Therefore, patterns of intron retention in this species can potentially provide insights into the forces governing intron evolution.

Methodology/Principal Findings

Here we studied intron features and levels of expression in O. lucimarinus using expressed sequence tags (ESTs) to annotate the current genome assembly. ESTs were assembled into unigene clusters that were mapped back to the O. lucimarinus Build 2.0 assembly using BLAST and the level of gene expression was inferred from the number of ESTs in each cluster. We find a positive correlation between expression levels and both intron number (R = +0.0893, p = <0.0005) and intron density (number of introns/kb of CDS; R = +0.0753, p = <0.005).

Conclusions/Significance

In a species with a genome that has been recently subjected to a great reduction of non-coding DNA, these results imply the existence of selective/functional roles for introns that are principally detectable in highly expressed genes. In these cases, introns are likely maintained by balancing the selective forces favoring their maintenance with other mutational and/or selective forces acting on genome size.     

Genetic Diversity,Heteroplasmy, and Recombination in Mitochondrial Genomes of Daphnia pulex,Daphnia pulicaria,and Daphnia obtusa

Genetic variants of mitochondrial DNA at the individual (heteroplasmy) and population (polymorphism) levels provide insight into their roles in multiple cellular and evolutionary processes. However, owing to the paucity of genome-wide data at the within-individual and population levels, the broad patterns of these two forms of variation remain poorly understood. Here, we analyze 1,804 complete mitochondrial genome sequences from Daphnia pulex, Daphnia pulicaria, and Daphnia obtusa. Extensive heteroplasmy is observed in D. obtusa, where the high level of intraclonal divergence must have resulted from a biparental-inheritance event, and recombination in the mitochondrial genome is apparent, although perhaps not widespread. Global samples of D. pulex reveal remarkably low mitochondrial effective population sizes, <3% of those for the nuclear genome. In addition, levels of population diversity in mitochondrial and nuclear genomes are uncorrelated across populations, suggesting an idiosyncratic evolutionary history of mitochondria in D. pulex. These population-genetic features appear to be a consequence of background selection associated with highly deleterious mutations arising in the strongly linked mitochondrial genome, which is consistent with polymorphism and divergence data suggesting a predominance of strong purifying selection. Nonetheless, the fixation of mildly deleterious mutations in the mitochondrial genome also appears to be driving positive selection on genes encoded in the nuclear genome whose products are deployed in the mitochondrion.     

Genome size estimation of an endoparasitoid wasp,Cotesia plutellae,using quantitative real-time polymerase chain reaction

An endoparasitoid wasp, Cotesia plutellae, is a natural enemy against the diamondback moth, Plutella xylostella, which is the most destructive insect pest of cruciferous crop plants. The wasp genome contains genetic information of several parasitic factors, such as its symbiotic virus (C. plutellae bracovirus), venom, teratocyte as well as the parasitoid itself. These parasitic factors interfere with physiological processes of the immature stages of P. xylostella and need to be analyzed concerning their genetic components. A full genome sequence of C. plutellae would be highly informative to determine functional genes associated with these parasitic factors. Before a full genome sequence analysis of C. plutellae can be undertaken, an estimate of genome size is needed. In this study, we used a strategy using a quantitative real-time polymerase chain reaction (qPCR) to measure the wasp genome size. qPCR determined the number of a single copy gene hexokinase in a total mass of genomic DNA. The resulting molecular weight of the DNA sample was used to calculate the genome size in base pair (bp). This estimation approach indicated a genome size of C. plutellae of 186.06 ± 1.21 Mb.     

The genome sequence of Xanthomonas oryzae pathovar oryzae KACC10331, the bacterial blight pathogen of rice

The nucleotide sequence was determined for the genome of Xanthomonas oryzae pathovar oryzae (Xoo) KACC10331, a bacterium that causes bacterial blight in rice (Oryza sativa L.). The genome is comprised of a single, 4 941 439 bp, circular chromosome that is G + C rich (63.7%). The genome includes 4637 open reading frames (ORFs) of which 3340 (72.0%) could be assigned putative function. Orthologs for 80% of the predicted Xoo genes were found in the previously reported X.axonopodis pv. citri (Xac) and X.campestris pv. campestris (Xcc) genomes, but 245 genes apparently specific to Xoo were identified. Xoo genes likely to be associated with pathogenesis include eight with similarity to Xanthomonas avirulence (avr) genes, a set of hypersensitive reaction and pathogenicity (hrp) genes, genes for exopolysaccharide production, and genes encoding extracellular plant cell wall-degrading enzymes. The presence of these genes provides insights into the interactions of this pathogen with its gramineous host.     

One Bacterial Cell,One Complete Genome

While the bulk of the finished microbial genomes sequenced to date are derived from cultured bacterial and archaeal representatives, the vast majority of microorganisms elude current culturing attempts, severely limiting the ability to recover complete or even partial genomes from these environmental species. Single cell genomics is a novel culture-independent approach, which enables access to the genetic material of an individual cell. No single cell genome has to our knowledge been closed and finished to date. Here we report the completed genome from an uncultured single cell of Candidatus Sulcia muelleri DMIN. Digital PCR on single symbiont cells isolated from the bacteriome of the green sharpshooter Draeculacephala minerva bacteriome allowed us to assess that this bacteria is polyploid with genome copies ranging from approximately 200–900 per cell, making it a most suitable target for single cell finishing efforts. For single cell shotgun sequencing, an individual Sulcia cell was isolated and whole genome amplified by multiple displacement amplification (MDA). Sanger-based finishing methods allowed us to close the genome. To verify the correctness of our single cell genome and exclude MDA-derived artifacts, we independently shotgun sequenced and assembled the Sulcia genome from pooled bacteriomes using a metagenomic approach, yielding a nearly identical genome. Four variations we detected appear to be genuine biological differences between the two samples. Comparison of the single cell genome with bacteriome metagenomic sequence data detected two single nucleotide polymorphisms (SNPs), indicating extremely low genetic diversity within a Sulcia population. This study demonstrates the power of single cell genomics to generate a complete, high quality, non-composite reference genome within an environmental sample, which can be used for population genetic analyzes.     

The Chlorella variabilis NC64A Genome Reveals Adaptation to Photosymbiosis,Coevolution with Viruses,and Cryptic Sex

Chlorella variabilis NC64A, a unicellular photosynthetic green alga (Trebouxiophyceae), is an intracellular photobiont of Paramecium bursaria and a model system for studying virus/algal interactions. We sequenced its 46-Mb nuclear genome, revealing an expansion of protein families that could have participated in adaptation to symbiosis. NC64A exhibits variations in GC content across its genome that correlate with global expression level, average intron size, and codon usage bias. Although Chlorella species have been assumed to be asexual and nonmotile, the NC64A genome encodes all the known meiosis-specific proteins and a subset of proteins found in flagella. We hypothesize that Chlorella might have retained a flagella-derived structure that could be involved in sexual reproduction. Furthermore, a survey of phytohormone pathways in chlorophyte algae identified algal orthologs of Arabidopsis thaliana genes involved in hormone biosynthesis and signaling, suggesting that these functions were established prior to the evolution of land plants. We show that the ability of Chlorella to produce chitinous cell walls likely resulted from the capture of metabolic genes by horizontal gene transfer from algal viruses, prokaryotes, or fungi. Analysis of the NC64A genome substantially advances our understanding of the green lineage evolution, including the genomic interplay with viruses and symbiosis between eukaryotes.     

Down,then up: non-parallel genome size changes and a descending chromosome series in a recent radiation of the Australian allotetraploid plant species,Nicotiana section Suaveolentes (Solanaceae)

Background and AimsThe extent to which genome size and chromosome numbers evolve in concert is little understood, particularly after polyploidy (whole-genome duplication), when a genome returns to a diploid-like condition (diploidization). We study this phenomenon in 46 species of allotetraploid Nicotiana section Suaveolentes (Solanaceae), which formed <6 million years ago and radiated in the arid centre of Australia.MethodsWe analysed newly assessed genome sizes and chromosome numbers within the context of a restriction site-associated nuclear DNA (RADseq) phylogenetic framework.Key ResultsRADseq generated a well-supported phylogenetic tree, in which multiple accessions from each species formed unique genetic clusters. Chromosome numbers and genome sizes vary from n = 2x = 15 to 24 and 2.7 to 5.8 pg/1C nucleus, respectively. Decreases in both genome size and chromosome number occur, although neither consistently nor in parallel. Species with the lowest chromosome numbers (n = 15–18) do not possess the smallest genome sizes and, although N. heterantha has retained the ancestral chromosome complement, n = 2x = 24, it nonetheless has the smallest genome size, even smaller than that of the modern representatives of ancestral diploids.ConclusionsThe results indicate that decreases in genome size and chromosome number occur in parallel down to a chromosome number threshold, n = 20, below which genome size increases, a phenomenon potentially explained by decreasing rates of recombination over fewer chromosomes. We hypothesize that, more generally in plants, major decreases in genome size post-polyploidization take place while chromosome numbers are still high because in these stages elimination of retrotransposons and other repetitive elements is more efficient. Once such major genome size change has been accomplished, then dysploid chromosome reductions take place to reorganize these smaller genomes, producing species with small genomes and low chromosome numbers such as those observed in many annual angiosperms, including Arabidopsis.     

Evolutionary history of mitogen-activated protein kinase (MAPK) genes in Lotus,Medicago, and Phaseolus

Mitogen-Activated Protein Kinase (MAPK) genes encode proteins that mediate various signaling pathways associated with biotic and abiotic stress responses in eukaryotes. The MAPK genes form a 3-tier signal transduction cascade between cellular stimuli and physiological responses. Recent identification of soybean MAPKs and availability of genome sequences from other legume species allowed us to identify their MAPK genes. The main objectives of this study were to identify MAPKs in 3 legume species, Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, and to assess their phylogenetic relationships. We used approaches in comparative genomics for MAPK gene identification and named the newly identified genes following Arabidopsis MAPK nomenclature model. We identified 19, 18, and 15 MAPKs and 7, 4, and 9 MAPKKs in the genome of Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, respectively. Within clade placement of MAPKs and MAPKKs in the 3 legume species were consistent with those in soybean and Arabidopsis. Among 5 clades of MAPKs, 4 founder clades were consistent to MAPKs of other plant species and orthologs of MAPK genes in the fifth clade-"Clade E" were consistent with those in soybean. Our results also indicated that some gene duplication events might have occurred prior to eudicot-monocot divergence. Highly diversified MAPKs in soybean relative to those in 3 other legume species are attributable to the polyploidization events in soybean. The identification of the MAPK genes in the legume species is important for the legume crop improvement; and evolutionary relationships and functional divergence of these gene members provide insights into plant genome evolution.     

Genome sequence and comparative analysis of DRQ-2, the type strain of Nonomuraea indica

Strain DRQ-2^T (type strain of Nonomuraea indica) is worthy for genome sequencing, due to its ability to produce a wide variety of industrially important enzymes such as amylase, asparaginase, cellulase, gelatinase, glutaminase, and protease. Genome sequencing and comparison of strain DRQ-2^T is described in the present work. The genome size was estimated to be 8,288,417 (bp) that consisted of 59 contigs. The G + C content of the genome was 72.4%. A total of 7730 genes were predicted with two rRNAs and 64 tRNAs. The genome analysis of the strain DRQ-2^T showed the presence of a wide range of secondary metabolite gene clusters. Pan-Genomes Analysis Pipeline (PGAP) indicated that strain DRQ-2^T had large numbers of unique genes. The majority of N. indica DRQ-2^T genes encode for hypothetical proteins, indicating the functions of these ortholog clusters were still remain to be determined.     

DcSto: carrot Stowaway-like elements are abundant, diverse, and polymorphic

We investigated nine families of Stowaway-like miniature inverted-repeat transposable elements (MITEs) in the carrot genome, named DcSto1 to DcSto9. All of them were AT-rich and shared a highly conserved 6 bp-long TIR typical for Stowaways. The copy number of DcSto1 elements was estimated as ca. 5,000 per diploid genome. We observed preference for clustered insertions of DcSto and other MITEs. Distribution of DcSto1 hybridization signals revealed presence of DcSto1 clusters within euchromatic regions along all chromosomes. An arrangement of eight regions encompassing DcSto insertion sites, studied in detail, was highly variable among plants representing different populations of Daucus carota. All of these insertions were polymorphic which most likely suggests a very recent mobilization of those elements. Insertions of DcSto near carrot genes and presence of putative promoters, regulatory motifs, and polyA signals within their sequences might suggest a possible involvement of DcSto in the regulation of gene expression.     

Complete chloroplast genome sequences of Dipterygium glaucum and Cleome chrysantha and other Cleomaceae Species,comparative analysis and phylogenetic relationships

This current study presents, for the first time, the complete chloroplast genome of two Cleomaceae species: Dipterygium glaucum and Cleome chrysantha in order to evaluate the evolutionary relationship. The cp genome is 158,576 bp in length with 35.74% GC content in D. glaucum and 158,111 bp with 35.96% GC in C. chrysantha. Inverted repeats IR 26,209 bp, 26,251 bp each, LSC of 87,738 bp, 87,184 bp and SSC of 18,420 bp, 18,425 bp respectively. There are 136 genes in the genome, which includes 80 protein coding genes, 31 tRNA genes and four rRNA genes were observed in both chloroplast genomes. 117 genes are unique while the remaining 19 genes are duplicated in IR regions. The analysis of repeats shows that the cp genome includes all types of repeats with more frequent occurrences of palindromic; Also, this analysis indicates that the total number of simple sequence repeats (SSR) were 323 in D. glaucum, and 313 in C. chrysantha, of which the majority of the SSRs in these plastid genomes were mononucleotide repeats A/T which are located in the intergenic spacer. Moreover, the comparative analysis of the four cp sequences revealed four hotspot genes (atpF, rpoC2, rps19, and ycf1), these variable regions could be used as molecular makers for the species authentication as well as resources for inferring phylogenetic relationships of the species. All the relationships in the phylogenetic tree are with high support, this indicate that the complete chloroplast genome is a useful data for inferring phylogenetic relationship within the Cleomaceae and other families. The simple sequence repeats identified will be useful for identification, genetic diversity, and other evolutionary studies of the species. This study reported the first cp genome of the genus Dipterygium and Cleome. The finding of this study will be beneficial for biological disciplines such as evolutionary and genetic diversity studies of the species within the core Cleomaceae.     

ISD1, an Insertion Element from the Sulfate-Reducing Bacterium Desulfovibrio vulgaris Hildenborough: Structure, Transposition, and Distribution

Insertion element ISD1, discovered when its transposition caused the insertional inactivation of an introduced sacB gene, is present in two copies in the genome of Desulfovibrio vulgaris Hildenborough. Southern blot analysis indicated at least two insertion sites in the sacB gene. Cloning and sequencing of a transposed copy of ISD1 indicated a length of 1,200 bp with a pair of 44-bp imperfect inverted repeats at the ends, flanked by a direct repeat of the 4-bp target sequence. AAGG and AATT were found to function as target sequences. ISD1 encodes a transposase from two overlapping open reading frames by programmed translational frameshifting at an A₆G shifty codon motif. Sequence comparison showed that ISD1 belongs to the IS3 family. Isolation and analysis of the chromosomal copies, ISD1-A and ISD1-B, by PCR and sequencing indicated that these are not flanked by direct repeats. ISD1-A is inserted in a region of the chromosome containing the gapdh-pgk genes (encoding glyceraldehyde-3-phosphate dehydrogenase and phosphoglycerate kinase). Active transposition to other loci in the genome was demonstrated, offering the potential of a new tool for gene cloning and mutagenesis. ISD1 is the first transposable element described for the sulfate reducers, a large and environmentally important group of bacteria. The distribution of ISD1 in genomes of sulfate-reducing bacteria is limited. A single copy is present in the genome of D. desulfuricans Norway.     

Genomic organization and genomic structural rearrangements of Sphingobium japonicum UT26, an archetypal γ-hexachlorocyclohexane-degrading bacterium

The complete genome sequencing of a γ-hexachlorocyclohexane-degrading strain, Sphingobium japonicum UT26, revealed that the genome consists of two circular chromosomes [with sizes of 3.5 Mb (Chr1) and 682 kb (Chr2)], a 191-kb large plasmid (pCHQ1), and two small plasmids with sizes of 32 and 5 kb. The lin genes are dispersed on Chr1, Chr2, and pCHQ1. Comparison of the UT26 genome with those of other sphingomonad strains demonstrated that the “specific” lin genes for conversion of γ-HCH to β-ketoadipate (linA, linB, linC, linRED, and linF) are located on the DNA regions unique to the UT26 genome, suggesting the acquisition of these lin genes by horizontal transfer events. On the other hand, linGHIJ and linKLMN are located on the regions conserved in the genomes of sphingomonads, suggesting that the linGHIJ-encoded β-ketoadipate pathway and the LinKLMN-type ABC transporter system are involved in core functions of sphingomonads. Based on these results, we propose a hypothesis that UT26 was created by recruiting the specific lin genes into a strain having core functions of sphingomonads. Most of the specific lin genes in UT26 are associated with IS6100. Our analysis of spontaneous linA-, linC-, and linRED-deletion mutants of UT26 revealed the involvement of IS6100 in their deduced genome rearrangements. These facts strongly suggest that IS6100 plays important roles both in the dissemination of the specific lin genes and in the genome rearrangements.     

Complete chloroplast genome of the medicinal plant Evolvulus alsinoides: comparative analysis,identification of mutational hotspots and evolutionary dynamics with species of Solanales

Evolvulus alsinoides, belonging to the family Convolvulaceae, is an important medicinal plant widely used as a nootropic in the Indian traditional medicine system. In the genus Evolvulus, no research on the chloroplast genome has been published. Hence, the present study focuses on annotation, characterization, identification of mutational hotspots, and phylogenetic analysis in the complete chloroplast genome (cp) of E. alsinoides. Genome comparison and evolutionary dynamics were performed with the species of Solanales. The cp genome has 114 genes (80 protein-coding genes, 30 transfer RNA, and 4 ribosomal RNA genes) that were unique with total genome size of 157,015 bp. The cp genome possesses 69 RNA editing sites and 44 simple sequence repeats (SSRs). Predicted SSRs were randomly selected and validated experimentally. Six divergent hotspots such as trnQ-UUG, trnF-GAA, psaI, clpP, ndhF, and ycf1 were discovered from the cp genome. These microsatellites and divergent hot spot sequences of the Taxa ‘Evolvulus’ could be employed as molecular markers for species identification and genetic divergence investigations. The LSC area was found to be more conserved than the SSC and IR region in genome comparison. The IR contraction and expansion studies show that nine genes rpl2, rpl23, ycf1, ycf2, ycf1, ndhF, ndhA, matK, and psbK were present in the IR-LSC and IR-SSC boundaries of the cp genome. Fifty-four protein-coding genes in the cp genome were under negative selection pressure, indicating that they were well conserved and were undergoing purifying selection. The phylogenetic analysis reveals that E. alsinoides is closely related to the genus Cressa with some divergence from the genus Ipomoea. This is the first time the chloroplast genome of the genus Evolvulus has been published. The findings of the present study and chloroplast genome data could be a valuable resource for future studies in population genetics, genetic diversity, and evolutionary relationship of the family Convolvulaceae.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12298-021-01051-w.