Identification and Expression Analysis of BURP Domain-Containing Genes in Medicago truncatula

Genome-wide identification and expression analysis of the BURP domain-containing genes in Malus domestica

The conserved BURP-containing proteins are specific to plants and play a crucial role in plant growth, development, and response to abiotic stresses. However, less is known about the systematic characterization of BURP-containing proteins in apple. This study aimed to identify and analyze all BURP-containing genes in the apple genome, as well as to examine their expression patterns through various bioinformatics methods. Eighteen members of BURP-containing genes were identified in apple, six members lacked signal peptides, and the secondary structure was mainly a Random coil of BURP-containing genes. Gene structure and Motif analysis showed that proteins have similar structures and are conserved at the C-terminal. Cis-acting element analysis revealed that the proteins contain phytohormone and stress response elements, and chromosomal localization revealed that the family is unevenly distributed across eight chromosomes, with duplication of fragments leading to the expansion of family proteins. Tissue expression showed that MdPG3 and MdPG4 were expressed in different tissues and different varieties, MdRD2 and MdRD7 were highly expressed in 'M74' fruits and MdRD7 in 'M49' leaves, while MdUSP1 was highly expressed in 'GD' roots. The quantitative real-time PCR analysis showed that the expressions of six and seven genes were significantly up-regulated under NaCl and PEG treatments, respectively, whereas MdRD7 was significantly up-regulated under NaCl and PEG treatment over time. This study offers a comprehensive identification and expression analysis of BURP-containing proteins in apple. The findings provide a theoretical foundation for further exploration of the functions of this protein family.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01393-7.

Identification and expression analysis of BURP domain-containing genes in jujube and their involvement in low temperature and drought response

Background

Plant-specific BURP domain-containing genes are involved in plant development and stress responses. However, the role of BURP family in jujube (Ziziphus jujuba Mill.) has not been investigated.

Results

In this study, 17 BURP genes belonging to four subfamilies were identified in jujube based on homology analysis, gene structures, and conserved motif confirmation. Gene duplication analysis indicated both tandem duplication and segmental duplication had contributed to ZjBURP expansion. The ZjBURPs were extensively expressed in flowers, young fruits, and jujube leaves. Transcriptomic data and qRT-PCR analysis further revealed that ZjBURPs also significantly influence fruit development, and most genes could be induced by low temperature, salinity, and drought stresses. Notably, several BURP genes significantly altered expression in response to low temperature (ZjPG1) and drought stresses (ZjBNM7, ZjBNM8, and ZjBNM9).

Conclusions

These results provided insights into the possible roles of ZjBURPs in jujube development and stress response. These findings would help selecting candidate ZjBURP genes for cold- and drought-tolerant jujube breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08907-9.

Deciphering genome-wide transcriptomic changes in grapevines heavily infested by spotted lanternflies

The spotted lanternfly, a newly invasive insect in the U.S. that is a great concern for the grapevine industry, produces damage on its host plants through aggressive feeding, using a piercing and sucking method to feed on the phloem of plants. In the eastern US, adult SLF can invade vineyards through fruit ripening until the end of the growing season; however, it is still unclear how prolonged late-season SLF feeding can affect the health of grapevines, as well as the host responses to this extensive damage. Thus, we have performed a comprehensive genome-wide transcriptome analysis in grapevines heavily infested by the spotted lanternfly, as it occurs in Pennsylvania vineyards, and compared it to other relevant transcriptomes in grapes with different degrees to susceptibility to similar pests. Among a variety of plant responses, we highlight here a subset of relevant biological pathways that distinguish or are common to the spotted lanternfly and other phloem feeders in grapevine. The molecular interaction between spotted lanternfly and the vine begins with activation of signal transduction cascades mediated mainly by protein kinase genes. It also induces the expression of transcription factors in the nucleus, of other signaling molecules like phytohormones and secondary metabolites, and their downstream target genes responsible for defense and physiological functions, such as detoxification and photosynthesis. Grapevine responses furthermore include the activation of genes for cell wall strengthening via biosynthesis of major structural components. With this study, we hope to provide the regulatory network to explain effects that the invasive spotted lanternfly has on grapevine health with the goal to improve its susceptibility.

Overexpressing PpBURP2 in Rice Increases Plant Defense to Abiotic Stress and Bacterial Leaf Blight

Mosses are one of the earliest diverging land plants that adapted to living on land. The BURP domain-containing proteins (BURP proteins) are plant-specific proteins that appeared when plants shifted from aquatic environments to land. Phylogenetic analysis revealed that the BURP domain of higher plants is originated from lower land plants and divergent because of motif conversion. To discover the function of BURP protein in moss, rice transgenics with ectopic expression of PpBURP2 were subjected to different abiotic stresses treatments. The results revealed that the ectopic expression of PpBURP2 enhanced the tolerance to osmotic and saline stresses at the seedling stage and drought stress at the adult stage. Further ectopic expression of PpBURP2 improved the cadmium (2+) (Cd2+) tolerance and reduced Cd2+ accumulation in rice leaves. Transcriptomic analysis of the transgenic PpBURP2 plants showed that the differentially expressed genes were involved in the metabolism of secondary metabolites, energy, oxidation-reduction process, and defense-related genes. Further experiments showed that the photosynthetic efficiency and resistance against bacterial leaf blight were obviously improved in transgenic plants. Yeast two-hybrid and bimolecular fluorescence complementation (BiFC) assays revealed the physical interaction of BURP domain protein from rice and moss with mitogen-activated protein kinase kinase (MKK) from rice. Therefore, our findings demonstrate that overexpressing PpBURP2 in rice confers resistance to abiotic stresses and bacterial leaf blight. They also suggested that the regulatory role of BURP-like proteins across lower and higher plants was evolutionary conservation of responses of different classes of plants to different environmental challenges.

De novo genome assembly and in natura epigenomics reveal salinity-induced DNA methylation in the mangrove tree Bruguiera gymnorhiza

Mangroves are adapted to harsh environments, such as high ultraviolet (UV) light, low nutrition, and fluctuating salinity in coastal zones. However, little is known about the transcriptomic and epigenomic basis of the resilience of mangroves due to limited available genome resources. We performed a de novo genome assembly and in natura epigenome analyses of the mangrove Bruguiera gymnorhiza, one of the dominant mangrove species. We also performed the first genome-guided transcriptome assembly for mangrove species. The 309 Mb of the genome is predicted to encode 34 403 genes and has a repeat content of 48%. Depending on its growing environment, the natural B. gymnorhiza population showed drastic morphological changes associated with expression changes in thousands of genes. Moreover, high-salinity environments induced genome-wide DNA hypermethylation of transposable elements (TEs) in the B. gymnorhiza. DNA hypermethylation was concurrent with the transcriptional regulation of chromatin modifier genes, suggesting robust epigenome regulation of TEs in the B. gymnorhiza genome under high-salinity environments. The genome and epigenome data in this study provide novel insights into the epigenome regulation of mangroves and a better understanding of the adaptation of plants to fluctuating, harsh natural environments.

Genome-wide analysis of BURP genes and identification of a BURP-V gene RcBURP4 in Rosa chinensis

Nine RcBURPs have been identified in Rosa chinensis, and overexpression of RcBURP4 increased ABA, NaCl sensitivity, and drought tolerance in transgenic Arabidopsis. BURP proteins are unique to plants and may contribute greatly to growth, development, and stress responses of plants. Despite the vital role of BURP proteins, little is known about these proteins in rose (Rosa spp.). In the present study, nine genes belonging to the BURP family in R. chinensis were identified using multiple bioinformatic approaches against the rose genome database. The nine RcBURPs, with diverse structures, were located on all chromosomes of the rose genome, except for Chr2 and Chr3. Phylogenic analysis revealed that these RcBURPs can be classified into eight subfamilies, including BNM2-like, PG1β-like, USP-like, RD22-like, BURP-V, BURP-VI, BURP-VII, and BURP-VIII. Conserved motif and exon-intron analyses indicated a conserved pattern within the same subfamily. The presumed cis-regulatory elements (CREs) within the promoter region of each RcBURP were analyzed and the results showed that all RcBURPs contained different types of CREs, including abiotic stress-, light response-, phytohormones response-, and plant growth and development-related CREs. Transcriptomic analysis revealed that a BURP-V member, RcBURP4, was induced in rose leaves and roots under mild and severe drought treatments. We then overexpressed RcBURP4 in Arabidopsis and examined its role under abscisic acid (ABA), NaCl, polyethylene glycol (PEG), and drought treatments. Nine stress-responsive genes expression were changed in RcBURP4-overexpressing leaves and roots. Furthermore, RcBURP4-silenced rose plants exhibited decreased tolerance to dehydration. The results obtained from this study provide the first comprehensive overview of RcBURPs and highlight the importance of RcBURP4 in rose plant.

Genome-wide identification of the BURP domain-containing genes in Phaseolus vulgaris

Plant-specific BURP domain-containing proteins have an essential role in the plant's development and stress responses. Although BURP domain-containing proteins have been identified in several plant species, genome-wide analysis of the BURP gene family has not been investigated in the common bean. In the present study, we identified 11 BURP family members in the common bean (Phaseolus vulgaris) genome with a comprehensive in silico analysis. Pairwise alignment and phylogenetic analyses grouped PvBURP members into four subfamilies [RD-22 like (3), PG1β-like (4), BNM2-like (3), and USP-like (1)] according to their amino acid motifs, protein domains and intron-exon structure. The physical and biochemical characteristics of amino acids, motif and intron-exon structure, and cis-regulatory elements of BURPs members were determined. Promoter regions of BURP members included stress, light, and hormone response-related cis-elements. Therefore, expression profiles of PvBURP genes were identified with in silico tools and qRT-PCR analyses under stress (salt and drought) and hormone treatment (ABA, IAA) in the current study. While significant activity changes were not observed in BURP genes in RNA-seq data sets related to salt stress, it was determined that some BURP genes were expressed differently in those with drought stress. We identified 12 different miRNA, including miRNA395, miRNA156, miRNA169, miRNA171, miRNA319, and miRNA390, targeting the nine PvBURP genes using two different in silico tools based on perfect or near-perfect complementarity to their targets. Here we present the first study to identify and characterize the BURP genes in common bean using whole-genome analysis, and the findings may serve as a reference for future functional research in common bean.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01052-9.

The ankyrin repeat gene family in Capsicum spp: Genome-wide survey, characterization and gene expression profile

The ankyrin (ANK) repeat protein family is largely distributed across plants and has been found to participate in multiple processes such as plant growth and development, hormone response, response to biotic and abiotic stresses. It is considered as one of the major markers of capsaicin content in pepper fruits. In this study, we performed a genome-wide identification and expression analysis of genes encoding ANK proteins in three Capsicum species: Capsicum baccatum, Capsicum annuum and Capsicum chinense. We identified a total of 87, 85 and 96 ANK genes in C. baccatum, C. annuum and C. chinense genomes, respectively. Next, we performed a comprehensive bioinformatics analysis of the Capsicum ANK gene family including gene chromosomal localization, Cis-elements, conserved motif identification, intron/exon structural patterns and gene ontology classification as well as profile expression. Phylogenetic and domain organization analysis grouped the Capsicum ANK gene family into ten subfamilies distributed across all 12 pepper chromosomes at different densities. Analysis of the expression of ANK genes in leaf and pepper fruits suggested that the ANKs have specific expression patterns at various developmental stages in placenta tissue. Our results provide valuable information for further studies of the evolution, classification and putative functions of ANK genes in pepper.

Comparative analysis of the accelerated aged seed transcriptome profiles of two maize chromosome segment substitution lines

Seed longevity is one of the most essential characteristics of seed quality. Two chromosome segment substitution lines, I178 and X178, which show significant differences in seed longevity, were subjected to transcriptome sequencing before and after five days of accelerated aging (AA) treatments. Compared to the non-aging treatment, 286 and 220 differentially expressed genes (DEGs) were identified after 5 days of aging treatment in I178 and X178, respectively. Of these DEGs, 98 were detected in both I178 and X178, which were enriched in Gene Ontology (GO) terms of the cellular component of the nuclear part, intracellular part, organelle and membrane. Only 86 commonly downregulated genes were enriched in GO terms of the carbohydrate derivative catabolic process. Additionally, transcriptome analysis of alternative splicing (AS) events in I178 and X178 showed that 63.6% of transcript isoforms occurred AS in all samples, and only 1.6% of transcript isoforms contained 169 genes that exhibited aging-specific AS arising after aging treatment. Combined with the reported QTL mapping result, 7 DEGs exhibited AS after aging treatment, and 13 DEGs in mapping interval were potential candidates that were directly or indirectly related to seed longevity.

RNAseq analysis reveals drought-responsive molecular pathways with candidate genes and putative molecular markers in root tissue of wheat

Drought is one of the major impediments in wheat productivity. Traditional breeding and marker assisted QTL introgression had limited success. Available wheat genomic and RNA-seq data can decipher novel drought tolerance mechanisms with putative candidate gene and marker discovery. Drought is first sensed by root tissue but limited information is available about how roots respond to drought stress. In this view, two contrasting genotypes, namely, NI5439 41 (drought tolerant) and WL711 (drought susceptible) were used to generate ~78.2 GB data for the responses of wheat roots to drought. A total of 45139 DEGs, 13820 TF, 288 miRNAs, 640 pathways and 435829 putative markers were obtained. Study reveals use of such data in QTL to QTN refinement by analysis on two model drought-responsive QTLs on chromosome 3B in wheat roots possessing 18 differentially regulated genes with 190 sequence variants (173 SNPs and 17 InDels). Gene regulatory networks showed 69 hub-genes integrating ABA dependent and independent pathways controlling sensing of drought, root growth, uptake regulation, purine metabolism, thiamine metabolism and antibiotics pathways, stomatal closure and senescence. Eleven SSR markers were validated in a panel of 18 diverse wheat varieties. For effective future use of findings, web genomic resources were developed. We report RNA-Seq approach on wheat roots describing the drought response mechanisms under field drought conditions along with genomic resources, warranted in endeavour of wheat productivity.

Unraveling the complex genome of Saccharum spontaneum using Polyploid Gene Assembler

The Polyploid Gene Assembler (PGA), developed and tested in this study, represents a new strategy to perform gene-space assembly from complex genomes using low coverage DNA sequencing. The pipeline integrates reference-assisted loci and de novo assembly strategies to construct high-quality sequences focused on gene content. Pipeline validation was conducted with wheat (Triticum aestivum), a hexaploid species, using barley (Hordeum vulgare) as reference, that resulted in the identification of more than 90% of genes and several new genes. Moreover, PGA was used to assemble gene content in Saccharum spontaneum species, a parental lineage for hybrid sugarcane cultivars. Saccharum spontaneum gene sequence obtained was used to reference-guided transcriptome analysis of six different tissues. A total of 39,234 genes were identified, 60.4% clustered into known grass gene families. Thirty-seven gene families were expanded when compared with other grasses, three of them highlighted by the number of gene copies potentially involved in initial development and stress response. In addition, 3,108 promoters (many showing tissue specificity) were identified in this work. In summary, PGA can reconstruct high-quality gene sequences from polyploid genomes, as shown for wheat and S. spontaneum species, and it is more efficient than conventional genome assemblers using low coverage DNA sequencing.

Differential gene expression and gene ontologies associated with increasing water-stress in leaf and root transcriptomes of perennial ryegrass (Lolium perenne)

Perennial ryegrass (Lolium perenne) is a forage and amenity grass species widely cultivated in temperate regions worldwide. As such, perennial ryegrass populations are exposed to a range of environmental conditions and stresses on a seasonal basis and from year to year. One source of potential stress is limitation on water availability. The ability of these perennial grasses to be able to withstand and recover after periods of water limitation or drought can be a key component of grassland performance. Thus, we were interested in looking at changes in patterns of gene expression associated with increasing water stress. Clones of a single genotype of perennial ryegrass were grown under non-flowering growth room conditions in vermiculite supplemented with nutrient solution. Leaf and root tissue was sampled at 4 times in quadruplicate relating to estimated water contents of 35%, 15%, 5% and 1%. RNA was extracted and RNAseq used to generate transcriptome profiles at each sampling point. Transcriptomes were assembled using the published reference genome sequence and differential gene expression analysed using 3 different programmes, DESeq2, edgeR and limma (with the voom transformation), individually and in combination, deriving Early, Middle and Late stage comparisons. Identified differentially expressed genes were then associated with enriched GO terms using BLAST2GO. For the leaf, up-regulated differentially expressed genes were strongly associated with GO terms only during the Early stage and the majority of GO terms were associated with only down-regulated genes at the Middle or Late stages. For the roots, few differentially expressed genes were identified at either Early or Middle stages. Only one replicate at 1% estimated water content produced high quality data for the root, however, this indicated a high level of differential expression. Again the majority of enriched GO terms were associated with down-regulated genes. The performance of the different analysis programmes and the annotations associated with identified differentially expressed genes is discussed.

Genome-wide identification and expression analysis of the BURP domain-containing genes in Gossypium hirsutum

Background

Many BURP domain-containing proteins, which are unique to plants, have been identified. They performed diverse functions in plant development and the stress response. To date, only a few BURP domain-containing genes have been studied, and no comprehensive analysis of the gene family in cotton has been reported.

Results

In this study, 18, 17 and 30 putative BURP genes were identified in G. raimondii (D₅), G. arboreum (A₂) and G. hirsutum (AD₁), respectively. These BURP genes were phylogenetically classified into eight subfamilies, which were confirmed by analyses of gene structures, motifs and protein domains. The uneven distribution of BURPs in chromosomes and gene duplication analysis indicated that segmental duplication might be the main driving force of the GhBURP family expansion. Promoter regions of all GhBURPs contained at least one putative stress-related cis-elements. Analysis of transcriptomic data and qRT-PCR showed that GhBURPs showed different expression patterns in different organs, and all of them, especially the members of the RD22-like subfamily, could be induced by different stresses, such as abscisic acid (ABA) and salicylic acid (SA), which indicated that the GhBURPs may performed important functions in cotton’s responses to various abiotic stresses.

Conclusions

Our study comprehensively analyzed BURP genes in G. hirsutum, providing insight into the functions of GhBURPs in cotton development and adaptation to stresses.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5948-y) contains supplementary material, which is available to authorized users.

Genome-wide identification and gene expression pattern of ABC transporter gene family in Capsicum spp

ATP-binding cassette (ABC) transporter genes act as transporters for different molecules across biological membranes and are involved in a diverse range of biological processes. In this study, we performed a genome-wide identification and expression analysis of genes encoding ABC transporter proteins in three Capsicum species, i.e., Capsicum annuum, Capsicum baccatum and Capsicum chinense. Capsicum is a valuable horticultural crop worldwide as an important constituent of many foods while containing several medicinal compounds including capsaicin and dihydrocapsaicin. Our results identified the presence of a total of 200, 185 and 187 ABC transporter genes in C. annuum, C. baccatum and C. chinense genomes, respectively. Capsaicin and dihydrocapsaicin content were determined in green pepper fruits (16 dpa). Additionally, we conducted different bioinformatics analyses including ABC genes classification, gene chromosomal location, Cis elements, conserved motifs identification and gene ontology classification, as well as profile expression of selected genes. Based on phylogenetic analysis and domain organization, the Capsicum ABC gene family was grouped into eight subfamilies. Among them, members within the ABCG, ABCB and ABCC subfamilies were the most abundant, while ABCD and ABCE subfamilies were less abundant throughout all species. ABC members within the same subfamily showed similar motif composition. Furthermore, common cis-elements involved in the transcriptional regulation were also identified in the promoter regions of all Capsicum ABC genes. Gene expression data from RNAseq and reverse transcription-semi-quantitative PCR analysis revealed development-specific stage expression profiles in placenta tissues. It suggests that ABC transporters, specifically the ABCC and ABCG subfamilies, may be playing important roles in the transport of secondary metabolites such as capsaicin and dihydrocapsaicin to the placenta vacuoles, effecting on their content in pepper fruits. Our results provide a more comprehensive understanding of ABC transporter gene family in different Capsicum species while allowing the identification of important candidate genes related to capsaicin content for subsequent functional validation.

Transcriptome sequencing and expression profiling of genes involved in the response to abiotic stress in Medicago ruthenica

Medicago ruthenica is a perennial forage legume with the remarkable ability to survive under unfavorable environmental conditions. It has been identified as an excellent species of Medicago that can adapt to various environmental stresses including low temperature, drought, and salinity. To investigate its potential as a genetic resource, we performed transcriptome sequencing and analysis in M. ruthenica under abiotic stresses. We generated >120 million reads from six cDNA libraries, resulting in 79,249 unique transcripts, most of which were highly similar to transcripts from M. truncatula (44,608, 56.3%) and alfalfa (M. sativa, 48,023, 60.6%). Based on gene expression profiles, 2,721 transcripts were identified as abiotic stress responsive genes which were predicted to be mainly involved in phytohormone signaling pathways, transcriptional regulation, and ROS-scavenging. These results suggest that they play critical roles in the response to abiotic stress. In summary, we identified genes in our transcriptome dataset involved in the regulation of the abiotic stress response in M. ruthenica which will provide a valuable resource for the future identification and functional analysis of candidate genes for adaption to unfavorable conditions. The genes identified here could be also useful for improving stress tolerance traits in alfalfa through molecular breeding in the future.

Genome-wide DNA Methylation analysis in response to salinity in the model plant caliph medic (Medicago truncatula)

Background

DNA methylation has a potential role in controlling gene expression and may, therefore, contribute to salinity adaptation in plants. Caliph medic (Medicago truncatula) is a model legume of moderate salinity tolerance capacity; however, a base-resolution DNA methylome map is not yet available for this plant.

Results

In this report, a differential whole-genome bisulfite sequencing (WGBS) was carried out using DNA samples extracted from root tissues exposed to either control or saline conditions. Around 50 million differentially methylated sites (DMSs) were recognized, 7% of which were significantly (p < 0.05, FDR < 0.05) altered in response to salinity. This analysis showed that 77.0% of the contexts of DMSs were mCHH, while only 9.1% and 13.9% were mCHG and mCG, respectively. The average change in methylation level was increased in all sequence contexts, ranging from 3.8 to 10.2% due to salinity stress. However, collectively, the level of the DNA methylation in the gene body slightly decreased in response to salinity treatment. The global increase in DNA methylation due to salinity was confirmed by mass spectrometry analysis. Gene expression analysis using qPCR did not reveal a constant relationship between the level of mCG methylation and the transcription abundance of some genes of potential importance in salinity tolerance, such as the potassium channel KAT3, the vacuolar H⁺-pyrophosphatase (V-PPase), and the AP2/ERF and bZIP transcription factors, implying the involvement of other epigenetic gene expression controllers. Computational functional prediction of the annotated genes that embrace DMSs revealed the presence of enzymes with potential cellular functions in biological processes associated with salinity tolerance mechanisms.

Conclusions

The information obtained from this study illustrates the effect of salinity on DNA methylation and shows how plants can remodel the landscape of 5-methylcytosine nucleotide (5-mC) in the DNA across gene structures, in response to salinity. This remodeling varies between gene regions and between 5-mC sequence contexts. The mCG has a vague impact on the expression levels of a few selected potentially important genes in salt tolerant mechanisms.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4484-5) contains supplementary material, which is available to authorized users.

Novel Alleles of Two Tightly Linked Genes Encoding Polygalacturonase-Inhibiting Proteins (VrPGIP1 and VrPGIP2) Associated with the Br Locus That Confer Bruchid (Callosobruchus spp.) Resistance to Mungbean (Vigna radiata) Accession V2709

Nearly all mungbean cultivars are completely susceptible to seed bruchids (Callosobruchus chinensis and Callosobruchus maculatus). Breeding bruchid-resistant mungbean is a major goal in mungbean breeding programs. Recently, we demonstrated in mungbean (Vigna radiata) accession V2802 that VrPGIP2, which encodes a polygalacturonase inhibiting protein (PGIP), is the Br locus responsible for resistance to C. chinensis and C. maculatus. In this study, mapping in mungbean accession V2709 using a BC₁₁F₂ population of 355 individuals revealed that a single major quantitative trait locus, which controlled resistance to both C. chinensis and C. maculatus, was located in a 237.35 Kb region of mungbean chromosome 5 that contained eight annotated genes, including VrPGIP1 (LOC106760236) and VrPGIP2 (LOC106760237). VrPGIP1 and VrPGIP2 are located next to each other and are only 27.56 Kb apart. Sequencing VrPGIP1 and VrPGIP2 in "V2709" revealed new alleles for both VrPGIP1 and VrPGIP2, named VrPGIP1-1 and VrPGIP2-2, respectively. VrPGIP2-2 has one single nucleotide polymorphism (SNP) at position 554 of wild type VrPGIP2. This SNP is a guanine to cystine substitution and causes a proline to arginine change at residue 185 in the VrPGIP2 of "V2709". VrPGIP1-1 has 43 SNPs compared with wild type and "V2802", and 20 cause amino acid changes in VrPGIP1. One change is threonine to proline at residue 185 in VrPGIP1, which is the same as in VrPGIP2. Sequence alignments of VrPGIP2 and VrPGIP1 from "V2709" with common bean (Phaseolus vulgaris) PGIP2 revealed that residue 185 in VrPGIP2 and VrPGIP1 contributes to the secondary structures of proteins that affect interactions between PGIP and polygalacturonase, and that some amino acid changes in VrPGIP1 also affect interactions between PGIP and polygalacturonase. Thus, tightly linked VrPGIP1 and VrPGIP2 are the likely genes at the Br locus that confer bruchid resistance in mungbean "V2709".

Genome-Wide Identification, Evolution, and Expression Analysis of the ATP-Binding Cassette Transporter Gene Family in Brassica rapa

ATP-binding cassette (ABC) proteins can act as transporters of different substrates across biological membranes by hydrolyzing ATP. However, little information is available about ABC transporters in Brassica rapa, an important leafy vegetable. In the present study, we carried out genome-wide identification, characterization and molecular evolution analyses of ABC gene family in B. rapa and 9 other plant species. A total of 179 B. rapa ABC genes (BraABCs) were identified. Among them, 173 BraABCs were identified on 10 chromosomes. Based on phylogenetic analysis and domain organization, the BraABC family could be grouped into eight subfamilies. BraABCs in the same subfamily showed similar motif composition and exon-intron organization. Common and unique cis-elements involved in the transcriptional regulation were also identified in the promoter regions of BraABCs. Tissue-expression analysis of BraABCs demonstrated their diverse spatiotemporal expression profiles. Influences of the whole genome triplication (WGT) on the evolution of BraABCs were studied in detail. BraABCs were preferentially retained compared with their neighboring genes during diploidization after WGT. Synteny analysis identified 76 pairs of syntenic BraABC paralogs among the three subgenomes of B. rapa, and 10 paralog pairs underwent positive selection with ω (= Ka/Ks) ratios greater than 1. Analyses of the expression patterns of syntenic BraABC paralogs pairs across five tissues and under stress treatments revealed their functional conservation, sub-functionalization, neo-functionalization and pseudogenization during evolution. Our study presents a comprehensive overview of the ABC gene family in B. rapa and will be helpful for the further functional study of BraABCs in plant growth, development, and stress responses.