Genome-wide identification and characterization of circRNAs in wheat tiller.

KEY MESSAGE: This study found that the intergenic circRNAs of wheat were more abundant than those of other plants. More importantly, a circRNA-mediated network associated with tillering was constructed for the first time. Circular RNAs (circRNAs) are a class of endogenous non-coding RNAs with covalently closed circular structures, which play an important role in transcriptional and post-transcriptional regulation. Tiller is an important agronomic trait that determines plant morphological architecture and affects spike number in wheat. However, no studies on the characteristics and functions of circRNAs involved in the regulation of wheat tiller. Here, we performed a genome-wide identification of circRNAs using ribosomal-depleted RNA-seq from wheat tiller of two pairs near-isogenic lines. A total of 686 circRNAs were identified and distributed on 21 chromosomes of wheat, of which 537 were novel circRNAs. Unlike other plants, the majority of these circRNAs (61.8%) were derived from intergenic regions. One circRNA-mediated network associated with tillering was constructed through weighted gene co-expression network analysis, including 323 circRNAs, 117 miRNAs, and 968 mRNAs. GO and pathway enrichment analysis of mRNAs suggested that these circRNAs are involved in cell cycle, ncRNA export from nucleus, developmental process, plant hormone signal transduction, MAPK signaling pathway, RNA degradation. Of these circRNAs, ten circRNAs are associated with known tillering/branching genes in rice or Arabidopsis thaliana, including OsCesA7, EBR1, DTE1, CRD1, LPA1, PAY1, LRK1, OsNR2, OsCCA1, OsBZR1. In summary, we present the first study of the identification and characterization of circRNAs in wheat tiller, and the results suggest these circRNAs associated with tillering could play an important role in wheat tiller formation and development.

PPVED: A machine learning tool for predicting the effect of single amino acid substitution on protein function in plants.

Single amino acid substitution (SAAS) produces the most common variant of protein function change under physiological conditions. As the number of SAAS events in plants has increased exponentially, an effective prediction tool is required to help identify and distinguish functional SAASs from the whole genome as either potentially causal traits or as variants. Here, we constructed a plant SAAS database that stores 12 865 SAASs in 6172 proteins and developed a tool called Plant Protein Variation Effect Detector (PPVED) that predicts the effect of SAASs on protein function in plants. PPVED achieved an 87% predictive accuracy when applied to plant SAASs, an accuracy that was much higher than those from six human database software: SIFT, PROVEAN, PANTHER-PSEP, PhD-SNP, PolyPhen-2, and MutPred2. The predictive effect of six SAASs from three proteins in Arabidopsis and maize was validated with wet lab experiments, of which five substitution sites were accurately predicted. PPVED could facilitate the identification and characterization of genetic variants that explain observed phenotype variations in plants, contributing to solutions for challenges in functional genomics and systems biology. PPVED can be accessed under a CC-BY (4.0) license via http://www.ppved.org.cn.

Teosinte confers specific alleles and yield potential to maize improvement.

KEY MESSAGE: Teosinte improves maize grain yield and broadens the maize germplasm. Seventy-one quantitative trait loci associated with 24 differential traits between maize and teosinte were identified. Maize is a major cereal crop with a narrow germplasm that has limited its production and breeding progress. Teosinte, an ancestor of maize, provides valuable genetic resources for maize breeding. To identify the favorable alien alleles in teosinte and its yield potential for maize breeding, 4 backcrossed maize-teosinte recombinant inbred line (RIL) populations were cultivated under five conditions. A North Carolina mating design II experiment was conducted on inbred lines with B73 and Mo17 pedigree backgrounds to analyze their combining ability. Abundant phenotypic variation on 26 traits of four RIL populations were found, of which barren tip length, kernel height, and test weight showed positive genetic improvement potential. The hybrid FM132 (BD138/MP116) showed a superior grain yield to that of the check, with an average yield gain of 4.86%. Moreover, inbred lines BD138 and MP048 showed a higher general grain yield combining ability than those of their corresponding checks. We screened 4,964,439 high-quality single-nucleotide polymorphisms in the BD (B73/Zea diploperennis) RIL population for bin construction and used 2322 bin markers for genetic map construction and quantitative trait loci (QTL) mapping. Via inclusive composite interval mapping, 71 QTL associated with 24 differential traits were identified. Gene annotation and transcriptional expression suggested that Zm00001eb352570 and Zm00001eb352580, both annotated as ethylene-responsive transcription factors, were key candidate genes that regulate ear height and the ratio of ear to plant height. Our results indicate that teosinte could broaden the narrow maize germplasm, improve yield potential, and provide desirable alleles for maize breeding.

Identification of lncRNAs involved in wheat tillering development in two pairs of near-isogenic lines.

Emerging evidence demonstrates that lncRNAs participate in various developmental processes in plants via post-transcription regulation. However, few lncRNAs have been identified as regulators of tiller development in wheat (Triticum aestivum L.). In this study, high-throughput ribosomal depleted RNA sequencing was performed on the tillering nodes of two pairs of near-isogenic lines that differed only in the tillering trait. We identified 5399 lncRNA transcripts using bioinformational analyses. KEGG pathway analysis revealed 74 common differentially expressed lncRNAs substantially enriched in photosynthesis-related, phenylpropanoid biosynthesis, phosphatidylinositol signaling, brassinosteroid biosynthesis, zeatin biosynthesis, and carotenoid biosynthesis pathways. Detailed functional annotations of target genes were used to identify 27 tillering-associated lncRNAs. Among these, 10 were in photosynthesis-related pathways; 15 were in secondary metabolite pathways; and 8 were in plant hormone pathways, with 6 enriched in two kinds of pathways. These findings contribute to identifying tillering-associated lncRNAs in wheat and enable further investigation into the functions and roles of key candidate lncRNAs, and more experimental evidence was also needed if breeders wanted to utilize these candidate lncRNAs in wheat crop yield improvement in the future.

Evolutionary patterns of DNA base composition at polymorphic sites highlight the role of the environment in shaping barley and rice genomes.

Insights into changes in genome base composition underlying crop domestication can be gained by using comparative genomics. With this approach, previous studies have reported that crop genomes during domestication accumulate more nucleotides adenine (A) and thymine (T) (termed as [AT]-increase) across polymorphic sites. However, the potential influence of the environment or its factors, for example, solar ultraviolet (UV) radiation and temperature, on the [AT]-increase has not been well elucidated. Here, we investigated the [AT]-increase in barley (Hordeum vulgare L.) and rice (Oryza sativa L.) and the association with natural environments, where accessions are distributed. With 12,798,376 and 2,861,535 single-nucleotide polymorphisms from 368 barley and 1375 rice accessions, respectively, we discovered that [AT] increases from wild accessions to improved cultivars, and genomic regions with larger [AT]-increase tend to have higher UV-related motif frequencies, suggesting solar UV radiation as a potential factor in driving genome variation. To link [AT] change with the geographic distribution, we gathered georeferenced accessions and examined their local environments. Interestingly, negative correlations between [AT] and environmental factors were observed (r = -0.39 ∼ -0.75) and modern accessions with higher [AT] values, as compared with wild relatives, are from the environments with lower solar UV radiation or lower temperature. With [AT] and environmental factors as phenotypes, genome-wide association mapping identified three candidate genes that have the potential to contribute to [AT] variation under the effect of environmental conditions. Our findings provide genomic and environmental insights into evolutionary pattern of DNA base composition and underlying mechanisms.

SSRMMD: A Rapid and Accurate Algorithm for Mining SSR Feature Loci and Candidate Polymorphic SSRs Based on Assembled Sequences.

Microsatellites or simple sequence repeats (SSRs) are short tandem repeats of DNA widespread in genomes and transcriptomes of diverse organisms and are used in various genetic studies. Few software programs that mine SSRs can be further used to mine polymorphic SSRs, and these programs have poor portability, have slow computational speed, are highly dependent on other programs, and have low marker development rates. In this study, we develop an algorithm named Simple Sequence Repeat Molecular Marker Developer (SSRMMD), which uses improved regular expressions to rapidly and exhaustively mine perfect SSR loci from any size of assembled sequence. To mine polymorphic SSRs, SSRMMD uses a novel three-stage method to assess the conservativeness of SSR flanking sequences and then uses the sliding window method to fragment each assembled sequence to assess its uniqueness. Furthermore, molecular biology assays support the polymorphic SSRs identified by SSRMMD. SSRMMD is implemented using the Perl programming language and can be downloaded from https://github.com/GouXiangJian/SSRMMD.