Comparative genomic analysis of pyrene-degrading <Emphasis Type="Italic">Mycobacterium</Emphasis> species: Genomic islands and ring-hydroxylating dioxygenases involved in pyrene degradation

The genome sequences of two pyrene-degrading bacterial strains of Mycobacterium spp. PYR10 and PYR15, isolated from the estuarine wetland of the Han river, South Korea, were determined using the PacBio RS II sequencing platform. The complete genome of strain PYR15 was 6,037,017 bp in length with a GC content of 66.5%, and contained 5,933 protein-coding genes. The genome of strain PYR10 was 5,999,427 bp in length with a GC content of 67.7%, and contained 5,767 protein-coding genes. Based on the average nucleotide identity values, these strains were designated as M. gilvum PYR10 and M. pallens PYR15. A genomic comparison of these pyrene-degrading Mycobacterium strains with pyrene-non-degrading strains revealed that the genomes of pyrene-degrading strains possessed similar repertoires of ringhydroxylating dioxygenases (RHDs), including the pyrenehydroxylating dioxygenases encoded by nidA and nidA3, which could be readily distinguished from those of pyrenenon-degraders. Furthermore, genomic islands, containing catabolic gene clusters, were shared only among the pyrenedegrading Mycobacterium strains and these gene clusters contained RHD genes, including nidAB and nidA3B3. Our genome data should facilitate further studies on the evolution of the polycyclic aromatic hydrocarbon-degradation pathways in the genus Mycobacterium.     

The complete chloroplast genome sequence and phylogenetic analysis of <Emphasis Type="Italic">Chuanminshen</Emphasis> (<Emphasis Type="Italic">Chuanminshenviolaceum</Emphasis> Sheh et Shan)

Chloroplast genome sequences are very useful for species identification and phylogenetics. Chuanminshen (Chuanminshen violaceum Sheh et Shan) is an important traditional Chinese medicinal plant, for which the phylogenetic position is still controversial. In this study, the complete chloroplast genome of Chuanminshen violaceum Sheh et Shan was determined. The total size of Chuanminshen chloroplast genome was 154,529 bp with 37.8% GC content. It has the typical quadripartite structure, a large single copy (17,800 bp) and a small single copy (84,171 bp) and a pair of inverted repeats (26,279 bp). The whole genome harbors 132 genes, which includes 85 protein coding genes, 37 tRNA genes, eight rRNA genes, and two pseudogenes. Thirty-nine SSR loci, 32 tandem repeats and 49 dispersed repeats were found. Phylogenetic analyses results with the help of MEGA showed a new insight for the Chuanminshen phylogenetic relationship with the reported chloroplast genomes in Apiales plants.     

Molecular variation of the testes-specific β<Emphasis Type="Italic">NACtes</Emphasis> genes in the <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> genome

The βNACtes gene family of the Drosophila melanogaster genome provides a model for investigating the mechanisms of the molecular evolution of recently evolved genes. The βNACtes genes code for proteins that are homologous to the subunit of the nascent polypeptide-associated complex (NAC), are expressed exclusively in the testis, and are localized on the X chromosome as two-gene clusters and one separate copy. Population polymorphism of the βNACtes genes was studied using several wild-type D. melanogaster stocks, and βNACtes paralogs were compared with each other. A heterogeneous pattern was observed for βNACtes polymorphism: the 3′ genes of the two-gene clusters were low polymorphic, whereas, separate, the βNACtes1 gene was the most variable. The 5′ βNACtes copies of the two-gene tandems were practically identical, whereas the 3′ βNACtes copies were highly diverged. Hence, local gene conversion was assumed to provide for the selective homogenization of the 5′ genes. A comparison of the βNACtes paralogs showed that the majority of amino acid differences were in the N-terminal region, containing the βNAC domain. The McDonald-Kreitman test was used to analyze the divergence of βNACtes paralogs and implicated positive selection in the evolution of the βNACtes gene family.     

Genome sequence and comparative analysis of <Emphasis Type="Italic">Bacillus cereus</Emphasis> BC04, reveals genetic diversity and alterations for antimicrobial resistance

In this study, we delineated the genome sequence of a Bacillus cereus strain BC04 isolated from a stool sample in India. The draft genome is 5.1 Mb in size and consists of total 109 scaffolds, GC content is 35.2% with 5182 coding genes. The comparative analysis with other completely sequenced genomes highlights the unique presence of genomic islands, hemolysin, capsular synthetic protein, modifying enzymes accC7 and catA15, regulators of antibiotic resistance MarR and LysR with annotated functions related to virulence, stress response, and antimicrobial resistance. Overall, this study not only signifies the genetic diversity in gut isolate BC04 in particular, but also pinpoints the presence of unique genes possessed by B. cereus which can be pertinently exploited to design novel drugs and intervention strategies for the treatment of food borne diseases.     

High quality reference genome of drumstick tree (<Emphasis Type="Italic">Moringa oleifera</Emphasis> Lam.), a potential perennial crop

The drumstick tree (Moringa oleifera Lam.) is a perennial crop that has gained popularity in certain developing countries for its high-nutrition content and adaptability to arid and semi-arid environments. Here we report a high-quality draft genome sequence of M. oleifera. This assembly represents 91.78% of the estimated genome size and contains 19,465 protein-coding genes. Comparative genomic analysis between M. oleifera and related woody plant genomes helps clarify the general evolution of this species, while the identification of several species-specific gene families and positively selected genes in M. oleifera may help identify genes related to M. oleifera’s high protein content, fast-growth, heat and stress tolerance. This reference genome greatly extends the basic research on M. oleifera, and may further promote applying genomics to enhanced breeding and improvement of M. oleifera.     

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

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

Identification and functional analysis of anthocyanin biosynthesis genes in <Emphasis Type="Italic">Phalaenopsis</Emphasis> hybrids

Phalaenopsis species are among the most popular potted flowers for their fascinating flowers. When their whole-genome sequencing was completed, they have become useful for studying the molecular mechanism of anthocyanin biosynthesis. Here, we identified 49 candidate anthocyanin synthetic genes in the Phalaenopsis genome. Our results showed that duplication events might contribute to the expansion of some gene families, such as the genes encoding chalcone synthase (PeCHS), flavonoid 3′-hydroxylase (PeF3′H), and myeloblastosis (PeMYB). To elucidate their functions in anthocyanin biosynthesis, we conducted a global expression analysis. We found that anthocyanin synthesis occurred during the very early flower development stage and that the flavanone 3-hydroxylase (F3H), F3′H, and dihydroflavonol 4-reductase (DFR) genes played key roles in this process. Over-expression of Phalaenopsis flavonoid 3′,5′-hydroxylase (F3′5′H) in petunia showed that it had no function in anthocyanin production. Furthermore, global analysis of sequences and expression patterns show that the regulatory genes are relatively conserved and might be important in regulating anthocyanin synthesis through different combined expression patterns. To determine the functions of MYB2, 11, and 12, we over-expressed them in petunia and performed yeast two-hybrid analysis with anthocyanin (AN)1 and AN11. The MYB2 protein had strong activity in regulating anthocyanin biosynthesis and induced significant pigment accumulation in transgenic plant petals, whereas MYB11 and MYB12 had lower activities. Our work provided important improvement in the understanding of anthocyanin biosynthesis and established a foundation for floral colour breeding in Phalaenopsis through genetic engineering.     

The chloroplast genome sequence from <Emphasis Type="Italic">Eugenia uniflora</Emphasis>, a Myrtaceae from Neotropics

Eugenia uniflora is a plant native to tropical America that holds great ecological and economic importance. The complete chloroplast (cp) genome sequence of Eugenia uniflora, a member of the Neotropical Myrtaceae family, is reported here. The genome is 158,445 bp in length and exhibits a typical quadripartite structure of the large (LSC, 87,459 bp) and small (SSC, 18,318 bp) single-copy regions, separated by a pair of inverted repeats (IRs, 26,334 bp). It contains 111 unique genes, including 77 protein-coding genes, 30 tRNAs and 4 rRNAs. The genome structure, gene order, GC content and codon usage are similar to the typical angiosperm cp genomes. Comparison of the entire cp genomes of E. uniflora L. and three other Myrtaceae revealed an expansion of 43 bp in the intergenic spacer located between the IRA/large single-copy (LSC) border and the first gene of LSC region. Simple sequence repeat (SSR) analysis revealed that most SSRs are AT rich, which contribute to the overall AT richness of the cp genome. Additionally, fewer SSRs are distributed in the protein-coding sequences compared to the noncoding regions. Phylogenetic analysis among 58 species based on 57 cp genes demonstrated a closer relationship between E. uniflora L. and Syzygium cumini (L). Skeels compared to the Eucalyptus clade in the Myrtaceae family. The complete cp genome sequence of E. uniflora reported here has importance for population genetics, as well as phylogenetic and evolutionary studies in this species and other Myrtaceae species from Neotropical regions.     

Genome sequence of <Emphasis Type="Italic">Candida versatilis</Emphasis> and comparative analysis with other yeast

The genome of Candida versatilis was sequenced to understand its characteristics in soy sauce fermentation. The genome size of C. versatilis was 9.7 Mb, the content of G + C was 39.74 %, scaffolds of N50 were 1,229,640 bp in length, containing 4711 gene. There were predicted 269 tRNA genes and 2201 proteins with clear function. Moreover, the genome information of C. versatilis was compared with another salt-tolerant yeast Zygosaccharomyces rouxii and the model organism Saccharomyces cerevisiae. C. versatilis and Z. rouxii genome size was close and both smaller than 12.1 for the Mb of S. cerevisiae. Using the OrthoMCL protein, three genomes were divided into 4663 groups. There were about 3326 homologous proteins in C. versatilis, Z. rouxii and S. cerevisiae.     

The mitochondrial DNA of <Emphasis Type="Italic">Dictyostelium discoideum</Emphasis>: complete sequence,gene content and genome organization

We present an overview of the gene content and organization of the mitochondrial genome of Dictyostelium discoideum. The mitochondria genome consists of 55,564?bp with an A + T content of 72.6%. The identified genes include those for two ribosomal RNAs (rnl and rns), 18 tRNAs, ten subunits of the NADH dehydrogenase complex (nad1, 2, 3, 4, 4L, 5, 6, 7, 9 and 11), apocytochrome b (cytb), three subunits of the cytochrome oxidase (cox1/2 and 3), four subunits of the ATP synthase complex (atp1, 6, 8 and 9), 15 ribosomal proteins, and five other ORFs, excluding intronic ORFs. Notable features of D. discoideum mtDNA include the following. (1) All genes are encoded on the same strand of the DNA and a universal genetic code is used. (2) The cox1 gene has no termination codon and is fused to the downstream cox2 gene. The 13 genes for ribosomal proteins and four ORF genes form a cluster 15.4?kb long with several gene overlaps. (3) The number of tRNAs encoded in the genome is not sufficient to support the synthesis of mitochondrial protein. (4) In total, five group I introns reside in rnl and cox1/2, and three of those in cox1/2 contain four free-standing ORFs. We compare the genome to other sequenced mitochondrial genomes, particularly that of Acanthamoeba castellanii.     

Complete chloroplast genome sequence of <Emphasis Type="Italic">Capsicum</Emphasis><Emphasis Type="Italic">baccatum</Emphasis> var. <Emphasis Type="Italic">baccatum</Emphasis>

Molecular markers derived from the complete chloroplast genome can provide effective tools for species identification and phylogenetic resolution. Complete chloroplast (cp) genome sequences of Capsicum species have been reported. We herein report the complete chloroplast genome sequence of Capsicum baccatum var. baccatum, a wild Capsicum species. The total length of the chloroplast genome is 157,145 bp with 37.7 % overall GC content. One pair of inverted repeats, 25,910 bp in length, was separated by a small single-copy region (17,974 bp) and large single-copy region (87,351 bp). This region contains 86 protein-coding genes, 30 tRNA genes, 4 rRNA genes, and 11 genes contain one or two introns. Pair-wise alignments of chloroplast genome were performed for genome-wide comparison. Analysis revealed a total of 134 simple sequence repeat (SSR) motifs and 282 insertions or deletions variants in the C. baccatum var. baccatum cp genome. The types and abundances of repeat units in Capsicum species were relatively conserved, and these loci could be used in future studies to investigate and conserve the genetic diversity of the Capsicum species.     

Comparative chloroplast genomes of <Emphasis Type="Italic">Paris</Emphasis> Sect. <Emphasis Type="Italic">Marmorata</Emphasis>: insights into repeat regions and evolutionary implications

Background

Species of Paris Sect. Marmorata are valuable medicinal plants to synthesize steroidal saponins with effective pharmacological therapy. However, the wild resources of the species are threatened by plundering exploitation before the molecular genetics studies uncover the genomes and evolutionary significance. Thus, the availability of complete chloroplast genome sequences of Sect. Marmorata is necessary and crucial to the understanding the plastome evolution of this section and facilitating future population genetics studies. Here, we determined chloroplast genomes of Sect. Marmorata, and conducted the whole chloroplast genome comparison.

Results

This study presented detailed sequences and structural variations of chloroplast genomes of Sect. Marmorata. Over 40 large repeats and approximately 130 simple sequence repeats as well as a group of genomic hotspots were detected. Inverted repeat contraction of this section was inferred via comparing the chloroplast genomes with the one of P. verticillata. Additionally, almost all the plastid protein coding genes were found to prefer ending with A/U. Mutation bias and selection pressure predominately shaped the codon bias of most genes. And most of the genes underwent purifying selection, whereas photosynthetic genes experienced a relatively relaxed purifying selection.

Conclusions

Repeat sequences and hotspot regions can be scanned to detect the intraspecific and interspecific variability, and selected to infer the phylogenetic relationships of Sect. Marmorata and other species in subgenus Daiswa. Mutation and natural selection were the main forces to drive the codon bias pattern of most plastid protein coding genes. Therefore, this study enhances the understanding about evolution of Sect. Marmorata from the chloroplast genome, and provide genomic insights into genetic analyses of Sect. Marmorata.

     

Phylogeny and Molecular Evolution of miR820 and miR396 microRNA Families in <Emphasis Type="Italic">Oryza</Emphasis> AA Genomes

The phylogeny and evolution of the microRNA families, miR820 and miR396, was analysed across the AA genomes of the Oryza species, the close relatives of domesticated rice. A highly dynamic evolution of the miR820 family was revealed. The number of copies of MIR820 genes, their chromosomal location and the mature microRNA sequence varied greatly with a total of 16 novel miR820 variants being identified. The phylogeny of pre-MIR820 sequences revealed that MIR820 genes of recently evolved Oryza AA genomes may have derived from sequence divergence of one or a few ancestral genes found in wild Australian perennial rice populations, Taxon B (jpn2)-MIR820 genes. Genomic scale duplication played an important role in the evolution of some of the miR396 family genes in AA genome Oryza species. miR396 family contained a MIR396 gene cluster (MIR396a and MIR396c) which was conserved across the cereal genomes. Nucleotide diversity analysis at these two MIR396 loci revealed that domesticated rice has retained less than 10% of the total diversity present in wild species. In contrast, the nucleotide sequence of four MIR396 loci remained almost conserved across domesticated and wild rices, indicating that they were under extreme functional constraint and may be involved in regulating some fundamental processes which are important both for wild and domesticated rices. Expression analysis demonstrated that miR820 variants were expressed in O. glaberrima O. barthi and O. longistaminata genome. These findings pose new challenges to explain the possible role of miR820 variants identified.     

Complete mitochondrial DNA sequence and phylogenetic analysis of Zhikong scallop <Emphasis Type="Italic">Chlamys farreri</Emphasis> (Bivalvia: Pectinidae)

The complete mitochondrial genome of Zhikong scallop Chlamys farreri is 21,695 bp in length and contains 12 protein-coding genes (the atp8 gene is absent, as in most bivalves), 2 ribosomal RNA genes, and 22 transfer RNA genes. The heavy strand has an overall A+T content of 58.7%. GC and AT skews for the mt genome of C. farreri are 0.337 and ?0.184, respectively, indicating the nucleotide bias against C and A. The mitochondrial gene order of C. farreri differs drastically from the scallops Argopecten irradians, Mimachlamys nobilis and Placopecten magellanicus, which belong to the same family Pectinidae. 6623 bp non-coding nucleotides exist intergenically in the mitogenome of C. farreri, with a large continuous sequence (4763 bp) between tRNA ^Val and tRNA ^Asn . Two repeat families are found in the large continuous sequence, which seems to be a common feature of scallops. Phylogenetic analysis based on 12 concatenated amino acid sequences of protein-coding genes supports the monophyly of Pectinidae and paraphyletic Pteriomorphia with respect to Heteroconchia.     

Genome-wide identification and expression analysis of sulphate transporter (<Emphasis Type="Italic">SULTR</Emphasis>) genes under sulfur deficiency in <Emphasis Type="Italic">Brachypodium distachyon</Emphasis>

Sulphur is an important mineral element for plant growth and development. It involves in a number of metabolic processes with crucial functions. This study has performed a genome-wide analysis of sulfate transporter (SULTR) genes in Brachypodium distachyon. Ten putative SULTR genes were identified in Brachypodium genome. BdSULTR genes included 6–17 exons encoding a protein of 647–693 residues with basic nature. BdSULTR proteins included both sulfate_transp (PF00916) and STAS (PF01740) domains. BdSULTRs were classified into 4 groups based on the phylogenetic distribution. Promoter regions of all BdSULTR genes, except for BdSULTR3;3 and 3;5 included the SURECOREATSULTR11 elements. A considerable structural overlap was identified between superimposed SULTR1;3 and 3;1 proteins, indicating that SULTR1 members may also involve in plant stress response/tolerance like SULTR3 members. Microarray and RNA-Seq analyses also revealed the differential expression of SULTR 1 and 3 genes under different biotic/abiotic stresses. Protein–protein interaction partners of BdSULTRs were mainly related with adenylyl-sulfate kinases, 5′-adenylylsulfate reductases, ATP sulfurylases, and acyl carrier proteins. Moreover, expression profiles of identified BdSULTR genes under S-deficiency were analyzed using RT-qPCR. It was revealed that BdSULTR1;1 and 3;1 are highly expressed in plant roots as ~tenfold and ~fivefold, respectively, while BdSULTR2 (~15-fold) and 3;1 (~twofold) are abundantly expressed in leaf tissues.