Genome-wide transcriptome analysis of the orphan crop tef (Eragrostis tef (Zucc.) Trotter) under long-term low calcium stress.

Calcium (Ca2+) is one of the essential mineral nutrients for plant growth and development. However, the effects of long-term Ca2+ deficiency in orphan crops such as Tef [(Eragrostis tef) (Zucc.) Trotter], which accumulate high levels of Ca in the grains, remained unknown. Tef is a staple crop for nearly 70 million people in East Africa, particularly in Ethiopia and Eritrea. It is one of the most nutrient-dense grains, and is also more resistant to marginal soils and climatic conditions than main cereals like corn, wheat, and rice. In this study, tef plants were grown in a hydroponic solution containing optimum (1 mM) or low (0.01 mM) Ca2+, and plant growth parameters and whole-genome transcriptome were analyzed. Ca+2-deficient plants exhibited leaf necrosis, leaf curling, and growth stunting symptoms. Ca2+ deficiency significantly decreased root and shoot Ca, potassium (K), and copper content in both root and shoots. At the same time, it greatly increased root iron (Fe) content, suggesting the role of Ca2+ in the uptake and/or translocation of these minerals. Transcriptomic analysis using RNA-seq revealed that members of Ca2+ channels, including the cyclic nucleotide-gated channels and glutamate receptor-like channels, Ca2+-transporters, Ca2+-binding proteins and Ca2+-dependent protein kinases were differentially regulated by Ca+2 treatment. Moreover, several Fe/metal transporters, including members of vacuolar Fe transporters, yellow stripe-like, natural resistance-associated macrophage protein, and oligo-peptide transporters, were differentially regulated between shoot and root in response to Ca2+ treatment. Taken together, our findings suggest that Ca2+ deficiency affects plant growth and mineral accumulation by regulating the transcriptomes of several transporters and signaling genes.

Analysis of miRNAs responsive to long-term calcium deficiency in tef (Eragrostis tef (Zucc.) Trotter).

MicroRNAs (miRNAs) play an important role in growth, development, stress resilience, and epigenetic modifications of plants. However, the effect of calcium (Ca2+) deficiency on miRNA expression in the orphan crop tef (Eragrostis tef) remains unknown. In this study, we analyzed expression of miRNAs in roots and shoots of tef in response to Ca2+ treatment. miRNA-seq followed by bioinformatic analysis allowed us to identify a large number of small RNAs (sRNAs) ranging from 17 to 35 nt in length. A total of 1380 miRNAs were identified in tef experiencing long-term Ca2+ deficiency while 1495 miRNAs were detected in control plants. Among the miRNAs identified in this study, 161 miRNAs were similar with those previously characterized in other plant species and 348 miRNAs were novel, while the remaining miRNAs were uncharacterized. Putative target genes and their functions were predicted for all the known and novel miRNAs that we identified. Based on gene ontology (GO) analysis, the predicted target genes are known to have various biological and molecular functions including calcium uptake and transport. Pairwise comparison of differentially expressed miRNAs revealed that some miRNAs were specifically enriched in roots or shoots of low Ca2+-treated plants. Further characterization of the miRNAs and their targets identified in this study may help in understanding Ca2+ deficiency responses in tef and related orphan crops.

Complete mitochondrial genome of Cyclograpsus intermedius Ortmann, 1894 (Crustacea: Decapoda: Grapsoidea) specimen collected in South Korea.

In this study, the complete 16,184 bp mitochondrial genome of Cyclograpsus intermedius was determined from a specimen collected in South Korea. It consists of 13 protein-coding, 22 tRNA, 2 rRNA genes, and a non-coding A + T rich region. The base composition of the heavy strand in the mitochondrial genome was 34.7% A, 10.7% G, 18.7% C, and 35.9% T, resulting in a G + C content of 29.4%. A maximum-likelihood phylogenetic tree based on the 13 mitochondrial protein-coding genes showed that C. intermedius clustered together with the Varunidae. These molecular data will be useful for studying the evolutionary relationships among crab species.

A chloroplast genome of Forsythia saxatilis (Nakai) Nakai, an endemic species in Korea.

Forsythia saxatilis (Nakai) Nakai is an endemic species in Korea. We present the second complete chloroplast genome sequence of F. saxatilis showing that the chloroplast genome is 156,376 bp in length containing four subregions: a large single-copy (LSC) region of 87,097 bp, a small single-copy (SSC) region of 17,859 bp, and a pair of inverted repeat regions (IRs) of 25,710 bp. The genome contains 131 genes including 88 protein-coding genes, 35 tRNAs, and eight rRNAs. We found 40 base differences in 17 sites between our chloroplast genome and the previously reported chloroplast genome sequences of F. saxatilis. All of the differences were found in a region of 832 bp in the LSC region except for one site, which indicated potential sequencing errors in building the previously reported genome because of little substitutions in the homologous region of the other Forsythia species and even of Abeliophyllum distichum, a sister to Forsythia. The new chloroplast genomes of F. saxatilis and those of F. x intermedia and F. viridissima are identical, which suggests that F. saxatilis is maternally closely related to the other two species.

Silicon Enhances Biomass and Grain Yield in an Ancient Crop Tef [Eragrostis tef (Zucc.) Trotter].

Silicon (Si) is one of the beneficial plant mineral nutrients which is known to improve biotic and abiotic stress resilience and productivity in several crops. However, its beneficial role in underutilized or "orphan" crop such as tef [Eragrostis tef (Zucc.) Trotter] has never been studied before. In this study, we investigated the effect of Si application on tef plant performance. Plants were grown in soil with or without exogenous application of Na2SiO3 (0, 1.0, 2.0, 3.0, 4.0, and 5.0 mM), and biomass and grain yield, mineral content, chlorophyll content, plant height, and expression patterns of putative Si transporter genes were studied. Silicon application significantly increased grain yield (100%) at 3.0 mM Si, and aboveground biomass yield by 45% at 5.0 mM Si, while it had no effect on plant height. The observed increase in grain yield appears to be due to enhanced stress resilience and increased total chlorophyll content. Increasing the level of Si increased shoot Si and Na content while it significantly decreased the content of other minerals including K, Ca, Mg, P, S, Fe, and Mn in the shoot, which is likely due to the use of Na containing Si amendment. A slight decrease in grain Ca, P, S, and Mn was also observed with increasing Si treatment. The increase in Si content with increasing Si levels prompted us to analyze the expression of Si transporter genes. The tef genome contains seven putative Si transporters which showed high homology with influx and efflux Lsi transporters reported in various plant species including rice. The tef Lsi homologs were deferentially expressed between tissues (roots, leaves, nodes, and inflorescences) and in response to Si, suggesting that they may play a role in Si uptake and/or translocation. Taken together, these results show that Si application improves stress resilience and yield and regulates the expression of putative Si transporter genes. However, further study is needed to determine the physiological function of the putative Si transporters, and to study the effect of field application of Si on tef productivity.

Lipid Composition of Latex and Rubber Particles in Hevea brasiliensis and Taraxacum kok-saghyz.

Natural rubber is usually synthesized in the rubber particles present in the latex of rubber-producing plants such as the Pará rubber tree (Hevea brasiliensis) and rubber dandelion (Taraxacum kok-saghyz). Since the detailed lipid compositions of fresh latex and rubber particles of the plants are poorly known, the present study reports detailed compound lipid composition, focusing on phospholipids and galactolipids in the latex and rubber particles of the plants. In the fresh latex and rubber particles of both plants, phospholipids were much more dominant (85-99%) compared to galactolipids. Among the nine classes of phospholipids, phosphatidylcholines (PCs) were most abundant, at ~80%, in both plants. Among PCs, PC (36:4) and PC (34:2) were most abundant in the rubber tree and rubber dandelion, respectively. Two classes of galactolipids, monogalactosyl diacylglycerol and digalactosyl diacylglycerol, were detected as 12% and 1%, respectively, of total compound lipids in rubber tree, whereas their percentages in the rubber dandelion were negligible (< 1%). Overall, the compound lipid composition differed only slightly between the fresh latex and the rubber particles of both rubber plants. These results provide fundamental data on the lipid composition of rubber particles in two rubber-producing plants, which can serve as a basis for artificial rubber particle production in the future.

Prediction of Dhurrin Metabolism by Transcriptome and Metabolome Analyses in Sorghum.

Sorghum (Sorghum bicolor (L.)) Moench is an important food for humans and feed for livestock. Sorghum contains dhurrin which can be degraded into toxic hydrogen cyanide. Here, we report the expression patterns of 14 candidate genes related to dhurrin ((S)-4-Hydroxymandelnitrile-ß-D-glucopyranoside) metabolism and the effects of the gene expression on specific metabolite content in selected sorghum accessions. Dhurrin-related metabolism is vigorous in the early stages of development of sorghum. The dhurrin contents of most accessions tested were in the range of approximately 6-22 µg mg-1 fresh leaf tissue throughout growth. The p-hydroxybenzaldehyde (pHB) contents were high at seedling stages, but almost nonexistent at adult stages. The contents of p-hydroxyphenylacetic acid (pHPAAc) were relatively low throughout growth compared to those of dhurrin or pHB. Generally, the expression of the candidate genes was higher at seedling stage than at other stages and decreased gradually as plants grew. In addition, we identified significant SNPs, and six of them were potentially associated with non-synonymous changes in CAS1. Our results may provide the basis for choosing breeding materials to regulate cyanide contents in sorghum varieties to prevent HCN toxicity of livestock or to promote drought tolerance or pathogen resistance.

Heterologous activation of the Hevea PEP16 promoter in the rubber-producing laticiferous tissues of Taraxacum kok-saghyz.

Hevea brasiliensis, the most abundant rubber crop, is used widely for the commercial production of natural rubber. To reduce the risk of a shortage in the supply of natural rubber that may arise from a single major rubber crop, rubber dandelion (Taraxacum kok-saghyz) has been developed as an alternative rubber-producing crop by using a transgenic approach. However, it is necessary to identify a suitable promoter for the transfer of rubber biosynthesis-related genes to the species. In this study, the promoter region of H. brasiliensis PEP16, which was isolated as a potentially important component in rubber biosynthesis, was sequenced and a pPEP16::GUS fusion construct was introduced into T. kok-saghyz. Histological and fluorometric studies using transgenic T. kok-saghyz plants indicated that the HbPEP16 promoter was highly activated in a laticiferous tissue-specific manner under normal growth conditions and that promoter activation was tightly regulated by various hormones and external signals. These findings suggested that the HbPEP16 promoter may be a useful molecular tool for the manipulation of gene expression in the laticiferous tissues of T. kok-saghyz.

Illustrations in the Journal of Craniofacial Surgery.

The aim of this study is to analyze the medical illustrations that appeared in the Journal of Craniofacial Surgery (JCFS) during the last 3 years (2016-2018).All illustrations in 2676 articles from 24 issues were analyzed.In total, 363 articles (13.6%) contained illustrations (for total of 589 illustrations), and 38 articles (1.4%) cited illustrations from other articles. The number of illustrations increased over time (152 in 2016, 181 in 2017, and 256 in 2018; P = 0.007). The signature of the illustrator appeared in 16.5% of the illustrations, and the illustrator's name was acknowledged in 7.6%. The most frequent contents of the 589 illustrations were surgical procedures (55.4%), followed by anatomy (34.1%). Among the 555 illustrations showing regions of the body, the most frequently illustrated region was the oral cavity (41.6%), followed by the eye and periorbital region (18.4%), the nose (12.4%), and the head (11.4%). Of the 555 illustrations showing tissue, the tissue most frequently depicted was bone (52.8%), followed by soft tissue (42.5%). The mean score for clarity of the intended message was 3.5 ±â€Š1.2. The mean score for artistic skill was 3.2 ±â€Š1.0.As the number of illustrations in JCFS increases, it is increasingly important that information on the illustrators be given. The journal editor should recommend that authors provide the illustrators' information. From the author's perspective, the best illustration is one that tells the story the author wants to convey to readers. Therefore, the crucial point is the expression of the author's goal, rather than the esthetics of the illustration.

Constructing a Reference Genome in a Single Lab: The Possibility to Use Oxford Nanopore Technology.

The whole genome sequencing (WGS) has become a crucial tool in understanding genome structure and genetic variation. The MinION sequencing of Oxford Nanopore Technologies (ONT) is an excellent approach for performing WGS and it has advantages in comparison with other Next-Generation Sequencing (NGS): It is relatively inexpensive, portable, has simple library preparation, can be monitored in real-time, and has no theoretical limits on reading length. Sorghum bicolor (L.) Moench is diploid (2n = 2x = 20) with a genome size of about 730 Mb, and its genome sequence information is released in the Phytozome database. Therefore, sorghum can be used as a good reference. However, plant species have complex and large genomes when compared to animals or microorganisms. As a result, complete genome sequencing is difficult for plant species. MinION sequencing that produces long-reads can be an excellent tool for overcoming the weak assembly of short-reads generated from NGS by minimizing the generation of gaps or covering the repetitive sequence that appears on the plant genome. Here, we conducted the genome sequencing for S. bicolor cv. BTx623 while using the MinION platform and obtained 895,678 reads and 17.9 gigabytes (Gb) (ca. 25× coverage of reference) from long-read sequence data. A total of 6124 contigs (covering 45.9%) were generated from Canu, and a total of 2661 contigs (covering 50%) were generated from Minimap and Miniasm with a Racon through a de novo assembly using two different tools and mapped assembled contigs against the sorghum reference genome. Our results provide an optimal series of long-read sequencing analysis for plant species while using the MinION platform and a clue to determine the total sequencing scale for optimal coverage that is based on various genome sizes.

A Transcriptome Approach Toward Understanding Fruit Softening in Persimmon.

Persimmon (Diospyros kaki Thunb.), which is a climacteric fruit, softens in 3-5 weeks after harvest. However, little is known regarding the transcriptional changes that underlie persimmon ripening. In this study, high-throughput de novo RNA sequencing was performed to examine differential expression between freshly harvested (FH) and softened (ST) persimmon fruit peels. Using the Illumina HiSeq platform, we obtained 259,483,704 high quality reads and 94,856 transcripts. After the removal of redundant sequences, a total of 31,258 unigenes were predicted, 1,790 of which were differentially expressed between FH and ST persimmon (1,284 up-regulated and 506 down-regulated in ST compared with FH). The differentially expressed genes (DEGs) were further subjected to KEGG pathway analysis. Several pathways were found to be up-regulated in ST persimmon, including "amino sugar and nucleotide sugar metabolism." Pathways down-regulated in ST persimmon included "photosynthesis" and "carbon fixation in photosynthetic organisms." Expression patterns of genes in these pathways were further confirmed using quantitative real-time RT-PCR. Ethylene gas production during persimmon softening was monitored with gas chromatography and found to be correlated with the fruit softening. Transcription involved in ethylene biosynthesis, perception and signaling was up-regulated. On the whole, this study investigated the key genes involved in metabolic pathways of persimmon fruit softening, especially implicated in increased sugar metabolism, decreased photosynthetic capability, and increased ethylene production and other ethylene-related functions. This transcriptome analysis provides baseline information on the identity and modulation of genes involved in softening of persimmon fruits and can underpin the future development of technologies to delay softening in persimmon.

On the study of microbial transcriptomes using second- and third-generation sequencing technologies.

Second-generation sequencing technologies transformed the study of microbial transcriptomes. They helped reveal the transcription start sites and antisense transcripts of microbial species, improving the microbial genome annotation. Quantification of genome-wide gene expression levels allowed for functional studies of microbial research. Ever-evolving sequencing technologies are reshaping approaches to studying microbial transcriptomes. Recently, Oxford Nanopore Technologies delivered a sequencing platform called MinION, a third-generation sequencing technology, to the research community. We expect it to be the next sequencing technology that enables breakthroughs in life science fields. The studies of microbial transcriptomes will be no exception. In this paper, we review microbial transcriptomics studies using second- generation sequencing technology. We also discuss the prospect of microbial transcriptomics studies with thirdgeneration sequencing.

Functional Characterization of Phalaenopsis aphrodite Flowering Genes PaFT1 and PaFD.

We show that the key flowering regulators encoded by Phalaenopsis aphrodite FLOWERING LOCUS T1 (PaFT1) and PaFD share high sequence homologies to these from long-day flowering Arabidopsis and short-day flowering rice. Interestingly, PaFT1 is specifically up-regulated during flowering inductive cooling treatment but is not subjected to control by photoperiod in P. aphrodite. Phloem or shoot apex-specific expression of PaFT1 restores the late flowering of Arabidopsis ft mutants. Moreover, PaFT1 can suppress the delayed flowering caused by SHORT VEGATATIVE PHASE (SVP) overexpression as well as an active FRIGIDA (FRI) allele, indicating the functional conservation of flowering regulatory circuit in different plant species. PaFT1 promoter:GUS in Arabidopsis showed similar staining pattern to that of Arabidopsis FT in the leaves and guard cells but different in the shoot apex. A genomic clone or heat shock-inducible expression of PaFT1 is sufficient to the partial complementation of the ft mutants. Remarkably, ectopic PaFT1 expression also triggers precocious heading in rice. To further demonstrate the functional conservation of the flowering regulators, we show that PaFD, a bZIP transcription factor involved in flowering promotion, interacts with PaFT1, and PaFD partially complemented Arabidopsis fd mutants. Transgenic rice expressing PaFD also flowered early with increased expression of rice homologues of APETALA1 (AP1). Consistently, PaFT1 knock-down Phalaenopsis plants generated by virus-induced gene silencing exhibit delayed spiking. These studies suggest functional conservation of FT and FD genes, which may have evolved and integrated into distinct regulatory circuits in monopodial orchids, Arabidopsis and rice that promote flowering under their own inductive conditions.

Draft Genome Sequence of a Taxonomically Unique Neisseria Strain Isolated from a Greater White-Fronted Goose (Anser albifrons) Egg on the North Slope of Alaska.

We report here the draft genome sequence of a unique Neisseria strain that was isolated from a greater white-fronted goose (Anser albifrons) egg. The sequencing was performed with an Illumina MiSeq system, and the sequence consists of 275 contigs. The total genome is 2,397,978 bp long and has a G+C content of 46.4%.

A unique open reading frame within the comX gene of Streptococcus mutans regulates genetic competence and oxidative stress tolerance.

Streptococcus mutans displays complex regulation of genetic competence, with ComX controlling late competence gene transcription. The rcrRPQ operon has been shown to link oxidative stress tolerance, (p)ppGpp metabolism and competence in S. mutans. Importantly, an rcrR polar (ΔrcrR-P) mutant is hyper-transformable, but an rcrR non-polar (ΔrcrR-NP) mutant cannot be transformed. Transcriptome comparisons of the rcrR mutants using RNA-Seq and quantitative real-time polymerase chain reaction revealed little expression in the 5' region of comX in ΔrcrR-NP, but high level expression in the 3' region. Northern blotting with comX probes revealed two distinct transcripts in the ΔrcrR-P and ΔrcrR-NP strains, and 5' Rapid Amplification of cDNA Ends mapped the 5' terminus of the shorter transcript to nt +140 of the comX structural gene, where a unique 69-aa open reading frame, termed XrpA, was encoded in a different reading frame than ComX. Two single-nucleotide substitution mutants (comX::T162C; comX::T210A) were introduced to disrupt XrpA without affecting the sequence of ComX. When the mutations were in the ΔrcrR-NP genetic background, ComX production and transformation were restored. Overexpression of xrpA led to impaired growth in aerobic conditions and decreased transformability. These results reveal an unprecedented mechanism for competence regulation and stress tolerance by a gene product encoded within the comX gene that appears unique to S. mutans.

Draft Genome Sequence of an Atypical Strain of Streptococcus pneumoniae Isolated from a Respiratory Infection.

Next-generation sequencing was used to investigate an unknown clinical respiratory infection. This new strain of Streptococcus pneumoniae, ASVL_JC_0001, was isolated from a clinical specimen from a patient with bronchitis and pulmonary inflammation. The draft genome sequence, obtained with an Illumina MiSeq sequencing system, consists of 83 large contigs, a total of 2,092,532 bp long, and has a GC content of 40.3%.

Trithorax group protein Oryza sativa Trithorax1 controls flowering time in rice via interaction with early heading date3.

Trithorax group proteins are chromatin-remodeling factors that activate target gene expression by antagonistically functioning against the Polycomb group. In Arabidopsis (Arabidopsis thaliana), Arabidopsis Trithorax protein1 (ATX1) regulates flowering time and floral organ identity. Here, we observed that suppression of Oryza sativa Trithorax1 (OsTrx1), an ortholog of ATX1, delayed flowering time in rice (Oryza sativa). Because the delay occurred only under long-day conditions, we evaluated the flowering signal pathways that specifically function under long-day conditions. Among them, the OsMADS50 and Heading date1 pathways were not affected by the mutation. However, the Grain number, plant height, and heading date7 (Ghd7) pathway was altered in ostrx1. Transcript levels of OsGI, phytochrome genes, and Early heading date3 (Ehd3), which function upstream of Ghd7, were unchanged in the mutant. Because Trx group proteins form a complex with other proteins to modify the chromatin structure of target genes, we investigated whether OsTrx1 interacts with a previously identified protein that functions upstream of Ghd7. We demonstrated that the plant homeodomain motif of OsTrx1 binds to native histone H3 from the calf thymus and that OsTrx1 binds to Ehd3 through the region between the plant homeodomain and SET domains. Finally, we showed that the SET domain at the C-terminal end of OsTrx1 has histone H3 methyltransferase activity when incubated with oligonucleosomes. Our results suggest that OsTrx1 plays an important role in regulating flowering time in rice by modulating chromatin structure.

Transcriptomic and genomic evidence for Streptococcus agalactiae adaptation to the bovine environment.

BACKGROUND: Streptococcus agalactiae is a major cause of bovine mastitis, which is the dominant health disorder affecting milk production within the dairy industry and is responsible for substantial financial losses to the industry worldwide. However, there is considerable evidence for host adaptation (ecotypes) within S. agalactiae, with both bovine and human sourced isolates showing a high degree of distinctiveness, suggesting differing ability to cause mastitis. Here, we (i) generate RNAseq data from three S. agalactiae isolates (two putative bovine adapted and one human) and (ii) compare publicly available whole genome shotgun sequence data from an additional 202 isolates, obtained from six host species, to elucidate possible genetic factors/adaptations likely important for S. agalactiae growth and survival in the bovine mammary gland. RESULTS: Tests for differential expression showed distinct expression profiles for the three isolates when grown in bovine milk. A key finding for the two putatively bovine adapted isolates was the up regulation of a lactose metabolism operon (Lac.2) that was strongly correlated with the bovine environment (all 36 bovine sourced isolates on GenBank possessed the operon, in contrast to only 8/151 human sourced isolates). Multi locus sequence typing of all genome sequences and phylogenetic analysis using conserved operon genes from 44 S. agalactiae isolates and 16 additional Streptococcus species provided strong evidence for acquisition of the operon via multiple lateral gene transfer events, with all Streptococcus species known to be major causes of mastitis, identified as possible donors. Furthermore, lactose fermentation tests were only positive for isolates possessing Lac.2. Combined, these findings suggest that lactose metabolism is likely an important adaptation to the bovine environment. Additional up regulation in the bovine adapted isolates included genes involved in copper homeostasis, metabolism of purine, pyrimidine, glycerol and glucose, and possibly aminoglycoside antibiotic resistance. CONCLUSION: We detected several genetic factors likely important in S. agalactiae's adaptation to the bovine environment, in particular lactose metabolism. Of concern is the up regulation of a putative antibiotic resistance gene (GCN5-related N-acetyltransferase) that might reflect an adaptation to the use of aminoglycoside antibiotics within this environment.

Gene regulation by CcpA and catabolite repression explored by RNA-Seq in Streptococcus mutans.

A bacterial transcriptome of the primary etiological agent of human dental caries, Streptococcus mutans, is described here using deep RNA sequencing. Differential expression profiles of the transcriptome in the context of carbohydrate source, and of the presence or absence of the catabolite control protein CcpA, revealed good agreement with previously-published DNA microarrays. In addition, RNA-seq considerably expanded the repertoire of DNA sequences that showed statistically-significant changes in expression as a function of the presence of CcpA and growth carbohydrate. Novel mRNAs and small RNAs were identified, some of which were differentially expressed in conditions tested in this study, suggesting that the function of the S. mutans CcpA protein and the influence of carbohydrate sources has a more substantial impact on gene regulation than previously appreciated. Likewise, the data reveal that the mechanisms underlying prioritization of carbohydrate utilization are more diverse than what is currently understood. Collectively, this study demonstrates the validity of RNA-seq as a potentially more-powerful alternative to DNA microarrays in studying gene regulation in S. mutans because of the capacity of this approach to yield a more precise landscape of transcriptomic changes in response to specific mutations and growth conditions.