Genome-wide identification of U-box gene family and expression analysis in response to saline-alkali stress in foxtail millet (Setaria italica L. Beauv).

E3 ubiquitin ligases are central modifiers of plant signaling pathways that regulate protein function, localization, degradation, and other biological processes by linking ubiquitin to target proteins. E3 ubiquitin ligases include proteins with the U-box domain. However, there has been no report about the foxtail millet (Setaria italica L. Beauv) U-box gene family (SiPUB) to date. To explore the function of SiPUBs, this study performed genome-wide identification of SiPUBs and expression analysis of them in response to saline-alkali stress. A total of 70 SiPUBs were identified, which were unevenly distributed on eight chromosomes. Phylogenetic and conserved motif analysis demonstrated that SiPUBs could be clustered into six subfamilies (I-VI), and most SiPUBs were closely related to the homologues in rice. Twenty-eight types of cis-acting elements were identified in SiPUBs, most of which contained many light-responsive elements and plant hormone-responsive elements. Foxtail millet had 19, 78, 85, 18, and 89 collinear U-box gene pairs with Arabidopsis, rice, sorghum, tomato, and maize, respectively. Tissue specific expression analysis revealed great variations in SiPUB expression among different tissues, and most SiPUBs were relatively highly expressed in roots, indicating that SiPUBs may play important roles in root development or other growth and development processes of foxtail millet. Furthermore, the responses of 15 SiPUBs to saline-alkali stress were detected by qRT-PCR. The results showed that saline-alkali stress led to significantly differential expression of these 15 SiPUBs, and SiPUB20/48/70 may play important roles in the response mechanism against saline-alkali stress. Overall, this study provides important information for further exploration of the biological function of U-box genes.

Genome-wide identification of AAAP gene family and expression analysis in response to saline-alkali stress in foxtail millet (Setaria italica L.).

Amino acid/auxin permease (AAAP) genes encode a large family of protein transporters that play important roles in various aspects of plant growth and development. Here, we performed genome-wide identification of members in the foxtail millet (Setaria italica L.) AAAP family (SiAAAP) and their saline-alkali stress-induced expression patterns, resulting in the identification of 65 SiAAAP genes, which could be divided into eight subfamilies. Except for SiAAAP65, the remaining 64 genes were located on nine chromosomes of foxtail millet. Gene structure and conserved motif analyses indicated that the members in the same subfamily are highly conserved. Gene duplication event analysis suggested that tandem duplication may be the main factor driving the expansion of this gene family, and Ka/Ks analysis indicated that all the duplicated genes have undergone purifying selection. Transcriptome analysis showed differential expression of SiAAAPs in roots, stems, leaves, and tassel inflorescence. Analysis of cis-acting elements in the promoter indicated that SiAAAPs contain stress-responsive cis-acting elements. Under saline-alkali stress, qRT-PCR analysis showed that SiAAP3, SiLHT2, and SiAAP16 were differentially expressed between salt-alkali tolerant millet variety JK3 and salt-alkali sensitive millet variety B175. These results suggest that these genes may be involved in or regulate the response to saline-alkali stress, providing a theoretical basis for further studying the function of SiAAAPs.

[Identification and expression analysis of TCP family members in tobacco (Nicotiana tabacum L.)].

TCP family as plant specific transcription factor, plays an important role in different aspects of plant development. In order to screen TCP family members in tobacco, the homologous sequences of tobacco and Arabidopsis TCP family were identified by genome-wide homologous alignment. The physicochemical properties, phylogenetic relationships and cis-acting elements were analyzed by bioinformatics. The homologous genes of AtTCP3/AtTCP4 were screened, and RT-qPCR was used to detect the changes of gene expression upon 20% PEG6000 treatment. The results show that tobacco contains 63 TCP family members. Their amino acid sequence length ranged from 89 aa to 596 aa, and their protein hydropathicity grand average of hydropathicity (GRAVY) ranged from -1.147 to 0.125. The isoelectric point (pI) ranges from 4.42 to 9.94, the number of introns is 0 to 3, and the subcellular location is all located in the nucleus. The results of conserved domain and phylogenetic relationship analysis showed that the tobacco TCP family can be divided into PCF, CIN and CYC/TB1 subfamilies, and each subfamily has a stable sequence. The results of cis-acting elements in gene promoter region showed that TCP family genes contain low docile acting elements (LTR) and a variety of stress and metabolic regulation related elements (MYB, MYC). Analysis of gene expression patterns showed that AtTCP3/AtTCP4 homologous genes (NtTCP6, NtTCP28, NtTCP30, NtTCP33, NtTCP42, NtTCP57, NtTCP63) accounted for 20% PEG6000 treatment significantly up-regulated/down-regulated expression, and NtTCP30 and NtTCP57 genes were selected as candidate genes in response to drought. The results of this study analyzed the TCP family in the tobacco genome and provided candidate genes for the study of drought-resistance gene function and variety breeding in tobacco.

Chloroplast Genes Are Involved in The Male-Sterility of K-Type CMS in Wheat.

The utilization of crop heterosis can greatly improve crop yield. The sterile line is vital for the heterosis utilization of wheat (Triticum aestivum L.). The chloroplast genomes of two sterile lines and one maintainer were sequenced using second-generation high-throughput technology and assembled. The nonsynonymous mutated genes among the three varieties were identified, the expressed difference was further analyzed by qPCR, and finally, the function of the differentially expressed genes was analyzed by the barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) method. A total of 16 genes containing 31 nonsynonymous mutations between K519A and 519B were identified. There were no base mutations in the protein-encoding genes between K519A and YS3038. The chloroplast genomes of 519B and K519A were closely related to the Triticum genus and Aegilops genus, respectively. The gene expression levels of the six selected genes with nonsynonymous mutation sites for K519A compared to 519B were mostly downregulated at the binucleate and trinucleate stages of pollen development. The seed setting rates of atpB-silenced or ndhH-silenced 519B plants by BSMV-VIGS method were significantly reduced. It can be concluded that atpB and the ndhH are likely to be involved in the reproductive transformation of 519B.

Dynamic expression of miRNAs and functional analysis of target genes involved in the response to male sterility of the wheat line YS3038.

Thermosensitive cytoplasmic male sterile (TCMS) lines play an important role in wheat breeding, heterosis utilization, and germplasm innovation. MicroRNAs (miRNAs) can regulate the expression level of target genes by inhibiting the translation of these genes. YS3038 is a wheat TCMS line. In this study, the fertility conversion mechanism of YS3038 was studied by examining the abortion characteristics of YS3038, the regulation pattern of miRNAs and the target genes of miRNAs in YS3038. MiRNA-seq was performed on three important stages of YS3038 under sterile and fertile conditions. Then, the clean reads were aligned with some databases to filter other ncRNAs and repeats. The known miRNAs and novel miRNAs were predicted by sequence comparison with known miRNAs from miRbase. Differential expression of miRNAs between different stages and between different fertile conditions was analyzed, and functional analysis of target genes with opposite expression patterns as those of the miRNAs was conducted. The Ubisch bodies and microspores of sterile anthers were covered with filamentous materials. The degradation of the tapetum cells, the chloroplast structure of endothecium cells, and the microspore structure were abnormal. Microspore development was hindered from the late uninucleate stage to the binucleate stage. Twenty, 52, and 68 differentially expressed miRNAs (DEmiRs) were identified at the early uninucleate, late uninucleate, and binucleate stages, respectively, and there were 0, 7, and 72 differentially expressed target genes (DETGs), respectively, at these three stages. At the binucleate stage, 29 DEmiRs had 41 target mRNAs in total, and the expression patterns of the 41 target mRNAs were opposite to those of the 29 miRNAs. Fifteen significantly enriched KEGG pathways were associated with the 41 target mRNAs. Leucine-rich repeat receptor-like kinases (LRR-RLKs) play important roles in plant developmental and physiological processes. Some studies have shown that the expression of LRR-RLKs is related to the differentiation of microsporocytes and tapetum cells and to male sterility. An LRR-RLK (TaeRPK) gene was silenced by the barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) method, and the seed setting rates of the TaeRPK-silenced plants (3.51%) were significantly lower than those of the negative control plants (88.78%) (P < 0.01). Thus, the TaeRPK gene is likely to be involved in the fertility conversion of YS3038.

Mitochondrial genes are involved in the fertility transformation of the thermosensitive male-sterile line YS3038 in wheat.

Heterosis can improve the stress resistance, quality, and yield of crops, and the male sterility of wheat can be utilized to accelerate the breeding process of hybrid. To determine whether mitochondrial genes are involved in the fertility of K-type cytoplasmic male-sterile (CMS) line and the YS-type thermosensitive male-sterile (TMS) line in wheat, we sequenced and assembled the mitochondrial genomes of K519A, 519B, and YS3038 by next-generation sequencing (NGS). The non-synonymous mutations were analyzed, and the first-generation sequencing was conducted to verify the non-synonymous mutation sites. Furthermore, the expression patterns of genes with non-synonymous mutations were analyzed. Finally, the candidate genes were silenced by barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) to test the functions of the candidate genes. The results revealed that the mitochondrial genomes of K519A, 519B, and YS3038 were 420,543, 433,560, and 452,567 bp in length, respectively. Besides, 33, 31, and 37 protein-coding genes were identified in K519A, 519B, and YS3038, respectively. There were 14 protein-coding genes and 83 open reading frame (ORF) sequences that differed between K519A and 519B and 10 protein-coding genes and 122 ORF sequences that differed between K519A and YS3038. At the binucleate stage, seven genes (nad6, ORF256, ORF216, ORF138, atp6, nad3, and cox1) were downregulated in K519A compared with 519B, and 10 genes (nad6, atp6, cox3, atp8, nad3, cox1, rps3, ORF216, ORF138, and ORF224) were downregulated in YS3038 compared with K519A. Besides, six genes (nad6, ORF138, cox3, cox1, rps3, and ORF224) were downregulated under fertile conditions relative to sterile conditions in YS3038. Gene silencing analysis showed that the silencing of cox1 significantly reduced the seed setting rate of YS3038, indicating that the cox1 gene may be involved in the fertility transformation of YS3038. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-021-01252-x.

Identification of ceRNA and candidate genes related to fertility conversion of TCMS line YS3038 in wheat.

Previous studies have indicated that noncoding RNAs are important factors in gene functions. To explore the mechanism of male sterility of YS3038, the sterile genes were mapped, and based on previous work, the expression of long noncoding RNAs (lncRNAs), circular RNAs (circRNAs), and their target genes was studied. Weighted gene coexpression network analysis (WGCNA) and competitive endogenous RNA (ceRNA) analysis were further performed for differentially expressed noncoding RNAs and target genes. At last, the candidate genes were silenced by barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) to prove their function. The sterile genes were mapped on chromosomes 1B and 6B based on chip mix pool analysis, and one major effect QTL (27.3190% variation) was found based on SSR primers. The WGCNA analysis revealed that the dark turquoise and steel blue modules were highly correlated with anther development and fertility conversion, respectively. The ceRNA analysis showed that a total of 184 RNAs interacted with each other, including 115 mRNAs, 55 microRNAs (miRNAs), eight circRNAs, and six lncRNAs. Finally, the seed setting rate of the plant was significantly decreased after fatty acyl-CoA reductase 5 silencing. This study provides breeders with a new option for the development of thermosensitive cytoplasmic male-sterile (TCMS) wheat lines, which will favor the sustainable development of two-line hybrid wheat.

Genome-wide identification, phylogeny and expression analysis of the SPL gene family in wheat.

BACKGROUND: Members of the plant-specific SPL gene family (squamosa promoter-binding protein -like) contain the SBP conserved domain and are involved in the regulation of plant growth and development, including the development of plant flowers and plant epidermal hair, the plant stress response, and the synthesis of secondary metabolites. This family has been identified in various plants. However, there is no systematic analysis of the SPL gene family at the genome-wide level of wheat. RESULTS: In this study, 56 putative TaSPL genes were identified using the comparative genomics method; we renamed them TaSPL001 - TaSPL056 on their chromosomal distribution. According to the un-rooted neighbor joining phylogenetic tree, gene structure and motif analyses, the 56 TaSPL genes were divided into 8 subgroups. A total of 81 TaSPL gene pairs were designated as arising from duplication events and 64 interacting protein branches were identified as involve in the protein interaction network. The expression patterns of 21 randomly selected TaSPL genes in different tissues (roots, stems, leaves and inflorescence) and under 4 treatments (abscisic acid, gibberellin, drought and salt) were detected by quantitative real-time polymerase chain reaction (qRT-PCR). CONCLUSIONS: The wheat genome contains 56 TaSPL genes and those in same subfamily share similar gene structure and motifs. TaSPL gene expansion occurred through segmental duplication events. Combining the results of transcriptional and qRT-PCR analyses, most of these TaSPL genes were found to regulate inflorescence and spike development. Additionally, we found that 13 TaSPLs were upregulated by abscisic acid, indicating that TaSPL genes play a positive role in the abscisic acid-mediated pathway of the seedling stage. This study provides comprehensive information on the SPL gene family of wheat and lays a solid foundation for elucidating the biological functions of TaSPLs and improvement of wheat yield.

Differential expression pattern of the proteome in response to cadmium stress based on proteomics analysis of wheat roots.

BACKGROUND: Heavy metal cadmium (Cd) is a common environmental pollutant in soils, which has an negative impacts on crop growth and development. At present, cadmium has become a major soil and water heavy metal pollutant, which not only causes permanent and irreversible health problems for humans, but also causes a significant reduction in crop yields. RESULTS: This study examined the chemical forms of Cd in the roots of two wheat varieties (M1019 and Xinong20) by continuous extraction and analyzed differences in distribution characteristics of Cd in the root cell wall, cytoplasm, and organelles by elemental content determination and subcellular separation. Furthermore, we conducted proteomics analysis of the roots of the two varieties under Cd pollution using mass spectrometry quantitative proteomics techniques. A total of 11,651 proteins were identified, of which 10,532 proteins contained quantitative information. In addition, the differentially expressed proteins in the two varieties were related to DNA replication and repair, protein metabolism, and the glutathione metabolism pathway. CONCLUSION: The results of this study improve our understanding of the mechanism of plant responses to Cd stress.

QTL Mapping of Kernel Traits and Validation of a Major QTL for Kernel Length-Width Ratio Using SNP and Bulked Segregant Analysis in Wheat.

One RIL population derived from the cross between Dalibao and BYL8 was used to examine the phenotypes of kernel-related traits in four different environments. Six important kernel traits, kernel length (KL), kernel width (KW), kernel perimeter (KP), kernel area (KA), kernel length/width ratio (KLW), and thousand-kernel weight (TKW) were evaluated in Yangling, Shaanxi Province, China (2016 and 2017), Nanyang, Henan Province, China (2017) and Suqian, Jiangsu Province, China (2017). A genetic linkage map was constructed using 205 SSR markers, and a total of 21 significant QTLs for KL, KW, KP, KA, KLW and TKW were located on 10 of the 21 wheat chromosomes, including 1A, 1B, 2A, 2B, 2D, 3D, 4D, 5A, 5B, and 7D, with a single QTL in different environments explaining 3.495-30.130% of the phenotypic variation. There were four loci for KLW, five for KA, five for KL, three for KP, two for KW, and two for TKW among the detected QTLs. We used BSA + 660 K gene chip technology to reveal the positions of major novel QTLs for KLW. A total of 670 out of 5285 polymorphic SNPs were detected on chromosome 2A. The SNPs in 2A are most likely related to the major QTL, and there may be minor QTLs on 5B, 7A, 3A and 4B. SSR markers were developed to verify the chromosome region associated with KLW. A linkage map was constructed with 7 SSR markers, and a major effect QTL was identified within a 21.55 cM interval, corresponding to a physical interval of 10.8 Mb in the Chinese Spring RefSeq v1.0 sequence. This study can provide useful information for subsequent construction of fine mapping and marker-assisted selection breeding.

Energy metabolism involved in fertility of the wheat TCMS line YS3038.

MAIN CONCLUSION: In the wheat TCMS line YS3038, the anther development is inhibited from late uninucleate stage to the binucleate stage. The disruption of energy metabolism pathways by aberrant transcriptional regulation causes the male sterility under low temperatures. The utilization of thermosensitive male sterile (TMS) lines provides a basis for two-line breeding. Previous work, including morphological and cytological observations, has shown that the development process of the TMS line YS3038 is inhibited from the late uninucleate stage to the binucleate stage. Transcriptomics studies could now help to elucidate the overall expression of related genes in a specific reproductive process, revealing the metabolic network and its regulatory mechanism of the reproductive process from the transcription level. Considering the fertility characteristics of YS3038, three important stages for transcriptome analysis were determined to be the early uninucleate, late uninucleate and binucleate stages. The number of differentially expressed genes (DEGs) was found to be highest in the binucleate stage, and most were related to energy metabolism. Quantitative PCR analysis of selected genes related to energy metabolism revealed that their expression patterns were consistent with the sequencing results. Analysis of the fertility mechanism of YS3038 showed that although the tapetum of anthers was degraded in advance of the tetrad stage, the development of microspores did not result in obvious abnormalities until the binucleate stage, because the genes involved in energy metabolism pathways, including starch and sucrose metabolism (SSM), glycolysis, the tricarboxylic acid (TCA) cycle, oxidative phosphorylation, and respiration electron transport chain are differentially expressed under sterile and fertile conditions. Therefore, the pollen in YS3038 was sterile.

Genome-wide identification, phylogeny and expression profiling of class III peroxidases gene family in Brachypodium distachyon.

Class III peroxidases are classical secretory plant peroxidases belonging to a large multi-gene family. Class III peroxidases are involved in various physical processes and the response to biotic and abiotic stress to protect plants from environmental adversities. In this study, 151 BdPrx genes were identified using HMM and Blastp program. According to their physical location, the 151 BdPrx genes were mapped on five chromosomes. The results of Gene Structure Display Serve and MEME revealed that BdPrxs in the same subgroup shared similar gene structure, and their protein sequences were highly conserved. Based on the analysis of evolutionary relationships and Ka/Ks, 151 BdPrx genes were divided into 15 subgroups, they have undergone purifying selection. In addition, the result of GO annotation showed that 100% of the BdPrxs participated in antioxidant. The protein-protein interaction network was constructed using the orthology-based method, found that 66 BdPrxs were involved in the regulatory network and 183 network branches were identified. Furthermore, analysis of the transcriptome data indicated that the BdPrx genes responded to low concentration of exogenous phytohormones and exhibited different levels of expression in the different tissues. Subsequently, 19 genes were selected for quantitative real-time PCR and found to be mainly expressed in the roots, might preferentially respond to hydrogen peroxide and gibberellin. Our results provide a foundation for further evolutionary and functional study of Prx genes in B. distachyon.

Construction of a high-density genetic map using specific-locus amplified fragments in sorghum.

BACKGROUND: Sorghum is mainly used as a human food and beverage source, playing an important role in the production of ethanol and other bio-industrial products. Thus it is regarded as a model crop for energy plants. Genetic map construction is the foundation for marker-assisted selection and gene cloning. So far several sorghum linkage maps have been reported using different kinds of molecular markers. However marker numbers and chromosome coverage are limited. As a result, it is difficult to get consistent results and the maps are hard to unify. In the present study, the genomes of 130 individuals consisting an F2 population together with their parents were surveyed using a high-throughput sequencing technique. A high-density linkage map was constructed using specific-locus amplified fragments (SLAF) markers. This map can provide information and serve as a reference for effective gene exploration, and for marker assisted-breeding program. RESULTS: A high-throughput sequencing method was adopted to screen SLAF markers with 130 F2 individuals from a cross between a grain sorghum variety, J204, and a sweet sorghum variety, Keter. In the present study, 52,928 suitable SLAF markers out of 43,528,021 pair-end reads were chosen to conduct genetic map construction, 12.0% of which were polymorphic. Among the 6353 polymorphic SLAF markers, 5829 (91.8%) were successfully genotyped in the F2 mapping population. Finally 2246 SLAF markers were obtained to construct a high-density genetic linkage map. The total distance of linkage map covering all 10 chromosomes was 2158.1 cM. The largest gap on each chromosome was 10.2 cM on average. The proportion of gaps less than and/or equal to 5.0 cM was averagely 98.1%. The markers on each chromosome ranged from 123 (chromosome 9) to 315 (chromosome 4) with a mean value of 224.6, the distance between adjacent markers ranged from 0.6 (chromosome 10) to 1.3 cM (chromosome 9) with an average distance of only 0.98 cM. CONCLUSION: A high density sorghum genetic map was constructed in this study. The total length was 2158.1 cM covering all 10 chromosomes with a total number of 2246 SLAF markers. The construction of this map can provide detailed information for accurate gene localization and cloning and application of marker-assisted breeding.

Combining Next Generation Sequencing with Bulked Segregant Analysis to Fine Map a Stem Moisture Locus in Sorghum (Sorghum bicolor L. Moench).

Sorghum is one of the most promising bioenergy crops. Stem juice yield, together with stem sugar concentration, determines sugar yield in sweet sorghum. Bulked segregant analysis (BSA) is a gene mapping technique for identifying genomic regions containing genetic loci affecting a trait of interest that when combined with deep sequencing could effectively accelerate the gene mapping process. In this study, a dry stem sorghum landrace was characterized and the stem water controlling locus, qSW6, was fine mapped using QTL analysis and the combined BSA and deep sequencing technologies. Results showed that: (i) In sorghum variety Jiliang 2, stem water content was around 80% before flowering stage. It dropped to 75% during grain filling with little difference between different internodes. In landrace G21, stem water content keeps dropping after the flag leaf stage. The drop from 71% at flowering time progressed to 60% at grain filling time. Large differences exist between different internodes with the lowest (51%) at the 7th and 8th internodes at dough stage. (ii) A quantitative trait locus (QTL) controlling stem water content mapped on chromosome 6 between SSR markers Ch6-2 and gpsb069 explained about 34.7-56.9% of the phenotypic variation for the 5th to 10th internodes, respectively. (iii) BSA and deep sequencing analysis narrowed the associated region to 339 kb containing 38 putative genes. The results could help reveal molecular mechanisms underlying juice yield of sorghum and thus to improve total sugar yield.

Molecular mapping of the brace root traits in sorghum (Sorghum bicolor L. Moench).

The presence and morphology of plant brace roots are important root architecture traits. Brace roots contribute significantly to effective anchorage and water and nutrient uptake during late growth and development, and more importantly, have a substantial influence on grain yield under soil flooding or water limited conditions. However, little is known about the genetic mechanisms that underlie brace root traits. In this study, quantitative trait loci (QTLs) for presence of brace roots from the sorghum landrace "Sansui" were mapped and associated molecular markers were identified. A linkage map was constructed with 109 assigned simple sequence repeat markers using a F2 mapping population derived from the cross Sansui/Jiliang 2. Two QTLs associated with presence of brace roots were localized on chromosomes 6 and 7. The major QTL on chromosome 7 between markers Dsenhsbm7 and Xcup 70 explained about 52.5% of the phenotypic variation, and the minor QTL on chromosome 6 was flanked by Xtxp127 and Xtxp6 and accounted for 7.0% of phenotypic variation. These results will provide information for the improvement of sorghum root architecture associated with brace roots.

The preparation of 3,5-dihydroxy-4-isopropylstilbene nanoemulsion and in vitro release.

We have reported a novel procedure to prepare 3,5-dihydroxy-4-isopropylstilbene (DHPS) nanoemulsion, using a low-energy emulsification method. Based on the phase diagram, the optimum prescription of nanoemulsion preparation was screened. With polyoxyethylenated castor oil (EL-40) as the surfactant, ethanol as the co-surfactant, and isopropyl myristate (IPM) as the oil phase, the DHPS nanoemulsion was obtained with a transparent appearance, little viscosity, and spherically uniform distribution verified by transmission electron microscopy and laser scattering analyzer. The nanoemulsion was also determined by FT-Raman spectroscopy. The DHPS nanoemulsion demonstrated good stability and stable physical and chemical properties. The nanoemulsion dramatically improved the transdermal release of DHPS (from 8.02 µg · cm(-2) to 273.15 µg · cm(-2)) and could become a favorable new dosage form for DHPS.