Genome-wide identification, characterization, and evolutionary analysis of the barley TALE gene family and its expression profiles in response to exogenous hormones.

Three-amino-loop-extension (TALE) family belongs to the homeobox gene superfamily and occurs widely in plants, playing a crucial role in regulating their growth and development. Currently, genome-wide analysis of the TALE family has been completed in many plants. However, the systematic identification and hormone response analysis of the TALE gene family in barley are still lacking. In this study, 21 TALE candidate genes were identified in barley, which can be divided into KNOX and BELL subfamilies. Barley TALE members in the same subfamily of the phylogenetic tree have analogically conserved motifs and gene structures, and segmental duplications are largely responsible for the expansion of the HvTALE family. Analysis of TALE orthologous and homologous gene pairs indicated that the HvTALE family has mainly undergone purifying selective pressure. Through spatial structure simulation, HvKNOX5-HvKNOX6 and HvKNOX5-HvBELL11 complexes are all formed through hydrogen bonding sites on both the KNOX2 and homeodomain (HD) domains of HvKNOX5, which may be essential for protein interactions among the HvTALE family members. Expression pattern analyses reveal the potential involvement of most HvTALE genes in responses to exogenous hormones. These results will lay the foundation for regulation and function analyses of the barley TALE gene family in plant growth and development by hormone regulation.

Genome-Wide Identification, Characterization, Evolutionary Analysis, and Expression Pattern of the GPAT Gene Family in Barley and Functional Analysis of HvGPAT18 under Abiotic Stress.

Glycerol-3-phosphoacyltransferase (GPAT) is an important rate-limiting enzyme in the biosynthesis of triacylglycerol (TAG), which is of great significance for plant growth, development, and response to abiotic stress. Although the characteristics of GPAT have been studied in many model plants, little is known about its expression profile and function in barley, especially under abiotic stress. In this study, 22 GPAT genes were identified in the barley genome and divided into three groups (I, II, III), with the latter Group III subdivided further into three subgroups based on the phylogenetic analysis. The analyses of conserved motifs, gene structures, and the three-dimensional structure of HvGPAT proteins also support this classification. Through evolutionary analysis, we determined that HvGPATs in Group I were the earliest to diverge during 268.65 MYA, and the differentiation of other HvGPATs emerged during 86.83-169.84 MYA. The tissue expression profile showed that 22 HvGPAT genes were almost not expressed in INF1 (inflorescence 1). Many functional elements related to stress responses and hormones in cis-element analysis, as well as qRT-PCR results, confirm that these HvGPAT genes were involved in abiotic stress responses. The expression level of HvGPAT18 was significantly increased under abiotic stress and its subcellular localization indicated its function in the endoplasmic reticulum. Various physiological traits under abiotic stress were evaluated using transgenic Arabidopsis to gain further insight into the role of HvGPAT18, and it was found that transgenic seedlings have stronger resistance under abiotic stress than to the wild-type (WT) plants. Overall, our results provide new insights into the evolution and function of the barley GPAT gene family and enable us to explore the molecular mechanism of functional diversity behind the evolutionary history of these genes.

Two new anti-inflammatory trachylobane diterpenoids from Euphorbia atoto.

Two new rare trachylobane euphoratones A-B (1-2), together with five known diterpenoids (compounds 3-7), were isolated from the aerial parts of Euphorbia atoto. Their structures were unambiguously elucidated through HRESIMS, 1D and 2D NMR spectral analysis. Compounds 1, 3, 4 and 7 showed weak anti-inflammatory activities (IC50 77.49 ± 6.34, 41.61 ± 14.49, 16.00 ± 1.71 and 33.41 ± 4.52 µM, respectively), compared to the positive control quercetin (IC50 15.23 ± 0.65 µM).

How do US corporations communicate interculturally with their Chinese stakeholders: Analysis of GM Company's social media posts from the cultural value perspective.

Social Media is an important means of communication with audiences around the world. The purpose of this study was to explore whether GM-a famous US auto company adapts its US Cultural values to suit the prevalent cultural values of its Chinese stakeholders on Chinese social media. Content analysis was used to evaluate the cultural content of GM Company's posts on Weibo and Twitter. Although influenced by the special features of the car industry, there is still enough evidence that the communication style of the US auto Company makes cultural adaption on Chinese social media, reflecting more Chinese prevalent cultural values.

Cytotoxic steroids from the stems of Strophanthus divaricatus.

Phytochemical investigation on the stems of Strophanthus divaricatus led to the isolation of four undescribed cardiac glycosides and one undescribed C21 pregnane, together with eleven known steroids. Their structures were elucidated by a comprehensive analysis of HRESIMS, 1D and 2D NMR spectra. The absolute configuration of 16 was determined by comparison of the experimental and computed ECD spectra. Compounds 1-13 and 15 displayed potent to significant cytotoxicity against human cancer cell lines K562, SGC-7901, A549 and HeLa with IC50 values of 0.02-16.08, 0.04-23.13, 0.06-22.31 and 0.06-15.13 µM, respectively.

Transcriptomic analysis of rapeseed (Brassica napus. L.) seed development in Xiangride, Qinghai Plateau, reveals how its special eco-environment results in high yield in high-altitude areas.

As one of the most important oil crops, rapeseed (Brassica napus) is cultivated worldwide to produce vegetable oil, animal feed, and biodiesel. As the population grows and the need for renewable energy increases, the breeding and cultivation of high-yield rapeseed varieties have become top priorities. The formation of a high rapeseed yield is so complex because it is influenced not only by genetic mechanisms but also by many environmental conditions, such as climatic conditions and different farming practices. Interestingly, many high-yield areas are located in special eco-environments, for example, in the high-altitude Xiangride area of the Qinghai Plateau. However, the molecular mechanisms underlying the formation of high yields in such a special eco-environment area remain largely unknown. Here, we conducted field yield analysis and transcriptome analysis in the Xiangride area. Compared with the yield and environmental factors in the Xinning area (a low-yielding area), we found that the relatively longer daylight length is the key to high rapeseed yield in the Xiangride area, which leads up to a 52.1% increase in rapeseed yield, especially the increase in thousand seed weight and silique number (SN). Combined with transcriptome H-cluster analysis and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analyses, we can assume that the grain development of rapeseed in the Xiangride area is ahead of schedule and lasts for a long time, leading to the high-yield results in the Xiangride area, confirmed by the expression analysis by quantitative real-time polymerase chain reaction (qRT-PCR) of yield-related genes. Our results provide valuable information for further exploring the molecular mechanism underlying high yield in special ecological environments and provide a helpful reference for studying seed development characteristics in special-producing regions for Brassica napus.

Genome-wide identification and expression profile under abiotic stress of the barley non-specific lipid transfer protein gene family and its Qingke Orthologues.

BACKGROUND: Plant non-specific lipid transfer proteins (nsLTPs), a group of small, basic ubiquitous proteins to participate in lipid transfer, cuticle formation and stress response, are involved in the regulation of plant growth and development. To date, although the nsLTP gene family of barley (Hordeum vulgare L.) has been preliminarily identified, it is still unclear in the recently completed genome database of barley and Qingke, and its transcriptional profiling under abiotic stress has not been elucidated as well. RESULTS: We identified 40 barley nsLTP (HvLTP) genes through a strict screening strategy based on the latest barley genome and 35 Qingke nsLTP (HtLTP) orthologues using blastp, and these LTP genes were divided into four types (1, 2, D and G). At the same time, a comprehensive analysis of the physical and chemical characteristics, homology alignment, conserved motifs, gene structure and evolution of HvLTPs and HtLTPs further supported their similar nsLTP characteristics and classification. The genomic location of HvLTPs and HtLTPs showed that these genes were unevenly distributed, and obvious HvLTP and HtLTP gene clusters were found on the 7 chromosomes including six pairs of tandem repeats and one pair of segment repeats in the barley genome, indicating that these genes may be co-evolutionary and co-regulated. A spatial expression analysis showed that most HvLTPs and HtLTPs had different tissue-specific expression patterns. Moreover, the upstream cis-element analysis of HvLTPs and HtLTPs showed that there were many different stress-related transcriptional regulatory elements, and the expression pattern of HvLTPs and HtLTPs under abiotic stress also indicated that numerous HvLTP and HtLTP genes were related to the abiotic stress response. Taken together, these results may be due to the differences in promoters rather than by genes themselves resulting in different expression patterns under abiotic stress. CONCLUSION: Due to a stringent screening and comprehensive analysis of the nsLTP gene family in barley and Qingke and its expression profile under abiotic stress, this study can be considered a useful source for the future studies of nsLTP genes in either barley or Qingke or for comparisons of different plant species.

The complete chloroplast genome of Salsola abrotanoides (Chenopodiaceae), a desert halophyte shrub in China.

Salsola abrotanoides, one of the dominant plant species of desert vegetation, adapts well to the arid, saline, and alkaline environment in the Qinghai-Tibetan Plateau. Here, we reported the complete chloroplast sequence and characters of S. abrotanoides based on the Illumina NovaSeq Platform. The chloroplast genome is 151,622 bp in length, containing a pair of inverted repeated (IR) regions of 23,701 bp, a large single copy (LSC) region of 84,658 bp, and a small single copy (SSC) region of 19,562 bp. And the chloroplast genome sequence encodes 130 genes totally, including 85 mRNA genes, 37 tRNA genes, and 8 rRNA genes. S. abrotanoides is the first species of Genus Salsola and the chloroplast sequence will provide a valuable resource for the phylogenetic studies of Chenopodiaceae.

Genome-wide characterization and expression profiling of the relation of the HD-Zip gene family to abiotic stress in barley (Hordeum vulgare L.).

The homeodomain-leucine zipper (HD-Zip) gene family plays an important role in plant growth and environmental responses. At present, research on the HD-Zip gene family of barley is incomplete. In this study, 32 HD-Zip genes (HvHD-Zip 1-32) were identified from the barley genome and were subsequently divided into four subfamilies according to conserved structure and motif analysis. Whole genome replication events in barley and Arabidopsis, rice, and wheat HD-Zip gene families were analyzed, yielding 3, 14 and 25 gene pairs, respectively, but no segmental or tandem duplication events were identified in the barley HD-Zip gene family. Subsequently, quantitative real-time PCR (qRT-PCR) analysis revealed that the HvHD-Zip gene is sensitive to drought stress and that members of the HD-Zip I and HD-Zip IV subfamilies are generally more sensitive to abiotic stresses. Our results suggest a relationship between barley resistance and the potential key HvHD-Zip gene, which lay the foundation for further functional studies.

Flow cytometic analysis of Penaeus monodon haemocyte responses to poly I:C.

This study was aimed at investigating the cellular responses of Penaeus monodon haemocytes to poly I:C stimulation using flow cytometric assay. Total haemocyte count (THC), percentages of different haemocyte subpopulations [hyaline cells (HC), semigranular cells (SGC) and granular cells (GC)], non-specific esterase activity (EA), total reactive oxygen species/reactive nitrogen species (ROS/RNS) production, nitric oxide (NO) production, apoptotic haemocyte ratio and plasmic phenoloxidase (PO) activity were determined in poly I:C-injected shrimp. Results showed that poly I:C at a low dose (5 µg shrimp-1) caused obvious increases in THC, GC proportion, ROS/RNS production and NO production, but had no significant effect on EA, apoptosis and PO activity. In the early stage of poly I:C injection at a high dose (20 µg shrimp-1), THC and GC proportion improvements could also be observed, suggesting that GC might be induced to release from hemocytopoietic or other tissues to participate in immune response, and this subpopulation might be the main cell type involved in the cellular defence against virus. In the later period, proportions of both GC and SGC reduced paralleled by THC reduction, indicating that depletion of GC and SGC was mainly contributed to the reduced count of circulating haemocyte. Obvious increases in ROS/RNS production and NO production were induced in haemocyte of shrimp under a high dose of poly I:C stimulation, but only slight rise of EA and suppression of PO activity could be observed in poly I:C-stimulated shrimp, suggesting that ROS/RNS-dependent system was vital in the immune defence of shrimp against virus. On the other hand, increase of apoptotic haemocyte ratio and THC reduction were presented after the drastic increases of ROS/RNS and NO productions, implying that the stimulated ROS/RNS might be excess and harmful, and was the major factor for the haemocyte apoptosis and depletion. THC recovered after 48 h injection, while haemocyte apoptosis also returned to the control level, suggesting that apoptosis might be contributed to eliminate damaged, weak or infected haemocytes to renew the circulating haemocytes, and it could be considered as an important defending strategy against virus.

Isolation and Characterization of Ftsz Genes in Cassava.

The filamenting temperature-sensitive Z proteins (FtsZs) play an important role in plastid division. In this study, three FtsZ genes were isolated from the cassava genome, and named MeFtsZ1, MeFtsZ2-1, and MeFtsZ2-2, respectively. Based on phylogeny, the MeFtsZs were classified into two groups (FtsZ1 and FtsZ2). MeFtsZ1 with a putative signal peptide at N-terminal, has six exons, and is classed to FtsZ1 clade. MeFtsZ2-1 and MeFtsZ2-2 without a putative signal peptide, have seven exons, and are classed to FtsZ2 clade. Subcellular localization found that all the three MeFtsZs could locate in chloroplasts and form a ring in chloroplastids. Structure analysis found that all MeFtsZ proteins contain a conserved guanosine triphosphatase (GTPase) domain in favor of generate contractile force for cassava plastid division. The expression profiles of MeFtsZ genes by quantitative reverse transcription-PCR (qRT-PCR) analysis in photosynthetic and non-photosynthetic tissues found that all of the MeFtsZ genes had higher expression levels in photosynthetic tissues, especially in younger leaves, and lower expression levels in the non-photosynthetic tissues. During cassava storage root development, the expressions of MeFtsZ2-1 and MeFtsZ2-2 were comparatively higher than MeFtsZ1. The transformed Arabidopsis of MeFtsZ2-1 and MeFtsZ2-2 contained abnormally shape, fewer number, and larger volume chloroplasts. Phytohormones were involved in regulating the expressions of MeFtsZ genes. Therefore, we deduced that all of the MeFtsZs play an important role in chloroplast division, and that MeFtsZ2 (2-1, 2-2) might be involved in amyloplast division and regulated by phytohormones during cassava storage root development.

Comparative transcriptome profile of the leaf elongation zone of wild barley (Hordeum spontaneum) eibi1 mutant and its isogenic wild type

Abstract The naturally occurring wild barley mutant eibi1/hvabcg31 suffers from severe water loss due to the permeable leaf cuticle. Eibi1/HvABCG31 encodes a full ATP-binding cassette (ABC) transporter, HvABCG31, playing a role in cutin deposition in the elongation zone of growing barley leaves. The eibi1 allele has pleiotropic effects on the appearance of leaves, plant stature, fertility, spike and grain size, and rate of germination. Comparative transcriptome profile of the leaf elongation zone of the eibi1 mutant as well as its isogenic wild type showed that various pathogenesis-related genes were up-regulated in the eibi1 mutant. The known cuticle-related genes that we analyzed did not show significant expression difference between the mutant and wild type. These results suggest that the pleiotropic effects may be a compensatory consequence of the activation of defense genes in the eibi1 mutation. Furthermore, we were able to find the mutation of the eibi1/hvabcg31 allele by comparing transcript sequences, which indicated that the RNA-Seq is useful not only for researches on general molecular mechanism but also for the identification of possible mutant genes.

Identification, Expression, and Functional Analysis of the Fructokinase Gene Family in Cassava.

Fructokinase (FRK) proteins play important roles in catalyzing fructose phosphorylation and participate in the carbohydrate metabolism of storage organs in plants. To investigate the roles of FRKs in cassava tuber root development, seven FRK genes (MeFRK1-7) were identified, and MeFRK1-6 were isolated. Phylogenetic analysis revealed that the MeFRK family genes can be divided into α (MeFRK1, 2, 6, 7) and ß (MeFRK3, 4, 5) groups. All the MeFRK proteins have typical conserved regions and substrate binding residues similar to those of the FRKs. The overall predicted three-dimensional structures of MeFRK1-6 were similar, folding into a catalytic domain and a ß-sheet ''lid" region, forming a substrate binding cleft, which contains many residues involved in the binding to fructose. The gene and the predicted three-dimensional structures of MeFRK3 and MeFRK4 were the most similar. MeFRK1-6 displayed different expression patterns across different tissues, including leaves, stems, tuber roots, flowers, and fruits. In tuber roots, the expressions of MeFRK3 and MeFRK4 were much higher compared to those of the other genes. Notably, the expression of MeFRK3 and MeFRK4 as well as the enzymatic activity of FRK were higher at the initial and early expanding tuber stages and were lower at the later expanding and mature tuber stages. The FRK activity of MeFRK3 and MeFRK4 was identified by the functional complementation of triple mutant yeast cells that were unable to phosphorylate either glucose or fructose. The gene expression and enzymatic activity of MeFRK3 and MeFRK4 suggest that they might be the main enzymes in fructose phosphorylation for regulating the formation of tuber roots and starch accumulation at the tuber root initial and expanding stages.

Comparative transcriptome profile of the leaf elongation zone of wild barley (Hordeum spontaneum) eibi1 mutant and its isogenic wild type.

The naturally occurring wild barley mutant eibi1/hvabcg31 suffers from severe water loss due to the permeable leaf cuticle. Eibi1/HvABCG31 encodes a full ATP-binding cassette (ABC) transporter, HvABCG31, playing a role in cutin deposition in the elongation zone of growing barley leaves. The eibi1 allele has pleiotropic effects on the appearance of leaves, plant stature, fertility, spike and grain size, and rate of germination. Comparative transcriptome profile of the leaf elongation zone of the eibi1 mutant as well as its isogenic wild type showed that various pathogenesis-related genes were up-regulated in the eibi1 mutant. The known cuticle-related genes that we analyzed did not show significant expression difference between the mutant and wild type. These results suggest that the pleiotropic effects may be a compensatory consequence of the activation of defense genes in the eibi1 mutation. Furthermore, we were able to find the mutation of the eibi1/hvabcg31 allele by comparing transcript sequences, which indicated that the RNA-Seq is useful not only for researches on general molecular mechanism but also for the identification of possible mutant genes.

Structure, Expression, and Functional Analysis of the Hexokinase Gene Family in Cassava.

Hexokinase (HXK) proteins play important roles in catalyzing hexose phosphorylation and sugar sensing and signaling. To investigate the roles of HXKs in cassava tuber root development, seven HXK genes (MeHXK1-7) were isolated and analyzed. A phylogenetic analysis revealed that the MeHXK family can be divided into five subfamilies of plant HXKs. MeHXKs were clearly divided into type A (MeHXK1) and type B (MeHXK2-7) based on their N-terminal sequences. MeHXK1-5 all had typical conserved regions and similar protein structures to the HXKs of other plants; while MeHXK6-7 lacked some of the conserved regions. An expression analysis of the MeHXK genes in cassava organs or tissues demonstrated that MeHXK2 is the dominant HXK in all the examined tissues (leaves, stems, fruits, tuber phloems, and tuber xylems). Notably, the expression of MeHXK2 and the enzymatic activity of HXK were higher at the initial and expanding tuber stages, and lower at the mature tuber stage. Furthermore, the HXK activity of MeHXK2 was identified by functional complementation of the HXK-deficient yeast strain YSH7.4-3C (hxk1, hxk2, glk1). The gene expression and enzymatic activity of MeHXK2 suggest that it might be the main enzyme for hexose phosphorylation during cassava tuber root development, which is involved in sucrose metabolism to regulate the accumulation of starch.

Characterization and genetic mapping of the ß-diketone deficient eceriferum-b barley mutant.

KEY MESSAGE: The barley eceriferum-b.2 (cer-b.2) mutant produces glossy leaf sheaths and is deficient in the cuticular wax component 14,16-hentriacontanedione. The mutated gene maps to a 1.3-cM interval on chromosome 3HL flanked by the genes MLOC_10972 and MLOC_69561. The cuticular wax coating of leaves and stems in many grass species is responsible for the plants' glaucous appearance. A major component of the wax is a group of ß-diketone compounds. The barley eceriferum-b.2 (cer-b.2) mutant produces glossy leaf sheaths and is deficient for the compound 14,16-hentriacontanedione. A linkage analysis based on 708 gametes allowed the gene responsible for the mutant phenotype to be mapped to a 1.3-cM interval on chromosome 3HL flanked by the two genes MLOC_10972 and _69561. The product of the wild type allele may represent a step in the ß-diketone synthesis pathway.

SpBADH of the halophyte Sesuvium portulacastrum strongly confers drought tolerance through ROS scavenging in transgenic Arabidopsis.

Glycine betaine (GB) accumulation is involved in abiotic stress. However, it is not known whether BADH, the key enzyme of GB synthesis, utilizes the antioxidant system to confer drought stress tolerance. In this study, a novel member of the ALDH10 gene family, SpBADH, was isolated from Sesuvium portulacastrum. The expression of this gene was up-regulated by NaCl, PEG6000, H2O2, ABA and high temperature in S. portulacastrum. SpBADH overexpression in Arabidopsis resulted in higher BADH activity and GB content and might increase tolerance to drought/osmotic stresses, specifically strong tolerance to drought stress. Transgenic lines exhibited lower MDA and H2O2 contents but higher proline, POD, SOD and CAT contents than the wild type under drought and osmotic stresses. SpBADH overexpression in Arabidopsis also enhanced the expression of ROS-related genes including AtSOD, AtPOD, AtCAT, AtAPX and Atpsb under drought and osmotic stresses. Thus, SpBADH increases plant tolerance to drought or osmotic stresses by reducing H2O2, increasing proline, and activating antioxidative enzymes to improve ROS scavenging.

SpAHA1 and SpSOS1 Coordinate in Transgenic Yeast to Improve Salt Tolerance.

In plant cells, the plasma membrane Na+/H+ antiporter SOS1 (salt overly sensitive 1) mediates Na+ extrusion using the proton gradient generated by plasma membrane H+-ATPases, and these two proteins are key plant halotolerance factors. In the present study, two genes from Sesuvium portulacastrum, encoding plasma membrane Na+/H+ antiporter (SpSOS1) and H+-ATPase (SpAHA1), were cloned. Localization of each protein was studied in tobacco cells, and their functions were analyzed in yeast cells. Both SpSOS1 and SpAHA1 are plasma membrane-bound proteins. Real-time polymerase chain reaction (PCR) analyses showed that SpSOS1 and SpAHA1 were induced by salinity, and their expression patterns in roots under salinity were similar. Compared with untransformed yeast cells, SpSOS1 increased the salt tolerance of transgenic yeast by decreasing the Na+ content. The Na+/H+ exchange activity at plasma membrane vesicles was higher in SpSOS1-transgenic yeast than in the untransformed strain. No change was observed in the salt tolerance of yeast cells expressing SpAHA1 alone; however, in yeast transformed with both SpSOS1 and SpAHA1, SpAHA1 generated an increased proton gradient that stimulated the Na+/H+ exchange activity of SpSOS1. In this scenario, more Na+ ions were transported out of cells, and the yeast cells co-expressing SpSOS1 and SpAHA1 grew better than the cells transformed with only SpSOS1 or SpAHA1. These findings demonstrate that the plasma membrane Na+/H+ antiporter SpSOS1 and H+-ATPase SpAHA1 can function in coordination. These results provide a reference for developing more salt-tolerant crops via co-transformation with the plasma membrane Na+/H+ antiporter and H+-ATPase.

Characterization and genetic mapping of eceriferum-ym (cer-ym), a cutin deficient barley mutant with impaired leaf water retention capacity.

The cuticle covers the aerial parts of land plants, where it serves many important functions, including water retention. Here, a recessive cuticle mutant, eceriferum-ym (cer-ym), of Hordeum vulgare L. (barley) showed abnormally glossy spikes, sheaths, and leaves. The cer-ym mutant plant detached from its root system was hypersensitive to desiccation treatment compared with wild type plants, and detached leaves of mutant lost 41.8% of their initial weight after 1 h of dehydration under laboratory conditions, while that of the wild type plants lost only 7.1%. Stomata function was not affected by the mutation, but the mutant leaves showed increased cuticular permeability to water, suggesting a defective leaf cuticle, which was confirmed by toluidine blue staining. The mutant leaves showed a substantial reduction in the amounts of the major cutin monomers and a slight increase in the main wax component, suggesting that the enhanced cuticle permeability was a consequence of cutin deficiency. cer-ym was mapped within a 0.8 cM interval between EST marker AK370363 and AK251484, a pericentromeric region on chromosome 4H. The results indicate that the desiccation sensitivity of cer-ym is caused by a defect in leaf cutin, and that cer-ym is located in a chromosome 4H pericentromeric region.

Molecular Mechanisms Underlying Hull-Caryopsis Adhesion/Separation Revealed by Comparative Transcriptomic Analysis of Covered/Naked Barley (Hordeum vulgare L.).

The covered/naked caryopsis trait of barley is an important agronomic trait because it is directly linked to dietary use. The formation of covered/naked caryopsis is controlled by an NUD transcription factor, which is involved in pericarp cuticle development. However, the molecular mechanism underlying this trait remains so far largely unknown. In this study, comparative transcriptomes of grains three weeks after anthesis of Tibetan Hulless barley landrace Dulihuang and covered barley Morex were analyzed using RNA-seq technique. A total of 4031 differentially expressed genes (DEGs) were identified. The Nud gene was overexpressed in Morex, with trace expression in Dulihuang. Among seventeen cuticle related DEGs, sixteen were down regulated and one up regulated in Morex. These results suggest that the Nud gene in covered caryopsis might down regulate cuticle related genes, which may cause a permeable cuticle over pericarp, leading to a hull-caryopsis organ fusion. A functional cuticle covering the pericarp of naked caryopsis might be the result of deletion or low expression level of the Nud gene. The functional cuticle defines a perfect boundary to separate the caryopsis from the hull in naked barley.