Xinan Zhou

Elucidation of miRNAs-Mediated Responses to Low Nitrogen Stress by Deep Sequencing of Two Soybean Genotypes

Nitrogen (N) is a major limiting factor in crop production, and plant adaptive responses to low N are involved in many post-transcriptional regulation. Recent studies indicate that miRNAs play important roles in adaptive responses. However, miRNAs in soybean adaptive responses to N limitation have been not reported. We constructed sixteen libraries to identify low N-responsive miRNAs on a genome-wide scale using samples from 2 different genotypes (low N sensitive and low N tolerant) subjected to various periods of low nitrogen stress. Using high-throughput sequencing technology (Illumina-Solexa), we identified 362 known miRNAs variants belonging to 158 families and 90 new miRNAs belonging to 55 families. Among these known miRNAs variants, almost 50% were not different from annotated miRNAs in miRBase. Analyses of their expression patterns showed 150 known miRNAs variants as well as 2 novel miRNAs with differential expressions. These differentially expressed miRNAs between the two soybean genotypes were compared and classified into three groups based on their expression patterns. Predicted targets of these miRNAs were involved in various metabolic and regulatory pathways such as protein degradation, carbohydrate metabolism, hormone signaling pathway, and cellular transport. These findings suggest that miRNAs play important roles in soybean response to low N and contribute to the understanding of the genetic basis of differences in adaptive responses to N limitation between the two soybean genotypes. Our study provides basis for expounding the complex gene regulatory network of these miRNAs.

Transcript Profile of the Response of Two Soybean Genotypes to Potassium Deficiency

The macronutrient potassium (K) is essential to plant growth and development. Crop yield potential is often affected by lack of soluble K. The molecular regulation mechanism of physiological and biochemical responses to K starvation in soybean roots and shoots is not fully understood. In the present study, two soybean varieties were subjected to low-K stress conditions: a low-K-tolerant variety (You06-71) and a low-K-sensitive variety (HengChun04-11). Eight libraries were generated for analysis: 2 genotypes ×2 tissues (roots and shoots) ×2 time periods [short term (0.5 to 12 h) and long term (3 to 12 d)]. RNA derived from the roots and shoots of these two varieties across two periods (short term and long term) were sequenced and the transcriptomes were compared using high-throughput tag-sequencing. To this end, a large number of clean tags (tags used for analysis after removal of dirty tags) corresponding to distinct tags (all types of clean tags) were identified in eight libraries (L1, You06-71-root short term; L2, HengChun04-11-root short term; L3, You06-71-shoot short term; L4, HengChun04-11-shoot short term; L5, You06-71-root long term; L6, HengChun04-11-root long term; L7, You06-71-shoot long term; L8, HengChun04-11-shoot long term). All clean tags were mapped to the available soybean (Glycine max) transcript database (http://www.soybase.org). Many genes showed substantial differences in expression across the libraries. In total, 5,440 transcripts involved in 118 KEGG pathways were either up- or down-regulated. Fifteen genes were randomly selected and their expression levels were confirmed using quantitative RT-PCR. Our results provide preliminary information on the molecular mechanism of potassium absorption and transport under low-K stress conditions in different soybean tissues.

Identification of Glycine Max MicroRNAs in response to phosphorus deficiency

MicroRNAs (miRNAs) are endogenous small RNAs regulating plant development and stress responses. In addition, phosphorus (P) is an important macronutrient for plant growth and development. More than two hundred miRNAs have been identified in Glycine Max and a few of miRNAs have been shown to respond to P deficiency, however, whether there are other miRNAs involved in P deficiency response is largely unknown. In this study, we used high-throughput small RNA sequencing and whole-genome-wide mining to identify the potential miRNAs in response to P deficiency. After sequencing, we deduced 183 known, 99 conserved and 126 novel miRNAs in Glycine Max. Among them, in response to P deficiency, the expressions of 27 known, 16 conserved and 12 novel miRNAs showed significant changes in roots, whereas the expressions of 34 known, 14 conserved and 7 novel miRNAs were significantly different in shoots. Furthermore, we validated the predicated novel miRNAs and found that three miRNAs in roots and five miRNAs in shoots responded to P deficiency. Some miRNAs were P-induced whereas some were P-suppressed. Together these results indicated that the miRNAs identified might play important roles in regulating P signaling pathway.

Genetic analysis and molecular mapping of resistance gene to Phakopsora pachyrhizi in soybean germplasm SX6907

Key messageIn this study,Rpp6907,a novel resistance gene/allele toPhakopsora pachyrhiziin soybean, was mapped in a 111.9-kb region, including three NBS-LRR type predicted genes, on chromosome 18.AbstractSoybean rust caused by Phakopsora pachyrhizi Sydow has been reported in numerous soybean-growing regions worldwide. The development of rust-resistant varieties is the most economical and environmentally safe method to control the disease. The Chinese soybean germplasm SX6907 is resistant to P. pachyrhizi and exhibits immune reaction compared with the known Rpp genes. These characteristics suggest that SX6907 may carry at least one novel Rpp gene/allele. Three F2 populations from the crosses of SX6907 (resistant) and Tianlong 1, Zhongdou40, and Pudou11 (susceptible) were used to map the Rpp gene. Three resistance responses (immune, red-brown, and tan-colored lesion) were observed from the F2 individuals. The segregation follows a ratio of 1(resistance):2(heterozygous):1(susceptible), indicating that the resistance in SX6907 is controlled by a single incomplete dominant gene (designated as Rpp6907). Results showed that Rpp6907 was mapped on soybean chromosome 18 (molecular linkage group G, MLG G) flanked by simple sequence repeat (SSR) markers SSR24 and SSR40 at a distance of 111.9 kb. Among the ten genes marked within this 111.9-kb region between the two markers, three genes (Glyma18g51930, Glyma18g51950, and Glyma18g51960) are nucleotide-binding site and leucine-rich repeat-type genes. These genes may be involved in recognizing the presence of pathogens and ultimately conferring resistance. Based on resistance spectrum analysis and mapping results, we inferred that Rpp6907 is a novel gene different from Rpp1 in PI 200492, PI 561356, PI 587880A, PI 587886, and PI 594538A, or a new Rpp1-b allele.

Large-scale sequencing of normalized full-length cDNA library of soybean seed at different developmental stages and analysis of the gene expression profiles based on ESTs

Although GenBank has now covered over 1,400,000 expressed sequence tags (ESTs) from soybean, most ESTs available to the public have been derived from tissues or environmental conditions rather than developing seeds. It is absolutely necessary for annotating the molecular mechanisms of soybean seed development to analyze completely the gene expression profiles of its immature seed at various stages. Here we have constructed a full-length-enriched cDNA library comprised of a total of 45,408 cDNA clones which cover various stages of soybean seed development. Furthermore, we have sequenced from 5′ ends of these clones, 36,656 ESTs were obtained in the present study. These EST sequences could be categorized into 27,982 unigenes, including 22,867 contigs and 5,115 singletons, among which 27,931 could be mapped onto soybean 20 chromosome sequences. Comparative genomic analysis with other plants has revealed that these unigenes include lots of candidate genes specific to dicot, legume and soybean. Approximately 1,789 of these unigenes currently show no homology to known soybean sequences, suggesting that many represent mRNAs specifically expressed in seeds. Novel abundant genes involved in the oil synthesis have been found in this study, may serve as a valuable resource for soybean seed improvement.

RNA-Seq Analysis of Differential Gene Expression Responding to Different Rhizobium Strains in Soybean (Glycine max) Roots

The root nodule symbiosis (RNS) between legume plants and rhizobia is the most efficient and productive source of nitrogen fixation, and has critical importance in agriculture and mesology. Soybean (Glycine max), one of the most important legume crops in the world, establishes a nitrogen-fixing symbiosis with different types of rhizobia, and the efficiency of symbiotic nitrogen fixation in soybean greatly depends on the symbiotic host-specificity. Although, it has been reported that rhizobia use surface polysaccharides, secretion proteins of the type-three secretion systems and nod factors to modulate host range, the host control of nodulation specificity remains poorly understood. In this report, the soybean roots of two symbiotic systems (Bradyrhizobium japonicum strain 113-2-soybean and Sinorhizobium fredii USDA205-soybean)with notable different nodulation phenotypes and the control were studied at five different post-inoculation time points (0.5, 7–24 h, 5, 16, and 21 day) by RNA-seq (Quantification). The results of qPCR analysis of 11 randomly-selected genes agreed with transcriptional profile data for 136 out of 165 (82.42%) data points and quality assessment showed that the sequencing library is of quality and reliable. Three comparisons (control vs. 113-2, control vs. USDA205 and USDA205 vs. 113-2) were made and the differentially expressed genes (DEGs) between them were analyzed. The number of DEGs at 16 days post-inoculation (dpi) was the highest in the three comparisons, and most of the DEGs in USDA205 vs. 113-2 were found at 16 dpi and 21 dpi. 44 go function terms in USDA205 vs. 113-2 were analyzed to evaluate the potential functions of the DEGs, and 10 important KEGG pathway enrichment terms were analyzed in the three comparisons. Some important genes induced in response to different strains (113-2 and USDA205) were identified and analyzed, and these genes primarily encoded soybean resistance proteins, NF-related proteins, nodulins and immunity defense proteins, as well as proteins involving flavonoids/flavone/flavonol biosynthesis and plant-pathogen interaction. Besides, 189 candidate genes are largely expressed in roots and\or nodules. The DEGs uncovered in this study provides molecular candidates for better understanding the mechanisms of symbiotic host-specificity and explaining the different symbiotic effects between soybean roots inoculated with different strains (113-2 and USDA205).

Genome-Wide Survey of the Soybean GATA Transcription Factor Gene Family and Expression Analysis under Low Nitrogen Stress

GATA transcription factors are transcriptional regulatory proteins that contain a characteristic type-IV zinc finger DNA-binding domain and recognize the conserved GATA motif in the promoter sequence of target genes. Previous studies demonstrated that plant GATA factors possess critical functions in developmental control and responses to the environment. To date, the GATA factors in soybean (Glycine max) have yet to be characterized. Thus, this study identified 64 putative GATA factors from the entire soybean genomic sequence. The chromosomal distributions, gene structures, duplication patterns, phylogenetic tree, tissue expression patterns, and response to low nitrogen stress of the 64 GATA factors in soybean were analyzed to further investigate the functions of these factors. Results indicated that segmental duplication predominantly contributed to the expansion of the GATA factor gene family in soybean. These GATA proteins were phylogenetically clustered into four distinct subfamilies, wherein their gene structure and motif compositions were considerably conserved. A comparative phylogenetic analysis of the GATA factor zinc finger domain sequences in soybean, Arabidopsis (Arabidopsis thaliana), and rice (Oryza sativa) revealed four major classes. The GATA factors in soybean exhibited expression diversity among different tissues; some of these factors showed tissue-specific expression patterns. Numerous GATA factors displayed upregulation or downregulation in soybean leaf in response to low nitrogen stress, and two GATA factors GATA44 and GATA58 were likely to be involved in the regulation of nitrogen metabolism in soybean. Overexpression of GmGATA44 complemented the reduced chlorophyll phenotype of the Arabidopsis ortholog AtGATA21 mutant, implying that GmGATA44 played an important role in modulating chlorophyll biosynthesis. Overall, our study provides useful information for the further analysis of the biological functions of GATA factors in soybean and other crops.

Analysis of QTLs for Root Traits at Seedling Stage in Soybean

Abstract Very few quantitative trait loci (QTLs) for root traits have been identified in soybean [ Glycine max (L.) Merr.], although root system plays important roles in the growth and development of a plant. The aim of this study was to localize QTLs for root traits at seedling stage in soybean. The mapping population contained 165 recombinant inbred lines (RILs), which was derived from the cross between Zhongdou 29 and Zhongdou 32. The seedlings were hydroponically cultured, and the root traits, including length of the main root, number of lateral roots, root weight, root volume, plant weight, shoot weight per plant, and ratio of root weight to shoot weight, were measured at V2 stage. Using composite interval mapping (CIM) method, 20 QTLs associated with roots and aerial part traits were mapped on 8 chromosomes, which included 9 major-effect QTLs on chromosomes 11 and 14 with the phenotypic contributions of 10.5-26.1%. On chromosomes 11 and 14, some QTLs for root traits and aerial part traits were distributed in clusters. This result is consistent with the correlation analysis of phenotypic traits, indicating that the root traits are associated with the shoot traits.

Search for Nodulation and Nodule Development-Related Cystatin Genes in the Genome of Soybean (Glycine max)

Nodulation, nodule development and senescence directly affects nitrogen fixation efficiency, and previous studies have shown that inhibition of some cysteine proteases delay nodule senescence, so their nature inhibitors, cystatin genes, are very important in nodulation, nodule development, and senescence. Although several cystatins are actively transcribed in soybean nodules, their exact roles and functional diversities in legume have not been well explored in genome-wide survey studies. In this report, we performed a genome-wide survey of cystatin family genes to explore their relationship to nodulation and nodule development in soybean and identified 20 cystatin genes that encode peptides with 97–245 amino acid residues, different isoelectric points (pI) and structure characteristics, and various putative plant regulatory elements in 3000 bp putative promoter fragments upstream of the 20 soybean cystatins in response to different abiotic/biotic stresses, hormone signals, and symbiosis signals. The expression profiles of these cystatin genes in soybean symbiosis with rhizobium strain Bradyrhizobium japonicum strain 113-2 revealed that 7 cystatin family genes play different roles in nodulation as well as nodule development and senescence. However, these genes were not root nodule symbiosis (RNS)—specific and did not encode special clade cystatin protein with structures related to nodulation and nodule development. Besides, only two of these soybean cystatins were not upregulated in symbiosis after ABA treatment. The functional analysis showed that a candidate gene Glyma.15G227500 (GmCYS16) was likely to play a positive role in soybean nodulation. Besides, evolutionary relationships analysis divided the cystatin genes from Arabidopsis thaliana, Nicotiana tabacum, rice, barley and four legume plants into three groups. Interestingly, Group A cystatins are special in legume plants, but only include one of the above-mentioned 7 cystatin genes related to nodulation and nodule development. Overall, our results provide useful information or clues for our understanding of the functional diversity of legume cystatin family proteins in soybean nodulation and nodule development and for finding nodule-specific cysteine proteases in soybean.

Identification and Comparative Analysis of CBS Domain-Containing Proteins in Soybean (Glycine max) and the Primary Function of GmCBS21 in Enhanced Tolerance to Low Nitrogen Stress

Nitrogen is an important macronutrient required for plant growth, and is a limiting factor for crop productivity. Improving the nitrogen use efficiency (NUE) is therefore crucial. At present, the NUE mechanism is unclear and information on the genes associated with NUE in soybeans is lacking. cystathionine beta synthase (CBS) domain-containing proteins (CDCPs) may be implicated in abiotic stress tolerance in plants. We identified and classified a CBS domain–containing protein superfamily in soybean. A candidate gene for NUE, GmCBS21, was identified. GmCBS21 gene characteristics, the temporal expression pattern of the GmCBS21 gene, and the phenotype of GmCBS21 overexpression in transgenic Arabidopsis thaliana under low nitrogen stress were analyzed. The phenotypes suggested that the transgenic Arabidopsis thaliana seedlings performed better under the nitrogen-deficient condition. GmCBS21-overexpressing transgenic plants exhibit higher low nitrogen stress tolerance than WT plants, and this suggests its role in low nitrogen stress tolerance in plants. We conclude that GmCBS21 may serve as an excellent candidate for breeding crops with enhanced NUE and better yield.