The fucosyltransferase gene family: an amazing summary of the underlying mechanisms of gene evolution

The fucosyltransferase gene family encodes enzymes that transfer fucose in 1,2, 1,3/4 and 1,6 linkages on a large variety of glycans. The most ancient genes harbour a split coding sequence, and encode enzyme that transfer fucose at or near O- and N-peptidic sites (serine, threonine or chitobiose unit). Conversely, the more recent genes have a monoexonic coding sequence, and encode enzymes that transfer fucose at the glycan periphery. All basic mechanisms of gene evolution contribute to this amazing scenario: exon shuffling, transposition, point mutations, and duplication. As typical examples: (i) exon shuffling leads to the ancestral organization of the 1,6 fucosyltransferase gene; (ii) the ancestor of 1,2 fucosyltransferase genes is reshaped by retrotransposition at the same locus; (iii) duplication associated to point mutations leads to the most recent 1,3/4 fucosyltransferase genes.     

A diffusion model for the fate of tandem gene duplicates in diploids

Suppose one chromosome in one member of a population somehow acquires a duplicate copy of the gene, fully linked to the original gene's locus. Preservation is the event that eventually every chromosome in the population is a descendant of the one which initially carried the duplicate. For a haploid population in which the absence of all copies of the gene is lethal, the probability of preservation has recently been estimated via a diffusion approximation. That approximation is shown to carry over to the case of diploids and arbitrary strong selection against the absence of the gene. The techniques used lead to some new results. In the large population limit, it is shown that the relative probability that descendants of a small number of individuals carrying multiple copies of the gene fix in the population is proportional to the number of copies carried. The probability of preservation is approximated when chromosomes carrying two copies of the gene are subject to additional, fully non-functionalizing mutations, thereby modelling either an additional cost of replicating a longer genome, or a partial duplication of the gene. In the latter case the preservation probability depends only on the mutation rate to null for the duplicated portion of the gene.     

Pigmentary function and evolution of tyrp1 gene duplicates in fish

The function of the tyrosinase‐related protein 1 (Tyrp1) has not yet been investigated in vertebrates basal to tetrapods. Teleost fishes have two duplicates of the tyrp1 gene. Here, we show that the teleost tyrp1 duplicates have distributed the ancestral gene expression in the retinal pigment epithelium (RPE) and melanophores in a species‐specific manner. In medaka embryos, tyrp1a expression is found in the RPE and in melanophores while tyrp1b is only expressed in melanophores. In zebrafish embryos, expression of tyrp1 paralogs overlaps in the RPE and in melanophores. Knockdown of each zebrafish tyrp1 duplicate alone does not show pigmentary defects, but simultaneous knockdown of both tyrp1 genes results in the formation of brown instead of black eumelanin accompanied by severe melanosome defects. Our study suggests that the brown melanosome color in Tyrp1‐deficient vertebrates is an effect of altered eumelanin synthesis. Black eumelanin formation essentially relies on the presence of Tyrp1 and some of its function is most likely conserved from the common ancestor of bony vertebrates.     

The role of gene flow asymmetry along an environmental gradient in constraining local adaptation and range expansion

Theoretically, asymmetric gene flow along an environmental gradient can limit species range expansion by keeping peripheral populations from locally adapting. However, few empirical studies have examined this potentially fundamental evolutionary mechanism. We address this possibility in the cricket Allonemobius socius, which exist along a season‐length gradient where the probability of producing a single generation per year (univoltinism) increases with latitude. As the probability of univoltinism increases northwards, populations are expected to hedge their bets by producing a greater proportion of diapause eggs when exposed to a mild diapause cue. However, gene flow from southern populations may disrupt local adaptation in the north by reducing the proportion of diapause eggs (expected to be 100% in pure univoltine environments). This may limit range expansion along the northern periphery where A. socius compete with A. fasciatus, a sister species that exhibits an invariant diapause‐only egg‐laying strategy. To assess the potential for range limitation, we examined diapause incidence (the proportion of diapause eggs produced under diapause conditions), gene flow symmetry and population structure across nine A. socius populations. We found that gene flow was asymmetric and biased northwards towards the periphery. Furthermore, peripheral populations that inhabited pure univoltine environments produced numerous nondiapause eggs (a southern, bivoltine diapause phenotype), which we assume to be a suboptimal phenotype. These patterns suggest that asymmetric gene flow along the gradient constrains adaptation in peripheral populations, potentially constraining species range expansion.     

Not all anthocyanins are born equal: distinct patterns induced by stress in Arabidopsis

Main Conclusion

Different abiotic stress conditions induce distinct sets of anthocyanins, indicating that anthocyanins have different biological functions, or that decoration patterns of each anthocyanin are used for unique purposes during stress. The induction of anthocyanin accumulation in vegetative tissues is often considered to be a response of plants to biotic or abiotic stress conditions. Arabidopsis thaliana (Arabidopsis) accumulates over 20 anthocyanins derived from the anthocyanidin cyanidin in an organ-specific manner during development, but the anthocyanin chemical diversity for their alleged stress protective functions remains unclear. We show here that, when grown in various abiotic stress conditions, Arabidopsis not only often accumulates significantly higher levels of total anthocyanins, but different stress conditions also favor the accumulation of different sets of anthocyanins. For example, the anthocyanin patterns of seedlings grown at pH 3.3 or in media lacking phosphate are very similar and characterized by relatively high levels of the anthocyanins A8 and A11. In contrast, anthocyanin inductive conditions (AIC) provided by high sucrose media are characterized by high accumulation of A9* and A5 relative to other stress conditions. The modifications present in each condition correlate reasonably well with the induction of the respective anthocyanin modification enzymes. Taken together, our results suggest that Arabidopsis anthocyanin profiles provide ‘fingerprints’ that reflect the stress status of the plants.     

Characterizing gene movements between chromosomes in Drosophila

Drosophila have a variety of innate immune strategies for defending itself from infection, including humoral and cell mediated responses to invading microorganisms. At the front lines of these responses, are a diverse group of pattern recognition receptors that recognize pathogen associated molecular patterns. These patterns include bacterial lipopolysaccharides, peptidoglycans, and fungal β?1,3 glucans. Some of the receptors catalytically modify the pathogenic determinant, but all are responsible for directly facilitating a signaling event that results in an immune response. Some of these events require multiple pattern recognition receptors acting sequentially to activate a pathway. In some cases, a signaling pathway may be activated by a variety of different pathogens, through parallel receptors detecting different pathogenic determinants. In this chapter, we review what is known about pattern recognition receptors in Drosophila, and how those lessons may be applied towards a broader understanding of immunity.     

Mixture analysis of asymmetry: modelling directional asymmetry, antisymmetry and heterogeneity in fluctuating asymmetry

The occurrence of different forms of asymmetry complicates the analysis and interpretation of patterns in asymmetry. Furthermore, between-individual heterogeneity in developmental stability (DS) and thus fluctuating asymmetry (FA), is required to find relationships between DS and other factors. Separating directional asymmetry (DA) and antisymmetry (AS) from real FA and understanding between-individual heterogeneity in FA is therefore crucial in the analysis and interpretation of patterns in asymmetry. In this paper we introduce and explore mixture analysis to (i) identify FA, DA and AS from the distribution of the signed asymmetry, and (ii) to model and quantify between-individual heterogeneity in developmental stability and FA. In addition, we expand mixtures to the estimation of the proportion of variation in the unsigned FA that can be attributed to between-individual heterogeneity in the presumed underlying developmental stability (the so-called hypothetical repeatability). Finally, we construct weighted normal probability plots to investigate the assumption of underlying normality of the different components. We specifically show that (i) model selection based on the likelihood ratio test has the potential to yield models that incorporate nearly all heterogeneity in FA; (ii) mixtures appear to be a powerful and sensitive statistical technique to identify the different forms of asymmetry; (iii) restricted measurement accuracy and the occurrence of many zero observations results in an overestimation of the hypothetical repeatability on the basis of the model parameters; and (iv) as judged from the high correlation coefficients of the normal probability plots, the underlying normality assumption appears to hold for the empirical data we analysed. In conclusion, mixtures provide a useful statistical tool to study patterns in asymmetry.     

Correlated asymmetry of sequence and functional divergence between duplicate proteins of Saccharomyces cerevisiae

The role of sequence divergence in functional divergence of duplicate genes is a topic of great interest. In this study, we compare the numbers of amino acid substitutions in each sequence since two yeast duplicates diverged, using a preduplication ancestral outgroup. Using this strategy, we explored the relationship between sequence divergence and functional divergence between duplicate partners. We show that the degree of relative functional asymmetry between duplicate proteins is proportional to the relative sequence divergence between them. Furthermore, of the two duplicates, the copy closer to their ancestral sequence (fewer number of amino acid substitutions) interacts with more proteins and affects fitness more severely when deleted. Therefore, asymmetric sequence divergence between duplicates is correlated with asymmetric functional divergence and may underlie the duplicate's role in genetic robustness against mutations. Among the functional traits considered, protein abundance appears to have the strongest correlation with the nonsynonymous divergence between duplicates. Taken together with the results from whole-genome analyses, our results indicate that within-species duplicates are subject to the same evolutionary force that acts on interspecific sequence and functional divergence. In particular, we detect signs of purifying selection on the more slowly evolving duplicate.     

Prenatal stress and increased fluctuating asymmetry in the parietal bones of neonatal rats

Recent studies have reported increased fluctuating dental and long bone asymmetry in neonates as a function of prenatal stress. The present study was designed to assess the effects of prenatal stress on a third calcium-dependent system, membranous bone. Pregnant rats were exposed to cold, heat, or noise from conception through parturition. Bilateral parietal bone lengths were measured in the term neonates. Levels of fluctuating asymmetry were found to be significantly increased (p less than .001) in all three stressed groups compared to unstressed controls. Results support the concept of a generalized stress response, and suggestions are offered for human osteological application.     

A novel gene family NBPF: intricate structure generated by gene duplications during primate evolution

Partial and complete genome duplications occurred during evolution and resulted in the creation of new genes and gene families. We identified a novel and intricate human gene family located primarily in regions of segmental duplications on human chromosome 1. We named it NBPF, for neuroblastoma breakpoint family, because one of its members is disrupted by a chromosomal translocation in a neuroblastoma patient. The NBPF genes have a repetitive structure with high intragenic and intergenic sequence similarity in both coding and noncoding regions. These similarities might expose these genomic regions to illegitimate recombination, resulting in structural variation in the NBPF genes. The encoded proteins contain a highly conserved domain of unknown function, which we have named the NBPF repeat. In silico analysis combined with the isolation of multiple full-length cDNA clones showed that several members of this gene family are abundantly expressed in a large variety of tissues and cell lines. Strikingly, no discernable orthologues could be identified in the completed genomes of fruit fly, nematode, mouse, or rat, but sequences with low homology could be isolated from the draft canine and bovine genomes. Interestingly, this gene family shows primate-specific duplications that result in species-specific arrays of NBPF homologous sequences. Overall, this novel NBPF family reflects the continuous evolution of primate genomes that resulted in large physiological differences, and its potential role in this process is discussed.     

Duplication and expression analysis of multicopy miRNA gene family members in Arabidopsis and rice

To understand the expansion ofmulticopy microRNA （miRNA） families in plants, we localized the reported miRNA genes from Arabidopsis and rice to their chromosomes, respectively, and observed that 37% of 117 miRNA genes from Arabidopsis and 35% of 173 miRNA genes from rice were segmental duplications in the genome. In order to characterize whether the expression diversification has occurred among plant multicopy miRNA family members, we designed PCR primers targeting 48 predicted miRNA precursors from 10 families in Arabidopsis and rice. Results from RT-PCR data suggest that the transcribed precursors of members within the same miRNA family were present at different expression levels. In addition, although miRl60 and miR162 sequences were conserved in Arabidopsis and rice, we found that the expression patterns of these genes differed between the two species. These data suggested that expression diversification has occurred in multicopy miRNA families, increasing our understanding of the expression regulation of miRNAs in plants.     

Not all scale-free networks are born equal: the role of the seed graph in PPI network evolution

The (asymptotic) degree distributions of the best-known “scale-free” network models are all similar and are independent of the seed graph used; hence, it has been tempting to assume that networks generated by these models are generally similar. In this paper, we observe that several key topological features of such networks depend heavily on the specific model and the seed graph used. Furthermore, we show that starting with the “right” seed graph (typically a dense subgraph of the protein–protein interaction network analyzed), the duplication model captures many topological features of publicly available protein–protein interaction networks very well.     

Left-right asymmetry in embryonic development: a comprehensive review

Embryonic morphogenesis occurs along three orthogonal axes. While the patterning of the anterior-posterior and dorsal-ventral axes has been increasingly well characterized, the left-right (LR) axis has only recently begun to be understood at the molecular level. The mechanisms which ensure invariant LR asymmetry of the heart, viscera, and brain represent a thread connecting biomolecular chirality to human cognition, along the way involving fundamental aspects of cell biology, biophysics, and evolutionary biology. An understanding of LR asymmetry is important not only for basic science, but also for the biomedicine of a wide range of birth defects and human genetic syndromes. This review summarizes the current knowledge regarding LR patterning in a number of vertebrate and invertebrate species, discusses several poorly understood but important phenomena, and highlights some important open questions about the evolutionary origin and conservation of mechanisms underlying embryonic asymmetry.     

Null alleles at two lactate dehydrogenase loci in rainbow trout are associated with decreased developmental stability

Previous studies with rainbow trout (Oncorhynchus mykiss) have shown that allozymic heterozygotes have increased developmental stability, as measured by reduced fluctuating bilateral asymmetry. In this paper, we examine the phenotypic effects of null alleles at two lactate dehydrogenase (LDH) loci. If the association between allozymic heterozygosity and developmental stability is due largely to linked chromosomal segments, then we would expect null allele heterozygotes to have increased developmental stability. In contrast, heterozygotes for LDH null alleles in three populations have reduced developmental stability. This suggests that the reduction in enzyme activity at these loci is having a deleterious effect on development that is strong enough to mask any beneficial effects that may be associated with heterozygosity for these chromosomal segments. The LDH loci examined in this study are members of two different paralogous pairs of duplicate genes produced by the polyploidization of the ancestral salmonid genome. The apparent deleterious effects of these null alleles in heterozygotes could retard the possible loss of duplicate gene expression.     

Recent gene conversions between duplicated glutamate decarboxylase genes (gadA and gadB) in pathogenic Escherichia coli

Escherichia coli have evolved adaptive systems to resist strongly acidic habitats in part through the production of 2 biochemically identical isoforms of glutamate decarboxylase (GAD), encoded by the gadA and gadB genes. These genes occur in E. coli and other members of the genospecies (e.g., Shigella spp.) and originated as part of a genomic fitness island acquired early in Escherichia evolution. The present duplicated gad loci are widely spaced on the E. coli chromosome, and the 2 genes are 97% similar in sequence. Comparison of the nucleotide sequences of the gadA and gadB in 16 strains of pathogenic E. coli revealed 3.8% and 5.0% polymorphism in the 2 genes, respectively. Alignment of the homologous genes identified a total of 120 variable sites, including 21 fixed nucleotide differences between the loci within the first 82 codons of the genes. Twenty-three phylogenetically informative sites were polymorphic for the same nucleotides in both genes suggesting recent gene conversions or intergenic recombination. Phylogenetic analysis based on the synonymous substitutions per synonymous site indicated 2 cases in which specific gadA and gadB alleles were more closely related to one another than to other alleles at the corresponding locus. The results indicate that at least 3 gene conversion events have occurred after the gad gene duplication in the evolution of E. coli. Despite multiple gene conversion events, the upstream regulatory regions and the 5' end of each gene remains distinct, suggesting that maintaining functionally different gad genes is important in this acid-resistance mechanism in pathogenic E. coli.     

Loss of protein interactions and regulatory divergence in yeast whole-genome duplicates

Whole-genome duplications are important for the growth of genome complexity. We investigated various factors involved in the evolution of yeast whole-genome duplicates (ohnologs) making emphasis on the analysis of protein interactions. We found that ohnologs have a lower number of protein interactions compared with small-scale duplicates and singletons (by about − 40%). The loss of interactions was proportional to their initial number and independent of ohnolog position in the protein interaction network. A faster evolving member of an ohnolog pair has a lower number of interactions compared to its counterpart. The Gene Ontology mapping of non-overlapping and overlapping interactants of paired ohnologs reveals a sharp asymmetry in GO terms related to regulation. The fraction of these terms is much higher in non-overlapping interactants (compared to overlapping interactants and total dataset). Network clustering coefficient is lower in ohnologs, yet they show an increased density of protein interactions restricted within the whole ohnologs set. These facts suggest that subfunctionalization (or subneofunctionalization) reflected in the loss of protein interactions was a prevailing process in the divergence of ohnologs, which distinguishes them from small-scale duplicates. The loss of protein interactions was associated with the regulatory divergence between the members of an ohnolog pair. A small-scale modularity (reflected in clustering coefficient) probably was not important for ohnologs retention, yet a larger-scale modularity could be involved in their evolution.     

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.