Pyropia yezoensis is a model organism often used to investigate the mechanisms underlying stress tolerance in intertidal zones. The digital gene expression (DGE) approach was used to characterize a genome-wide comparative analysis of differentially expressed genes (DEGs) that influence the physiological, developmental or biochemical processes in samples subjected to 4 treatments: high-temperature stress (HT), chilling stress (CS), freezing stress (FS) and normal temperature (NT). Equal amounts of total RNAs collected from 8 samples (two biological replicates per treatment) were sequenced using the Illumina/Solexa platform. Compared with NT, a total of 2202, 1334 and 592 differentially expressed unigenes were detected in HT, CS and FS respectively. Clustering analysis suggested P. yezoensis acclimates to low and high-temperature stress condition using different mechanisms: In heat stress, the unigenes related to replication and repair of DNA and protein processing in endoplasmic reticulum were active; however at low temperature stresses, unigenes related to carbohydrate metabolism and energy metabolism were active. Analysis of gene differential expression showed that four categories of DEGs functioning as temperature sensors were found, including heat shock proteins, H2A, histone deacetylase complex and transcription factors. Heat stress caused chloroplast genes down-regulated and unigenes encoding metacaspases up-regulated, which is an important regulator of PCD. Cold stress caused an increase in the expression of FAD to improve the proportion of polyunsaturated fatty acids. An up-regulated unigene encoding farnesyl pyrophosphate synthase was found in cold stress, indicating that the plant hormone ABA also played an important role in responding to temperature stress in P. yezoensis. The variation of amount of unigenes and different gene expression pattern under different temperature stresses indicated the complicated and diverse regulation mechanism in response to temperature stress in P. yezoensis. Several common metabolism pathways were found both in P. yezoensis and in higher plants, such as FAD in low-temperature stress and HSP in heat stress. Meanwhile, many chloroplast genes and unigene related to the synthesis of abscisic acid were detected, revealing its unique temperature-regulation mechanism in this intertidal species. This sequencing dataset and analysis may serve as a valuable resource to study the mechanisms involved in abiotic stress tolerance in intertidal seaweeds.

Pyropia yezoensis, a marine red alga, is an ideal research model for studying the mechanisms of abiotic stress tolerance in intertidal seaweed. Real-time quantitative polymerase chain reaction (RT-qPCR) is the most commonly used method to analyze gene expression levels. To accurately quantify gene expression, selection and validation of stable reference genes is required. We used transcriptome profiling data from different abiotic stress treatments to identify six genes with relatively stable expression levels: MAP, ATPase, CGS1, PPK, DPE2, and FHP. These six genes and three conventional reference genes, UBC, EF1-α, and eif4A, were chosen as candidates for optimal reference gene selection. Five common statistical approaches (geNorm, ΔCt method, NormFinder, BestKeeper, and ReFinder) were used to identify the stability of each reference gene. Our results show that: MAP, UBC, and FHP are stably expressed in all analyzed conditions; CGS1 and UBC are stably expressed under conditions of dehydration stress; and MAP, UBC, and CGS1 are stably expressed under conditions of temperature stress. We have identified appropriate reference genes for RT-qPCR in P. yezoensis under different abiotic stress conditions which will facilitate studies of gene expression under these conditions.

In this study, we determined the mitochondrial genome of Rhynchophorus ferrugineus. The mitochondrial genome is 15 924 bp in length (GC content: 25.6%), encodes 2 ribosomal RNA genes, 13 protein-coding genes, 21 transfer RNA genes, and 1 D-loop region. nad6 and cob are overlapped by 30 bp and atp8 and atp6 are overlapped by 12 bp. The phylogenetic tree involving 29 available closely related species further validated the new determined sequences and phylogeny of R. ferrugineus.

Buddleja colvilei Hook.f. & Thomson (Scrophulariaceae) is a threatened alpine plant with a distribution throughout the Himalayas, also used as an ornamental plant. The name Buddleja sessilifolia B.S. Sun ex S.Y. Pao was assigned in 1983 to a plant distributed throughout the Gaoligong Mountains, but the name was later placed in synonymy with B. colvilei in the Flora of China. In this study we sequenced the complete chloroplast (cp) genomes of two individuals of B. colvilei and three individuals of B. sessilifolia from across the range. Both molecular and morphological analysis support the revision of B. sessilifolia. The phylogenetic analysis constructed with the whole cp genomes, the large single-copy regions (LSC), small single-copy regions (SSC), inverted repeat (IR) and the nuclear genes 18S/ITS1/5.8S/ITS2/28S all supported B. sessilifolia as a distinct species. Additionally, coalescence-based species delimitation methods (bGMYC, bPTP) using the whole chloroplast datasets also supported B. sessilifolia as a distinct species. The results suggest that the B. sessilifolia lineage was early diverging among the Asian Buddleja species. Overall gene contents were similar and gene arrangements were found to be highly conserved in the two species, however, fixed differences were found between the two species. A total of 474 single nucleotide polymorphisms (SNPs) were identified between the two species. The Principal Coordinate Analysis of the morphological characters resolved two groups and supported B. sessilifolia as a distinct species. Discrimination of B. colvilei and B. sessilifolia using morphological characters and the redescription of B. sessilifolia are detailed here.

In this study, complete mitochondrial genome of the large earth bumblebee, Bombus terrestris, was first determined. The mitochondrial genome is 174 000 bp in length, and encodes 1 D-loop region, 2 ribosomal RNA genes, 13 protein-coding genes and 20 transfer RNA genes. Average GC content of this genome is 15%. nad6 and cob are overlapped with 12 bp, while atp8 and atp6 are overlapped with 18 bp. Phylogeny analysis indicates Bombus species exhibit closer genetic distance with Melipona species than with Apis species.