Evidence For Functional Divergence Research Articles

Tandemly duplicated MYB genes are functionally diverged in the regulation of anthocyanin biosynthesis in soybean.

Gene duplications have long been recognized as a driving force in the evolution of genes, giving rise to novel functions. The soybean (Glycine max) genome is characterized by a large number of duplicated genes. However, the extent and mechanisms of functional divergence among these duplicated genes in soybean remain poorly understood. In this study, we revealed that 4 MYB genes (GmMYBA5, GmMYBA2, GmMYBA1, and Glyma.09g235000)-presumably generated by tandem duplication specifically in the Phaseoleae lineage-exhibited a stronger purifying selection in soybean compared to common bean (Phaseolus vulgaris). To gain insights into the diverse functions of these tandemly duplicated MYB genes in anthocyanin biosynthesis, we examined the expression, transcriptional activity, induced metabolites, and evolutionary history of these 4 MYB genes. Our data revealed that Glyma.09g235000 is a pseudogene, while the remaining 3 MYB genes exhibit strong transcriptional activation activity, promoting anthocyanin biosynthesis in different soybean tissues. GmMYBA5, GmMYBA2, and GmMYBA1 induced anthocyanin accumulation by upregulating the expression of anthocyanin pathway-related genes. Notably, GmMYBA5 showed a lower capacity for gene induction compared to GmMYBA2 and GmMYBA1. Metabolomics analysis further demonstrated that GmMYBA5 induced distinct anthocyanin accumulation in Nicotiana benthamiana leaves and soybean hairy roots compared to GmMYBA2 and GmMYBA1, suggesting their functional divergence leading to the accumulation of different metabolites accumulation following gene duplication. Together, our data provide evidence of functional divergence within the MYB gene cluster following tandem duplication, which sheds light on the potential evolutionary directions of gene duplications during legume evolution.

Nuku, a family of primate retrocopies derived from KU70.

The gene encoding the ubiquitous DNA repair protein, Ku70p, has undergone extensive copy number expansion during primate evolution. Gene duplications of KU70 have the hallmark of long interspersed element-1 mediated retrotransposition with evidence of target-site duplications, the poly-A tails, and the absence of introns. Evolutionary analysis of this expanded family of KU70-derived “NUKU” retrocopies reveals that these genes are both ancient and also actively being created in extant primate species. NUKU retrocopies show evidence of functional divergence away from KU70, as evinced by their altered pattern of tissue expression and possible tissue-specific translation. Molecular modeling predicted that amino acid changes in Nuku2p at the interaction interface with Ku80p would prevent the assembly of the Ku heterodimer. The lack of Nuku2p-Ku80p interaction was confirmed by yeast two-hybrid assay, which contrasts the robust interaction of Ku70p-Ku80p. While several NUKU retrocopies appear to have been degraded by mutation, NUKU2 shows evidence of positive natural selection, suggesting that this retrocopy is undergoing neofunctionalization. Although Nuku proteins do not appear to antagonize retrovirus transduction in cell culture, the observed expansion and rapid evolution of NUKUs could be being driven by alternative selective pressures related to infectious disease or an undefined role in primate physiology.

Molecular evolution of Pr1 proteases depicts ongoing diversification in Metarhizium spp.

The Pr1 family of serine endopeptidases plays an important role in pathogenicity and virulence of entomopathogens such as Metarhizium anisopliae (Ascomycota: Hypocreales). These virulence factors allow for the penetration of the host cuticle, a vital step in the infective process of this fungus, which possesses 11 Pr1 isoforms (Pr1A through Pr1K). The family is divided into two classes with Class II (proteinase K-like) comprising 10 isoforms further split into three subfamilies. It is believed that these isoforms act synergistically and with other virulence factors, allowing pathogenicity to multiple hosts. As virulence coevolves through reciprocal selection with hosts, positive selection may lead to the evolution of new protease families or isoforms of extant ones that can withstand host defenses. This work tests this hypothesis in Class II Pr1 proteins, focusing on M. anisopliae, employing different methods for phylogenetic inference in amino acid and nucleotide datasets in multiple arrangements for Metarhizium spp. and related species. Phylogenies depict groups that match the taxonomy of their respective organisms with high statistical support, with minor discrepancies. Positively selected sites were identified in six out of ten Pr1 isoforms, most of them located in the proteolytic domain and spatially close to the catalytic residues. Moreover, there was evidence of functional divergence in the majority of pairwise comparisons. These results imply the existence of differential selective pressure acting on Pr1 proteins and a potential new isoform, likely affecting host specificities, virulence, or even adapting the organism to different host-independent lifestyles.

Structural organization and functional divergence of high isoelectric point \u03b1-amylase genes in bread wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.)

BackgroundHigh isoelectric point α-amylase genes (Amy1) play major roles during cereal seed germination, and are associated with unacceptable high residual α-amylase activities in ripe wheat grains. However, in wheat and barley, due to extremely high homology of duplicated copies, and large and complex genome background, the knowledge on this multigene family is limited.ResultsIn the present work, we identified a total of 41 Amy1 genes among 13 investigated grasses. By using genomic resources and experimental validation, the exact copy numbers and chromosomal locations in wheat and barley were determined. Phylogenetic and syntenic analyses revealed tandem gene duplication and chromosomal rearrangement leading to separation of Amy1 into two distinct loci, Amy1θ and Amy1λ. The divergence of Amy1λ from Amy1θ was driven by adaptive selection pressures performed on two amino acids, Arg97 and Asn233 (P > 0.95*). The predicted protein structural alteration caused by substitution of Asp233Asn in the conserved starch binding surface site, and significantly expressional differentiation during seed germination and grain development provided evidence of functional divergence between Amy1θ and Amy1λ genes. We screened out candidate copies (TaAmy1-A1/A2 and TaAmy1-D1) associated with high residual α-amylase activities in ripe grains. Furthermore, we proposed an evolutionary model for expansion dynamics of Amy1 genes.ConclusionsOur study provides comprehensive analyses of the Amy1 multigene family, and defines the fixation of two spatially structural Amy1 loci in wheat and barley. Potential functional divergence between them is reflected by their sequence features and expressional patterns, and driven by gene duplication, chromosome rearrangement and natural selections during gene family evolution. Furthermore, the discrimination of differentially effective copies during seed germination and/or grain development will provide guidance to manipulation of α-amylase activity in wheat and barley breeding for better yield and processing properties.

Distinct Domains of CheA Confer Unique Functions in Chemotaxis and Cell Length in Azospirillum brasilense Sp7.

Chemotaxis is the movement of cells in response to gradients of diverse chemical cues. Motile bacteria utilize a conserved chemotaxis signal transduction system to bias their motility and navigate through a gradient. A central regulator of chemotaxis is the histidine kinase CheA. This cytoplasmic protein interacts with membrane-bound receptors, which assemble into large polar arrays, to propagate the signal. In the alphaproteobacterium Azospirillum brasilense, Che1 controls transient increases in swimming speed during chemotaxis, but it also biases the cell length at division. However, the exact underlying molecular mechanisms for Che1-dependent control of multiple cellular behaviors are not known. Here, we identify specific domains of the CheA1 histidine kinase implicated in modulating each of these functions. We show that CheA1 is produced in two isoforms: a membrane-anchored isoform produced as a fusion with a conserved seven-transmembrane domain of unknown function (TMX) at the N terminus and a soluble isoform similar to prototypical CheA. Site-directed and deletion mutagenesis combined with behavioral assays confirm the role of CheA1 in chemotaxis and implicate the TMX domain in mediating changes in cell length. Fluorescence microscopy further reveals that the membrane-anchored isoform is distributed around the cell surface while the soluble isoform localizes at the cell poles. Together, the data provide a mechanism for the role of Che1 in controlling multiple unrelated cellular behaviors via acquisition of a new domain in CheA1 and production of distinct functional isoforms.IMPORTANCE Chemotaxis provides a significant competitive advantage to bacteria in the environment, and this function has been transferred laterally multiple times, with evidence of functional divergence in different genomic contexts. The molecular principles that underlie functional diversification of chemotaxis in various genomic contexts are unknown. Here, we provide a molecular mechanism by which a single CheA protein controls two unrelated functions: chemotaxis and cell length. Acquisition of this multifunctionality is seemingly a recent evolutionary event. The findings illustrate a mechanism by which chemotaxis function may be co-opted to regulate additional cellular functions.

Evidence for functional divergence in AM fungal communities from different montane altitudes

Abstract A considerable amount of symbiotic functional variability has been demonstrated among arbuscular mycorrhizal fungal (AMF) species. However, little is known about the functional divergence in AMF communities from contrasting climates. Ecotypes of Pennisetum flaccidum were grown in growth chambers with AMF communities originating from two elevations (3 105 and 4 176 m asl) on Mount Segrila in Tibet. The combinations were grown under two temperature regimes (7–12 °C and 14–19 °C), equivalent to the mean temperatures during the growing season at the two elevations. AMF communities from different elevations showed different fungal performance at both temperatures. In most cases, AMF from cold habitats developed higher intraradical and extraradical colonization at 7–12 °C than AMF from the warmer habitat, and vice versa. For mycorrhizal growth responses (MGR), there was no consistent local versus foreign effect, but AMF inoculum from colder habitats consistently led to a higher MGR than the other inoculum. These results provide strong evidence for functional differences among AMF communities from contrasting climates.

Differential selection within the Drosophila retinal determination network and evidence for functional divergence between paralog pairs

The retinal determination (RD) network in Drosophila comprises 14 known nuclear proteins that include DNA-binding proteins, transcriptional coactivators, kinases, and phosphatases. The composition of the network varies considerably throughout the animal kingdom, with the network in several basal insects having fewer members and with vertebrates having potentially significantly higher numbers of RD genes. One important contributing factor for the variation in gene number within the network is gene duplication. For example, 10 members of the RD network in Drosophila are derived from duplication events. Here we present an analysis of the coding regions of the five pairs of duplicate genes from within the RD network of several different Drosophila species. We demonstrate that there is differential selection across the coding regions of all RD genes. Additionally, some of the most significant differences in ratios of non-silent-to-silent site substitutions (d(N)/d(S)) between paralog pairs are found within regions that have no ascribed function. Previous structure/function analyses of several duplicate genes have identified areas within one gene that contain novel activities when compared with its paralog. The evolutionary analysis presented here identifies these same areas in the paralogs as being under high levels of relaxed selection. We suggest that sequence divergence between paralogs and selection signatures can be used as a reasonable predictor of functional changes in rapidly evolving motifs.

A comparative genomic analysis of the small heat shock proteins in Caenorhabditis elegans and briggsae

The small heat shock proteins (sHSPs) are a ubiquitous family of molecular chaperones. We have identified 18 sHSPs in the Caenorhabditis elegans genome and 20 sHSPs in the Caenorhabditis briggsae genome. Analysis of phylogenetic relationships and evolutionary dynamics of the sHSPs in these two genomes reveals a very complex pattern of evolution. The sHSPs in C. elegans and C. briggsae do not display clear orthologous relationships with other invertebrate sHSPs. But many sHSPs in C. elegans have orthologs in C. briggsae. One group of sHSPs, the HSP16s, has a very unusual evolutionary history. Although there are a number of HSP16s in both the C. elegans and C. briggsae genomes, none of the HSP16s display orthologous relationships across these two species. The HSP16s have an unusual gene pair structure and a complex evolutionary history shaped by gene duplication, gene conversion, and purifying selection. We found no evidence of recent positive selection acting on any of the sHSPs in C. elegans or in C. briggsae. There is also no evidence of functional divergence within the pairs of orthologous C. elegans and C. briggsae sHSPs. However, the evolutionary patterns do suggest that functional divergence has occurred between the sHSPs in C. elegans and C. briggsae and the sHSPs in more distantly related invertebrates.

Full-length sequence and expression analysis of Toll-like receptor 9 in the gilthead seabream ( Sparus aurata L.)

Toll-Like Receptors (TLRs) have recently emerged as key sensors of invading microbes, acting through recognition of pathogen-associated molecular patterns. It has been demonstrated that TLR9 is involved in the recognition of unmethylated CpG motifs in mice, humans, and pigs. We report here the full-length sequence of TLR9 cDNA in the gilthead sea bream ( Sparus aurata L.). The predicted protein (1063 amino acids) was similar to mammalian TLR9s, showing 21 leucine-rich repeats in the extracellular region and a typical Toll/IL-1R (TIR) domain in the intracellular region. Comparative analysis of TLR9 sequences indicated that critical residues for ligand-binding are conserved across vertebrate lineages, although evidence of functional divergence was observed. Analysis of the genomic structure of sea bream TLR9 gene revealed the presence of two intervening sequences. Retention of the second intron produced an alternatively spliced mRNA ( TLR9B) showing differential expression among tissues or developmental stages compared to the wild-type isoform ( TLR9A). RT-PCR analysis indicated a broad expression of TLR9A, especially in immune-related organs (spleen, head-kidney) and mucosal–epithelial barriers (gills, gut, skin). Using quantitative Real-Time RT-PCR, no statistically significant variation was observed for TLR9 mRNAs expression in the spleen of experimentally infected animals compared to healthy controls. Comparing sequence and expression profile of sea bream TLR9 with mammalian TLR9s suggested that the main function of TLR9 might be conserved across vertebrates, although species-specific features are present (modulation of ligand-binding specificity, alternative splicing).

Paralogy and orthology in the MALVACEAE rpb2 gene family: investigation of gene duplication in hibiscus.

A sample of the second largest subunit of low-copy nuclear RNA polymerase II (rpb2) sequences from Malvaceae subfamily Malvoideae suggests that rpb2 has been duplicated early in the subfamily's history. Hibiscus and related taxa possess two rpb2 genes, both of which produce congruent phylogenetic patterns that are largely concordant with cpDNA topologies. No evidence of functional divergence or disruption was found among duplicated copies, suggesting that long-term maintenance of duplicated copies of rpb2 is usual in this lineage. Therefore, this gene may be suitable for the potential diagnosis of relatively old polyploid events. One probable pseudogene was found in Radyera farragei and a single chimeric sequence was recovered from Howittia trilocularis, suggesting that the rpb2 locus is not as prone to evolutionary processes that can confound phylogenetic inferences based on nDNA sequences. The pattern of relationships among rpb2 sequences, coupled with chromosome number information and Southern hybridization data, suggests that an early polyploid event was not the cause of the duplication, despite independent evidence of paleopolyploidy in some members of Malvoideae. Rpb2 exons and introns together are suitable for phylogenetic analysis, producing well-resolved and well-supported results that were robust to model permutation and congruent with previous studies of subfamily Malvoideae using cpDNA characters.