Proteome evaluation of homolog abundance patterns in Arachis hypogaea cv. Tifrunner

Background Cultivated peanut (Arachis hypogaea, AABB genome), an allotetraploid from a cross between A. duranensis (AA genome) and A. ipaensis (BB genome), is an important oil and protein crop with released genome and RNA-seq sequence datasets. These datasets provide the molecular foundation for studying gene expression and evolutionary patterns. However, there are no reports on the proteomic data of A. hypogaea cv. Tifrunner, which limits understanding of its gene function and protein level evolution. Results This study sequenced the A. hypogaea cv. Tifrunner leaf and root proteome using the tandem mass tag technology. A total of 4803 abundant proteins were identified. The 364 differentially abundant proteins were estimated by comparing protein abundances between leaf and root proteomes. The differentially abundant proteins enriched the photosystem process. The number of biased abundant homeologs between the two sub-genomes A (87 homeologs in leaf and root) and B (69 and 68 homeologs in leaf and root, respectively) was not significantly different. However, homeologous proteins with biased abundances in different sub-genomes enriched different biological processes. In the leaf, homeologs biased to sub-genome A enriched biosynthetic and metabolic process, while homeologs biased to sub-genome B enriched iron ion homeostasis process. In the root, homeologs with biased abundance in sub-genome A enriched inorganic biosynthesis and metabolism process, while homeologs with biased abundance in sub-genome B enriched organic biosynthesis and metabolism process. Purifying selection mainly acted on paralogs and homeologs. The selective pressure values were negatively correlated with paralogous protein abundance. About 77.42% (24/31) homeologous and 80% (48/60) paralogous protein pairs had asymmetric abundance, and several protein pairs had conserved abundances in the leaf and root tissues. Conclusions This study sequenced the proteome of A. hypogaea cv. Tifrunner using the leaf and root tissues. Differentially abundant proteins were identified, and revealed functions. Paralog abundance divergence and homeolog bias abundance was elucidated. These results indicate that divergent abundance caused retention of homologs in A. hypogaea cv. Tifrunner.

Development and validation of diagnostic SNP markers for quality control genotyping in a collection of four rice (Oryza) species

Morphological identification of closely related rice species, particularly those in the Oryza AA genome group, presents major challenges and often results in cases of misidentification. Recent work by this group identified diagnostic single nucleotide polymorphic (SNP) markers specific for several rice species and subspecies based on DArTseq next-generation sequencing technology (“DArTseq”). These SNPs can be used for quality control (QC) analysis in rice breeding and germplasm maintenance programs. Here, we present the DArTseq-based diagnostic SNPs converted into Kompetitive allele-specific PCR (KASPar or KASP) assays and validation data for a subset of them; these can be used for low-cost routine genotyping quality control (QC) analysis. Of the 224 species/subspecies’ diagnostic SNPs tested, 158 of them produced working KASP assays, a conversion success rate of 70%. Two validation experiments were run with 87 of the 158 SNP markers to ensure that the assays amplified, were polymorphic, and distinguished the five species/subspecies tested. Based on these validation test results, we recommend a panel of 36 SNP markers that clearly delineate O. barthii, O. glaberrima, O. longistaminata, O. sativa spp. indica and japonica. The KASP assays provide a flexible, rapid turnaround and cost-effective tool to facilitate germplasm curation and management of these four Oryza AA genome species across multiple genebanks.

Analysis of Genomic Variation in Different Brassica Napus Synthetic Allopolyploids

Allopolyploidy is an evolutionary and mechanisticaly intriguing process involving the reconciliation of two or more sets of diverged genomes and regulatory interactions, resulting in new phenotypes. In this study, we explored the genomic variation of eight F2 synthetic B. napus using whole-genome sequencing. We found that there was a genetic variation in the F2 generation. Part of the variation was consistent in the F2 generation, and a small number of mutations only appeared in a single plant of the F2 generation. The analysis of copy number variation (CNV) found that most of the AA genome was lost, and most of the CC genome was obtained. In addition, there was inter-chromosomal translocation (CTX) in the F2 generation and the number of each plant was different. The above results indicate that the F2 generation showed genetic variation and there was a difference between eight plants, which may lay a molecular basis for the unique field performance of the offspring. It provides a new perspective of genomic variation and trait separation in the early stages of allopolyploid polyploid formation.

Development of a genome-wide InDel marker set for allele discrimination between rice (Oryza sativa) and the other seven AA-genome Oryza species

Wild relatives of rice in the genus Oryza (composed of 24 species with 11 different genome types) have been significantly contributing to the varietal improvement of rice (Oryza sativa). More than 4000 accessions of wild rice species are available and they are regarded as a “genetic reservoir” for further rice improvement. DNA markers are essential tools in genetic analysis and breeding. To date, genome-wide marker sets for wild rice species have not been well established and this is one of the major difficulties for the efficient use of wild germplasm. Here, we developed 541 genome-wide InDel markers for the discrimination of alleles between the cultivated species O. sativa and the other seven AA-genome species by positional multiple sequence alignments among five AA-genome species with four rice varieties. The newly developed markers were tested by PCR-agarose gel analysis of 24 accessions from eight AA genome species (three accessions per species) along with two representative cultivars (O. sativa subsp. indica cv. IR24 and subsp. japonica cv. Nipponbare). Marker polymorphism was validated for 475 markers. The number of polymorphic markers between IR24 and each species (three accessions) ranged from 338 (versus O. rufipogon) to 416 (versus O. longistaminata) and the values in comparison with Nipponbare ranged from 179 (versus O. glaberrima) to 323 (versus O. glumaepatula). These marker sets will be useful for genetic studies and use of the AA-genome wild rice species.

Exploring the Loci Responsible for Awn Development in Rice through Comparative Analysis of All AA Genome Species

Wild rice species have long awns at their seed tips, but this trait has been lost through rice domestication. Awn loss mitigates harvest and seed storage; further, awnlessness increases the grain number and, subsequently, improves grain yield in Asian cultivated rice, highlighting the contribution of the loss of awn to modern rice agriculture. Therefore, identifying the genes regulating awn development would facilitate the elucidation of a part of the domestication process in rice and increase our understanding of the complex mechanism in awn morphogenesis. To identify the novel loci regulating awn development and understand the conservation of genes in other wild rice relatives belonging to the AA genome group, we analyzed the chromosome segment substitution lines (CSSL). In this study, we compared a number of CSSL sets derived by crossing wild rice species in the AA genome group with the cultivated species Oryza sativa ssp. japonica. Two loci on chromosomes 7 and 11 were newly discovered to be responsible for awn development. We also found wild relatives that were used as donor parents of the CSSLs carrying the functional alleles responsible for awn elongation, REGULATOR OF AWN ELONGATION 1 (RAE1) and RAE2. To understand the conserveness of RAE1 and RAE2 in wild rice relatives, we analyzed RAE1 and RAE2 sequences of 175 accessions among diverse AA genome species retrieved from the sequence read archive (SRA) database. Comparative sequence analysis demonstrated that most wild rice AA genome species maintained functional RAE1 and RAE2, whereas most Asian rice cultivars have lost either or both functions. In addition, some different loss-of-function alleles of RAE1 and RAE2 were found in Asian cultivated species. These findings suggest that different combinations of dysfunctional alleles of RAE1 and RAE2 were selected after the speciation of O. sativa, and that two-step loss of function in RAE1 and RAE2 contributed to awnlessness in Asian cultivated rice.

Genome-wide analysis of long noncoding RNAs, 24-nt siRNAs, DNA methylation and H3K27me3 marks in Brassica rapa

Long noncoding RNAs (lncRNAs) are RNA fragments that generally do not code for a protein but are involved in epigenetic gene regulation. In this study, lncRNAs of Brassica rapa were classified into long intergenic noncoding RNAs, natural antisense RNAs, and intronic noncoding RNAs and their expression analyzed in relation to genome-wide 24-nt small interfering RNAs (siRNAs), DNA methylation, and histone H3 lysine 27 trimethylation marks (H3K27me3). More than 65% of the lncRNAs analyzed consisted of one exon, and more than 55% overlapped with inverted repeat regions (IRRs). Overlap of lncRNAs with IRRs or genomic regions encoding for 24-nt siRNAs resulted in increased DNA methylation levels when both were present. LncRNA did not overlap greatly with H3K27me3 marks, but the expression level of intronic noncoding RNAs that did coincide with H3K27me3 marks was higher than without H3K27me3 marks. The Brassica genus comprises important vegetables and oil seed crops grown across the world. B. rapa is a diploid (AA genome) thought to be one of the ancestral species of both B. juncea (AABB genome) and B. napus (AACC) through genome merging (allotetrapolyploidization). Complex genome restructuring and epigenetic alterations are thought to be involved in these allotetrapolyploidization events. Comparison of lncRNAs between B. rapa and B. nigra, B. oleracea, B. juncea, and B. napus showed the highest conservation with B. oleracea. This study presents a comprehensive analysis of the epigenome structure of B. rapa at multi-epigenetic levels (siRNAs, DNA methylation, H3K27me3, and lncRNAs) and identified a suite of candidate lncRNAs that may be epigenetically regulated in the Brassica genus.

Identification of MaWRKY40 and MaDLO1 as effective marker genes for tracking the salicylic acid-mediated immune response in bananas

Bananas lie among the world’s most important cash and staple crops but are threatened by various devastating pathogens. The phytohormone salicylic acid (SA) plays a key role in the regulation of plant immune response. Tracking the expression of SA-responsive marker genes under pathogen infection is important in pathogenesis elucidation. However, the common SA-responsive marker genes are not consistently induced in different banana cultivars or different organs. Here, we conducted transcriptome analysis for SA response of a banana cultivar, ‘Pei-Chiao’ (Cavendish, AAA genome), and identified three genes, MaWRKY40, MaWRKY70, and Downy Mildew Resistant 6 (DMR6)-Like Oxygenase 1 (MaDLO1) that are robustly induced upon SA treatment in both the leaves and roots. Consistent induction of these three genes by SA treatment was also detected in both the leaves and roots of bananas belonging to different genome types such as ‘Tai-Chiao No. 7’ (Cavendish, AAA genome), ‘Pisang Awak’ (ABB genome), and ‘Lady Finger’ (AA genome). Furthermore, the biotrophic pathogen cucumber mosaic virus elicited the expression of MaWRKY40 and MaDLO1 in infected-leaves of susceptible cultivars. The hemi-biotrophic fungal pathogen Fusarium oxysporum f. sp. cubense tropical race 4 (Foc TR4) also consistently induced the expression of MaWRKY40 and MaDLO1 in the infected-roots of the Foc TR4-resistant cultivar. These results indicate that MaWRKY40 and MaDLO1 can serve as reliable SA-responsive marker genes for the study of plant immunity in banana. Revealing SA-responsive marker genes provides a stepping-stone for further studies in banana resistance to pathogens.

A Genetic Resource for Rice Improvement：introgression Library of Agronomic Traits for All AA Genome Species in Genus Oryza

Background: Rice improvement depends on the availability of genetic variation, and AA genome Oryza species are the natural reservoir of favorable genes which are useful for rice breeding. Developing the introgression library using multiple AA genome species was rarely reported.Results: In this study, to systematically evaluate and utilize potentially valuable QTLs/genes or allelic variations, based on the evaluation and selection of agronomic traits, 6372 introgression lines (ILs) were raised by crossing 330 accessions of 7 AA genome species as the donor parents, with three elite cultivars of O. sativa, Dianjingyou 1, Yundao 1 and RD23 as the recurrent parents, respectively. Further, twenty-six, twenty-six and nineteen loci were detected in the multiple donors using 1,401 ILs in the Dianjingyou 1 background for grain length, grain width, and the ratio of grain length to grain width, respectively. Interestingly, ten loci had opposite effect on grain length in the different donors, so did grain width. Moreover, one locus for grain width, qGW3.1, was validated using the segregation population derived from the donor of O. glumaepatula. Conclusions: This introgression library provided the powerful resource for future rice improvement and genetic dissection of allelic variations. Selections of favorable alleles that are present in wild relatives proceed the driving force of the rice domestication.

Analysis of Transcriptional Changes in Different Brassica napus Synthetic Allopolyploids

Allopolyploidy is an evolutionary and mechanistically intriguing process involving the reconciliation of two or more sets of diverged genomes and regulatory interactions, resulting in new phenotypes. In this study, we explored the gene expression patterns of eight F2 synthetic Brassica napus using RNA sequencing. We found that B. napus allopolyploid formation was accompanied by extensive changes in gene expression. A comparison between F2 and the parent shows a certain proportion of differentially expressed genes (DEG) and activation\silent gene, and the two genomes (female parent (AA)\male parent (CC) genomes) showed significant differences in response to whole-genome duplication (WGD); non-additively expressed genes represented a small portion, while Gene Ontology (GO) enrichment analysis showed that it played an important role in responding to WGD. Besides, genome-wide expression level dominance (ELD) was biased toward the AA genome, and the parental expression pattern of most genes showed a high degree of conservation. Moreover, gene expression showed differences among eight individuals and was consistent with the results of a cluster analysis of traits. Furthermore, the differential expression of waxy synthetic pathways and flowering pathway genes could explain the performance of traits. Collectively, gene expression of the newly formed allopolyploid changed dramatically, and this was different among the selfing offspring, which could be a prominent cause of the trait separation. Our data provide novel insights into the relationship between the expression of differentially expressed genes and trait segregation and provide clues into the evolution of allopolyploids.