De novo assembly and characterization of the complete chloroplast genome of radish (Raphanus sativus L.)

Radish (Raphanus sativus L.) is an edible root vegetable crop that is cultivated worldwide and whose genome has been sequenced. Here we report the complete nucleotide sequence of the radish cultivar WK10039 chloroplast (cp) genome, along with a de novo assembly strategy using whole genome shotgun sequence reads obtained by next generation sequencing. The radish cp genome is 153,368 bp in length and has a typical quadripartite structure, composed of a pair of inverted repeat regions (26,217 bp each), a large single copy region (83,170 bp), and a small single copy region (17,764 bp). The radish cp genome contains 87 predicted protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Sequence analysis revealed the presence of 91 simple sequence repeats (SSRs) in the radish cp genome.     

The complete chloroplast genome sequence of Mahonia bealei (Berberidaceae) reveals a significant expansion of the inverted repeat and phylogenetic relationship with other angiosperms

Mahonia bealei (Berberidaceae) is a frequently-used traditional Chinese medicinal plant with efficient anti-inflammatory ability. This plant is one of the sources of berberine, a new cholesterol-lowering drug with anti-diabetic activity. We have sequenced the complete nucleotide sequence of the chloroplast (cp) genome of M. bealei. The complete cp genome of M. bealei is 164,792 bp in length, and has a typical structure with large (LSC 73,052 bp) and small (SSC 18,591 bp) single-copy regions separated by a pair of inverted repeats (IRs 36,501 bp) of large size. The Mahonia cp genome contains 111 unique genes and 39 genes are duplicated in the IR regions. The gene order and content of M. bealei are almost unarranged which is consistent with the hypothesis that large IRs stabilize cp genome and reduce gene loss-and-gain probabilities during evolutionary process. A large IR expansion of over 12 kb has occurred in M. bealei, 15 genes (rps19, rpl22, rps3, rpl16, rpl14, rps8, infA, rpl36, rps11, petD, petB, psbH, psbN, psbT and psbB) have expanded to have an additional copy in the IRs. The IR expansion rearrangement occurred via a double-strand DNA break and subsequence repair, which is different from the ordinary gene conversion mechanism. Repeat analysis identified 39 direct/inverted repeats 30 bp or longer with a sequence identity ≥ 90%. Analysis also revealed 75 simple sequence repeat (SSR) loci and almost all are composed of A or T, contributing to a distinct bias in base composition. Comparison of protein-coding sequences with ESTs reveals 9 putative RNA edits and 5 of them resulted in non-synonymous modifications in rpoC1, rps2, rps19 and ycf1. Phylogenetic analysis using maximum parsimony (MP) and maximum likelihood (ML) was performed on a dataset composed of 65 protein-coding genes from 25 taxa, which yields an identical tree topology as previous plastid-based trees, and provides strong support for the sister relationship between Ranunculaceae and Berberidaceae. Molecular dating analyses suggest that Ranunculaceae and Berberidaceae diverged between 90 and 84 mya, which is congruent with the fossil records and with recent estimates of the divergence time of these two taxa.     

The complete chloroplast genome sequence of Taxus chinensis var. mairei (Taxaceae): loss of an inverted repeat region and comparative analysis with related species

Taxus chinensis var. mairei (Taxaceae) is a domestic variety of yew species in local China. This plant is one of the sources for paclitaxel, which is a promising antineoplastic chemotherapy drugs during the last decade. We have sequenced the complete nucleotide sequence of the chloroplast (cp) genome of T. chinensis var. mairei. The T. chinensis var. mairei cp genome is 129,513 bp in length, with 113 single copy genes and two duplicated genes (trnI-CAU, trnQ-UUG). Among the 113 single copy genes, 9 are intron-containing. Compared to other land plant cp genomes, the T. chinensis var. mairei cp genome has lost one of the large inverted repeats (IRs) found in angiosperms, fern, liverwort, and gymnosperm such as Cycas revoluta and Ginkgo biloba L. Compared to related species, the gene order of T. chinensis var. mairei has a large inversion of ~ 110 kb including 91 genes (from rps18 to accD) with gene contents unarranged. Repeat analysis identified 48 direct and 2 inverted repeats 30 bp long or longer with a sequence identity greater than 90%. Repeated short segments were found in genes rps18, rps19 and clpP. Analysis also revealed 22 simple sequence repeat (SSR) loci and almost all are composed of A or T.     

Complete chloroplast genome sequences of Praxelis (Eupatorium catarium Veldkamp), an important invasive species

Praxelis (Eupatorium catarium Veldkamp) is a new hazardous invasive plant species that has caused serious economic losses and environmental damage in the Northern hemisphere tropical and subtropical regions. Although previous studies focused on detecting the biological characteristics of this plant to prevent its expansion, little effort has been made to understand the impact of Praxelis on the ecosystem in an evolutionary process. The genetic information of Praxelis is required for further phylogenetic identification and evolutionary studies. Here, we report the complete Praxelis chloroplast (cp) genome sequence. The Praxelis chloroplast genome is 151,410 bp in length including a small single-copy region (18,547 bp) and a large single-copy region (85,311 bp) separated by a pair of inverted repeats (IRs; 23,776 bp). The genome contains 85 unique and 18 duplicated genes in the IR region. The gene content and organization are similar to other Asteraceae tribe cp genomes. We also analyzed the whole cp genome sequence, repeat structure, codon usage, contraction of the IR and gene structure/organization features between native and invasive Asteraceae plants, in order to understand the evolution of organelle genomes between native and invasive Asteraceae. Comparative analysis identified the 14 markers containing greater than 2% parsimony-informative characters, indicating that they are potential informative markers for barcoding and phylogenetic analysis. Moreover, a sister relationship between Praxelis and seven other species in Asteraceae was found based on phylogenetic analysis of 28 protein-coding sequences. Complete cp genome information is useful for plant phylogenetic and evolutionary studies within this invasive species and also within the Asteraceae family.     

Comparative genomics of four Liliales families inferred from the complete chloroplast genome sequence of Veratrum patulum O. Loes. (Melanthiaceae)

The sequence of the chloroplast genome, which is inherited maternally, contains useful information for many scientific fields such as plant systematics, biogeography and biotechnology because its characteristics are highly conserved among species. There is an increase in chloroplast genomes of angiosperms that have been sequenced in recent years. In this study, the nucleotide sequence of the chloroplast genome (cpDNA) of Veratrum patulum Loes. (Melanthiaceae, Liliales) was analyzed completely. The circular double-stranded DNA of 153,699 bp consists of two inverted repeat (IR) regions of 26,360 bp each, a large single copy of 83,372 bp, and a small single copy of 17,607 bp. This plastome contains 81 protein-coding genes, 30 distinct tRNA and four genes of rRNA. In addition, there are six hypothetical coding regions (ycf1, ycf2, ycf3, ycf4, ycf15 and ycf68) and two open reading frames (ORF42 and ORF56), which are also found in the chloroplast genomes of the other species. The gene orders and gene contents of the V. patulum plastid genome are similar to that of Smilax china, Lilium longiflorum and Alstroemeria aurea, members of the Smilacaceae, Liliaceae and Alstroemeriaceae (Liliales), respectively. However, the loss rps16 exon 2 in V. patulum results in the difference in the large single copy regions in comparison with other species. The base substitution rate is quite similar among genes of these species. Additionally, the base substitution rate of inverted repeat region was smaller than that of single copy regions in all observed species of Liliales. The IR regions were expanded to trnH_GUG in V. patulum, a part of rps19 in L. longiflorum and A. aurea, and whole sequence of rps19 in S. china. Furthermore, the IGS lengths of rbcL-accD-psaI region were variable among Liliales species, suggesting that this region might be a hotspot of indel events and the informative site for phylogenetic studies in Liliales. In general, the whole chloroplast genome of V. patulum, a potential medicinal plant, will contribute to research on the genetic applications of this genus.     

The complete mitochondrial genome of Sasakia funebris (Leech) (Lepidoptera: Nymphalidae) and comparison with other Apaturinae insects

Sasakia funebris, a member of the lepidopteran family, Nymphalidae (superfamily Papilionoidea) is a rare species and is found only in some areas of South China. In this study, the 15,233 bp long complete mitochondrial genome of S. funebris was determined, and harbors the gene arrangement identical to all other sequenced lepidopteran insects. The nucleotide composition of the genome is highly A + T biased, accounting for 81.2%. All protein-coding genes (PCGs) start with typical ATN codons, except for COI which begins with the CGA codon. All tRNAs have a typical clover-leaf secondary structure, except for tRNA^Ser(AGN), the dihydrouridine (DHU) arm of which forms a simple loop. The S. funebris A + T-rich region of 370 bp contains several features common to the Lepidoptera insects, including the motif ATAGA followed by a 19 bp poly-T stretch, and two tandem repeats consisting of 18 bp repeat units and 14 bp repeat units. The phylogenetic analyses of Apaturinae based on mitogenome sequences showed: (S. funebris + Sasakia charonda) + (Apatura metis + Apatura ilia). This result is consistent with the morphological classification.     

The co-existence of two growth hormone receptors and their differential expression profiles between female and male tongue sole (Cynoglossus semilaevis)

Growth hormone receptor (Ghr) is a single-transmembrane pass protein which is important in initiating the ability of growth hormone (Gh) to regulate development and somatic growth in vertebrates. In this study, molecular cloning, expression analysis of two different ghr genes (ghr1 and ghr2) in the tongue sole (Cynoglossus semilaevis) was conducted. As a result, the ghr1 and ghr2 cDNA sequences are 2364 bp and 3125 bp, each of which encodes a transmembrane protein of 633 and 561 amino acids (aa), respectively. Besides, the ghr1 gene includes nine exons and eight introns. The sex-specific tissue expression was analyzed by using 14 tissues from females, normal males and extra-large male adults. Both the ghr1 and ghr2 were predominantly expressed in the liver, and the ghr1 expression level in normal males was 1.6 and 1.4 times as much as those in females and extra-large males, while the ghr2 mRNA expression level in normal males was 1.1 and 1.2 times as much as those in females and extra-large males, respectively. Ontogenetic expression analysis at early life stages indicated that the ghr1 and ghr2 mRNAs were detected at all of the 35 sampling points (from oosphere to 410 days-old). Furthermore, the sex differences in ghr mRNA expressions were also examined by using a full-sib family of C. semilaevis. Significantly higher levels of ghr1 mRNA were observed in males than in females at most stages of the sampling period (P < 0.01). The ghr2 mRNA expression at most stages exhibited a significant sexual difference at each sampling point (P < 0.01) without any variation trend related with the sexes during the whole sampling period.     

Genomic structure, polymorphism and expression analysis of the growth hormone (GH) gene in female and male Half-smooth tongue sole (Cynoglossus semilaevis)

Growth hormone (GH) is a polypeptide which is an important regulator of development and somatic growth in teleosts, and may be associated with the mechanisms which drive sexual growth dimorphism in the Half-smooth tongue sole (Cynoglossus semilaevis). In this study, the full length gh cDNA was cloned from C. semilaevis by homology cloning and the rapid amplification of cDNA ends polymerase chain reaction (RACE-PCR). The full-length gh cDNA is 826 bp and contains an open reading frame (ORF) of 603 bp encoding a protein of 200 amino acids (AA). The precursor of gh consists of a 17 amino-acid signal peptide followed by a 183 amino-acid mature polypeptide. GH gene sequences obtained from female and male adults consist of 3428 bp and 3371 bp, respectively, each of which includes six exons and five introns, and the difference in the GH gene size was mainly caused by the microsatellites. When 14 tissues from females, normal males and extra-large male adults were analyzed for sex-specific tissue expression, the gh mRNA was found to be predominantly expressed in the pituitary, and the expression levels in females were 3.6 times as much as those in normal males, while the mRNA expression in extra-large males was 1.7 times as much as those in normal males. Sex differences in gh mRNA expression during development were also examined by using a full-sib family of C. semilaevis, and the gh mRNA was detected at all of the 12 time points sampled from 10 to 380 days-old. A significant increase in gh mRNA was detected starting in 80 day old fish and was then followed by a drop to very low levels starting at 230 day old fish. Differential expression indicated that the gh expression level in females was significantly higher than males (P < 0.01) at all of the stages except for 10 days-old. Two microsatellite loci were identified in the second intron of the GH gene. Using these two polymorphic markers to genotype 224 individuals, there was no significant difference between the females and males in the Bohai Sea, the Yellow Sea and the hatchery samples.     

Molecular cloning and expression in vitro of a carboxylesterase gene from the Glanville fritillary butterfly (Melitaea cinxia)

Carboxylesterase (EC 3.1.1.1) is a member of the carboxyl/cholinesterase (CCE) superfamily, which is widely distributed in animals, plants and microorganisms. This enzyme has been known to be associated with insecticide resistance and detoxification. Although CCEs have been extensively studied in insects, including lepidopterans, the research on butterflies, a major subgroup in Lepidoptera, is still poor. In the present study, we cloned a CCE gene (McCCE1) from the Glanville fritillary butterfly (Melitaea cinxia, Lepidoptera: Nymphalidae). The full-length cDNA encoding McCCE1 was 1786 bp, containing a 1641 bp open reading frame encoding 546 amino acids, a 38 bp 5′-untranslated region (5′-UTR), and a 107 bp 3′-UTR with a poly(A) tail. The functionally conserved amino acids in McCCE1 shared the 55% identity with the cytoplasmic esterase CCE017a in Helicoverpa armigera (Lepidoptera: Noctuidae), which has been associated with detoxification. Assays in vitro showed that the recombinant McCCE1 could hydrolyze α- and β-naphthyl acetate. Thus, the present study adds to the body of knowledge concerning the detoxification of pesticides by lepidopterans.     

Characterization of the chloroplast genome sequence of oil palm (Elaeis guineensis Jacq.)

Oil palm (Elaeis guineensis Jacq.) is an economically important crop, which is grown for oil production. To better understand the molecular basis of oil palm chloroplasts, we characterized the complete chloroplast (cp) genome sequence obtained from 454 pyrosequencing. The oil palm cp genome is 156,973 bp in length consisting of a large single-copy region of?85,192 bp flanked on each side by inverted repeats of 27,071 bp with a small single-copy region of 17,639 bp joining the?repeats. The genome contains 112 unique genes: 79 protein-coding genes, 4 ribosomal RNA genes and 29 tRNA genes. By aligning the cp?genome sequence with oil palm cDNA sequences, we observed 18 non-silent and 10 silent RNA editing events among 19 cp protein-coding genes. Creation of an initiation codon by RNA editing in rpl2 has been reported in several monocots and was also found in the oil palm cp genome. Fifty common chloroplast protein-coding genes from 33 plant taxa were used to construct ML and MP?phylogenetic trees. Their topologies are similar and strongly support for the position of E. guineensis as the sister of closely related species Phoenix dactylifera in Arecaceae (palm families) of monocot subtrees.     

Molecular cloning and mRNA expression analysis of ornithine decarboxylase antizyme 2 in ovarian follicles of the Sichuan white goose (Anser cygnoides)

The ornithine decarboxylase antizyme 2 (OAZ2) gene is a member of the antizyme gene family. Antizymes play pivotal roles in various cellular pathways, including polyamine anabolism and apoptosis. The molecular structure and expression profile of the OAZ2 in goose ovarian follicles have not been reported. In this study, the OAZ2 cDNA sequence of the Sichuan white goose was cloned (Anser cygnoides), and phylogenetic and structural analyses of the OAZ2 were performed. The expression profiling of OAZ2 mRNA in goose ovarian follicles was examined using quantitative real-time PCR. The sequence analysis showed that the 756 bp OAZ2 sequence contained two overlapping open reading frames (ORF). ORF1 was 99 bp in length, and encoded a 32 aa polypeptide. ORF2 was 477 bp in length, and encoded a 158 aa polypeptide. The frameshift site that initiates the translation of ORF2 was located at nucleotide position 97 in the OAZ2. The analysis of OAZ2 mRNA expression in hierarchical follicles showed that the level of OAZ2 mRNA was higher in the SWF and F2 follicular stages than that in the ovarian stroma (P < 0.05). The lowest level of OAZ2 expression was detected in the ovarian stroma. These results suggest that the highly conserved frameshift region plays an important role in sustaining the function of OAZs. Furthermore, the significantly higher level of OAZ2 mRNA in the SWF stage indicates that OAZ2 may be involved in recruiting hierarchical follicles. Our results also suggest that OAZ2 may augment the effects of OAZ1 in follicle development.     

High level expression of a novel β-mannanase from Chaetomium sp. exhibiting efficient mannan hydrolysis

A novel β-mannanase gene (CsMan5A) was cloned from Chaetomium sp. CQ31 and expressed in Pichia pastoris. It had an open reading frame of 1251 bp encoding 416 amino acids and contained two introns. The deduced amino acid sequence shared the highest similarity (73%) with the β-mannanase from Emericella nidulans and belongs to glycosyl hydrolase family 5. The recombinant β-mannanase (CsMan5A) was secreted at extremely high levels of 50,030 U mL⁻¹ and 6.1 mg mL⁻¹ in high cell density fermentor. The purified enzyme was optimally active at pH 5.0 and 65 °C and displayed broad pH stability (pH 5.0-11.0) and exhibited specificity towards locust bean gum (K_m = 3.1 mg mL⁻¹), guar gum (K_m = 9.3 mg mL⁻¹) and konjac powder (K_m = 10.5 mg mL⁻¹). It efficiently degraded mannan polysaccharides into mannose and mannooligosacccharides, and also hydrolyzed mannotriose and mannotetraose. These properties make CsMan5A highly useful in food, feed and paper/pulp industries.