Purifying selection and birth-and-death evolution in the histone H4 gene family

Parallel gene size and isoform expansion of ancient neuronal genes

How nervous systems evolved is a central question in biology. A diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning hundreds of thousands of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. Although many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and, in many cases, the emergence of neurons as dedicated cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes potentially driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection, as demonstrated by low dN/dS ratios, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.

Arabidopsis thaliana: a powerful model organism to explore histone modifications and their upstream regulations

Histones are subjected to extensive covalent modifications that affect inter-nucleosomal interactions as well as alter chromatin structure and DNA accessibility. Through switching the corresponding histone modifications, the level of transcription and diverse downstream biological processes can be regulated. Although animal systems are widely used in studying histone modifications, the signalling processes that occur outside the nucleus prior to histone modifications have not been well understood due to the limitations including non viable mutants, partial lethality, and infertility of survivors. Here, we review the benefits of using Arabidopsis thaliana as the model organism to study histone modifications and their upstream regulations. Similarities among histones and key histone modifiers such as the Polycomb group (PcG) and Trithorax group (TrxG) in Drosophila, Human, and Arabidopsis are examined. Furthermore, prolonged cold-induced vernalization system has been well-studied and revealed the relationship between the controllable environment input (duration of vernalization), its chromatin modifications of FLOWERING LOCUS C (FLC), following gene expression, and the corresponding phenotypes. Such evidence suggests that research on Arabidopsis can bring insights into incomplete signalling pathways outside of the histone box, which can be achieved through viable reverse genetic screenings based on the phenotypes instead of direct monitoring of histone modifications among individual mutants. The potential upstream regulators in Arabidopsis can provide cues or directions for animal research based on the similarities between them.

Ancient origins of complex neuronal genes

How nervous systems evolved is a central question in biology. An increasing diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning millions of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. While many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and in many cases the emergence of dedicated neuronal cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection as demonstrated by low dN/dS scores, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.

Casting histone variants during mammalian reproduction

During mammalian reproduction, germ cell chromatin packaging is key to prepare parental genomes for fertilization and to initiate embryonic development. While chromatin modifications such as DNA methylation and histone post-translational modifications are well known to carry regulatory information, histone variants have received less attention in this context. Histone variants alter the stability, structure and function of nucleosomes and, as such, contribute to chromatin organization in germ cells. Here, we review histone variants expression dynamics during the production of male and female germ cells, and what is currently known about their parent-of-origin effects during reproduction. Finally, we discuss the apparent conundrum behind these important functions and their recent evolutionary diversification.

Varying strength of selection contributes to the intragenomic diversity of rRNA genes

Ribosome biogenesis in eukaryotes is supported by hundreds of ribosomal RNA (rRNA) gene copies that are encoded in the ribosomal DNA (rDNA). The multiple copies of rRNA genes are thought to have low sequence diversity within one species. Here, we present species-wide rDNA sequence analysis in Saccharomyces cerevisiae that challenges this view. We show that rDNA copies in this yeast are heterogeneous, both among and within isolates, and that many variants avoided fixation or elimination over evolutionary time. The sequence diversity landscape across the rDNA shows clear functional stratification, suggesting different copy-number thresholds for selection that contribute to rDNA diversity. Notably, nucleotide variants in the most conserved rDNA regions are sufficiently deleterious to exhibit signatures of purifying selection even when present in only a small fraction of rRNA gene copies. Our results portray a complex evolutionary landscape that shapes rDNA sequence diversity within a single species and reveal unexpectedly strong purifying selection of multi-copy genes.

Evolution of methyltransferase like (METTL) proteins in Metazoa: A complex gene family involved in epitranscriptomic regulation and other epigenetic processes

The methyltransferase-like (METTL) proteins constitute a family of seven-beta-strand methyltransferases with S-adenosyl methionine-binding domains that modify DNA, RNA, and proteins. Methylation by METTL proteins contributes to the epigenetic, and in the case of RNA modifications, epitranscriptomic regulation of a variety of biological processes. Despite their functional importance, most investigations of the substrates and functions of METTLs within metazoans have been restricted to model vertebrate taxa. In the present work, we explore the evolutionary mechanisms driving the diversification and functional differentiation of 33 individual METTL proteins across Metazoa. Our results show that METTLs are nearly ubiquitous across the animal kingdom, with most having arisen early in metazoan evolution (i.e., occur in basal metazoan phyla). Individual METTL lineages each originated from single independent ancestors, constituting monophyletic clades, which suggests that each METTL was subject to strong selective constraints driving its structural and/or functional specialization. Interestingly, a similar process did not extend to the differentiation of nucleoside-modifying and protein-modifying METTLs (i.e., each METTL type did not form a unique monophyletic clade). The members of these two types of METTLs also exhibited differences in their rates of evolution. Overall, we provide evidence that the long-term evolution of METTL family members was driven by strong purifying selection, which in combination with adaptive selection episodes, led to the functional specialization of individual METTL lineages. This work contributes useful information regarding the evolution of a gene family that fulfills a variety of epigenetic functions, and can have profound influences on molecular processes and phenotypic traits.

Evolutionary Dynamics of the Repetitive DNA in the Karyotypes of Pipa carvalhoi and Xenopus tropicalis (Anura, Pipidae)

The large amphibian genomes contain numerous repetitive DNA components that have played an important role in the karyotypic diversification of this vertebrate group. Hypotheses based on the presumable primitive karyotype (2n = 20) of the anurans of the family Pipidae suggest that they have evolved principally through intrachromosomal rearrangements. Pipa is the only South American pipid, while all the other genera are found in Africa. The divergence of the South American lineages from the African ones occurred at least 136 million years ago and is thought to have had a strong biogeographic component. Here, we tested the potential of the repetitive DNA to enable a better understanding of the differentiation of the karyotype among the family Pipidae and to expand our capacity to interpret the chromosomal evolution in this frog family. Our results indicate a long history of conservation in the chromosome bearing the H3 histone locus, corroborating inferences on the chromosomal homologies between the species in pairs 6, 8, and 9. The chromosomal distribution of the microsatellite motifs also provides useful markers for comparative genomics at the chromosome level between Pipa carvalhoi and Xenopus tropicalis, contributing new insights into the evolution of the karyotypes of these species. We detected similar patterns in the distribution and abundance of the microsatellite arrangements, which reflect the shared organization in the terminal/subterminal region of the chromosomes between these two species. By contrast, the microsatellite probes detected a differential arrangement of the repetitive DNA among the chromosomes of the two species, allowing longitudinal differentiation of pairs that are identical in size and morphology, such as pairs 1, 2, 4, and 5. We also found evidence of the distinctive composition of the repetitive motifs of the centromeric region between the species analyzed in the present study, with a clear enrichment of the (CA) and (GA) microsatellite motifs in P. carvalhoi. Finally, microsatellite enrichment in the pericentromeric region of chromosome pairs 6, 8, and 9 in the P. carvalhoi karyotype, together with interstitial telomeric sequences (ITS), validate the hypothesis that pericentromeric inversions occurred during the chromosomal evolution of P. carvalhoi and reinforce the role of the repetitive DNA in the remodeling of the karyotype architecture of the Pipidae.

Discovery of polypeptide ligand-receptor pairs based on their co-evolution

Sequencing diverse genomes allowed the tracing of orthologous and paralogous genes to understand the co-evolution of polypeptide ligands and receptors. This review documents the discovery of several polypeptide ligands and their cognate receptors mainly expressed in the reproductive tissue using evolutionary genomics. We discussed the sub-functionization of paralogs and co-evolution of ligand-receptor families. Based on the conserved signaling among paralogous receptors and common knock-out phenotypes of ligand-receptor pairs, relationships between relaxin family peptides and leucine-rich repeat-containing, G protein-coupled receptors (LGR) were revealed. We also described the identification of a novel paralogous glycoprotein hormone thyrostimulin and design of a long-acting FSH. Human stresscopin and stresscopin-related peptide, paralogous to CRH, were also identified based on the conserved signaling pathways. Recently, a novel ligand placensin expressed in human placenta was found based on the paralogous relationship with a metabolic hormone asprosin. Placensin likely contributes to stage-dependent increases in insulin resistance during human pregnancy and its elevated secretion was associated with gestational diabetes mellitus. Although many ligands were predicted based on sequence signatures, ligands of shorter sequences have not been identified, together with many "orphan" receptors without known ligands. Future development of tools for predicting ligands and high throughput assays to identify ligand-receptor pairs based on ligand binding and/or signal transduction could advance hormone-based physiology and pathophysiology.

The origin of chromosomal histones in a 30S ribosomal protein

Histones are genes that regulate chromatin structure. They are present in both eukaryotes and archaea, and form nucleosomes with DNA, but their exact evolutionary origins have hitherto remained a mystery. A longstanding hypothesis is that they have precursors in ribosomal proteins with whom they share much in common in terms of size and chemistry. By examining the proteome of the Asgard archaeon, Lokiarchaeum, the most conserved of all the histones, H4, is found to plausibly be homologous with one of its 30S ribosomal proteins, RPS6. This is based on both sequence identity and statistical analysis. The N-terminal tail, containing key sites involved in post-translational modifications, is notably present in the precursor gene. Although other archaeal groups possess similar homologs, these are not as close to H4 as the one found in Lokiarchaeum. The other core histones, H2A, H2B and H3, appear to have also evolved from the same ribosomal protein. Parts of H4 are also similar to another ribosomal protein, namely RPS15, suggesting that the ancestral precursor could have resembled both. Eukaryotic histones, in addition, appear to have an independent origin to that of their archaeal counterparts that evolved from similar, but still different, 30S subunit proteins, some of which are much more like histones in terms of their physical structure. The nucleosome may, therefore, be not only of archaeal but also of ribosomal origin.

Episodic positive diversifying selection on key immune system genes in major avian lineages

The major histocompatibility complex (MHC) of the adaptive immune system and the toll-like receptor (TLR) family of the innate immune system are involved in the detection of foreign invaders, and thus are subject to parasite-driven molecular evolution. Herein, we tested for macroevolutionary signatures of selection in these gene families within and among all three major clades of birds (Paleognathae, Galloanserae, and Neoaves). We characterized evolutionary relationships of representative immune genes (Mhc1 and Tlr2b) and a control gene (ubiquitin, Ubb), using a relatively large and phylogenetically diverse set of species with complete coding sequences (34 orthologous loci for Mhc1, 29 for Tlr2b, and 37 for Ubb). Episodic positive diversifying selection was found in the gene-wide phylogenies of the two immune genes, as well as at specific sites within each gene (8.5% of codon sites in Mhc1 and 2.7% in Tlr2b), but not in the control gene (Ubb). We found 20% of lineages under episodic diversifying selection in Mhc1 versus 9.1% in Tlr2b. For Mhc1, selection was relaxed in the Galloanserae and intensified in the Neoaves relative to the other clades, but no differences were detected among clades in the Tlr2b gene. In summary, we provide evidence of episodic positive diversifying selection in key immune genes and demonstrate differential strengths of selection within Class Aves, with the adaptive gene showing an increased divergence and evolutionary rate over the innate gene, contributing to the growing understanding of vertebrate immune gene evolution.

Evolution and Diversification of FRUITFULL Genes in Solanaceae

Ecologically and economically important fleshy edible fruits have evolved from dry fruit numerous times during angiosperm diversification. However, the molecular mechanisms that underlie these shifts are unknown. In the Solanaceae there has been a major shift to fleshy fruits in the subfamily Solanoideae. Evidence suggests that an ortholog of FRUITFULL (FUL), a transcription factor that regulates cell proliferation and limits the dehiscence zone in the silique of Arabidopsis, plays a similar role in dry-fruited Solanaceae. However, studies have shown that FUL orthologs have taken on new functions in fleshy fruit development, including regulating elements of tomato ripening such as pigment accumulation. FUL belongs to the core eudicot euFUL clade of the angiosperm AP1/FUL gene lineage. The euFUL genes fall into two paralogous clades, euFULI and euFULII. While most core eudicots have one gene in each clade, Solanaceae have two: FUL1 and FUL2 in the former, and MBP10 and MBP20 in the latter. We characterized the evolution of the euFUL genes to identify changes that might be correlated with the origin of fleshy fruit in Solanaceae. Our analyses revealed that the Solanaceae FUL1 and FUL2 clades probably originated through an early whole genome multiplication event. By contrast, the data suggest that the MBP10 and MBP20 clades are the result of a later tandem duplication event. MBP10 is expressed at weak to moderate levels, and its atypical short first intron lacks putative transcription factor binding sites, indicating possible pseudogenization. Consistent with this, our analyses show that MBP10 is evolving at a faster rate compared to MBP20. Our analyses found that Solanaceae euFUL gene duplications, evolutionary rates, and changes in protein residues and expression patterns are not correlated with the shift in fruit type. This suggests deeper analyses are needed to identify the mechanism underlying the change in FUL ortholog function.

Contrasting patterns of coding and flanking region evolution in mammalian keratin associated protein-1 genes

Mammalian genomes contain a number of duplicated genes, and sequence identity between these duplicates can be maintained by purifying selection. However, between-duplicate recombination can also maintain sequence identity between copies, resulting in a pattern known as concerted evolution where within-genome repeats are more similar to each other than to orthologous repeats in related species. Here we investigated the tandemly-repeated keratin-associated protein 1 (KAP1) gene family, KRTAP1, which encodes proteins that are important components of hair and wool in mammals. Comparison of eutherian mammal KRTAP1 gene repeats within and between species shows a strong pattern of concerted evolution. However, in striking contrast to the coding regions of these genes, we find that the flanking regions have a divergent pattern of evolution. This contrast in evolutionary pattern transitions abruptly near the start and stop codons of the KRTAP1 genes. We reveal that this difference in evolutionary patterns is not explained by conventional purifying selection, nor is it likely a consequence of codon adaptation or reverse transcription of KRTAP1-n mRNA. Instead, the evidence suggests that these contrasting patterns result from short-tract gene conversion events that are biased to the KRTAP1 coding region by selection and/or differential sequence divergence. This work demonstrates the power that gene conversion has to finely shape the evolution of repetitive genes, and provides another distinctive pattern of contrasting evolutionary outcomes that results from gene conversion. A greater emphasis on exploring the evolution of multi-gene eukaryotic families will reveal how common different contrasting evolutionary patterns are in gene duplicates.

Shedding light on lipid metabolism in Kinetoplastida: A phylogenetic analysis of phospholipase D protein homologs

Unicellular flagellates that make up the class Kinetoplastida include multiple parasites responsible for public health concerns, including Trypanosoma brucei and T. cruzi (agents of African sleeping sickness and Chagas disease, respectively), and various Leishmania species, which cause leishmaniasis. These diseases are generally difficult to eradicate, with treatments often having lethal side effects and/or being effective only during the acute phase of the diseases, when most patients are still asymptomatic. Phospholipid signaling and metabolism are important in the different life stages of Trypanosoma, including playing a role in transitions between stages and in immune system evasion, thus, making the responsible enzymes into potential therapeutic targets. However, relatively little is understood about how the pathways function in these pathogens. Thus, in this study we examined evolutionary history of proteins from one such signaling pathway, namely phospholipase D (PLD) homologs. PLD is an enzyme responsible for synthesizing phosphatidic acid (PA) from membrane phospholipids. PA is not only utilized for phospholipid synthesis, but is also involved in many other signaling pathways, including biotic and abiotic stress response. 37 different representative Kinetoplastida genomes were used for an exhaustive search to identify putative PLD homologs. The genome of Bodo saltans was the only one of surveyed Kinetoplastida genomes that encoded a protein that clustered with plant PLDs. The representatives from other Kinetoplastida species clustered together in two different clades, thought to be homologous to the PLD superfamily, but with shared sequence similarity with cardiolipin synthases (CLS), and phosphatidylserine synthases (PSS). The protein structure predictions showed that most Kinetoplastida sequences resemble CLS and PSS, with the exception of 5 sequences from Bodo saltans that shared significant structural similarities with the PLD sequences, suggesting the loss of PLD-like sequences during the evolution of parasitism in kinetoplastids. On the other hand, diacylglycerol kinase (DGK) homologs were identified for all species examined in this study, indicating that DGK could be the only pathway for the synthesis of PA involved in lipid signaling in these organisms due to genome streamlining during transition to parasitic lifestyle. Our findings offer insights for development of potential therapeutic and/or intervention approaches, particularly those focused on using PA, PLD and/or DGK related pathways, against trypanosomiasis, leishmaniasis, and Chagas disease.

Evolutionary origins and diversification of testis-specific short histone H2A variants in mammals

Eukaryotic genomes must accomplish both compact packaging for genome stability and inheritance, as well as accessibility for gene expression. They do so using post-translational modifications of four ancient canonical histone proteins (H2A, H2B, H3, and H4) and by deploying histone variants with specialized chromatin functions. Some histone variants are conserved across all eukaryotes, whereas others are lineage-specific. Here, we performed detailed phylogenomic analyses of “short H2A histone” variants found in mammalian genomes. We discovered a previously undescribed typically-sized H2A variant in monotremes and marsupials, H2A.R, which may represent the common ancestor of the short H2As. We also discovered a novel class of short H2A histone variants in eutherian mammals, H2A.Q. We show that short H2A variants arose on the X Chromosome in the common ancestor of all eutherian mammals and diverged into four evolutionarily distinct clades: H2A.B, H2A.L, H2A.P, and H2A.Q. However, the repertoires of short histone H2A variants vary extensively among eutherian mammals due to lineage-specific gains and losses. Finally, we show that all four short H2As are subject to accelerated rates of protein evolution relative to both canonical and other variant H2A proteins including H2A.R. Our analyses reveal that short H2As are a unique class of testis-restricted histone variants displaying an unprecedented evolutionary dynamism. Based on their X-Chromosomal localization, genetic turnover, and testis-specific expression, we hypothesize that short H2A variants may participate in genetic conflicts involving sex chromosomes during reproduction.

Inferring the evolution of the major histocompatibility complex of wild pigs and peccaries using hybridisation DNA capture-based sequencing

The major histocompatibility complex (MHC) is a key genomic model region for understanding the evolution of gene families and the co-evolution between host and pathogen. To date, MHC studies have mostly focused on species from major vertebrate lineages. The evolution of MHC classical (Ia) and non-classical (Ib) genes in pigs has attracted interest because of their antigen presentation roles as part of the adaptive immune system. The pig family Suidae comprises over 18 extant species (mostly wild), but only the domestic pig has been extensively sequenced and annotated. To address this, we used a DNA-capture approach, with probes designed from the domestic pig genome, to generate MHC data for 11 wild species of pigs and their closest living family, Tayassuidae. The approach showed good efficiency for wild pigs (~80% reads mapped, ~87× coverage), compared to tayassuids (~12% reads mapped, ~4× coverage). We retrieved 145 MHC loci across both families. Phylogenetic analyses show that the class Ia and Ib genes underwent multiple duplications and diversifications before suids and tayassuids diverged from their common ancestor. The histocompatibility genes mostly form orthologous groups and there is genetic differentiation for most of these genes between Eurasian and sub-Saharan African wild pigs. Tests of selection showed that the peptide-binding region of class Ib genes was under positive selection. These findings contribute to better understanding of the evolutionary history of the MHC, specifically, the class I genes, and provide useful data for investigating the immune response of wild populations against pathogens.

Dissection of genomic features and variations of three pathotypes of Puccinia striiformis through whole genome sequencing

Stripe rust of wheat, caused by Puccinia striiformis f. sp. tritici, is one of the important diseases of wheat. We used NGS technologies to generate a draft genome sequence of two highly virulent (46S 119 and 31) and a least virulent (K) pathotypes of P. striiformis from the Indian subcontinent. We generated ~24,000-32,000 sequence contigs (N50;7.4-9.2 kb), which accounted for ~86X-105X sequence depth coverage with an estimated genome size of these pathotypes ranging from 66.2-70.2 Mb. A genome-wide analysis revealed that pathotype 46S 119 might be highly evolved among the three pathotypes in terms of year of detection and prevalence. SNP analysis revealed that ~47% of the gene sets are affected by nonsynonymous mutations. The extracellular secreted (ES) proteins presumably are well conserved among the three pathotypes, and perhaps purifying selection has an important role in differentiating pathotype 46S 119 from pathotypes K and 31. In the present study, we decoded the genomes of three pathotypes, with 81% of the total annotated genes being successfully assigned functional roles. Besides the identification of secretory genes, genes essential for pathogen-host interactions shall prove this study as a huge genomic resource for the management of this disease using host resistance.

Disparate Modes of Evolution Shaped Modern Prion (PRNP) and Prion-Related Doppel (PRND) Variation in Domestic Cattle

Previous investigations aimed at determining whether the mammalian prion protein actually facilitates tangible molecular aspects of either a discrete or pleiotropic functional niche have been debated, especially given the apparent absence of overt behavioral or physiological phenotypes associated with several mammalian prion gene (PRNP) knockout experiments. Moreover, a previous evaluation of PRNP knockout cattle concluded that they were normal, suggesting that the bovine prion protein is physiologically dispensable. Herein, we examined the frequency and distribution of nucleotide sequence variation within the coding regions of bovine PRNP and the adjacent Doppel (PRND) gene, a proximal paralogue to PRNP on BTA13. Evaluation of PRND variation demonstrated that the gene does not depart from a strictly neutral model of molecular evolution, and would therefore not be expected to influence tests of selection within PRNP. Collectively, our analyses confirm that intense purifying selection is indeed occurring directly on bovine PRNP, which is indicative of a protein with an important role. These results suggest that the lack of observed fitness effects may not manifest in the controlled environmental conditions used to care for and raise PRNP knockout animals.

Patterns of Gene Conversion in Duplicated Yeast Histones Suggest Strong Selection on a Coadapted Macromolecular Complex

We find evidence for interlocus gene conversion in five duplicated histone genes from six yeast species. The sequences of these duplicated genes, surviving from the ancient genome duplication, show phylogenetic patterns inconsistent with the well-resolved orthology relationships inferred from a likelihood model of gene loss after the genome duplication. Instead, these paralogous genes are more closely related to each other than any is to its nearest ortholog. In addition to simulations supporting gene conversion, we also present evidence for elevated rates of radical amino acid substitutions along the branches implicated in the conversion events. As these patterns are similar to those seen in ribosomal proteins that have undergone gene conversion, we speculate that in cases where duplicated genes code for proteins that are a part of tightly interacting complexes, selection may favor the fixation of gene conversion events in order to maintain high protein identities between duplicated copies.

Cloning and Functional Analysis of Histones H3 and H4 in Nuclear Shaping during Spermatogenesis of the Chinese Mitten Crab, Eriocheir sinensis

During spermatogenesis in most animals, the basic proteins associated with DNA are continuously changing and somatic-typed histones are partly replaced by sperm-specific histones, which are then successively replaced by transition proteins and protamines. With the replacement of sperm nuclear basic proteins, nuclei progressively undergo chromatin condensation. The Chinese Mitten Crab (Eriocheir sinensis) is also known as the hairy crab or river crab (phylum Arthropoda, subphylum Crustacea, order Decapoda, and family Grapsidae). The spermatozoa of this species are aflagellate, and each has a spherical acrosome surrounded by a cup-shaped nucleus, peculiar to brachyurans. An interesting characteristic of the E. sinensis sperm nucleus is its lack of electron-dense chromatin. However, its formation is not clear. In this study, sequences encoding histones H3 and H4 were cloned by polymerase chain reaction amplification. Western blotting indicated that H3 and H4 existed in the sperm nuclei. Immunofluorescence and ultrastructural immunocytochemistry demonstrated that histones H3 and H4 were both present in the nuclei of spermatogonia, spermatocytes, spermatids and mature spermatozoa. The nuclear labeling density of histone H4 decreased in sperm nuclei, while histone H3 labeling was not changed significantly. Quantitative real-time PCR showed that the mRNA expression levels of histones H3 and H4 were higher at mitotic and meiotic stages than in later spermiogenesis. Our study demonstrates that the mature sperm nuclei of E. sinensis contain histones H3 and H4. This is the first report that the mature sperm nucleus of E. sinensis contains histones H3 and H4. This finding extends the study of sperm histones of E. sinensis and provides some basic data for exploring how decapod crustaceans form uncondensed sperm chromatin.

Structure and Organization of the Engraulidae Family U2 snRNA: An Evolutionary Model Gene?

The U2 snRNA multigene family has been analyzed in four species of the Engraulidae family--Engraulis encrasicolus, Engraulis mordax, Engraulis ringens, and Engraulis japonicas--with the object of understanding more about the structure of this multigene family in these pelagic species and studying their phylogenetic relationships. The results showed that the cluster of this gene family in the Engraulis genus is formed by the U2-U5 snRNA with highly conserved sequences of mini- and micro-satellites, such as (CTGT)n, embedded downstream of the transcription unit; findings indicate that this gene family evolved following the concerted model. The phylogenetic analysis of the non-transcribed spacer of cluster U2-U5 snDNA in the 4 species showed that the sequences of the species E. encrasicolus and E. japonicus are closely related; these two are genetically close to E. mordax and slightly more distant from E. ringens. The data obtained by molecular analysis of U2-U5 snDNA and their secondary structure, with the presence of the micro-satellite (CTGT)n and mini-satellites, show clearly that the species E. encrasicolus and E. japonicus are closely related and would be older than E. mordax and E. ringens.

Methyl jasmonate represses growth and affects cell cycle progression in cultured Taxus cells

Methyl jasmonate elicitation of Taxus cultures enhances paclitaxel accumulation, but represses growth by inhibition of cell cycle progression. Growth repression is evident both at the culture level and transcriptional level. Methyl jasmonate (MeJA) elicitation is an effective strategy to induce and enhance synthesis of the anticancer agent paclitaxel (Taxol(®)) in Taxus cell suspension cultures; however, concurrent decreases in growth are often observed, which is problematic for large-scale bioprocessing. Here, increased accumulation of paclitaxel in Taxus cuspidata suspension cultures with MeJA elicitation was accompanied by a concomitant decrease in cell growth, evident within the first 3 days post-elicitation. Both MeJA-elicited and mock-elicited cultures exhibited similar viability with no apoptosis up to day 16 and day 24 of the cell culture period, respectively, suggesting that growth repression is not attributable to cell death. Flow cytometric analyses demonstrated that MeJA perturbed cell cycle progression of asynchronously dividing Taxus cells. MeJA slowed down cell cycle progression, impaired the G1/S transition as observed by an increase in G0/G1 phase cells, and decreased the number of actively dividing cells. Through a combination of deep sequencing and gene expression analyses, the expression status of Taxus cell cycle-associated genes correlated with observations at the culture level. Results from this study provide valuable insight into the mechanisms governing MeJA perception and subsequent events leading to repression of Taxus cell growth.

Independent evolutionary origin of fem paralogous genes and complementary sex determination in hymenopteran insects

The primary signal of sex determination in the honeybee, the complementary sex determiner (csd) gene, evolved from a gene duplication event from an ancestral copy of the fem gene. Recently, other paralogs of the fem gene have been identified in several ant and bumblebee genomes. This discovery and the close phylogenetic relationship of the paralogous gene sequences led to the hypothesis of a single ancestry of the csd genetic system of complementary sex determination in the Hymenopteran insects, in which the fem and csd gene copies evolved as a unit in concert with the mutual transfers of sequences (concerted evolution). Here, we show that the paralogous gene copies evolved repeatedly through independent gene duplication events in the honeybee, bumblebee, and ant lineage. We detected no sequence tracts that would indicate a DNA transfer between the fem and the fem1/csd genes between different ant and bee species. Instead, we found tracts of duplication events in other genomic locations, suggesting that gene duplication was a frequent event in the evolution of these genes. These and other evidences suggest that the fem1/csd gene originated repeatedly through gene duplications in the bumblebee, honeybee, and ant lineages in the last 100 million years. Signatures of concerted evolution were not detectable, implicating that the gene tree based on neutral synonymous sites represents the phylogenetic relationships and origins of the fem and fem1/csd genes. Our results further imply that the fem1 and csd gene in bumblebees, honeybees, and ants are not orthologs, because they originated independently from the fem gene. Hence, the widely shared and conserved complementary sex determination mechanism in Hymenopteran insects is controlled by different genes and molecular processes. These findings highlight the limits of comparative genomics and emphasize the requirement to study gene functions in different species and major hymenopteran lineages.

Combined haploinsufficiency and purifying selection drive retention of RPL36a paralogs in Arabidopsis

Whole-genome duplication events have driven to a large degree the evolution of angiosperm genomes. Although the majority of redundant gene copies after a genome duplication are lost, subfunctionalization or gene balance account for the retention of gene copies. The Arabidopsis 80S ribosome represents an excellent model to test the gene balance hypothesis as it consists of 80 ribosomal proteins, all of them encoded by genes belonging to small gene families. Here, we present the isolation of mutant alleles of the APICULATA2 (API2) and RPL36aA paralogous genes, which encode identical ribosomal proteins but share a similarity of 89% in their coding sequences. RPL36aA was found expressed at a higher level than API2 in the wild type. The loss-of-function api2 and rpl36aa mutations are recessive and affect leaf development in a similar way. Their double mutant combinations with asymmetric leaves2-1 (as2-1) caused leaf polarity defects that were stronger in rpl36aa as2-1 than in api2 as2-1. Our results highlight the role of combined haploinsufficiency and purifying selection in the retention of these paralogous genes in the Arabidopsis genome.

Disparate molecular evolution of two types of repetitive DNAs in the genome of the grasshopper Eyprepocnemis plorans

Wide arrays of repetitive DNA sequences form an important part of eukaryotic genomes. These repeats appear to evolve as coherent families, where repeats within a family are more similar to each other than to other orthologous representatives in related species. The continuous homogenization of repeats, through selective and non-selective processes, is termed concerted evolution. Ascertaining the level of variation between repeats is crucial to determining which evolutionary model best explains the homogenization observed for these sequences. Here, for the grasshopper Eyprepocnemis plorans, we present the analysis of intragenomic diversity for two repetitive DNA sequences (a satellite DNA (satDNA) and the 45S rDNA) resulting from the independent microdissection of several chromosomes. Our results show different homogenization patterns for these two kinds of paralogous DNA sequences, with a high between-chromosome structure for rDNA but no structure at all for the satDNA. This difference is puzzling, considering the adjacent localization of the two repetitive DNAs on paracentromeric regions in most chromosomes. The disparate homogenization patterns detected for these two repetitive DNA sequences suggest that several processes participate in the concerted evolution in E. plorans, and that these mechanisms might not work as genome-wide processes but rather as sequence-specific ones.

Variations and evolution of polyubiquitin genes from ciliates

Polyubiquitin genes from seven ciliate species were amplified, cloned and sequenced. It is estimated that Strombidium sulcatum, Euplotes vannus, E. rariseta and Anteholosticha manca have a polyubiquitin gene of 3 repeats, and A. parawarreni, Paramecium caudatum and Pseudokeronopsis flava 4 repeats. The newly obtained ubiquitins mostly differ from that of humans by 1-5 residues in amino acid sequences. A neighbor-joining tree constructed based on monomeric ubiquitin genes supports the monophyly of an assemblage comprising the litostomateans and some oligohymenophoreans, but not the class Spirotrichea. The monomers from the same species are generally placed together and highly supported for the class Litostomatea, the genera Paramecium and Ichthyophthirius, but not for other species. The non-synonymous/synonymous rate ratio (dN/dS) at the protein level are less than 1, and the synonymous nucleotide differences per synonymous site (p(S)) from intraspecific comparisons are fairly high (0.02-0.72). These results indicate that ciliates have not only the conserved, but also some quite divergent, polyubiquitin genes and confirm that the polyubiquitin genes in ciliates evolve according to the birth-and-death mode of evolution under strong purifying selection.

Molecular evolution of the clustered MIC-3 multigene family of Gossypium species

The Gossypium MIC-3 (Meloidogyne Induced Cotton-3) gene family is of great interest for molecular evolutionary studies because of its uniqueness to Gossypium species, multi-gene content, clustered localization, and root-knot nematode resistance-associated features. Molecular evolution of the MIC-3 gene family was studied in 15 tetraploid and diploid Gossypium genotypes that collectively represent seven phylogenetically distinct genomes. Synonymous (d(S)) and non-synonymous (d(N)) nucleotide substitution rates suggest that the second of the two exons of the MIC-3 genes has been under strong positive selection pressure, while the first exon has been under strong purifying selection to preserve function. Based on nucleotide substitution rates, we conclude that MIC-3 genes are evolving by a birth-and-death process and that a 'gene amplification' mechanism has helped to retain all duplicate copies, which best fits with the "bait and switch" model of R-gene evolution. The data indicate MIC-3 gene duplication events occurred at various rates, once per 1 million years (MY) in the allotetraploids, once per ~2 MY in the A/F genome clade, and once per ~8 MY in the D-genome clade. Variations in the MIC-3 gene family seem to reflect evolutionary selection for increased functional stability, while also expanding the capacity to develop novel "switch" pockets for responding to diverse pests and pathogens. Such evolutionary roles are congruent with the hypothesis that members of this unique resistance gene family provide fitness advantages in Gossypium.

The 5S rDNA family evolves through concerted and birth-and-death evolution in fish genomes: an example from freshwater stingrays

Background

Ribosomal 5S genes are well known for the critical role they play in ribosome folding and functionality. These genes are thought to evolve in a concerted fashion, with high rates of homogenization of gene copies. However, the majority of previous analyses regarding the evolutionary process of rDNA repeats were conducted in invertebrates and plants. Studies have also been conducted on vertebrates, but these analyses were usually restricted to the 18S, 5.8S and 28S rRNA genes. The recent identification of divergent 5S rRNA gene paralogs in the genomes of elasmobranches and teleost fishes indicate that the eukaryotic 5S rRNA gene family has a more complex genomic organization than previously thought. The availability of new sequence data from lower vertebrates such as teleosts and elasmobranches enables an enhanced evolutionary characterization of 5S rDNA among vertebrates.

Results

We identified two variant classes of 5S rDNA sequences in the genomes of Potamotrygonidae stingrays, similar to the genomes of other vertebrates. One class of 5S rRNA genes was shared only by elasmobranches. A broad comparative survey among 100 vertebrate species suggests that the 5S rRNA gene variants in fishes originated from rounds of genome duplication. These variants were then maintained or eliminated by birth-and-death mechanisms, under intense purifying selection. Clustered multiple copies of 5S rDNA variants could have arisen due to unequal crossing over mechanisms. Simultaneously, the distinct genome clusters were independently homogenized, resulting in the maintenance of clusters of highly similar repeats through concerted evolution.

Conclusions

We believe that 5S rDNA molecular evolution in fish genomes is driven by a mixed mechanism that integrates birth-and-death and concerted evolution.

Cap-assisted internal initiation of translation of histone H4

In eukaryotes, a crucial step of translation initiation is the binding of the multifactor complex eIF4F to the 5' end of the mRNA, a prerequisite to recruitment of the activated small ribosomal 43S particle. Histone H4 mRNAs have short 5'UTRs, which do not conform to the conventional scanning-initiation model. Here we show that the ORF of histone mRNA contains two structural elements critical for translation initiation. One of the two structures binds eIF4E without the need of the cap. Ribosomal 43S particles become tethered to this site and directly loaded in the vicinity of the AUG. The other structure, 19 nucleotides downstream of the initiation codon, forms a three-way helix junction, which sequesters the m(7)G cap. This element facilitates direct positioning of the ribosome on the cognate start codon. This unusual translation initiation mode might be considered as a hybrid mechanism between the canonical and the IRES-driven translation initiation process.

Duplication, concerted evolution and purifying selection drive the evolution of mosquito vitellogenin genes

Background

Mosquito vitellogenin (Vtg) genes belong to a small multiple gene family that encodes the major yolk protein precursors required for egg production. Multiple Vtg genes have been cloned and characterized from several mosquito species, but their origin and molecular evolution are poorly understood.

Results

Here we used in silico and molecular cloning techniques to identify and characterize the evolution of the Vtg gene family from the genera Culex, Aedes/Ochlerotatus, and Anopheles. We identified the probable ancestral Vtg gene among different mosquito species by its conserved association with a novel gene approximately one kilobase upstream of the start codon. Phylogenetic analysis indicated that the Vtg gene family arose by duplication events, but that the pattern of duplication was different in each mosquito genera. Signatures of purifying selection were detected in Culex, Aedes and Anopheles. Gene conversion is a major driver of concerted evolution in Culex, while unequal crossover is likely the major driver of concerted evolution in Anopheles. In Aedes, smaller fragments have undergone gene conversion events.

Conclusions

The study shows concerted evolution and purifying selection shaped the evolution of mosquito Vtg genes following gene duplication. Additionally, similar evolutionary patterns were observed in the Vtg genes from other invertebrate and vertebrate organisms, suggesting that duplication, concerted evolution and purifying selection may be the major evolutionary forces driving Vtg gene evolution across highly divergent taxa.

Containment of extended length polymorphisms in silk proteins

The spider silk gene family to the current date has been developed by gene duplication and homogenization events as well as conservation of crucial sequence parts. These evolutionary processes have created an amazing diversity of silk types each associated with specific properties and functions. In addition, they have led to allelic and gene variants within a species as exemplified by the major ampullate spidroin 1 gene of Nephila clavipes. Due to limited numbers of individuals screened to date little is known about the extent of these heterogeneities and how they are finally manifested in the proteins. Using expanded sample sizes, we show that sequence variations expressed as deletions or insertions of tri-nucleotides lead to different sized and structured repetitive units throughout a silk protein. Moreover, major ampullate spidroins 1 can quite dramatically differ in their overall lengths; however, extreme variants do not spread widely in a spider population. This suggests that a certain size range stabilized by purifying selection is important for spidroin 1 gene integrity and protein function. More than one locus for spidroin 1 genes possibly exist within one individual genome, which are homogenized in size, are differentially expressed and give a spider a certain degree of adaptation on silk's composition and properties. Such mechanisms are shared to a lesser extent by the second major ampullate spidroin gene.

Birth-and-death long-term evolution promotes histone H2B variant diversification in the male germinal cell line

The rich diversity within each of the five histone families (H1, H2A, H2B, H3, and H4) can hardly be reconciled with the notion of homogenizing evolution. The prevalence of birth-and-death long-term evolution over concerted evolution has already been demonstrated in the linker histone H1 family as well as for the H2A, H3, and H4 core histone families. However, information about histone H2B is lacking. In the present work, we have analyzed the diversity of the members of this histone family across different eukaryotic genomes and have characterized the mechanisms involved in their long-term evolution. Our results reveal that, quite in contrast with other histones, H2B variants are subject to a very rapid process of diversification that primarily affects the male germinal cell lineage and involves their functional specialization probably as a consequence of neofunctionalization and subfunctionalization events after gene duplication. The overall parallelism observed between the molecular phylogenies and the relationships among the electrostatic potentials of the different variants suggests that the latter may have played a major structural selective constraint during H2B evolution. It thus seems that the reorganization of chromatin structure during spermiogenesis might have affected the evolutionary constraints driving histone H2B evolution, leading to an increase in diversity.

Probabilistic orthology analysis

Orthology analysis aims at identifying orthologous genes and gene products from different organisms and, therefore, is a powerful tool in modern computational and experimental biology. Although reconciliation-based orthology methods are generally considered more accurate than distance-based ones, the traditional parsimony-based implementation of reconciliation-based orthology analysis (most parsimonious reconciliation [MPR]) suffers from a number of shortcomings. For example, 1) it is limited to orthology predictions from the reconciliation that minimizes the number of gene duplication and loss events, 2) it cannot evaluate the support of this reconciliation in relation to the other reconciliations, and 3) it cannot make use of prior knowledge (e.g., about species divergence times) that provides auxiliary information for orthology predictions. We present a probabilistic approach to reconciliation-based orthology analysis that addresses all these issues by estimating orthology probabilities. The method is based on the gene evolution model, an explicit evolutionary model for gene duplication and gene loss inside a species tree, that generalizes the standard birth-death process. We describe the probabilistic approach to orthology analysis using 2 experimental data sets and show that the use of orthology probabilities allows a more informative analysis than MPR and, in particular, that it is less sensitive to taxon sampling problems. We generalize these anecdotal observations and show, using data generated under biologically realistic conditions, that MPR give false orthology predictions at a substantial frequency. Last, we provide a new orthology prediction method that allows an orthology and paralogy classification with any chosen sensitivity/specificity combination from the spectra of achievable combinations. We conclude that probabilistic orthology analysis is a strong and more advanced alternative to traditional orthology analysis and that it provides a framework for sophisticated comparative studies of processes in genome evolution.

Preliminary study of sperm chromatin characteristics of the brachyuran crab Maja brachydactyla. Histones and nucleosome-like structures in decapod crustacean sperm nuclei previously described without SNBPs

An interesting characteristic of decapod crustacean sperm nuclei is that they do not contain highly packaged chromatin. In the present study we re-examine the presence of DNA-interacting proteins in sperm nuclei of the brachyuran Maja brachydactyla. Although previous reports have indicated that, unlike the majority of sperm cells, DNA of decapod sperm is not organized by basic proteins, in this work we show that: (1) histones are present in sperm of M. brachydactyla; (2) histones are associated with sperm DNA; (3) histone H3 appears in lower proportions than the other core histones, while histone H2B appears in higher proportions; and (4) histone H3 in sperm nuclei is acetylated. This work complements a previous study of sperm histones of Cancer pagurus and supports the suggestion that decapod crustacean sperm chromatin deserves further attention.

The evolutionary differentiation of two histone H2A.Z variants in chordates (H2A.Z-1 and H2A.Z-2) is mediated by a stepwise mutation process that affects three amino acid residues

Background

The histone H2A family encompasses the greatest number of core histone variants of which the replacement variant H2A.Z is currently one of the most heavily studied. No clear mechanism for the functional variability that H2A.Z imparts to chromatin has yet been proposed. While most of the past studies have referred to H2A.Z generically as a single protein, in vertebrates it is a mixture of two protein forms H2A.Z-1 (previously H2A.Z) and H2A.Z-2 (previously H2A.F/Z or H2A.V) that differ by three amino acids.

Results

We have performed an extensive study on the long-term evolution of H2A.Z across metazoans with special emphasis on the possible selective mechanisms responsible for the differentiation between H2A.Z-1 and H2A.Z-2. Our results reveal a common origin of both forms early in chordate evolution. The evolutionary process responsible for the differentiation involves refined stepwise mutation change within the codons of the three differential residues. This eventually led to differences in the intensity of the selective constraints acting upon the different H2A.Z forms in vertebrates.

Conclusion

The results presented in this work definitively reveal that the existence of H2A.Z-1 and H2A.Z-2 is not a whim of random genetic drift. Our analyses demonstrate that H2A.Z-2 is not only subject to a strong purifying selection but it is significantly more evolutionarily constrained than H2A.Z-1. Whether or not the evolutionary drift between H2A.Z-1 and H2A.Z-2 has resulted in a functional diversification of these proteins awaits further research. Nevertheless, the present work suggests that in the process of their differently constrained evolutionary pathways, these two forms may have acquired new or complementary functions.

Histones and nucleosomes in Cancer sperm (Decapod: Crustacea) previously described as lacking basic DNA-associated proteins: a new model of sperm chromatin

To date several studies have been carried out which indicate that DNA of crustacean sperm is neither bound nor organized by basic proteins and, contrary to the rest of spermatozoa, do not contain highly packaged chromatin. Since this is the only known case of this type among metazoan cells, we have re-examined the composition, and partially the structure, of the mature sperm chromatin of Cancer pagurus, which has previously been described as lacking basic DNA-associated proteins. The results we present here show that: (a) sperm DNA of C. pagurus is bound by histones forming nucleosomes of 170 base pairs, (b) the ratio [histones/DNA] in sperm of two Cancer species is 0.5 and 0.6 (w/w). This ratio is quite lower than the proportion [proteins/DNA] that we found in other sperm nuclei with histones or protamines, whose value is from 1.0 to 1.2 (w/w), (c) histone H4 is highly acetylated in mature sperm chromatin of C. pagurus. Other histones (H3 and H2B) are also acetylated, though the level is much lower than that of histone H4. The low ratio of histones to DNA, along with the high level of acetylation of these proteins, explains the non-compact, decondensed state of the peculiar chromatin in the sperm studied here. In the final section we offer an explanation for the necessity of such decondensed chromatin during gamete fertilization of this species.

Biological applications of the theory of birth-and-death processes

In this review, we discuss applications of the theory of birth-and-death processes to problems in biology, primarily, those of evolutionary genomics. The mathematical principles of the theory of these processes are briefly described. Birth-and-death processes, with some straightforward additions such as innovation, are a simple, natural and formal framework for modeling a vast variety of biological processes such as population dynamics, speciation, genome evolution, including growth of paralogous gene families and horizontal gene transfer and somatic evolution of cancers. We further describe how empirical data, e.g. distributions of paralogous gene family size, can be used to choose the model that best reflects the actual course of evolution among different versions of birth-death-and-innovation models. We conclude that birth-and-death processes, thanks to their mathematical transparency, flexibility and relevance to fundamental biological processes, are going to be an indispensable mathematical tool for the burgeoning field of systems biology.

Rapid rates of lineage-specific gene duplication and deletion in the alpha-globin gene family

Phylogeny reconstructions of the globin gene families have revealed that paralogous genes within species are often more similar to one another than they are to their orthologous counterparts in closely related species. This pattern has been previously attributed to mechanisms of concerted evolution such as interparalog gene conversion that homogenize sequence variation between tandemly duplicated genes and therefore create the appearance of recent common ancestry. Here we report a comparative genomic analysis of the alpha-globin gene family in mammals that reveal a surprisingly high rate of lineage-specific gene duplication and deletion via unequal crossing-over. Results of our analysis reveal that patterns of sequence similarity between paralogous alpha-like globin genes from the same species are only partly explained by concerted evolution between preexisting gene duplicates. In a number of cases, sequence similarity between paralogous sequences from the same species is attributable to recent ancestry between the products of de novo gene duplications. As a result of this surprisingly rapid rate of gene gain and loss, many mammals possess alpha-like globin genes that have no orthologous counterparts in closely related species. The resultant variation in gene copy number among species may represent an important source of regulatory variation that affects physiologically important aspects of blood oxygen transport and aerobic energy metabolism.

Excystation of Eimeria tenella sporozoites impaired by antibody recognizing gametocyte/oocyst antigens GAM22 and GAM56

Eimeria tenella is the causative agent of coccidiosis in poultry. Infection of the chicken intestine begins with ingestion of sporulated oocysts releasing sporocysts, which in turn release invasive sporozoites. The monoclonal antibody E2E5 recognizes wall-forming body type II (WFBII) in gametocytes and the WFBII-derived inner wall of oocysts. Here we describe that this antibody also binds to the stieda body of sporocysts and significantly impairs in vitro excystation of sporozoites. Using affinity chromatography and protein sequence analysis, E2E5 is shown to recognize EtGAM56, the E. tenella ortholog of the Eimeria maxima gametocyte-specific GAM56 protein. In addition, this antibody was used to screen a genomic phage display library presenting E. tenella antigens as fusion proteins with the gene VIII product on the surfaces of phagemid particles and identified the novel 22-kDa histidine- and proline-rich protein EtGAM22. The Etgam22 mRNA is expressed predominantly at the gametocyte stage, as detected by Northern blotting. Southern blot analysis in combination with data from the E. tenella genome project revealed that Etgam22 is an intronless multicopy gene, with approximately 12 to 22 copies in head-to-tail arrangement. Conspicuously, Etgam56 is also intronless and is localized adjacent to another gam56-like gene, Etgam59. Our data suggest that amplification is common for genes encoding oocyst wall proteins.

Euglena light-harvesting complexes are encoded by multifarious polyprotein mRNAs that evolve in concert

Light-harvesting complexes (LHCs) are a superfamily of chlorophyll- and carotenoid-binding proteins that are responsible for the capture of light energy and its transfer to the photosynthetic reaction centers. Unlike those of most eukaryotes, the LHCs of Euglena gracilis are translated from large mRNAs, producing polyprotein precursors consisting of multiple concatenated LHC subunits that are separated by conserved decapeptide linkers. These precursors are posttranslationally targeted to the chloroplast and cleaved into individual proteins. We analyzed expressed sequence tags from Euglena to further characterize the structural features of the LHC polyprotein-coding genes and to examine the evolution of this multigene family. Of the 19 different LHC transcriptional units we detected, 17 encoded polyproteins composed of both tandem and nontandem repeats of LHC subunits; organizations that likely occurred through unequal crossing-over. Of the 2 nonpolyprotein-encoding LHC transcripts detected, 1 evolved from the truncation of a polyprotein-coding gene. Duplication of LHC polyprotein-coding genes was particularly important in the LHCI gene family where multiple paralogous sequences were detected. Intriguingly, several of the individual LHC-coding subunits both within and between transcriptional units appeared to be evolving in concert, suggesting that gene conversion has been a significant mechanism for LHC evolution in Euglena.

Using likelihood-free inference to compare evolutionary dynamics of the protein networks of H. pylori and P. falciparum

Gene duplication with subsequent interaction divergence is one of the primary driving forces in the evolution of genetic systems. Yet little is known about the precise mechanisms and the role of duplication divergence in the evolution of protein networks from the prokaryote and eukaryote domains. We developed a novel, model-based approach for Bayesian inference on biological network data that centres on approximate Bayesian computation, or likelihood-free inference. Instead of computing the intractable likelihood of the protein network topology, our method summarizes key features of the network and, based on these, uses a MCMC algorithm to approximate the posterior distribution of the model parameters. This allowed us to reliably fit a flexible mixture model that captures hallmarks of evolution by gene duplication and subfunctionalization to protein interaction network data of Helicobacter pylori and Plasmodium falciparum. The 80% credible intervals for the duplication–divergence component are [0.64, 0.98] for H. pylori and [0.87, 0.99] for P. falciparum. The remaining parameter estimates are not inconsistent with sequence data. An extensive sensitivity analysis showed that incompleteness of PIN data does not largely affect the analysis of models of protein network evolution, and that the degree sequence alone barely captures the evolutionary footprints of protein networks relative to other statistics. Our likelihood-free inference approach enables a fully Bayesian analysis of a complex and highly stochastic system that is otherwise intractable at present. Modelling the evolutionary history of PIN data, it transpires that only the simultaneous analysis of several global aspects of protein networks enables credible and consistent inference to be made from available datasets. Our results indicate that gene duplication has played a larger part in the network evolution of the eukaryote than in the prokaryote, and suggests that single gene duplications with immediate divergence alone may explain more than 60% of biological network data in both domains.

The importance of gene duplication to biological evolution has been recognized since the 1930s. For more than a decade, substantial evidence has been collected from genomic sequence data in order to elucidate the importance and the mechanisms of gene duplication; however, most biological characteristics arise from complex interactions between the cell's numerous constituents. Recently, preliminary descriptions of the protein interaction networks have become available for species of different domains. Adapting novel techniques in stochastic simulation, the authors demonstrate that evolutionary inferences can be drawn from large-scale, incomplete network data by fitting a stochastic model of network growth that captures hallmarks of evolution by duplication and divergence. They have also analyzed the effect of summarizing protein networks in different ways, and show that a reliable and consistent analysis requires many aspects of network data to be considered jointly; in contrast to what is commonly done in practice. Their results indicate that duplication and divergence has played a larger role in the network evolution of the eukaryote P. falciparum than in the prokaryote H. pylori, and emphasize at least for the eukaryote the potential importance of subfunctionalization in network evolution.

Unique nucleotide polymorphism of ankyrin gene cluster in Arabidopsis

The ankyrin (ANK) gene cluster is a part of a multigene family encoding ANK transmembrane proteins in Arabidopsis thaliana, and plays an important role in protein-protein interactions and in signal pathways. In contrast to other regions of a genome, the ANK gene cluster exhibits an extremely high level of DNA polymorphism in an approximately 5-kb region, without apparent decay. Phylogenetic analysis detects two clear, deeply differentiated haplotypes (dimorphism). The divergence between haplotypes of accession Col-0 and Ler-0 (Hap-C and Hap-L) is estimated to be 10.7%, approximately equal to the 10.5% average divergence between A. thaliana and A. lyrata. Sequence comparisons for the ANK gene cluster homologues in Col-0 indicate that the members evolve independently, and that the similarity among paralogues is lower than between alleles. Very little intralocus recombination or gene conversion is detected in ANK regions. All these characteristics of the ANK gene cluster are consistent with a tandem gene duplication and birth-and-death process. The possible mechanisms for and implications of this elevated nucleotide variation are also discussed, including the suggestion of balancing selection.

Spliceosomal small nuclear RNA genes in 11 insect genomes

The removal of introns from the primary transcripts of protein-coding genes is accomplished by the spliceosome, a large macromolecular complex of which small nuclear RNAs (snRNAs) are crucial components. Following the recent sequencing of the honeybee (Apis mellifera) genome, we used various computational methods, ranging from sequence similarity search to RNA secondary structure prediction, to search for putative snRNA genes (including their promoters) and to examine their pattern of conservation among 11 available insect genomes (A. mellifera, Tribolium castaneum, Bombyx mori, Anopheles gambiae, Aedes aegypti, and six Drosophila species). We identified candidates for all nine spliceosomal snRNA genes in all the analyzed genomes. All the species contain a similar number of snRNA genes, with the exception of A. aegypti, whose genome contains more U1, U2, and U5 genes, and A. mellifera, whose genome contains fewer U2 and U5 genes. We found that snRNA genes are generally more closely related to homologs within the same genus than to those in other genera. Promoter regions for all spliceosomal snRNA genes within each insect species share similar sequence motifs that are likely to correspond to the PSEA (proximal sequence element A), the binding site for snRNA activating protein complex, but these promoter elements vary in sequence among the five insect families surveyed here. In contrast to the other insect species investigated, Dipteran genomes are characterized by a rapid evolution (or loss) of components of the U12 spliceosome and a striking loss of U12-type introns.

Alternative starts of transcription, several paralogues, and almost-fixed interspecific differences of the gene fruitless in a hemimetabolous insect

In fruit flies, the gene fruitless (fru) governs the establishment of the potential for male sexual behavior. We partially cloned fru from a hemimetabolous insect for the first time and we compared fru among three closely related and acoustically communicating grasshopper species: Chorthippus biguttulus, C. brunneus, and C. mollis. The fru of grasshoppers is organized similarly to fru of holometabolous insects, with a BTB and Zn-finger domains separated by a nonconserved repetitive linker. As in Drosophila, several transcripts of fru are found in grasshoppers. We also present evidence for the coexistence of several copies of fru in the grasshopper genome. Within species these copies are almost identical and carry almost-fixed species-specific differences. This suggests that the paralogous copies of fru in grasshoppers do not evolve independently from each other.