Paralogous Pairs Research Articles

Genome-Wide Identification and Expression Profiling of the PDI Gene Family Reveals Their Probable Involvement in Abiotic Stress Tolerance in Tomato (Solanum Lycopersicum L.).

Protein disulfide isomerases (PDI) and PDI-like proteins catalyze the formation and isomerization of protein disulfide bonds in the endoplasmic reticulum and prevent the buildup of misfolded proteins under abiotic stress conditions. In the present study, we conducted the first comprehensive genome-wide exploration of the PDI gene family in tomato (Solanum lycopersicum L.). We identified 19 tomato PDI genes that were unevenly distributed on 8 of the 12 tomato chromosomes, with segmental duplications detected for 3 paralogous gene pairs. Expression profiling of the PDI genes revealed that most of them were differentially expressed across different organs and developmental stages of the fruit. Furthermore, most of the PDI genes were highly induced by heat, salt, and abscisic acid (ABA) treatments, while relatively few of the genes were induced by cold and nutrient and water deficit (NWD) stresses. The predominant expression of SlPDI1-1, SlPDI1-3, SlPDI1-4, SlPDI2-1, SlPDI4-1, and SlPDI5-1 in response to abiotic stress and ABA treatment suggested they play regulatory roles in abiotic stress tolerance in tomato in an ABA-dependent manner. Our results provide new insight into the structure and function of PDI genes and will be helpful for the selection of candidate genes involved in fruit development and abiotic stress tolerance in tomato.

Genome-Wide Analysis of Heat Shock Transcription Factors in Ziziphus jujuba Identifies Potential Candidates for Crop Improvement Under Abiotic Stress.

Plant heat shock transcription factors (Hsfs) play a significant role in adoption under abiotic stress conditions by modulating the expression of several stress-responsive genes. Analysis of the Hsf gene family will serve to understand the molecular mechanism which is involved in response to abiotic stress. The Ziziphus species grows in warm and dry regions and is inherently tolerant to abiotic stress conditions; thus, Ziziphus is a highly enriched source of genes conferring abiotic stress tolerance. Therefore, the present study provides a comprehensive genome-wide analysis of the Hsf gene family in Z. jujuba. Identified 21 non-redundant Hsf genes were grouped into three major classes (classes A, B, and C) based on the phylogenetic analysis. Promoter and gene ontology analysis suggested that ZjHsfs perform diverse functions in response to abiotic stress conditions. Two paralogous pairs resulting from tandem gene duplication events were identified. Also, physio-chemical properties of chromosomal locations, gene structure, motifs, and protein domain organization of Hsfs were analyzed. Real-time PCR expression analyses revealed that most of the Z. jujuba Hsf genes are differentially expressed in response to heat stress. The analysis suggested ZjHsf-2, ZjHsf-3, ZjHsf-5, ZjHsf-7, ZjHsf-8, ZjHsf-10, ZjHsf-12, ZjHsf-17, and ZjHsf-18 were the outstanding candidate genes for imparting heat stress tolerance and for future functional analysis. The present analysis laid the foundation for understanding the molecular mechanism of the Hsf gene family regulating Z. jujuba development and tolerance to abiotic stress conditions.

The small heat shock proteins, chaperonin 10, in plants: An evolutionary view and emerging functional diversity

Small heat shock proteins (sHSPs) constitute a class of molecular chaperones, which are evolutionarily conserved yet diverse group of molecules, rapidly produced in response to stress. In this study, we sought to identify plant sHSPs, especially chaperonin 10 (Cpn10) family members in major evolutionary lineages, and determine their biological significance. Multiple sequence alignment of Cpn10 domains revealed divergent amino acids as well as conserved sites. Phylogenetic tree depicted the diversification and expansion of Cpn10 gene family. During the process of evolution, the Ka/Ks ratio of orthologous and paralogous pairs was <1, suggesting their evolutionary convergence and biological relevance. Functional annotations demonstrated that Cpn10 are involved in protein folding, regulation of metabolic processes and abiotic stress responses. Furthermore, subcellular localization prediction revealed that Cpn10 proteins are localized in multiple compartments, indicating a critical cell-coordinated defense. In-silico gene expression analysis exhibited their expression in most tissues examined, implying functional redundancy. Interactome analysis illustrated their interaction with chloroplast and mitochondrial genes, which are majorly involved in protein folding and assembly. The transcriptional regulation revealed their stress-responsive and distinct physiological roles. Our findings would contribute to new insights on the evolutionary history of Cpn10 gene family and the distinct biological roles.

Genome-wide survey of soybean papain-like cysteine proteases and their expression analysis in root nodule symbiosis

BackgroundPlant papain-like cysteine proteases (PLCPs) are a large class of proteolytic enzymes and play important roles in root nodule symbiosis (RNS), while the whole-genome studies of PLCP family genes in legume are quite limited, and the roles of Glycine max PLCPs (GmPLCPs) in nodulation, nodule development and senescence are not fully understood.ResultsIn the present study, we identified 97 GmPLCPs and performed a genome-wide survey to explore the expansion of soybean PLCP family genes and their relationships to RNS. Nineteen paralogous pairs of genomic segments, consisting of 77 GmPLCPs, formed by whole-genome duplication (WGD) events were identified, showing a high degree of complexity in duplication. Phylogenetic analysis among different species showed that the lineage differentiation of GmPLCPs occurred after family expansion, and large tandem repeat segment were specifically in soybean. The expression patterns of GmPLCPs in symbiosis-related tissues and nodules identified RNS-related GmPLCPs and provided insights into their putative symbiotic functions in soybean. The symbiotic function analyses showed that a RNS-related GmPLCP gene (Glyma.04G190700) really participate in nodulation and nodule development.ConclusionsOur findings improved our understanding of the functional diversity of legume PLCP family genes, and provided insights into the putative roles of the legume PLCPs in nodulation, nodule development and senescence.

Expression Partitioning of Duplicate Genes at Single Cell Resolution in Arabidopsis Roots.

Gene duplication is a key evolutionary phenomenon, prevalent in all organisms but particularly so in plants, where whole genome duplication (WGD; polyploidy) is a major force in genome evolution. Much effort has been expended in attempting to understand the evolution of duplicate genes, addressing such questions as why some paralog pairs rapidly return to single copy status whereas, in other pairs, both paralogs are retained and may diverge in expression pattern or function. The effect of a gene – its site of expression and thus the initial locus of its function – occurs at the level of a cell comprising a single cell type at a given state of the cell’s development. Using Arabidopsis thaliana single cell transcriptomic data we categorized patterns of expression for 11,470 duplicate gene pairs across 36 cell clusters comprising nine cell types and their developmental states. Among these 11,470 pairs, 10,187 (88.8%) had at least one copy expressed in at least one of the 36 cell clusters. Pairs produced by WGD more often had both paralogs expressed in root cells than did pairs produced by small scale duplications. Three quarters of gene pairs expressed in the 36 cell clusters (7,608/10,187) showed extreme expression bias in at least one cluster, including 352 cases of reciprocal bias, a pattern consistent with expression subfunctionalization. More than twice as many pairs showed reciprocal expression bias between cell states than between cell types or between roots and leaves. A group of 33 gene pairs with reciprocal expression bias showed evidence of concerted divergence of gene networks in stele vs. epidermis. Pairs with both paralogs expressed without bias were less likely to have paralogs with divergent mutant phenotypes; such bias-free pairs showed evidence of preservation by maintenance of dosage balance. Overall, we found considerable evidence of shifts in gene expression following duplication, including in >80% of pairs encoding 7,653 genes expressed ubiquitously in all root cell types and states for which we inferred the polarity of change.

ZmSKS13, a cupredoxin domain-containing protein, is required for maize kernel development via modulation of redox homeostasis.

The SKU5 similar (SKS) genes encode a family of multi-copper-oxidase-like proteins with cupredoxin domains similar to those in laccase and ascorbate oxidase. Although SKS proteins are known to function in root growth and cotyledon vascular patterning in Arabidopsis, their role in plant reproductive processes is poorly understood. Here, we identified a seed mutant of maize (Zea mays), generated by ethyl methane sulfonate (EMS) mutagenesis, that we designated defective kernel-zk1 (dek-zk1). The mutant produced small, shriveled kernels with an aberrant basal endosperm transfer layer (BETL) and placento-chalazal (PC) layer and irregular starch granules. Map-based cloning revealed that Dek-zk1 encodes an SKU5 similar 13 (GenBank: ONM36900.1), so it was named ZmSKS13. ZmSKS13 comprises a paralogous pair with Zm00001d012524, but the transcript abundance of ZmSKS13 in developing kernels is 15 times higher than that of Zm00001d012524, resulting in dek-zk1 mutation conveying a distinct kernel phenotype. ZmSKS13 loss of function led to overaccumulation of reactive oxygen species (ROS) and severe DNA damage in the nucellus and BETL and PC layer cells, and exogenous antioxidants significantly alleviated the defects of the mutant kernels. Our results thus demonstrate that ZmSKS13 is a novel regulator that plays a crucial role in kernel development in maize through the modulation of ROS homeostasis.

Genome-wide analysis of HSP90 gene family in the Mediterranean olive (Olea europaea subsp. europaea) provides insight into structural patterns, evolution and functional diversity.

Plants regularly experience multiple abiotic and biotic pressures affecting their normal development. The 90-kDa heat shock protein (HSP90) plays a dynamic role in countering abiotic and biotic stresses via a plethora of functional mechanisms. The HSP90 has been investigated in many plant species. However, there is little information available about this gene family in the cultivated Mediterranean olive tree, Olea europaea subsp. europaea var. europaea. In thecurrentstudy, we systematicallyperformed genome-wide identification and characterization of the HSP90 gene family in O. europaea var. europaea (OeHSP90s). Twelve regular OeHSP90s were identified, which were phylogenetically grouped into two major clusters and four sub-clusters, showing five paralogous gene pairs evolving under purifying selection. All of the 12 proteins contained a Histidine kinase-like ATPase (HATPase_c) domain, justifying the role played by HSP90 proteins in ATP binding and hydrolysis. The predicted 3D structure of OeHSP90 proteins provided information to understand their functions at the biochemical level. Consistent with their phylogenetic relationships, OeHSP90 members were predicted to be localized in different cellular compartments, suggesting their involvement in various subcellular processes. In consonance with their spatial organization, olive HSP90 family members were found to share similar motif arrangements and similar number of exons. We found that OeHSP90 promoters contained various cis-acting elements associated with light responsiveness, hormone signaling pathways and reaction to various stress conditions. In addition, expression sequence tags (ESTs) analysis offered a view of OeHSP90 tissue- and developmental stage specific pattern of expression. Proteins interacting with OeHSP90s were predicted and their potential roles were discussed. Overall, our results offer premises for further investigation of the implication of HSP90 genes in the physiological processes of the olive and its adaptation to stresses.

Multiplex enCas12a screens detect functional buffering among paralogs otherwise masked in monogenic Cas9 knockout screens

BackgroundPooled library CRISPR/Cas9 knockout screening across hundreds of cell lines has identified genes whose disruption leads to fitness defects, a critical step in identifying candidate cancer targets. However, the number of essential genes detected from these monogenic knockout screens is low compared to the number of constitutively expressed genes in a cell.ResultsThrough a systematic analysis of screen data in cancer cell lines generated by the Cancer Dependency Map, we observe that half of all constitutively expressed genes are never detected in any CRISPR screen and that these never-essentials are highly enriched for paralogs. We investigated functional buffering among approximately 400 candidate paralog pairs using CRISPR/enCas12a dual-gene knockout screening in three cell lines. We observe 24 synthetic lethal paralog pairs that have escaped detection by monogenic knockout screens at stringent thresholds. Nineteen of 24 (79%) synthetic lethal interactions are present in at least two out of three cell lines and 14 of 24 (58%) are present in all three cell lines tested, including alternate subunits of stable protein complexes as well as functionally redundant enzymes.ConclusionsTogether, these observations strongly suggest that functionally redundant paralogs represent a targetable set of genetic dependencies that are systematically under-represented among cell-essential genes in monogenic CRISPR-based loss of function screens.

Expansion of the circadian transcriptome in Brassica rapa and genome-wide diversification of paralog expression patterns.

An important challenge of crop improvement strategies is assigning function to paralogs in polyploid crops. Here we describe the circadian transcriptome in the polyploid crop Brassica rapa. Strikingly, almost three-quarters of the expressed genes exhibited circadian rhythmicity. Genetic redundancy resulting from whole genome duplication is thought to facilitate evolutionary change through sub- and neo-functionalization among paralogous gene pairs. We observed genome-wide expansion of the circadian expression phase among retained paralogous pairs. Using gene regulatory network models, we compared transcription factor targets between B. rapa and Arabidopsis circadian networks to reveal evidence for divergence between B. rapa paralogs that may be driven in part by variation in conserved non-coding sequences (CNS). Additionally, differential drought response among retained paralogous pairs suggests further functional diversification. These findings support the rapid expansion and divergence of the transcriptional network in a polyploid crop and offer a new approach for assessing paralog activity at the transcript level.

BcMF30a and BcMF30c, Two Novel Non-Tandem CCCH Zinc-Finger Proteins, Function in Pollen Development and Pollen Germination in Brassica campestris ssp. chinensis

Chinese cabbage (Brassica campestris) is an economically important leaf vegetable crop worldwide. Mounting studies have shown that cysteine-cysteine-cysteine-histidine (CCCH) zinc-finger protein genes are involved in various plant growth and development processes. However, research on the involvement of these genes in male reproductive development is still in its infancy. Here, we identified 11 male fertility-related CCCH genes in Chinese cabbage. Among them, a pair of paralogs encoding novel non-tandem CCCH zinc-finger proteins, Brassica campestris Male Fertility 30a (BcMF30a) and BcMF30c, were further characterized. They were highly expressed in pollen during microgametogenesis and continued to express in germinated pollen. Further analyses demonstrated that both BcMF30a and BcMF30c may play a dual role as transcription factors and RNA-binding proteins in plant cells. Functional analysis showed that partial bcmf30a bcmf30c pollen grains were aborted due to the degradation of pollen inclusion at the microgametogenesis phase, and the germination rate of viable pollen was also greatly reduced, indicating that BcMF30a and BcMF30c are required for both pollen development and pollen germination. This research provided insights into the function of CCCH proteins in regulating male reproductive development and laid a theoretical basis for hybrid breeding of Chinese cabbage.

Determining the interaction status and evolutionary fate of duplicated homomeric proteins

Oligomeric proteins are central to life. Duplication and divergence of their genes is a key evolutionary driver, also because duplications can yield very different outcomes. Given a homomeric ancestor, duplication can yield two paralogs that form two distinct homomeric complexes, or a heteromeric complex comprising both paralogs. Alternatively, one paralog remains a homomer while the other acquires a new partner. However, so far, conflicting trends have been noted with respect to which fate dominates, primarily because different methods and criteria are being used to assign the interaction status of paralogs. Here, we systematically analyzed all Saccharomyces cerevisiae and Escherichia coli oligomeric complexes that include paralogous proteins. We found that the proportions of homo-hetero duplication fates strongly depend on a variety of factors, yet that nonetheless, rigorous filtering gives a consistent picture. In E. coli about 50%, of the paralogous pairs appear to have retained the ancestral homomeric interaction, whereas in S. cerevisiae only ~10% retained a homomeric state. This difference was also observed when unique complexes were counted instead of paralogous gene pairs. We further show that this difference is accounted for by multiple cases of heteromeric yeast complexes that share common ancestry with homomeric bacterial complexes. Our analysis settles contradicting trends and conflicting previous analyses, and provides a systematic and rigorous pipeline for delineating the fate of duplicated oligomers in any organism for which protein-protein interaction data are available.

Patterns of Expansion and Expression Divergence of the Polygalacturonase Gene Family in Brassica oleracea.

Plant polygalacturonases (PGs) are closely related to cell-separation events during plant growth and development by degrading pectin. Identifying and investigating their diversification of evolution and expression could shed light on research on their function. We conducted sequence, molecular evolution, and gene expression analyses of PG genes in Brassica oleracea. Ninety-nine B. oleracea PGs (BoPGs) were identified and divided into seven clades through phylogenetic analysis. The exon/intron structures and motifs were conserved within, but divergent between, clades. The second conserved domain (GDDC) may be more closely related to the identification of PGs. There were at least 79 common ancestor PGs between Arabidopsis thaliana and B. oleracea. The event of whole genome triplication and tandem duplication played important roles in the rapid expansion of the BoPG gene family, and gene loss may be an important mechanism in the generation of the diversity of BoPGs. By evaluating the expression in five tissues, we found that most of the expressed BoPGs in clades A, B, and E showed ubiquitous expression characteristics, and the expressed BoPGs in clades C, D, and F were mainly responsible for reproduction development. Most of the paralogous gene pairs (76.2%) exhibited divergent expression patterns, indicating that they may have experienced neofunctionalization or subfunctionalization. The cis-elements analysis showed that up to 96 BoPGs contained the hormone response elements in their promoters. In conclusion, our comparative analysis may provide a valuable data foundation for the further functional analysis of BoPGs during the development of B. oleracea.

Recurrent sequence evolution after independent gene duplication

BackgroundConvergent and parallel evolution provide unique insights into the mechanisms of natural selection. Some of the most striking convergent and parallel (collectively recurrent) amino acid substitutions in proteins are adaptive, but there are also many that are selectively neutral. Accordingly, genome-wide assessment has shown that recurrent sequence evolution in orthologs is chiefly explained by nearly neutral evolution. For paralogs, more frequent functional change is expected because additional copies are generally not retained if they do not acquire their own niche. Yet, it is unknown to what extent recurrent sequence differentiation is discernible after independent gene duplications in different eukaryotic taxa.ResultsWe develop a framework that detects patterns of recurrent sequence evolution in duplicated genes. This is used to analyze the genomes of 90 diverse eukaryotes. We find a remarkable number of families with a potentially predictable functional differentiation following gene duplication. In some protein families, more than ten independent duplications show a similar sequence-level differentiation between paralogs. Based on further analysis, the sequence divergence is found to be generally asymmetric. Moreover, about 6% of the recurrent sequence evolution between paralog pairs can be attributed to recurrent differentiation of subcellular localization. Finally, we reveal the specific recurrent patterns for the gene families Hint1/Hint2, Sco1/Sco2 and vma11/vma3.ConclusionsThe presented methodology provides a means to study the biochemical underpinning of functional differentiation between paralogs. For instance, two abundantly repeated substitutions are identified between independently derived Sco1 and Sco2 paralogs. Such identified substitutions allow direct experimental testing of the biological role of these residues for the repeated functional differentiation. We also uncover a diverse set of families with recurrent sequence evolution and reveal trends in the functional and evolutionary trajectories of this hitherto understudied phenomenon.

In Kluyveromyces lactis a Pair of Paralogous Isozymes Catalyze the First Committed Step of Leucine Biosynthesis in Either the Mitochondria or the Cytosol.

Divergence of paralogous pairs, resulting from gene duplication, plays an important role in the evolution of specialized or novel gene functions. Analysis of selected duplicated pairs has elucidated some of the mechanisms underlying the functional diversification of Saccharomyces cerevisiae (S. cerevisiae) paralogous genes. Similar studies of the orthologous pairs extant in pre-whole genome duplication yeast species, such as Kluyveromyces lactis (K. lactis) remain to be addressed. The genome of K. lactis, an aerobic yeast, includes gene pairs generated by sporadic duplications. The genome of this organism comprises the KlLEU4 and KlLEU4BIS paralogous pair, annotated as putative α-isopropylmalate synthases (α-IPMSs), considered to be the orthologs of the S. cerevisiae ScLEU4/ScLEU9 paralogous genes. The enzymes encoded by the latter two genes are mitochondrially located, differing in their sensitivity to leucine allosteric inhibition resulting in ScLeu4-ScLeu4 and ScLeu4-ScLeu9 sensitive dimers and ScLeu9-ScLeu9 relatively resistant homodimers. Previous work has shown that, in a Scleu4Δ mutant, ScLEU9 expression is increased and assembly of ScLeu9-ScLeu9 leucine resistant homodimers results in loss of feedback regulation of leucine biosynthesis, leading to leucine accumulation and decreased growth rate. Here we report that: (i) K. lactis harbors a sporadic gene duplication, comprising the KlLEU4, syntenic with S. cerevisiae ScLEU4 and ScLEU9, and the non-syntenic KlLEU4BIS, arising from a pre-WGD event. (ii) That both, KlLEU4 and KlLEU4BIS encode leucine sensitive α-IPMSs isozymes, located in the mitochondria (KlLeu4) and the cytosol (KlLeu4BIS), respectively. (iii) That both, KlLEU4 or KlLEU4BIS complement the Scleu4Δ Scleu9Δ leucine auxotrophic phenotype and revert the enhanced ScLEU9 transcription observed in a Scleu4Δ ScLEU9 mutant. The Scleu4Δ ScLEU9 growth mutant phenotype is only fully complemented when transformed with the syntenic KlLEU4 mitochondrial isoform. KlLEU4 and KlLEU4BIS underwent a different diversification pathways than that leading to ScLEU4/ScLEU9. KlLEU4 could be considered as the functional ortholog of ScLEU4, since its encoded isozyme can complement both the Scleu4Δ Scleu9Δ leucine auxotrophy and the Scleu4Δ ScLEU9 complex phenotype.

YkgM and YkgO maintain translation by replacing their paralogs, zinc-binding ribosomal proteins L31 and L36, with identical activities.

When a cell is zinc-deficient, ykgM and ykgO, which encode paralogs of the zinc-binding ribosomal proteins L31 and L36, are expressed from the ykgM operon, which is ordinarily held inactive by the Zur repressor. In ribosomes lacking L31, ribosomal subunit association is weakened, resulting in reduced in vitro translation and the deletion mutants of rpmE, the gene encoding L31, forming small colonies. We isolated four suppressor mutants of ∆rpmE that formed normal colonies. All four mutation sites were located in zur, and ribosomes of zur mutant cells contained one copy of YkgM and had translational activities equivalent to those of ribosomes containing L31. L36 is highly conserved among bacteria, chloroplast and mitochondria. Analysis of a deletion mutant of rpmJ, which encodes L36, suggested that L36 is involved in late assembly of the 50S particle, in vitro translation and cell growth. In zur mutant cells lacking rpmJ, the paralog YkgO was expressed and took over the functions of L36. zur mutant cells contained four types of ribosomes containing combinations of L31 or YkgM, and L36 or YkgO. Copy numbers of L31 and YkgM, and L36 and YkgO, summed to 1, indicating that each paralog pair shares a binding site.

Genome-wide identification of plasma membrane aquaporin gene family in Populus and functional identification of PIP1;1 involved in osmotic stress

Aquaporins are the membrane channels that facilitate the transportation of water and small neutral molecules in plants. As such, they play an important role in stress tolerance. This function may be more pronounced in woody plants, which must efficiently absorb water and transport it over long distances. Therefore, a functional analysis of plasma membrane aquaporins (PIP) under drought stress can provide insights into the roles of these proteins in woody plants. In this work, we obtained 15 members of the PtPIP gene family in poplar. They formed two groups in a phylogenetic analysis. Gene structure, including intron phases, was found to be remarkably well conserved among members of this gene family. Five paralogous gene pairs were identified in the PtPIP family, and four of those pairs arose via whole genome duplication. A selection pressure analysis showed that PtPIP genes have undergone positive selection. The transcripts of six PtPIP genes including PtPIP1;1 were predominantly detected in the root. Transgenic poplar overexpressing PtoPIP1;1 was sensitive to osmotic stress, which phenocopied PagMYB216-overexpressing poplar plants. Transactivation and yeast one-hybrid assays showed that PagMYB216, whose homologs in other plants are induced by abiotic stresses, is able to activate PtoPIP1;1 expression through binding to the MBS motif in its promoter. Our results provide new information about PtoPIPs and their roles in response to abiotic stresses, which will be helpful for better understanding their functions in Populus.

Molecular characterization of ascorbate peroxidase (APX) and APX-related (APX-R) genes in Triticum aestivum L.

Ascorbate peroxidases (APXs) are heme-dependent H2O2 scavenging enzymes involved in myriad biological processes. Herein, a total of 21 TaAPX and six TaAPX-R genes were identified from the A, B and D sub-genomes of Triticum aestivum L. The occurrence of three paralogous gene pairs with unequal evolutionary rate suggested functional divergence. The phylogenetic analysis formed four distinct clades having conserved gene and protein architecture, and sub-cellular localization. The tertiary structure analysis revealed the presence of helices and coils and residues involved in ligand binding. Transcriptional profiling of each TaAPX and TaAPX-R gene suggested their specific role during development and stress response. Modulated transcript expression and APX enzyme activity during various stress conditions indicated their role in stress response. Interaction analyses suggested their association with other genes, miRNAs and various legends. The present study reported numerous features of these genes, and may provide a platform for their detailed functional characterization in future studies.

A Chinese White Pear (Pyrus bretschneideri) BZR Gene PbBZR1 Act as a Transcriptional Repressor of Lignin Biosynthetic Genes in Fruits.

BZR transcription factors play essential roles in plant growth and environmental stimuli, and they are also the positive regulators of Brassinosteroid (BR) signal transduction in diverse plants. In addition, BZR TFs, as crucial regulators of BR synthesis, may have multiple stress-resistance functions and their related regulatory mechanisms have been well illustrated in model plants. Here, we carried out a genome-wide identification of BZR members in Chinese pear (Pyrus bretschneideri) and identified 13 members. By comparative analysis in five Rosaceae genomes, BZR members in the pear genome may have undergone large-scale duplication events during evolution. Purifying selection played an important role in almost all of the orthologous and paralogous gene pairs. According to the expression analysis of the PbBZRs during fruit development, three PbBZRs were selected for detailed analysis. Transcriptional activation assays presented that PbBZR1 repressed the promoters of P. bretschneideri lignin biosynthetic genes, such as PbCES9, PbCOMT3, and PbHCT6. Our study traces the evolution of BZR gene family members in Rosaceae genomes and illustrates that the rates of gene loss and gain are far from equilibrium in different species. At the same time, our results suggest that PbBZR1 may be involved in the negative regulation of lignin biosynthesis.

The differential expression patterns and co-expression networks of paralogs as an indicator of the TNM stages of lung adenocarcinoma and squamous cell carcinoma

Cancers constitute a severe threat to human health. Elucidating the association between the expression patterns of the paralogous genes and transcription factors (TF) and the progression of cancers by comprehensively investigating the expression patterns and co-expression networks will contribute to the in-depth understanding of the pathogenesis of cancers. Here, we identified the paralogous gene pairs and systematically analyzed the expression patterns of these paralogs and the known TFs to elucidate the associations with Tumor, Node, Metastasis (TNM) staging information across ten cancers. We found that the expression of ~60% paralogs was cancer-dependent, and more than 50% of the differentially expressed TFs pairs showed positive expression correlations. The down-regulation patterns of paralogs and TFs were closely associated with the M and N developmental stages of lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC). Our results will help to understand the roles of paralogs and TFs in cancer progression and to screen prognostic biomarkers for early cancer diagnosis.

Electrostatic Frustration Shapes Folding Mechanistic Differences in Paralogous Bacterial Stress Response Proteins

Paralogous proteins play a vital role in evolutionary adaptation of organisms and species divergence. One outstanding question is the molecular basis for how folding mechanisms differ in paralogs that not only exhibit similar topologies but also evolve under near-identical selection pressures. Here, we address this question by studying a paralogous protein pair from enterobacteria, Hha and Cnu, combining experiments, simulations and statistical modeling. We find that Hha is less stable and folds an order of magnitude slower than Cnu despite similar packing and topological features. Differences in surface charge–charge interactions, however, promote a N-terminal biased unfolding mechanism in Hha unlike Cnu that unfolds via the C terminus. Our work highlights how electrostatic frustration contributes to the population of heterogeneous native ensembles in paralogs and the avenues through which evolutionary topological constraints could be overcome by modulating charge–charge interactions.