Sequential phosphorylation of CST subunits by different cyclin-Cdk1 complexes orchestrate telomere replication

Telomeres are nucleoprotein structures that cap the ends of linear chromosomes. Telomere homeostasis is central to maintaining genomic integrity. In budding yeast, Cdk1 phosphorylates the telomere-specific binding protein, Cdc13, promoting the recruitment of telomerase to telomere and thereby telomere elongation. Cdc13 is also an integral part of the CST (Cdc13-Stn1-Ten1) complex that is essential for telomere capping and counteracting telomerase-dependent telomere elongation. Therefore, telomere length homeostasis is a balance between telomerase-extendable and CST-unextendable states. In our earlier work, we showed that Cdk1 also phosphorylates Stn1 which occurs sequentially following Cdc13 phosphorylation during cell cycle progression. This stabilizes the CST complex at the telomere and results in telomerase inhibition. Hence Cdk1-dependent phosphorylations of Stn1 acts like a molecular switch that drives Cdc13 to complex with Stn1-Ten1 rather than with telomerase. However, the underlying mechanism of how a single cyclin-dependent kinase phosphorylates Cdc13 and Stn1 in temporally distinct windows is largely unclear. Here, we show that S phase cyclins are necessary for telomere maintenance. The S phase and mitotic cyclins facilitate Cdc13 and Stn1 phosphorylation respectively, to exert opposing outcomes at the telomere. Thus, our results highlight a previously unappreciated role for cyclins in telomere replication.     

The evolution of disease resistance genes

Several common themes have shaped the evolution of plant disease resistance genes. These include duplication events of progenitor resistance genes and further expansion to create clustered gene families. Variation can arise from both intragenic and intergenic recombination and gene conversion. Recombination has also been implicated in the generation of novel resistance specificities. Resistance gene clusters appear to evolve more rapidly than other regions of the genome. In addition, domains believed to be involved in recognitional specificity, such as the leucine-rich repeat (LRR), are subject to adaptive selection. Transposable elements have been associated with some resistance gene clusters, and may generate further variation at these complexes.     

Lack of age-associated telomere shortening in long- and short-lived species of sea urchins

The red sea urchin, Strongylocentrotus franciscanus, can live in excess of 100 years while the sea urchin Lytechinus variegatus has an estimated lifespan of only 3-4 years. In an effort to understand the molecular mechanism underlying the difference in their longevity we characterized telomere biology in these species of sea urchins. Telomerase activity was found throughout early stages of development in L. variegatus and is maintained in adult tissues of L. variegatus and S. franciscanus. Terminal restriction fragment analysis indicated a lack of age-associated telomere shortening. These data suggest that long- and short-lived sea urchins do not utilize telomerase repression as a mechanism to suppress neoplastic transformation.     

Conservation and divergence of APETALA1/FRUITFULL-like gene function in grasses: evidence from gene expression analyses

Duplicated APETALA1/FRUITFULL (AP1/FUL) genes show distinct but overlapping patterns of expression within rice (Oryza sativa) and within ryegrass (Lolium temulentum), suggesting discrete functional roles in the transition to flowering, specification of spikelet meristem identity, and specification of floral organ identity. In this study, we analyzed the expression of the AP1/FUL paralogues FUL1 and FUL2 across phylogenetically disparate grasses to test hypotheses of gene function. In combination with other studies, our data support similar roles for both genes in spikelet meristem identity, a general role for FUL1 in floral organ identity, and a more specific role for FUL2 in outer floral whorl identity. In contrast to Arabidopsis AP1/FUL genes, expression of FUL1 and FUL2 is consistent with an early role in the transition to flowering. In general, FUL1 has a wider expression pattern in all spikelet organs than FUL2, but both genes are expressed in all spikelet organs in some cereals. FUL1 and FUL2 appear to have multiple redundant functions in early inflorescence development. We hypothesize that sub-functionalization of FUL2 and interaction of FUL2 with LHS1 could specify lemma and palea identity in the grass floret.     

Molecular evolution and population genetics of duplicated accessory gland protein genes in Drosophila

To investigate the potential importance of gene duplication in D. melanogaster accessory gland protein (Acp) gene evolution we carried out a computational analysis comparing annotated D. melanogaster Acp genes to the entire D. melanogaster genome. We found that two known Acp genes are actually members of small multigene families. Polymorphism and divergence data from these duplicated genes suggest that in at least four cases, protein divergence between D. melanogaster and D. simulans is a result of directional selection. One putative Acp revealed by our computational analysis shows evidence of a recent selective sweep in a non-African population (but not in an African population). These data support the idea that selection on reproduction-related genes may drive divergence of populations within species, and strengthen the conclusion that Acps may often be under directional selection in Drosophila.     

Evaluation of whether accelerated protein evolution in chordates has occurred before, after, or simultaneously with gene duplication

Gene duplication and loss are predicted to be at least of the order of the substitution rate and are key contributors to the development of novel gene function and overall genome evolution. Although it has been established that proteins evolve more rapidly after gene duplication, we were interested in testing to what extent this reflects causation or association. Therefore, we investigated the rate of evolution prior to gene duplication in chordates. Two patterns emerged; firstly, branches, which are both preceded by a duplication and followed by a duplication, display an elevated rate of amino acid replacement. This is reflected in the ratio of nonsynonymous to synonymous substitution (mean nonsynonymous to synonymous nucleotide substitution rate ratio [Ka:Ks]) of 0.44 compared with branches preceded by and followed by a speciation (mean Ka:Ks of 0.23). The observed patterns suggest that there can be simultaneous alteration in the selection pressures on both gene duplication and amino acid replacement, which may be consistent with co-occurring increases in positive selection, or alternatively with concurrent relaxation of purifying selection. The pattern is largely, but perhaps not completely, explained by the existence of certain families that have elevated rates of both gene duplication and amino acid replacement. Secondly, we observed accelerated amino acid replacement prior to duplication (mean Ka:Ks for postspeciation preduplication branches was 0.27). In some cases, this could reflect adaptive changes in protein function precipitating a gene duplication event. In conclusion, the circumstances surrounding the birth of new proteins may frequently involve a simultaneous change in selection pressures on both gene-copy number and amino acid replacement. More precise modeling of the relative importance of preduplication, postduplication, and simultaneous amino acid replacement will require larger and denser genomic data sets from multiple species, allowing simultaneous estimation of lineage-specific fluctuations in mutation rates and adaptive constraints.     

The complete chloroplast genome sequence of Pelargonium x hortorum: organization and evolution of the largest and most highly rearranged chloroplast genome of land plants

The chloroplast genome of Pelargonium x hortorum has been completely sequenced. It maps as a circular molecule of 217,942 bp and is both the largest and most rearranged land plant chloroplast genome yet sequenced. It features 2 copies of a greatly expanded inverted repeat (IR) of 75,741 bp each and, consequently, diminished single-copy regions of 59,710 and 6,750 bp. Despite the increase in size and complexity of the genome, the gene content is similar to that of other angiosperms, with the exceptions of a large number of pseudogenes, the recognition of 2 open reading frames (ORF56 and ORF42) in the trnA intron with similarities to previously identified mitochondrial products (ACRS and pvs-trnA), the losses of accD and trnT-ggu and, in particular, the presence of a highly divergent set of rpoA-like ORFs rather than a single, easily recognized gene for rpoA. The 3-fold expansion of the IR (relative to most angiosperms) accounts for most of the size increase of the genome, but an additional 10% of the size increase is related to the large number of repeats found. The Pelargonium genome contains 35 times as many 31 bp or larger repeats than the unrearranged genome of Spinacia. Most of these repeats occur near the rearrangement hotspots, and 2 different associations of repeats are localized in these regions. These associations are characterized by full or partial duplications of several genes, most of which appear to be nonfunctional copies or pseudogenes. These duplications may also be linked to the disruption of at least 1 but possibly 2 or 3 operons. We propose simple models that account for the major rearrangements with a minimum of 8 IR boundary changes and 12 inversions in addition to several insertions of duplicated sequence.     

The evolution of protein sequences by repetitious gene duplication: Clostridial flavodoxin

Summary Internal regularities of amino acid sequences of flavodoxins, FMN-containing, low molecular weight flavoproteins, were statistically examined using the minimum mutation method. The sequence ofClostridium pasteurianum flavodoxin shows statistically significant evidence of repetitious internal gene duplications at different levels of structure. Peptide pairs with a low chance probability of occurrence were frequently observed at a shift of 5 residues. The pairs with the lowest chance probabilites are a pair of heptapeptides at positions 39–45 vs. 44–50, a 5 residue shift (p = 9 × 10^–6). Most of the related pairs are consistent and could best be explained by the repeating pentapeptide sequence: (Lys-Gly-Ala-Asp-Val-)_n and appropriate gaps. Internal repetitions with longer shifts were also suggested for other flavodoxins. Repetitious gene duplication is proposed for the early stages of flavodoxin evolution.     

The evolution of mitochondrial genomes in modern frogs (Neobatrachia): nonadaptive evolution of mitochondrial genome reorganization

Background

Although mitochondrial (mt) gene order is highly conserved among vertebrates, widespread gene rearrangements occur in anurans, especially in neobatrachians. Protein coding genes in the mitogenome experience adaptive or purifying selection, yet the role that selection plays on genomic reorganization remains unclear. We sequence the mitogenomes of three species of Glandirana and hot spots of gene rearrangements of 20 frog species to investigate the diversity of mitogenomic reorganization in the Neobatrachia. By combing these data with other mitogenomes in GenBank, we evaluate if selective pressures or functional constraints act on mitogenomic reorganization in the Neobatrachia. We also look for correlations between tRNA positions and codon usage.

Results

Gene organization in Glandirana was typical of neobatrachian mitogenomes except for the presence of pseudogene trnS (AGY). Surveyed ranids largely exhibited gene arrangements typical of neobatrachian mtDNA although some gene rearrangements occurred. The correlation between codon usage and tRNA positions in neobatrachians was weak, and did not increase after identifying recurrent rearrangements as revealed by basal neobatrachians. Codon usage and tRNA positions were not significantly correlated when considering tRNA gene duplications or losses. Change in number of tRNA gene copies, which was driven by genomic reorganization, did not influence codon usage bias. Nucleotide substitution rates and d_N/d_S ratios were higher in neobatrachian mitogenomes than in archaeobatrachians, but the rates of mitogenomic reorganization and mt nucleotide diversity were not significantly correlated.

Conclusions

No evidence suggests that adaptive selection drove the reorganization of neobatrachian mitogenomes. In contrast, protein-coding genes that function in metabolism showed evidence for purifying selection, and some functional constraints appear to act on the organization of rRNA and tRNA genes. As important nonadaptive forces, genetic drift and mutation pressure may drive the fixation and evolution of mitogenomic reorganizations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-691) contains supplementary material, which is available to authorized users.     

Simulating evolution by gene duplication of protein features that require multiple amino acid residues

Gene duplication is thought to be a major source of evolutionary innovation because it allows one copy of a gene to mutate and explore genetic space while the other copy continues to fulfill the original function. Models of the process often implicitly assume that a single mutation to the duplicated gene can confer a new selectable property. Yet some protein features, such as disulfide bonds or ligand binding sites, require the participation of two or more amino acid residues, which could require several mutations. Here we model the evolution of such protein features by what we consider to be the conceptually simplest route-point mutation in duplicated genes. We show that for very large population sizes N, where at steady state in the absence of selection the population would be expected to contain one or more duplicated alleles coding for the feature, the time to fixation in the population hovers near the inverse of the point mutation rate, and varies sluggishly with the lambda(th) root of 1/N, where lambda is the number of nucleotide positions that must be mutated to produce the feature. At smaller population sizes, the time to fixation varies linearly with 1/N and exceeds the inverse of the point mutation rate. We conclude that, in general, to be fixed in 10(8) generations, the production of novel protein features that require the participation of two or more amino acid residues simply by multiple point mutations in duplicated genes would entail population sizes of no less than 10(9).     

Organization and structural evolution of four multigene families in Arabidopsis thaliana: AtLCAD,AtLGT, AtMYST and AtHD-GL2

The Arabidopsis Genome Initiative has released up to now more than 80% of the genome sequence of Arabidopsis thaliana. About 70% of the identified genes have at least one paralogue. In order to understand the biological function of individual genes, it is essential to study the structure, expression and organization of the entire multigene family. A systematic analysis of multigene families, made possible by the amount of genomic sequence data available, provides important clues for the understanding of genome evolution and plasticity. In this paper, four multigene families of A. thaliana are characterized, namely LCAD, HD-GL2, LGT and MYST. Members of HD-GL2 and LCAD have already been reported in plants. The LGT genes specify proteins containing motifs of glycosyl transferase. No plant genes similar to the LGT genes have been reported to date. The novel MYST family, most likely plant-specific, encodes proteins with no identified function. Sequencing and in silico analysis led to the characterization of 29 novel genes belonging to these four gene families. The organization, structure and evolution of all the members of the four families are discussed, as well as their chromosome location. Expression data of some of the paralogues of each family are also presented.     

The interface of protein structure, protein biophysics, and molecular evolution

Abstract The interface of protein structural biology, protein biophysics, molecular evolution, and molecular population genetics forms the foundations for a mechanistic understanding of many aspects of protein biochemistry. Current efforts in interdisciplinary protein modeling are in their infancy and the state-of-the art of such models is described. Beyond the relationship between amino acid substitution and static protein structure, protein function, and corresponding organismal fitness, other considerations are also discussed. More complex mutational processes such as insertion and deletion and domain rearrangements and even circular permutations should be evaluated. The role of intrinsically disordered proteins is still controversial, but may be increasingly important to consider. Protein geometry and protein dynamics as a deviation from static considerations of protein structure are also important. Protein expression level is known to be a major determinant of evolutionary rate and several considerations including selection at the mRNA level and the role of interaction specificity are discussed. Lastly, the relationship between modeling and needed high-throughput experimental data as well as experimental examination of protein evolution using ancestral sequence resurrection and in vitro biochemistry are presented, towards an aim of ultimately generating better models for biological inference and prediction.     

Divergence of duplicated genes in maize: evolution of contrasting targeting information for enzymes in the porphyrin pathway

The divergence of sequence and expression pattern of duplicated genes provides a means for genetic innovation to occur without sacrificing an essential function. The cpx1 and cpx2 genes of maize are a singular example of duplicated genes that have diverged by deletion and creation of protein targeting information. The cpx genes encode coproporphyrinogen III oxidase ('coprogen oxidase'), which catalyzes a step in the synthesis of chlorophyll and heme. In plants, this enzyme has been found exclusively in the plastids. The cpx1 and cpx2 genes encode almost identical, catalytically active enzymes with distinctive N-terminal peptide sequences. The cpx1 gene encodes the expected plastid transit peptide, but this region is deleted from the cpx2 gene. While the 5' regions of both messenger RNAs are highly similar, the cpx2 gene has an open-reading frame that could encode a new targeting signal. GFP fused with CPX1 localized to the plastids. In contrast, the GFP fusion with CPX2 did not target plastids and appeared to localize to mitochondria. Both cpx genes are expressed ubiquitously but, based on mutant phenotype, they seem to have discrete biological roles. Seedlings homozygous for a null mutation in the cpx1 gene completely lack chlorophyll and develop necrotic lesions in the light. However, the mutant seedlings and callus cultures will grow in tissue culture in the dark, implying that they retain a capacity to produce heme. We discuss models for the evolution of the cpx genes and possible roles of mitochondrion-localized coprogen oxidase activity in maize.     

Comparative methods for the analysis of gene-expression evolution: an example using yeast functional genomic data

Understanding the evolution of gene function is a primary challenge of modern evolutionary biology. Despite an expanding database from genomic and developmental studies, we are lacking quantitative methods for analyzing the evolution of some important measures of gene function, such as gene-expression patterns. Here, we introduce phylogenetic comparative methods to compare different models of gene-expression evolution in a maximum-likelihood framework. We find that expression of duplicated genes has evolved according to a nonphylogenetic model, where closely related genes are no more likely than more distantly related genes to share common expression patterns. These results are consistent with previous studies that found rapid evolution of gene expression during the history of yeast. The comparative methods presented here are general enough to test a wide range of evolutionary hypotheses using genomic-scale data from any organism.