Genome-Wide Identification of Common Bean PvLTP Family Genes and Expression Profiling Analysis in Response to Drought Stress.

Common bean is one of the most important legume crops for human consumption. Its yield is adversely affected by environmental stress. Plant non-specific lipid transfer proteins (nsLTPs) are essential for plant growth, development, and resistance to abiotic stress, such as salt, drought, and alkali. However, changes in nsLTP family genes responding to drought stress are less known. The PvLTP gene family in the common bean was identified by a comprehensive genome-wide analysis. Molecular weights, theoretical isoelectric points, phylogenetic tree, conserved motifs, gene structures, gene duplications, chromosome localization, and expression profiles were analyzed by SignalP 5.0, ExPASy, ClustalX 2.1, MEGA 7.0, NCBI-CDD, MEME, Weblogo, and TBtools 1.09876, respectively. Heatmap and qRT-PCR analyses were performed to validate the expression profiles of PvLTP genes in different organs. In addition, the expression patterns of nine PvLTP genes in common beans treated with drought stress were investigated by qRT-PCR. We obtained 58 putative PvLTP genes in the common bean genome via genome-wide analyses. Based on the diversity of the eight-cysteine motif (ECM), these genes were categorized into five types (I, II, IV, V, and VIII). The signal peptides of the PvLTP precursors were predicted to be from 16 to 42 amino acid residues. PvLTPs had a predicated theoretical isoelectric point of 3.94-10.34 and a molecular weight of 7.15-12.17 kDa. The phylogenetic analysis showed that PvLTPs were closer to AtLTPs than OsLTPs. Conserved motif and gene structure analyses indicated that PvLTPs were randomly distributed on all chromosomes except chromosome 9. In addition, 23 tandem duplicates of PvLTP genes were arranged in 10 gene clusters on chromosomes 1 and 2. The heatmap and qRT-PCR showed that PvLTP expression significantly varied in different tissues. Moreover, 9 PvLTP genes were up-regulated under drought treatment. Our results reveal that PvLTPs play potentially vital roles in plants and provide a comprehensive reference for studies on PvLTP genes and a theoretical basis for further analysis of regulatory mechanisms influencing drought tolerance in the common bean.

Construction of A GBS-Based High-Density Genetic Map and Flower Color-Related Loci Mapping in Grasspea (Lathyrus sativus L.).

Grasspea (Lathyrus sativus L.), a legume crop with excellent resistance to a broad array of environmental stressors, has, to this point, been poorly genetically characterized. High-density genetic linkage maps are critical for draft genome assembly, quantitative trait loci (QTLs) analysis, and gene mining. The lack of a high-density genetic linkage map has limited both genomic studies and selective breeding in grasspea. Here, we developed a high-density genetic linkage map of grasspea using genotyping-by-sequencing (GBS) to sequence 154 grasspea plants, comprising 2 parents and 152 F2 progeny. In all, 307.74 Gb of data was produced, including 2,108,910,938 paired-end reads, as well as 3536 SNPs mapped to seven linkage groups (LG1-LG7). With an average length of 996.52 cM per LG, the overall genetic distance was 6975.68 cM. Both the χ2 test and QTL analysis, based on the Kruskal-Wallis (KW) test and interval mapping (IM) analysis, revealed the monogenic inheritance of flower color in grasspea, with the responsible QTL located between 308.437 cM and 311.346 cM in LG4. The results can aid grasspea genome assembly and accelerate the selective breeding of new grasspea germplasm resources.

An RNA Sequencing Transcriptome Analysis of Grasspea (Lathyrus sativus L.) and Development of SSR and KASP Markers.

Grasspea (Lathyrus sativus L., 2n = 14) has great agronomic potential because of its ability to survive under extreme conditions, such as drought and flood. However, this legume is less investigated because of its sparse genomic resources and very slow breeding process. In this study, 570 million quality-filtered and trimmed cDNA sequence reads with total length of over 82 billion bp were obtained using the Illumina NextSeqTM 500 platform. Approximately two million contigs and 142,053 transcripts were assembled from our RNA-Seq data, which resulted in 27,431 unigenes with an average length of 1,250 bp and maximum length of 48,515 bp. The unigenes were of high-quality. For example, the stay-green (SGR) gene of grasspea was aligned with the SGR gene of pea with high similarity. Among these unigenes, 3,204 EST-SSR primers were designed, 284 of which were randomly chosen for validation. Of these validated unigenes, 87 (30.6%) EST-SSR primers produced polymorphic amplicons among 43 grasspea accessions selected from different geographical locations. Meanwhile, 146,406 SNPs were screened and 50 SNP loci were randomly chosen for the kompetitive allele-specific PCR (KASP) validation. Over 80% (42) SNP loci were successfully transformed to KASP markers. Comparison of the dendrograms according to the SSR and KASP markers showed that the different marker systems are partially consistent with the dendrogram constructed in our study.

Soil bacterial diversity changes in different broomcorn millet intercropping systems.

Plants growing in soil and the diverse microorganisms with which they are in direct contact have adapted to exploit their close association for mutual benefit. Various intercropping systems have been used to control plant disease and improve productivity in fields. Although high-throughput sequencing approaches have provided new insights into the soil bacterial community, current knowledge of intercropping of broomcorn millet with different leguminous plants is limited. In this study, characterization of different bacterial communities of monoculture and intercropping systems was achieved by deep sequencing. A total of 4684 operational taxonomic units were classified to the species level with good sampling depth and sequencing coverage. The abundance of Actinobacteria, Bacteroidetes, Proteobacteria, Acidobacteria, and Gemmatimonadetes varied at different growth stages and was related to growth of the intercropped plant. According to diversity analyses, Glycomyces, Aeromicrobium, Adhaeribacter, and Streptomyces were the dominant genera. In addition, we predicted functional gene composition based on bacterial OTUs present. Functional results showed that membrane transport and nutrient metabolism was highly abundant in all samples, although abundance varied at different growth stages, which indicated these pathways might be affected by the dominant categories of bacterial community. The dynamic changes observed during intercropping of broomcorn millet with different leguminous plants suggest that soil bacterial community structure exhibits a crop species-specific pattern. Further, agronomic trait data from different broomcorn millet intercropping systems were consistent with functional results and suggest that agronomic traits may be influenced by soil bacterial communities.

Large-scale microsatellite development in grasspea (Lathyrus sativus L.), an orphan legume of the arid areas.

BACKGROUND: Grasspea (Lathyrus sativus L., 2n = 14), a member of the family Leguminosae, holds great agronomic potential as grain and forage legume crop in the arid areas for its superb resilience to abiotic stresses such as drought, flood and salinity. The crop could not make much progress through conventional breeding in the past, and there are hardly any detailed molecular biology studies due to paucity of reliable molecular markers representative of the entire genome. RESULTS: Using the 454 FLX Titanium pyrosequencing technique, 651,827 simple sequence repeat (SSR) loci were identified and 50,144 nonredundant primer pairs were successfully designed, of which 288 were randomly selected for validation among 23 L. sativus and one L. cicera accessions of diverse provenance. 74 were polymorphic, 70 monomorphic, and 144 with no PCR product. The number of observed alleles ranged from two to five, the observed heterozygosity from 0 to 0.9545, and Shannon's information index ranged from 0.1013 to 1.0980, respectively. The dendrogram constructed by using unweighted pair group method with arithmetic mean (UPGMA) based on Nei's genetic distance, showed obvious distinctions and understandable relationships among the 24 accessions. CONCLUSIONS: The large number of SSR primer pairs developed in this study would make a significant contribution to genomics enabled improvement of grasspea.