Transitions between phases of genomic differentiation during stick-insect speciation.

Speciation can involve a transition from a few genetic loci that are resistant to gene flow to genome-wide differentiation. However, only limited data exist concerning this transition and the factors promoting it. Here, we study phases of speciation using data from >100 populations of 11 species of Timema stick insects. Consistent with early phases of genic speciation, adaptive colour-pattern loci reside in localized genetic regions of accentuated differentiation between populations experiencing gene flow. Transitions to genome-wide differentiation are also observed with gene flow, in association with differentiation in polygenic chemical traits affecting mate choice. Thus, intermediate phases of speciation are associated with genome-wide differentiation and mate choice, but not growth of a few genomic islands. We also find a gap in genomic differentiation between sympatric taxa that still exchange genes and those that do not, highlighting the association between differentiation and complete reproductive isolation. Our results suggest that substantial progress towards speciation may involve the alignment of multi-faceted aspects of differentiation.

Myristoylated CIL-7 regulates ciliary extracellular vesicle biogenesis.

The cilium both releases and binds to extracellular vesicles (EVs). EVs may be used by cells as a form of intercellular communication and mediate a broad range of physiological and pathological processes. The mammalian polycystins (PCs) localize to cilia, as well as to urinary EVs released from renal epithelial cells. PC ciliary trafficking defects may be an underlying cause of autosomal dominant polycystic kidney disease (PKD), and ciliary-EV interactions have been proposed to play a central role in the biology of PKD. In Caenorhabditis elegans and mammals, PC1 and PC2 act in the same genetic pathway, act in a sensory capacity, localize to cilia, and are contained in secreted EVs, suggesting ancient conservation. However, the relationship between cilia and EVs and the mechanisms generating PC-containing EVs remain an enigma. In a forward genetic screen for regulators of C. elegans PKD-2 ciliary localization, we identified CIL-7, a myristoylated protein that regulates EV biogenesis. Loss of CIL-7 results in male mating behavioral defects, excessive accumulation of EVs in the lumen of the cephalic sensory organ, and failure to release PKD-2::GFP-containing EVs to the environment. Fatty acylation, such as myristoylation and palmitoylation, targets proteins to cilia and flagella. The CIL-7 myristoylation motif is essential for CIL-7 function and for targeting CIL-7 to EVs. C. elegans is a powerful model with which to study ciliary EV biogenesis in vivo and identify cis-targeting motifs such as myristoylation that are necessary for EV-cargo association and function.

The role of DNA insertions in phenotypic differentiation between humans and other primates.

What makes us human is one of the most interesting and enduring questions in evolutionary biology. To assist in answering this question, we have identified insertions in the human genome which cannot be found in five comparison primate species: Chimpanzee, gorilla, orangutan, gibbon, and macaque. A total of 21,269 nonpolymorphic human-specific insertions were identified, of which only 372 were found in exons. Any function conferred by the remaining 20,897 is likely to be regulatory. Many of these insertions are likely to have been fitness neutral; however, a small number has been identified in genes showing signs of positive selection. Insertions found within positively selected genes show associations to neural phenotypes, which were also enriched in the whole data set. Other phenotypes that are found to be enriched in the data set include dental and sensory perception-related phenotypes, features which are known to differ between humans and other apes. The analysis provides several likely candidates, either genes or regulatory regions, which may be involved in the processes that differentiate humans from other apes.

Transposable element invasions.

Transposable elements have an ongoing, largely parasitic interaction with their hosts. We are interested in the timescale of this interaction. In a recent publication, we have examined the sequence divergence between class II DNA transposons from mammalian genomes. We asked whether these sequences undergo a continuing process of turnover, keeping a family as an integrated whole, as members of the family are continually created and lost. Alternatively, we envisaged that elements might have been involved in a burst of amplification, soon after they first occupied a mammalian genome, and the shared ancestry of present-day elements harks back to this initial amplification, a process that we termed a "life cycle." We resolved between these processes by estimating the time to common ancestry predicted from the genetic diversity of sequences found in a transposon family, and also estimating, from the mammalian orders that currently possess copies of the family, the time when the family first entered the mammalian genome. These times are approximately the same, supporting the "life cycle" model. This casts light on how far we can infer genetic changes in the past through the study of DNA sequences from the present.

Alu elements in primates are preferentially lost from areas of high GC content.

The currently-accepted dogma when analysing human Alu transposable elements is that 'young' Alu elements are found in low GC regions and 'old' Alus in high GC regions. The correlation between high GC regions and high gene frequency regions make this observation particularly difficult to explain. Although a number of studies have tackled the problem, no analysis has definitively explained the reason for this trend. These observations have been made by relying on the subfamily as a proxy for age of an element. In this study, we suggest that this is a misleading assumption and instead analyse the relationship between the taxonomic distribution of an individual element and its surrounding GC environment. An analysis of 103906 Alu elements across 6 human chromosomes was carried out, using the presence of orthologous Alu elements in other primate species as a proxy for age. We show that the previously-reported effect of GC content correlating with subfamily age is not reflected by the ages of the individual elements. Instead, elements are preferentially lost from areas of high GC content over time. The correlation between GC content and subfamily may be due to a change in insertion bias in the young subfamilies. The link between Alu subfamily age and GC region was made due to an over-simplification of the data and is incorrect. We suggest that use of subfamilies as a proxy for age is inappropriate and that the analysis of ortholog presence in other primate species provides a deeper insight into the data.

The diversity of class II transposable elements in mammalian genomes has arisen from ancestral phylogenetic splits during ancient waves of proliferation through the genome.

DNA transposons make up 3% of the human genome, approximately the same percentage as genes. However, because of their inactivity, they are often ignored in favor of the more abundant, active, retroelements. Despite this relative ignominy, there are a number of interesting questions to be asked of these transposon families. One particular question relates to the timing of proliferation and inactivation of elements in a family. Does an ongoing process of turnover occur, or is the process more akin to a life cycle for the family, with elements proliferating rapidly before deactivation at a later date? We answer this question by tracing back to the most recent common ancestor (MRCA) of each modern transposon family, using two different methods. The first method identifies the MRCA of the species in which a family of transposon fossils can still be found, which we assume will have existed soon after the true origin date of the transposon family. The second method uses molecular dating techniques to predict the age of the MRCA element from which all elements found in a modern genome are descended. Independent data from five pairs of species are used in the molecular dating analysis: human-chimpanzee, human-orangutan, dog-panda, dog-cat, and cow-pig. Orthologous pairs of elements from host species pairs are included, and the divergence dates of these species are used to constrain the analysis. We discover that, in general, the times to element common ancestry for a given family are the same for the different species pairs, suggesting that there has been no order-specific process of turnover. Furthermore, for most families, the ages of the common ancestor of the host species and of that of the elements are similar, suggesting a life cycle model for the proliferation of transposons. Where these two ages differ, in families found only in Primates and Rodentia, for example, we find that the host species date is later than that of the common ancestor of the elements, implying that there may be large deletions of elements from host species, examples of which were found in their ancestors.

Investigation of the origin and spread of a Mammalian transposable element based on current sequence diversity.

Almost half the human genome consists of mobile DNA elements, and their analysis is a vital part of understanding the human genome as a whole. Many of these elements are ancient and have persisted in the genome for tens or hundreds of millions of years, providing a window into the evolution of modern mammals. The Golem family have been used as model transposons to highlight computational analyses which can be used to investigate these elements, particularly the use of molecular dating with large transposon families. Whole-genome searches found Golem sequences in 20 mammalian species. Golem A and B subsequences were only found in primates and squirrel. Interestingly, the full-length Golem, found as a few copies in many mammalian genomes, was found abundantly in horse. A phylogenetic profile suggested that Golem originated after the eutherian-metatherian divergence and that the A and B subfamilies originated at a much later date. Molecular dating based on sequence diversity suggests an early age, of 175 Mya, for the origin of the family and that the A and B lineages originated much earlier than expected from their current taxonomic distribution and have subsequently been lost in some lineages. Using publically available data, it is possible to investigate the evolutionary history of transposon families. Determining in which organisms a transposon can be found is often used to date the origin and expansion of the families. However, in this analysis, molecular dating, commonly used for determining the age of gene sequences, has been used, reducing the likelihood of errors from deleted lineages.

Evaluation of a prediction protocol to identify potential targets of epigenetic reprogramming by the cancer associated Epstein Barr virus.

BACKGROUND: Epstein Barr virus (EBV) infects the majority of the human population, causing fatal diseases in a small proportion in conjunction with environmental factors. Following primary infection, EBV remains latent in the memory B cell population for life. Recurrent reactivation of the virus occurs, probably due to activation of the memory B-lymphocytes, resulting in viral replication and re-infection of B-lymphocytes. Methylation of the viral DNA at CpG motifs leads to silencing of viral gene expression during latency. Zta, the key viral protein that mediates the latency/reactivation balance, interacts with methylated DNA. Zta is a transcription factor for both viral and host genes. A sub-set of its DNA binding sites (ZREs) contains a CpG motif, which is recognised in its methylated form. Detailed analysis of the promoter of the viral gene BRLF1 revealed that interaction with a methylated CpG ZRE (RpZRE3) is key to overturning the epigenetic silencing of the gene. METHODOLOGY AND PRINCIPAL FINDINGS: Here we question whether we can use this information to identify which host genes contain promoters with similar response elements. A computational search of human gene promoters identified 274 targets containing the 7-nucleotide RpZRE3 core element. DNA binding analysis of Zta with 17 of these targets revealed that the flanking context of the core element does not have a profound effect on the ability of Zta to interact with the methylated sites. A second juxtaposed ZRE was observed for one promoter. Zta was able to interact with this site, although co-occupancy with the RpZRE3 core element was not observed. CONCLUSIONS/SIGNIFICANCE: This research demonstrates 274 human promoters have the potential to be regulated by Zta to overturn epigenetic silencing of gene expression during viral reactivation from latency.

Molecular characterisation of the SAND protein family: a study based on comparative genomics, structural bioinformatics and phylogeny.

The activities of vertebrate lysosomes are critical to many essential cellular processes. The yeast vacuole is analogous to the mammalian lysosome and is used as a tool to gain insights into vesicle mediated vacuolar/lysosome transport. The protein SAND, which does not contain a SAND domain (PFAM accession number PF01342), has recently been shown to function at the tethering/docking stage of vacuole fusion as a critical component of the vacuole SNARE complex. In this publication we have identified SAND in diverse eukaryotes, from single celled organisms such as the yeasts to complex multi-cellular chordates such as mammals. We have demonstrated subfamily divisions in the SAND proteins and show that in vertebrates, a duplication event gave rise to two SAND sequences. This duplication appears to have occurred during early vertebrate evolution and conceivably with the evolution of lysosomes. Using bioinformatics we predict a secondary structure, solvent accessibility profile and protein fold for the SAND proteins and determine conserved sequence motifs, present in all SAND proteins and those that are specific to subsets. A comprehensive evaluation of yeast and human functional studies in conjunction with our in silico analysis has identified potential roles for some of these motifs.