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

Integrative analysis indicates the potential values of ANKRD53 in stomach adenocarcinoma

BACKGROUND

Ankyrin repeat domain 53 (ANKRD53) plays an important role in maintaining chromosome integrity and stability, and chromosome instability is associated with cancer. Through integrative analysis, this study investigates the potential value of ANKRD53 in stomach adenocarcinoma (STAD).

METHODS

RNA-seq and scRNA-seq data were used for integrative analysis based on online databases. Expression of ANKRD53 was confirmed by RT-PCR after bioinformatic analysis. Kaplan-Meier and Cox regression analyses were performed to evaluate the prognostic value of ANKRD53 in STAD. Gene set enrichment analysis (GSEA) was performed to evaluate ANKRD53-related signaling pathways. In addition, the interaction of ANKRD53 with immunity was also investigated.

RESULTS

RT-PCR in STAD cell lines confirmed that ANKRD53 was downregulated in STAD samples compared to normal samples in the online databases. As an independent predictive biomarker, ANKRD53 was combined with other clinicopathological parameters to create a prognostic nomogram. Using GSEA, ANKRD53 was found to be involved in five pathways, including the TGF-β signaling pathway. Further investigation revealed that ANKRD53 was associated with immune checkpoint molecules, immunological pathways, and immunotherapy, in addition to MSI, TMB and neoantigens. In addition, scRNA-seq data revealed that ANKRD53 is mainly expressed in CD8+ T and dendritic cells.

CONCLUSIONS

ANKRD53 is an important biomarker for STAD that deserves further attention.

Identification and expression analyses of B3 genes reveal lineage-specific evolution and potential roles of REM genes in pepper

BACKGROUND

The B3 gene family, one of the largest plant-specific transcription factors, plays important roles in plant growth, seed development, and hormones. However, the B3 gene family, especially the REM subfamily, has not been systematically and functionally studied.

RESULTS

In this study, we performed genome-wide re-annotation of B3 genes in five Solanaceae plants, Arabidopsis thaliana, and Oryza sativa, and finally predicted 1,039 B3 genes, including 231 (22.2%) newly annotated genes. We found a striking abundance of REM genes in pepper species (Capsicum annuum, Capsicum baccatum, and Capsicum chinense). Comparative motif analysis revealed that REM and other subfamilies (ABI3/VP1, ARF, RAV, and HSI) consist of different amino acids. We verified that the large number of REM genes in pepper were included in the specific subgroup (G8) through the phylogenetic analysis. Chromosome location and evolutionary analyses suggested that the G8 subgroup genes evolved mainly via a pepper-specific recent tandem duplication on chromosomes 1 and 3 after speciation between pepper and other Solanaceae. RNA-seq analyses suggested the potential functions of REM genes under salt, heat, cold, and mannitol stress conditions in pepper (C. annuum).

CONCLUSIONS

Our study provides evolutionary and functional insights into the REM gene family in pepper.

De Novo Assembly of the Genome of the Sea Urchin Paracentrotus lividus (Lamarck 1816)

The Mediterranean purple sea urchin Paracentrotus lividus (Lamarck 1816) is a remarkable model system for molecular, evolutionary and cell biology studies, particularly in the field of developmental biology. We sequenced the genome, performed a de novo assembly, and analysed the assembly content. The genome of P. lividus was sequenced using Illumina NextSeq 500 System (Illumina) in a 2 × 150 paired-end format. More than 30,000 open reading frames (ORFs), (more than 8000 are unique), were identified and analysed to provide molecular tools accessible for the scientific community. In particular, several genes involved in complex innate immune responses, oxidative metabolism, signal transduction, and kinome, as well as genes regulating the membrane receptors, were identified in the P. lividus genome. In this way, the employment of the Mediterranean sea urchin for investigations and comparative analyses was empowered, leading to the explanation of cis-regulatory networks and their evolution in a key developmental model occupying an important evolutionary position with respect to vertebrates and humans.

Abundance of distal repetitive DNA sequences in Capsicum L. (Solanaceae) chromosomes

Chili peppers (Solanaceae family) have great commercial value. They are commercialized in natura and used as spices and for ornamental and medicinal purposes. Although three whole genomes have been published, limited information about satellite DNA sequences, their composition, and genomic distribution has been provided. Here, we exploited the noncoding repetitive fraction, represented by satellite sequences, that tends to accumulate in blocks along chromosomes, especially near the chromosome ends of peppers. Two satellite DNA sequences were identified (CDR-1 and CDR-2), characterized and mapped in silico in three Capsicum genomes (C. annuum, C. chinense, and C. baccatum) using data from the published high-coverage sequencing and repeats finding bioinformatic tools. Localization using FISH in the chromosomes of these species and in two others (C. frutescens and C. chacoense), totaling five species, showed signals adjacent to the rDNA sites. A sequence comparison with existing Solanaceae repeats showed that CDR-1 and CDR-2 have different origins but without homology to rDNA sequences. Satellites occupied subterminal chromosomal regions, sometimes collocated with or adjacent to 35S rDNA sequences. Our results expand knowledge about the diversity of subterminal regions of Capsicum chromosomes, showing different amounts and distributions within and between karyotypes. In addition, these sequences may be useful for future phylogenetic studies.

Genetic Analyses of Seed Longevity in Capsicum annuum L. in Cold Storage Conditions

Seed longevity is the most important trait in the genebank management system. No seed can remain infinitely viable. There are 1241 accessions of Capsicum annuum L. available at the German Federal ex situ genebank at IPK Gatersleben. C. annuum (Capsicum) is the most economically important species of the genus Capsicum. So far, there is no report that has addressed the genetic basis of seed longevity in Capsicum. Here, we convened a total of 1152 Capsicum accessions that were deposited in Gatersleben over forty years (from 1976 to 2017) and assessed their longevity by analyzing the standard germination percentage after 5-40 years of storage at -15/-18 °C. These data were used to determine the genetic causes of seed longevity, along with 23,462 single nucleotide polymorphism (SNP) markers covering all of the 12 Capsicum chromosomes. Using the association-mapping approach, we identified a total of 224 marker trait associations (MTAs) (34, 25, 31, 35, 39, 7, 21 and 32 MTAs after 5-, 10-, 15-, 20-, 25-, 30-, 35- and 40-year storage intervals) on all the Capsicum chromosomes. Several candidate genes were identified using the blast analysis of SNPs, and these candidate genes are discussed.

Comparative Transcriptome Analysis of CCCH Family in Roles of Flower Opening and Abiotic Stress in Osmanthus fragrans

CCCH is a zinc finger family with a typical CCCH-type motif which performs a variety of roles in plant growth and development and responses to environmental stressors. However, the information about this family has not been reported for Osmanthus fragrans. In this study, a total of 66 CCCH predicted genes were identified from the O. fragrans genome, the majority of which had multiple CCCH motifs. The 66 OfCCCHs were found to be unevenly distributed on 21 chromosomes and were clustered into nine groups based on their phylogenetic analysis. In each group, the gene structure and domain makeup were comparatively conserved. The expression profiles of the OfCCCH genes were examined in various tissues, the flower-opening processes, and under various abiotic stresses using transcriptome sequencing and qRT-PCR (quantitative real-time PCR). The results demonstrated the widespread expression of OfCCCHs in various tissues, the differential expression of 22 OfCCCHs during flower-opening stages, and the identification of 4, 5, and 13 OfCCCHs after ABA, salt, and drought stress treatment, respectively. Furthermore, characterization of the representative OfCCCHs (OfCCCH8, 23, 27, and 36) revealed that they were all localized in the nucleus and that the majority of them had transcriptional activation in the yeast system. Our research offers the first thorough examination of the OfCCCH family and lays the groundwork for future investigations regarding the functions of CCCH genes in O. fragrans.

Identification of genomic regions of wheat associated with grain Fe and Zn content under drought and heat stress using genome-wide association study

Wheat is the staple food crop of global importance for its grain nutrient quality. Grain iron and zinc content of the wheat grain is an important quantitatively inherited trait that is influenced by the environmental factors such as drought and heat stress. Phenotypic evaluation of 295 advanced breeding lines from the wheat stress breeding program of IARI was carried out under timely sown irrigated (IR), restricted irrigated, and late-sown conditions at New Delhi during the cropping season of 2020-21, and grain iron (GFeC) and zinc (GZnC) contents were estimated from both control and treatments. A statistically significant increase in GFeC and GZnC was observed under stress conditions compared to that of the control. Genotyping was carried out with the SNPs from the 35K Axiom Breeder's array, and marker-trait association was identified by GWAS analysis. Of the 23 MTAs identified, seven were linked with GFeC and sixteen were linked with GZnC. In silico analysis revealed a few important transcripts involved in various plant metabolism, growth, and development activities such as auxin response factor, root UVB sensitive proteins, potassium transporter, glycosyl transferase, COBRA, and F-box-like domain. The identified MTAs can be used for molecular breeding after validation and also for rapid development of micronutrient-rich varieties of wheat to mitigate hidden hunger.

Genetic Mapping to Detect Stringent QTLs Using 1k-RiCA SNP Genotyping Platform from the New Landrace Associated with Salt Tolerance at the Seedling Stage in Rice

Rice is the world’s most important food crop, providing the daily calorie intake for more than half of the world’s population. Rice breeding has always been preoccupied with maximizing yield potential. However, numerous abiotic factors, such as salt, cold, drought, and heat, significantly reduce rice productivity. Salinity, one of the major abiotic stresses, reduces rice yield worldwide. This study was conducted to determine new quantitative trait loci (QTLs) that regulate salt tolerance in rice seedlings. One F_2:3 mapping population was derived from a cross between BRRI dhan49 (a popular but sensitive rainfed rice variety) and Akundi (a salt-tolerant rice landrace in Bangladesh used as a donor parent). The 1k-Rice Custom Amplicon (1k-RiCA) single-nucleotide polymorphism (SNP) markers were used to genotype this mapping population. After removing segregation distortion and monomorphic markers, 884 SNPs generated a 1526.8 cM-long genetic linkage map with a mean marker density of 1.7 cM for the 12 linkage groups. By exploiting QGene and ICIM-ADD, a sum of 15 QTLs for nine traits was identified in salt stress on seven chromosomes. Four important genomic loci were identified (qSES1, qSL1, qSUR1 and qRL1) on chromosome 1. Out of these 15 QTLs, 14 QTLs are unique, as no other study has mapped in the same chromosomal location. We also detected 15 putative candidate genes and their functions. The ICIM-EPI approach identified 43 significant pairwise epistasis interactions between regions associated with and unassociated with QTLs. Apart from more well-known donors, Akundi serves as an important new donor source for global salt tolerance breeding initiatives, including Bangladesh. The introgression of the novel QTLs identified in this study will accelerate the development of new salt-tolerant varieties that are highly resistant to salt stress using marker-enabled breeding.

Magnesium- a Forgotten Element: Phenotypic Variation and Genome Wide Association Study in Turkish Common Bean Germplasm

Magnesium (Mg) is the fourth most abundant element in the human body and plays the role of cofactor for more than 300 enzymatic reactions. In plants, Mg is involved in various key physiological and biochemical processes like growth, development, photophosphorylation, chlorophyll formation, protein synthesis, and resistance to biotic and abiotic stresses. Keeping in view the importance of this element, the present investigation aimed to explore the Mg contents diversity in the seeds of Turkish common bean germplasm and to identify the genomic regions associated with this element. A total of 183 common bean accessions collected from 19 provinces of Turkey were used as plant material. Field experiments were conducted according to an augmented block design during 2018 in two provinces of Turkey, and six commercial cultivars were used as a control group. Analysis of variance depicted that Mg concentration among common bean accessions was statistically significant (p < 0.05) within each environment, however genotype × environment interaction was non-significant. A moderate level (0.60) of heritability was found in this study. Overall mean Mg contents for both environments varied from 0.33 for Nigde-Dermasyon to 1.52 mg kg⁻¹ for Nigde-Derinkuyu landraces, while gross mean Mg contents were 0.92 mg kg⁻¹. At the province level, landraces from Bolu were rich while the landraces from Bitlis were poor in seed Mg contents respectively. The cluster constellation plot divided the studied germplasm into two populations on the basis of their Mg contents. Marker-trait association was performed using a mixed linear model (Q + K) with a total of 7,900 DArTseq markers. A total of six markers present on various chromosomes (two at Pv01, and one marker at each chromosome i.e., Pv03, Pv07, Pv08, Pv11) showed statistically significant association for seed Mg contents. Among these identified markers, the DArT-3367607 marker present on chromosome Pv03 contributed to maximum phenotypic variation (7.5%). Additionally, this marker was found within a narrow region of previously reported markers. We are confident that the results of this study will contribute significantly to start common bean breeding activities using marker assisted selection regarding improved Mg contents.

PbANK facilitates the long-distance movement of the PbWoxT1-PbPTB3 RNP complex by degrading deposited callose

Grafting horticultural crops can result in phenotypic changes in the grafted materials due to the movement of macromolecular signals, including RNAs and proteins, across the graft union; however, little is known about the composition of trafficking ribonucleoprotein (RNP) complexes or how these macromolecules are transported. Here, we used the core of PbPTB3-PbWoxT1 RNP complex, PbPTB3, as bait to screen Pyrus betulaefolia cDNA library for its interaction partners. We identified an ankyrin protein, PbANK, that interacts with PbPTB3 to facilitate its transport through the phloem alongside PbWoxT1 mRNA. Heterografting experiments showed that silencing PbANK in rootstock prevented the transport of PbPTB3 and PbWoxT1 mRNA from the rootstock to the scion. Similarly, heterologous grafting experiments demonstrated that PbANK itself cannot be transported over long distances through a graft union. Fluorescence microscopy showed that silencing ANK affected the intercellular diffusion of PbPTB3 and increased callose deposition at plasmodesmata. Collectively, these findings demonstrate that PbANK mediates the long-distance movement of PbPTB3 and PbWoxT1 by degrading callose to increase the efficiency of cell-to-cell movement.

Genome-Wide Identification and Characterization of DnaJ Gene Family in Grape (Vitis vinifera L.)

Grape production in southern China suffers great loss due to various environmental stresses. To understand the mechanism of how the grape plants respond to these stresses is an active area of research in developing cultivation techniques. Plant stress resistance is known to rely on special proteins. Amongst them, DnaJ protein (HSP40) serves as co-chaperones of HSP70, playing crucial roles in various stress response. However, the DnaJ proteins encoded by the DnaJ gene family in Vitis vinifera L. have not been fully described yet. In this study, we identified 78 VvDnaJs in the grape genome that can be classified into three groups—namely, DJA, DJB, and DJC. To reveal the evolutionary and stress response mechanisms for the VvDnaJ gene family, their evolutionary and expression patterns were analyzed using the bioinformatic approach and qRT-PCR. We found that the members in the same group exhibited a similar gene structure and protein domain organization. Gene duplication analysis demonstrated that segmental and tandem duplication may not be the dominant pathway of gene expansion in the VvDnaJ gene family. Codon usage pattern analysis showed that the codon usage pattern of VvDnaJs differs obviously from the monocotyledon counterparts. Tissue-specific analysis revealed that 12 VvDnaJs present a distinct expression profile, implying their distinct roles in various tissues. Cis-acting element analysis showed that almost all VvDnaJs contained the elements responsive to either hormones or stresses. Therefore, the expression levels of VvDnaJs subjected to exogenous hormone applications and stress treatments were determined, and we found that VvDnaJs were sensitive to hormone treatments and shade, salt, and heat stresses, especially VIT_00s0324g00040. The findings of this study could provide comprehensive information for the further investigation on the genetics and protein functions of the DnaJ gene family in grape.

Study of Triticum aestivum Resistome in Response to Wheat dwarf India Virus Infection

Susceptible and resistant germplasm respond differently to pathogenic attack, including virus infections. We compared the transcriptome changes between a resistant wheat cultivar, Sonalika, and a susceptible cultivar, WL711, to understand this process in wheat against wheat dwarf India virus (WDIV) infection. A total of 2760 and 1853 genes were differentially expressed in virus-infected and mock-inoculated Sonalika, respectively, compared to WL711. The overrepresentation of genes involved in signaling, hormone metabolism, enzymes, secondary metabolites, proteolysis, and transcription factors was documented, including the overexpression of multiple PR proteins. We hypothesize that the virus resistance in Sonalika is likely due to strong intracellular surveillance via the action of multiple PR proteins (PR1, RAR1, and RPM1) and ChiB. Other genes such as PIP1, LIP1, DnaJ, defensins, oxalate oxidase, ankyrin repeat protein, serine-threonine kinase, SR proteins, beta-1,3-glucanases, and O-methyltransferases had a significant differential expression and play roles in stress tolerance, may also be contributing towards the virus resistance in Sonalika. In addition, we identified putative genes with unknown functions, which are only expressed in response to WDIV infection in Sonalika. The role of these genes could be further validated and utilized in engineering resistance in wheat and other crops.

Association Mapping of Verticillium Wilt Disease in a Worldwide Collection of Cotton (Gossypium hirsutum L.)

Cotton (Gossypium spp.) is the best plant fiber source in the world and provides the raw material for industry. Verticillium wilt caused by Verticillium dahliae Kleb. is accepted as a major disease of cotton production. The most practical way to deal with verticillium wilt is to develop resistant/tolerant varieties after cultural practices. One of the effective selections in plant breeding is the use of marker-assisted selection (MAS) via quantitative trait loci (QTL). Therefore, in this study, we aimed to discover the genetic markers associated with the disease. Through the association mapping analysis, common single nucleotide polymorphism (SNP) markers were obtained using 4730 SNP alleles. As a result, twenty-three markers were associated with defoliating (PYDV6 isolate) pathotype, twenty-one markers with non-defoliating (Vd11 isolate) pathotype, ten QTL with Disease Severity Index (DSI) of the leaves at the 50–60% boll opening period and eight markers were associated with DSI in the stem section. Some of the markers that show significant associations are located on protein coding genes such as protein Mpv17-like, 21 kDa protein-like, transcription factor MYB113-like, protein dehydration-induced 19 homolog 3-like, F-box protein CPR30-like, extracellular ribonuclease LE-like, putative E3 ubiquitin-protein ligase LIN, pentatricopeptide repeat-containing protein At3g62890-like, fructose-1,6-bisphosphatase, tubby-like F-box protein 8, endoglucanase 16-like, glucose-6-phosphate/phosphate translocator 2, metal tolerance protein 11-like, VAN3-binding protein-like, transformation/transcription domain-associated protein-like, pyruvate kinase isozyme A, ethylene-responsive transcription factor CRF2-like, molybdate transporter 2-like, IRK-interacting protein-like, glycosylphosphatidylinositol anchor attachment 1 protein, U3 small nucleolar RNA-associated protein 4-like, microtubule-associated protein futsch-like, transport and Golgi organization 2 homolog, splicing factor 3B subunit 3-like, mediator of RNA polymerase II transcription subunit 15a-like, putative ankyrin repeat protein, and protein networked 1D-like. It has been reported in previous studies that most of these genes are associated with biotic and abiotic stress factors. As a result, once validated, it would be possible to use the markers obtained in the study in Marker Assisted Selection (MAS) breeding.

The Ankyrin-Repeat Gene GmANK114 Confers Drought and Salt Tolerance in Arabidopsis and Soybean

Ankyrin repeat (ANK) proteins are essential in cell growth, development, and response to hormones and environmental stresses. In the present study, 226 ANK genes were identified and classified into nine subfamilies according to conserved domains in the soybean genome (Glycine max L.). Among them, the GmANK114 was highly induced by drought, salt, and abscisic acid. The GmANK114 encodes a protein that belongs to the ANK-RF subfamily containing a RING finger (RF) domain in addition to the ankyrin repeats. Heterologous overexpression of GmANK114 in transgenic Arabidopsis improved the germination rate under drought and salt treatments compared to wild-type. Homologous overexpression of GmANK114 improved the survival rate under drought and salt stresses in transgenic soybean hairy roots. In response to drought or salt stress, GmANK114 overexpression in soybean hairy root showed higher proline and lower malondialdehyde contents, and lower H2O2 and O2- contents compared control plants. Besides, GmANK114 activated transcription of several abiotic stress-related genes, including WRKY13, NAC11, DREB2, MYB84, and bZIP44 under drought and salt stresses in soybean. These results provide new insights for functional analysis of soybean ANK proteins and will be helpful for further understanding how ANK proteins in plants adapt to abiotic stress.