Paralogous Genes Research Articles

Two glyceraldehyde-3-phosphate dehydrogenases with distinctive roles in Pseudomonas syringae pv. tomato DC3000

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH or Gap) is a ubiquitously distributed enzyme that plays an essential role in the glycolytic and gluconeogenic pathways. However, additional roles have been described unrelated to its enzymatic function in diverse organisms, often linked to its presence in the cell surface or as a secreted protein. Despite being a paradigm among multifunctional/moonlighting proteins, little is known about its possible roles in phytopathogenic bacteria. In the present work we have studied three putative gap paralogous genes identified in the genome of Pseudomonas syringae pv. tomato (Pto) DC3000, an important model in molecular plant pathology, with the aim of determining their physiological and possible non-canonical roles in this bacterium and in the plant infection process. We have established that the Gap1 protein has a predominantly glycolytic activity, whereas the NADPH-dependent Gap2 main activity is gluconeogenic. The third paralogue lacks GAPDH activity in Pto but is indispensable for vitamin B6 metabolism and displays erythrose-4-phosphate dehydrogenase activity, thus referred as epd. Both Gap enzymes exhibit distinct functional characteristics depending on the bacterium physiological state, with Gap1 presenting a substantial role in motility, biosurfactant production and biofilm formation. On the other hand, solely Gap2 appears to be essential for growth on tomato plant. Furthermore, Gap1 and Gap2 present a distinctive transcriptional regulation and both have been identified exported outside the cells with different definite media compositions. This serves as compelling evidence of additional roles beyond their central metabolic functions.

The promoters of two CpMYB106-like paralog genes respond to abiotic stresses and phytohormones and drive differential expression in contrasting cultivars of Cucurbita pepo

Zucchini fruit (Cucurbita pepo L.) is susceptible to chilling injury (CI) during its postharvest life, with this response being cultivar-dependent. A previous transcriptomic analysis comparing cold-tolerant fruit from the cultivar ‘Natura’ with cold-sensitive fruit from the cultivar ‘Sinatra’, revealed the transcription factor (TF) CpMYB106-like as a relevant candidate gene for the acquisition of postharvest low-temperature tolerance. C. pepo has two CpMYB106-like paralogs, which are differentially present in the ‘Natura’ (CpMYB106-likeA) and ‘Sinatra’ (CpMYB106-likeB) cultivars due to gene and promoter deletions. The aim of this study was to investigate the transcriptional regulation of the CpMYB106-like paralogs to unravel their role in the postharvest life of zucchini fruit. For this, gene expression was analyzed and promoter sequences and their activities were studied through reporter β-glucuronidase (GUS) gene analysis in heterologous systems. An expression analysis showed that CpMYB106-like mRNA levels were induced only in the cold-stored fruit of ‘Natura’. The analysis of the promoter sequences showed numerous cis-regulatory elements (CREs), many of which have a differential distribution or frequency between the two paralogs. Promoter basal activity was determined by transient transformation of Nicotiana benthamiana leaves, which showed a higher expression driven by the promoter from the cold-tolerant cultivar ‘Natura’, which could be explained by the presence of three additional copies of the enhancer element EECCRCAH1 in proMYB106A. The tissue-specific expression pattern, and the response to abiotic stresses and phytohormones controlled by CpMYB106-like promoters was analyzed in transgenic Arabidopsis thaliana plants. GUS activity was mainly detected in the vascular system, leaves, and inflorescences, especially in siliques. With respect to abiotic stresses, low temperature decreased GUS expression in Arabidopsis, but increased it in plants carrying the ‘Natura’ promoter after a 24-h acclimation period at 22 °C. Salt treatment induced GUS activity driven by both promoters in Arabidopsis seedlings, but it was higher with the ‘Natura’ promoter. Furthermore, methyl jasmonate (MeJA) also induced a sharp increase of GUS controlled by proMYB106A. These results suggest that the regulation of CpMYB106-like TFs in zucchini is promoter-driven, and reveal the involvement of these genes in the acquisition of cold tolerance in fruit during postharvest stress conditions.

Identification, Evolution and Expression Analysis of GRF Family Reveals Their Involvement in Shoot Growth and Abiotic Stress Response in Moso Bamboo

Growth-regulating factors (GRFs) play an important role in regulating plant organ development, acting primarily as positive regulators of cell proliferation. However, research on the evolutionary history and expression patterns of the moso bamboo GRF family has been limited. In this study, a total of 24 GRFs have been identified in the Moso bamboo genome, and they have been categorized into four subfamilies. Estimation of the divergence time of paralogous gene pairs provided evidence supporting the significant contribution of recent whole-genome duplication events in the expansion of the GRF gene family. Sliding window analysis revealed that coding regions of a few PheGRFs, including the WRC and QLQ domains, may have undergone positive selection, possibly due to the redundant functions of paralogous genes. Coexpression network analysis further revealed the regulatory role of miR396 and various lncRNAs in controlling PheGRF expression. Based on the analysis of tissue-specific expression patterns using transcriptome sequencing, qRT-PCR results, and in situ hybridization, it was observed that most GRFs, particularly PheGRF6a and PheGRF9b, exhibited high levels of accumulation in the moso bamboo shoot. This suggests that the involvement of most PheGRF genes may be crucial for the growth and development of the bamboo shoot. A yeast two-hybrid screening revealed interactions between PheGRF9b and several proteins associated with plant growth and development, including PH02Gene11097.t1 (GIF3), PH02Gene37618.t (Phytochrome B), and PH02Gene01921.t3 (WD40). Based on transcriptome expression analysis, it was observed that a substantial number of PheGRFs exhibited significant variations under cold or drought stress treatments, and most of these genes were found to be downregulated, suggesting their role as abiotic stress-responsive genes. Our findings offer new insights into the GRF family of moso bamboo and provide some experimental evidence to support further gene functional validation research of PheGRF.

Comprehensive analyses of annexins in naked carp (Gymnocypris przewalskii) unveil their roles in saline-alkaline stress

Naked carp (Gymnocypris przewalskii) is a migratory fish endemic to the Qinghai Lake, the largest inland saline-alkaline lake in China. Based on our previous results of comparative transcriptome sequencing, we found that pathways related to regulation of actin cytoskeleton, plasma membrane damage, and inflammation-related signaling were significantly enriched during saline-alkaline stress, suggesting involvement of annexin (ANN) family proteins in stress adaption in naked carp. Herein, to explore the role of annexin family in naked carp under saline-alkaline stress, twenty GpANNs family members were identified at the genome-wide scale, and bioinformatics and expression analysis under saline-alkaline stress were performed. The predicted subcellular localization showed that most members of GpANNs were distributed in the cytoplasm, and the rest were localized in the cytoskeleton, nucleus, or mitochondria, respectively. GpANNs family proteins contained specific type II Ca2+-binding motifs, actin-binding motifs, GTP-binding motifs, etc. Phylogenetic tree analysis showed that the annexin homologues from naked carp, carp (Cyprinus carpio), zebrafish (Danio rerio), and chinese cavefish (Sinocyclocheilus rhinocerous) were divided into eight groups, i.e. ANN1 ∼ 6, ANN11 and ANN13. But homologues of ANN7 ∼ 10 groups were not found in these three fish species. Twenty-eight, nineteen and twenty-nine orthologous gene pairs were identified between naked carp and carp, zebrafish and chinese cavefish, respectively. Six paralogous gene pairs were identified in naked carp, all of which had a Ka/Ks ratio less than one, suggesting that the GpANNs family may have undergone purifying selection pressure during evolution. Expression of GpANN3 was significantly up-regulated in the kidney tissue under saline-alkaline stress, while GpANN2 was significantly down-regulated. Similar trends of protein expression were also observed by Western blot, suggesting involvement of GpANN3 and GpANN2 in stress response. Confocal microscopy revealed that GpANN2 was localized in both cytoplasm and nucleus, but GpANN3 was preferentially localized in the cytoplasm. Our study demonstrated the role of annexins in stress adaption, and provided the basic information required for further function analyses in naked carp.

Dissection of the role of a Src homology 3 domain in the evolution of binding preference of paralogous proteins.

Protein-protein interactions (PPIs) drive many cellular processes. Some interactions are directed by Src homology 3 (SH3) domains that bind proline-rich motifs on other proteins. The evolution of the binding specificity of SH3 domains is not completely understood, particularly following gene duplication. Paralogous genes accumulate mutations that can modify protein functions and, for SH3 domains, their binding preferences. Here, we examined how the binding of the SH3 domains of 2 paralogous yeast type I myosins, Myo3 and Myo5, evolved following duplication. We found that the paralogs have subtly different SH3-dependent interaction profiles. However, by swapping SH3 domains between the paralogs and characterizing the SH3 domains freed from their protein context, we find that very few of the differences in interactions, if any, depend on the SH3 domains themselves. We used ancestral sequence reconstruction to resurrect the preduplication SH3 domains and examined, moving back in time, how the binding preference changed. Although the most recent ancestor of the 2 domains had a very similar binding preference as the extant ones, older ancestral domains displayed a gradual loss of interaction with the modern interaction partners when inserted in the extant paralogs. Molecular docking and experimental characterization of the free ancestral domains showed that their affinity with the proline motifs is likely not the cause for this loss of binding. Taken together, our results suggest that a SH3 and its host protein could create intramolecular or allosteric interactions essential for the SH3-dependent PPIs, making domains not functionally equivalent even when they have the same binding specificity.

Subfunctionalisation of paralogous genes and evolution of differential codon usage preferences: The showcase of polypyrimidine tract binding proteins.

Gene paralogs are copies of an ancestral gene that appear after gene or full genome duplication. When two sister gene copies are maintained in the genome, redundancy may release certain evolutionary pressures, allowing one of them to access novel functions. Here, we focused our study on gene paralogs on the evolutionary history of the three polypyrimidine tract binding protein genes (PTBP) and their concurrent evolution of differential codon usage preferences (CUPrefs) in vertebrate species. PTBP1-3 show high identity at the amino acid level (up to 80%) but display strongly different nucleotide composition, divergent CUPrefs and, in humans and in many other vertebrates, distinct tissue-specific expression levels. Our phylogenetic inference results show that the duplication events leading to the three extant PTBP1-3 lineages predate the basal diversification within vertebrates, and genomic context analysis illustrates that local synteny has been well preserved over time for the three paralogs. We identify a distinct evolutionary pattern towards GC3-enriching substitutions in PTBP1, concurrent with enrichment in frequently used codons and with a tissue-wide expression. In contrast, PTBP2s are enriched in AT-ending, rare codons, and display tissue-restricted expression. As a result of this substitution trend, CUPrefs sharply differ between mammalian PTBP1s and the rest of PTBPs. Genomic context analysis suggests that GC3-rich nucleotide composition in PTBP1s is driven by local substitution processes, while the evidence in this direction is thinner for PTBP2-3. An actual lack of co-variation between the observed GC composition of PTBP2-3 and that of the surrounding non-coding genomic environment would raise an interrogation on the origin of CUPrefs, warranting further research on a putative tissue-specific translational selection. Finally, we communicate an intriguing trend for the use of the UUG-Leu codon, which matches the trends of AT-ending codons. Our results are compatible with a scenario in which a combination of directional mutation-selection processes would have differentially shaped CUPrefs of PTBPs in vertebrates: the observed GC-enrichment of PTBP1 in placental mammals may be linked to genomic location and to the strong and broad tissue-expression, while AT-enrichment of PTBP2 and PTBP3 would be associated with rare CUPrefs and thus, possibly to specialized spatio-temporal expression. Our interpretation is coherent with a gene subfunctionalisation process by differential expression regulation associated with the evolution of specific CUPrefs.

Efficient repair of human genetic defect by CRISPR/Cas9-mediated interlocus gene conversion

Abstract DNA double-strand breaks (DSBs) induced by gene-editing tools are primarily repaired through non-homologous end joining (NHEJ) or homology-directed repair (HDR) using synthetic DNA templates. However, error-prone NHEJ may result in unexpected indels at the targeted site. For most genetic disorders, precise HDR correction using exogenous homologous sequence is ideal. But, the therapeutic application of HDR might be especially challenging given the requirement for the codelivery of exogenous DNA templates with toxicity into cells, and the low efficiency of HDR could also limit its clinical application. In this study, we efficiently repair pathogenic mutations in HBB coding regions of hematopoietic stem cells (HSCs) using CRISPR/Cas9-mediated gene conversion (CRISPR/GC) using the paralog gene HBD as the internal template. After transplantation, these edited HSCs successfully repopulate the hematopoietic system and generate erythroid cells with significantly reduced thalassemia propensity. Moreover, a range of pathogenic gene mutations causing β-thalassemia in HBB coding regions were effectively converted to normal wild-type sequences without exogenous DNA templates using CRISPR/GC. This highlights the promising potential of CRISPR/GC, independent of synthetic DNA templates, for genetic disease gene therapy.

Genetic Channelization Mechanism of Four Chalcone Isomerase Homologous Genes for Synergistic Resistance to Fusarium wilt in Gossypium barbadense L.

Duplication events occur very frequently during plant evolution. The genes in the duplicated pathway or network can evolve new functions through neofunctionalization and subfunctionalization. Flavonoids are secondary metabolites involved in plant development and defense. Our previous transcriptomic analysis of F6 recombinant inbred lines (RILs) and the parent lines after Fusarium oxysporum f. sp. vasinfectum (Fov) infection showed that CHI genes have important functions in cotton. However, there are few reports on the possible neofunctionalization differences of CHI family paralogous genes involved in Fusarium wilt resistance in cotton. In this study, the resistance to Fusarium wilt, expression of metabolic pathway-related genes, metabolite content, endogenous hormone content, reactive oxygen species (ROS) content and subcellular localization of four paralogous CHI family genes in cotton were investigated. The results show that the four paralogous CHI family genes may play a synergistic role in Fusarium wilt resistance. These results revealed a genetic channelization mechanism that can regulate the metabolic flux homeostasis of flavonoids under the mediation of endogenous salicylic acid (SA) and methyl jasmonate (MeJA) via the four paralogous CHI genes, thereby achieving disease resistance. Our study provides a theoretical basis for studying the evolutionary patterns of homologous plant genes and using homologous genes for molecular breeding.

Expression Levels of Flowering Time Genes (CONZ1, GIGZ1A, GIGZ1B, FKF1A, and FKF1B) in Seedlings under Long-Day Conditions Differentiates Early and Late Zea mays L. Lines.

Phenophase durations, including the timing of flowering and ripeness, were characterized in 20inbred lines of the maize Zea mays L. Expression of key flowering initiation genes (CONZ1, GIGZ1a, GIGZ1b, ZmFKF1a, and ZmFKF1b) under long-photoperiod conditions was studied in seedlings of six maize lines that differed in ripeness time. Significantly lower transcription levels of all of the five genes were found in early-ripening lines compared with late-ripening lines. Similar expression patterns were observed for the GIGZ1a and GIGZ1b paralogous genes, while ZmFKF1a significantly predominated in expression over its paralog ZmFKF1b.

Wnt9 directs zebrafish heart tube assembly via a combination of canonical and non-canonical pathway signaling.

During zebrafish heart formation, cardiac progenitor cells converge at the embryonic midline where they form the cardiac cone. Subsequently, this structure transforms into a heart tube. Only little is known about the molecular mechanisms that control these morphogenetic processes. Here, we use light sheet microscopy and combine genetic, molecular biological and pharmacological tools to show that the paralogous genes wnt9a/b are required for the assembly of the nascent heart tube. In wnt9a/b double mutants, cardiomyocyte progenitor cells are delayed in their convergence towards the embryonic midline, the formation of the heart cone is impaired and the transformation into an elongated heart tube fails. The same cardiac phenotype occurs when both canonical and non-canonical Wnt signaling pathways are simultaneously blocked by pharmacological inhibition. This demonstrates that Wnt9a/b and canonical and non-canonical Wnt signaling regulate the migration of cardiomyocyte progenitor cells and control the formation of the cardiac tube. This can be partly attributed to their regulation of the timing of cardiac progenitor cell differentiation. Our study demonstrates how these morphogens activate a combination of downstream pathways to direct cardiac morphogenesis.

Multi-target genome editing reduces polyphenol oxidase activity in wheat (Triticum aestivum L.) grains.

Polyphenol oxidases (PPO) are dual activity metalloenzymes that catalyse the production of quinones. In plants, PPO activity may contribute to biotic stress resistance and secondary metabolism but is undesirable for food producers because it causes the discolouration and changes in flavour profiles of products during post-harvest processing. In wheat (Triticum aestivum L.), PPO released from the aleurone layer of the grain during milling results in the discolouration of flour, dough, and end-use products, reducing their value. Loss-of-function mutations in the PPO1 and PPO2 paralogous genes on homoeologous group 2 chromosomes confer reduced PPO activity in the wheat grain. However, limited natural variation and the proximity of these genes complicates the selection of extremely low-PPO wheat varieties by recombination. The goal of the current study was to edit all copies of PPO1 and PPO2 to drive extreme reductions in PPO grain activity in elite wheat varieties. A CRISPR/Cas9 construct with one single guide RNA (sgRNA) targeting a conserved copper binding domain was used to edit all seven PPO1 and PPO2 genes in the spring wheat cultivar 'Fielder'. Five of the seven edited T1 lines exhibited significant reductions in PPO activity, and T2 lines had PPO activity up to 86.7% lower than wild-type. The same construct was transformed into the elite winter wheat cultivars 'Guardian' and 'Steamboat', which have five PPO1 and PPO2 genes. In these varieties PPO activity was reduced by >90% in both T1 and T2 lines. In all three varieties, dough samples from edited lines exhibited reduced browning. This study demonstrates that multi-target editing at late stages of variety development could complement selection for beneficial alleles in crop breeding programs by inducing novel variation in loci inaccessible to recombination.

Population Genomics of Two Closely Related Anhydrobiotic Midges Reveals Differences in Adaptation to Extreme Desiccation.

The sleeping chironomid Polypedilum vanderplanki is capable of anhydrobiosis, a striking example of adaptation to extreme desiccation. Tolerance to complete desiccation in this species is associated with emergence of multiple paralogs of protective genes. One of the gene families highly expressed under anhydrobiosis and involved in this process is protein-L-isoaspartate (D-aspartate) O-methyltransferases (PIMTs). Recently, another closely related midge was discovered, Polypedilum pembai, which is able not only to tolerate desiccation but also to survive multiple desiccation-rehydration cycles. To investigate the evolution of anhydrobiosis in these species, we sequenced and assembled the genome of P. pembai and compared it with P. vanderplanki and also performed a population genomics analysis of several populations of P. vanderplanki and one population of P. pembai. We observe positive selection and radical changes in the genetic architecture of the PIMT locus between the two species, including its amplification in the P. pembai lineage. In particular, PIMT-4, the most highly expressed of these PIMTs, is present in six copies in the P. pembai; these copies differ in expression profiles, suggesting possible sub- or neofunctionalization. The nucleotide diversity of the genomic region carrying these new genes is decreased in P. pembai, but not in the orthologous region carrying the ancestral gene in P. vanderplanki, providing evidence for a selective sweep associated with postduplication adaptation in the former. Overall, our results suggest an extensive relatively recent and likely ongoing adaptation of the mechanisms of anhydrobiosis.

Genome-Wide Identification and Analysis of TCP Gene Family among Three Dendrobium Species.

Dendrobium orchids, which are among the most well-known species of orchids, are appreciated for their aesthetic appeal across the globe. Furthermore, due to their strict living conditions, they have accumulated high levels of active ingredients, resulting not only in their medicinal value but also in their strong ability to respond to harsh environments. The TCP gene family plays an important role in plant growth and development, and signal transduction. However, these genes have not been systematically investigated in Dendrobium species. In this study, we detected a total of 24, 23, and 14 candidate TCP members in the genome sequences of D. officinale, D. nobile, and D. chrysotoxum, respectively. These genes were classified into three clades on the basis of a phylogenetic analysis. The TCP gene numbers among Dendrobium species were still highly variable due to the independent loss of genes in the CIN clade. However, only three gene duplication events were detected, with only one tandem duplication event (DcTCP9/DcTCP10) in D. chrysotoxum and two pairs of paralogous DoTCP gene duplication events (DoTCP1/DoTCP23 and DoTCP16/DoTCP24) in D. officinale. A total of 25 cis-acting elements of TCPs related to hormone/stress and light responses were detected. Among them, the proportions of hormone response, light response, and stress response elements in D. officinale (100/421, 127/421, and 171/421) were similar to those in D. nobile (83/352, 87/352, and 161/352). Using qRT-PCR to determine their expression patterns under MeJA treatment, four DoTCPs (DoTCP2, DoTCP4, DoTCP6, and DoTCP14) were significantly upregulated under MeJA treatment, which indicates that TCP genes may play important roles in responding to stress. Under ABA treatment, seven DoTCPs (DoTCP3, DoTCP7, DoTCP9, DoTCP11, DoTCP14, DoTCP15, and DoTCP21) were significantly upregulated, indicating that TCP genes may also play an important role in hormone response. Therefore, these results can provide useful information for studying the evolution and function of TCP genes in Dendrobium species.

Global Brassicaceae phylogeny based on filtering of 1,000-gene dataset

The mustard family (Brassicaceae) is a scientifically and economically important family, containing the model plant Arabidopsis thaliana and numerous crop species that feed billions worldwide. Despite its relevance, most phylogenetic trees of the family are incompletely sampled and often contain poorly supported branches. Here, we present the most complete Brassicaceae genus-level family phylogenies to date (Brassicaceae Tree of Life or BrassiToL) based on nuclear (1,081 genes, 319 of the 349 genera; 57 of the 58 tribes) and plastome (60 genes, 265 genera; all tribes) data. We found cytonuclear discordance between the two, which is likely a result of rampant hybridization among closely and more distantly related lineages. To evaluate the impact of such hybridization on the nuclear phylogeny reconstruction, we performed five different gene sampling routines, which increasingly removed putatively paralog genes. Our cleaned subset of 297 genes revealed high support for the tribes, whereas support for the main lineages (supertribes) was moderate. Calibration based on the 20 most clock-like nuclear genes suggests a late Eocene to late Oligocene origin of the family. Finally, our results strongly support a recently published new family classification, dividing the family into two subfamilies (one with five supertribes), together representing 58 tribes. This includes five recently described or re-established tribes, including Arabidopsideae, a monogeneric tribe accommodating Arabidopsis without any close relatives. With a worldwide community of thousands of researchers working on Brassicaceae and its diverse members, our new genus-level family phylogeny will be an indispensable tool for studies on biodiversity and plant biology.

Dysregulation of core neurodevelopmental pathways-a common feature of cancers with perineural invasion.

Background: High nerve density in tumors and metastasis via nerves (perineural invasion-PNI) have been reported extensively in solid tumors throughout the body including pancreatic, head and neck, gastric, prostate, breast, and colorectal cancers. Ablation of tumor nerves results in improved disease outcomes, suggesting that blocking nerve-tumor communication could be a novel treatment strategy. However, the molecular mechanisms underlying this remain poorly understood. Thus, the aim here was to identify molecular pathways underlying nerve-tumor crosstalk and to determine common molecular features between PNI-associated cancers. Results: Analysis of head and neck (HNSCC), pancreatic, and gastric (STAD) cancer Gene Expression Omnibus datasets was used to identify differentially expressed genes (DEGs). This revealed extracellular matrix components as highly dysregulated. To enrich for pathways associated with PNI, genes previously correlated with PNI in STAD and in 2 HNSCC studies where tumor samples were segregated by PNI status were analyzed. Neurodevelopmental genes were found to be enriched with PNI. In datasets where tumor samples were not segregated by PNI, neurodevelopmental pathways accounted for 12%-16% of the DEGs. Further dysregulation of axon guidance genes was common to all cancers analyzed. By examining paralog genes, a clear pattern emerged where at least one family member from several axon guidance pathways was affected in all cancers examined. Overall 17 different axon guidance gene families were disrupted, including the ephrin-Eph, semaphorin-neuropilin/plexin, and slit-robo pathways. These findings were validated using The Cancer Genome Atlas and cross-referenced to other cancers with a high incidence of PNI including colon, cholangiocarcinoma, prostate, and breast cancers. Survival analysis revealed that the expression levels of neurodevelopmental gene families impacted disease survival. Conclusion: These data highlight the importance of the tumor as a source of signals for neural tropism and neural plasticity as a common feature of cancer. The analysis supports the hypothesis that dysregulation of neurodevelopmental programs is a common feature associated with PNI. Furthermore, the data suggested that different cancers may have evolved to employ alternative genetic strategies to disrupt the same pathways. Overall, these findings provide potential druggable targets for novel therapies of cancer management and provide multi-cancer molecular biomarkers.

Comprehensive Analysis of the DNA Methyltransferase Genes and Their Association with Salt Response in Pyrus betulaefolia

DNA methylation plays an indispensable role in genome stability, regulation of gene expression and plant stress response. It is mediated by DNA methyltransferases (MTases). Twelve putative MTases of P. betulaefolia were identified and were classified into MET1, CMT, DRM2 and Dnmt2 groups based on the organization of various characteristic domains. Three pairs of paralogous genes were identified with the Ka/Ks ratio varied from 0.232 for PbeMET1a and PbeMET1b to 0.251 for PbeCMT2 and PbeCMT3, respectively. In addition, the Ka/Ks ratio for nine pairs of orthologous gene pairs between P. betulaefolia and apple were varied from 0.053 for PbeDRM3 and MD17G1031900 to 0.278 for PbeDnmt2b and MD15G1120500, respectively. All the 12 members of MTase genes were located on nine chromosomes out of 17 P. betulaefolia chromosomes with highly conserved protein sequence structures. The isoelectric points (pI) of MTases ranged from 4.74 to 7.24, while molecular weight varied from 35.99 to 174.32. The expression profiles of MTase and other salt-responsive genes under salt stress treatment revealed their important roles involved in salt response in P. betulaefolia. Furthermore, three selected salt-responsive genes (PbeNHX2.1, PbeCBL2 and PbeAKT2) were found altered in methylation level of promoters (which contain CpG islands) under salt stress. Especially, the PbeAKT2 promoter regions showed high CHG and CHH methylation types. Our study provided a genome-wide survey of the MTase gene family and highlighted their roles in salt response. These results also provided an effective way for the breeding and improvement of salt-tolerant pear varieties.

Polymorphism in FAE1gene paralogs of Brassica juncea associated with erucic acid content

Breeding for low erucic acid is a major objective to improve the edible oil quality in Brassica juncea. FAE1, theintron-less, 1521-bp gene encoding fatty acid elongase, genetically control the biosynthesis of erucic acid in oilseed Brassicas. The Brassica juncea being an amphidiploid, contains two paralogs of this gene (FAE1.1 and FAE1.2). Allelic variations in the gene such as SNPs and InDels are predominantly responsible for variability in erucic acid content in the seed, among the Brassica genotypes. The coding sequence of FAE1 gene paralogs were sequenced from low (LEA) and high (HEA) erucic acid genotypes which had a few SNPs that differentiated LAE and HEA genotypes. Two SNPs in FAE1.1 at position 591 and 1265 and one in FAE1.2 at 237, which led to a change in the recognition site of Hpy99I, BglII and MnlI restriction enzymes, respectively, converted into CAPS markers to differentiate LEA and HEA genotypes. The efficacy of these CAPS markers was found 100 per cent when validated in Brassica juncea, B. rapa and B. nigra genotypes and used in back-cross breeding. However, owing to functional complexity of CAPS in routine analysis; nucleotide sequence variability in the promoter region of FAE1 gene paralogs were expedited to develop simple PCR based markers. The upstream regions of FAE were sequenced from LEA cultivar Pusa Mustard30 and HEA cultivar Pusa Bold. A28-bp deletion in the promoter of FAE1.1 and a 340- bp insertion of a transposonlike element in the FAE1.2 gene promoter was discovered in LEA genotype in comparison to HEA. Markers based on these sequence variabilities in the promoter regions of FAE1.1 and FAE1.2 were found to completely co-segregate with the seed erucic acid content. These markers would be effectively used in marker assisted selection for development of low erucic acid cultivars in B. juncea.

Berberine bridge enzyme-like oxidases of cellodextrins and mixed-linked β-glucans control seed coat formation.

Plants have evolved various resistance mechanisms to cope with biotic stresses that threaten their survival. The BBE23 member (At5g44360/BBE23) of the Arabidopsis berberine bridge enzyme-like (BBE-l) protein family (Arabidopsis thaliana) has been characterized in this paper in parallel with the closely related and previously described CELLOX (At4g20860/BBE22). In addition to cellodextrins, both enzymes, renamed here as CELLODEXTRIN OXIDASE 2 and 1 (CELLOX2 and CELLOX1), respectively, oxidize the mixed-linked β-1→3/β-1→4-glucans (MLGs), recently described as capable of activating plant immunity, reinforcing the view that the BBE-l family includes members that are devoted to the control of the homeostasis of potential cell wall-derived damage-associated molecular patterns (DAMPs). The 2 putatively paralogous genes display different expression profiles. Unlike CELLOX1, CELLOX2 is not expressed in seedlings or adult plants and is not involved in immunity against Botrytis cinerea. Both are instead expressed in a concerted manner in the seed coat during development. Whereas CELLOX2 is expressed mainly during the heart stage, CELLOX1 is expressed at the immediately later stage, when the expression of CELLOX2 decreases. Analysis of seeds of cellox1 and cellox2 knockout mutants shows alterations in the coat structure: the columella area is smaller in cellox1, radial cell walls are thicker in both cellox1 and cellox2, and the mucilage halo is reduced in cellox2. However, the coat monosaccharide composition is not significantly altered, suggesting an alteration of the organization of the cell wall, thus reinforcing the notion that the architecture of the cell wall in specific organs is determined not only by the dynamics of the synthesis/degradation of the main polysaccharides but also by its enzymatic oxidation.

Genome-Wide Syntenic and Evolutionary Analysis of 30 Key Genes Found in Ten Oryza Species

Rice is a vital staple food crop worldwide, providing nutrition and sustenance to a significant portion of the global population. The genetic diversity of cultivated rice species has been significantly reduced during domestication, resulting in the loss of favorable alleles. To overcome this limitation, wild rice species have been used in introgression breeding programs to introduce beneficial alleles. In this study, we performed syntenic and phylogenetic analyses for 10 Oryza species, comprising both cultivar and wild species. Pairwise syntenic analysis revealed 3885 synteny blocks containing 1,023,342 gene pairs among 10 species. O. nivara contained the most blocks that were syntenous with the other nine species. In total, 425 paralogous and orthologous genes were identified for 30 key genes involved in rice breeding. His1 (43), GS3 (28), and qSW5/GW5 (27) had the most paralogous and orthologous genes. For GS3 and qSW5/GW5, two gene transfer events were detected. These findings have implications for rice breeding strategies, particularly with respect to gene pyramiding and introgression breeding programs. This research will contribute to the development of elite cultivars with improved quality and yield to meet the growing global demand for high-quality rice.

Analysis of the Arabidopsis venosa4‐0 mutant supports the role of VENOSA4 in dNTP metabolism

Human Sterile alpha motif and histidine-aspartate domain containing protein 1 (SAMHD1) functions as a dNTPase to maintain dNTP pool balance. In eukaryotes, the limiting step in de novo dNTP biosynthesis is catalyzed by RIBONUCLEOTIDE REDUCTASE (RNR). In Arabidopsis, the RNR1 subunit of RNR is encoded by CRINKLED LEAVES 8 (CLS8), and RNR2 by three paralogous genes, including TSO MEANING 'UGLY' IN CHINESE 2 (TSO2). In plants, DIFFERENTIAL DEVELOPMENT OF VASCULAR ASSOCIATED CELLS 1 (DOV1) catalyzes the first step of the de novo biosynthesis of purines. Here, to explore the role of VENOSA4 (VEN4), the most likely Arabidopsis ortholog of human SAMHD1, we studied the ven4‐0 point mutation, whose leaf phenotype was stronger than those of its insertional alleles. Structural predictions suggested that the E249L substitution in the mutated VEN4-0 protein rigidifies its 3D structure. The morphological phenotypes of the ven4, cls8, and dov1 single mutants were similar, and those of the ven4 tso2 and ven4 dov1 double mutants were synergistic. The ven4‐0 mutant had reduced levels of four amino acids related to dNTP biosynthesis, including glutamine and glycine, which are precursors in the de novo purine biosynthesis. Our results reveal high functional conservation between VEN4 and SAMHD1 in dNTP metabolism.