Integration of multiple repeats of geminiviral DNA into the nuclear genome of tobacco during evolution.

Integration of viral DNA into the host nuclear genome, although not unusual in bacterial and animal systems, has surprisingly not been reported for plants. We have discovered geminvirus-related DNA (GRD) sequences, in the form of distinct sets of multiple direct repeats comprising three related repeat classes, situated in a unique locus in the Nicotiana tabacum (tobacco) nuclear genome. The organization of these sequences is similar or identical in eight different tobacco cultivars we have examined. DNA sequence analysis reveals that each repeat has sequences most resembling those of the New World geminiviral DNA replication origin plus the adjacent AL1 gene, encoding the viral replication protein. We believe these GRD sequences originated quite recently in Nicotiana evolution through integration of geminiviral DNA by some combination of the processes of illegitimate recombination, amplification, deletions, and rearrangements. These events must have occurred in plant tissue that was subsequently able to contribute to meristematic tissue yielding gametes. GRD may have been retained in tobacco by selection or by random fixation in a small evolving population. Although we cannot detect transcription of these sequences, this does not exclude the possibility that they may originally have been expressed.

Sequence similarity analysis of Escherichia coli proteins: functional and evolutionary implications.

A computer analysis of 2328 protein sequences comprising about 60% of the Escherichia coli gene products was performed using methods for database screening with individual sequences and alignment blocks. A high fraction of E. coli proteins--86%--shows significant sequence similarity to other proteins in current databases; about 70% show conservation at least at the level of distantly related bacteria, and about 40% contain ancient conserved regions (ACRs) shared with eukaryotic or Archaeal proteins. For > 90% of the E. coli proteins, either functional information or sequence similarity, or both, are available. Forty-six percent of the E. coli proteins belong to 299 clusters of paralogs (intraspecies homologs) defined on the basis of pairwise similarity. Another 10% could be included in 70 superclusters using motif detection methods. The majority of the clusters contain only two to four members. In contrast, nearly 25% of all E. coli proteins belong to the four largest superclusters--namely, permeases, ATPases and GTPases with the conserved "Walker-type" motif, helix-turn-helix regulatory proteins, and NAD(FAD)-binding proteins. We conclude that bacterial protein sequences generally are highly conserved in evolution, with about 50% of all ACR-containing protein families represented among the E. coli gene products. With the current sequence databases and methods of their screening, computer analysis yields useful information on the functions and evolutionary relationships of the vast majority of genes in a bacterial genome. Sequence similarity with E. coli proteins allows the prediction of functions for a number of important eukaryotic genes, including several whose products are implicated in human diseases.

One sequence, two folds: a metastable structure of CD2.

When expressed as part of a glutathione S-transferase fusion protein the NH2-terminal domain of the lymphocyte cell adhesion molecule CD2 is shown to adopt two different folds. The immunoglobulin superfamily structure of the major (85%) monomeric component has previously been determined by both x-ray crystallography and NMR spectroscopy. We now describe the structure of a second, dimeric, form present in about 15% of recombinant CD2 molecules. After denaturation and refolding in the absence of the fusion partner, dimeric CD2 is converted to monomer, illustrating that the dimeric form represents a metastable folded state. The crystal structure of this dimeric form, refined to 2.0-A resolution, reveals two domains with overall similarity to the IgSF fold found in the monomer. However, in the dimer each domain is formed by the intercalation of two polypeptide chains. Hence each domain represents a distinct folding unit that can assemble in two different ways. In the dimer the two domains fold around a hydrophilic interface believed to mimic the cell adhesion interaction at the cell surface, and the formation of dimer can be regulated by mutating single residues at this interface. This unusual misfolded form of the protein, which appears to result from inter- rather than intramolecular interactions being favored by an intermediate structure formed during the folding process, illustrates that evolution of protein oligomers is possible from the sequence for a single protein domain.

Documentation of reticulate evolution in peonies (Paeonia) using internal transcribed spacer sequences of nuclear ribosomal DNA: implications for biogeography and concerted evolution.

The internal transcribed spacers (ITS) of nuclear ribosomal DNA of 33 species of genus Paeonia (Paeoniaceae) were sequenced. In section Paeonia, different patterns of nucleotide additivity were detected in 14 diploid and tetraploid species at sites that are variable in the other 12 species of the section, suggesting that reticulate evolution has occurred. Phylogenetic relationships of species that do not show additivity, and thus ostensibly were not derived through hybridization, were reconstructed by parsimony analysis. The taxa presumably derived through reticulate evolution were then added to the phylogenetic tree according to additivity from putative parents. The study provides an example of successfully using ITS sequences to reconstruct reticulate evolution in plants and further demonstrates that the sequence data could be highly informative and accurate for detecting hybridization. Maintenance of parental sequences in the species of hybrid origin is likely due to slowing of concerted evolution caused by the long generation time of peonies. The partial and uneven homogenization of parental sequences displayed in nine species of putative hybrid origin may have resulted from gradients of gene conversion. The documented hybridizations may have occurred since the Pleistocene glaciations. The species of hybrid origin and their putative parents are now distantly allopatric. Reconstruction of reticulate evolution with sequence data, therefore, provides gene records for distributional histories of some of the parental species.

Gaia-ESO Survey: Analysis of pre-main sequence stellar spectra

Context. The Gaia-ESO Public Spectroscopic Survey is obtaining high-quality spectroscopy of some 100 000 Milky Way stars using the FLAMES spectrograph at the VLT, down to V = 19 mag, systematically covering all the main components of the Milky Way and providing the first homogeneous overview of the distributions of kinematics and chemical element abundances in the Galaxy. Observations of young open clusters, in particular, are giving new insights into their initial structure, kinematics, and their subsequent evolution. Aims. This paper describes the analysis of UVES and GIRAFFE spectra acquired in the fields of young clusters whose population includes pre-main sequence (PMS) stars. The analysis is applied to all stars in such fields, regardless of any prior information on membership, and provides fundamental stellar atmospheric parameters, elemental abundances, and PMS-specific parameters such as veiling, accretion, and chromospheric activity. Methods. When feasible, different methods were used to derive raw parameters (e.g. line equivalent widths) fundamental atmospheric parameters and derived parameters (e.g. abundances). To derive some of these parameters, we used methods that have been extensively used in the past and new ones developed in the context of the Gaia-ESO survey enterprise. The internal precision of these quantities was estimated by inter-comparing the results obtained by these different methods, while the accuracy was estimated by comparison with independent external data, such as effective temperature and surface gravity derived from angular diameter measurements, on a sample of benchmarks stars. A validation procedure based on these comparisons was applied to discard spurious or doubtful results and produce recommended parameters. Specific strategies were implemented to resolve problems of fast rotation, accretion signatures, chromospheric activity, and veiling. Results. The analysis carried out on spectra acquired in young cluster fields during the first 18 months of observations, up to June 2013, is presented in preparation of the first release of advanced data products. These include targets in the fields of the ρ Oph, Cha I, NGC 2264, γ Vel, and NGC 2547 clusters. Stellar parameters obtained with the higher resolution and larger wavelength coverage from UVES are reproduced with comparable accuracy and precision using the smaller wavelength range and lower resolution of the GIRAFFE setup adopted for young stars, which allows us to provide stellar parameters with confidence for the much larger GIRAFFE sample. Precisions are estimated to be ≈120 K rms in Teff, ≈0.3 dex rms in log g, and ≈0.15 dex rms in [Fe/H] for the UVES and GIRAFFE setups.

On the formation and evolution of the first Be star in a black hole binary MWC 656

We find that the formation of MWC 656 (the first Be binary containing a black hole) involves a common envelope phase and a supernova explosion. This result supports the idea that a rapidly rotating Be star can emerge out of a common envelope phase, which is very intriguing because this evolutionary stage is thought to be too fast to lead to significant accretion and spin up of the B star. We predict ∼10–100 of B-BH binaries to currently reside in the Galactic disc, among which around 1/3 contain a Be star, but there is only a small chance to observe a system with parameters resembling MWC 656. If MWC 656 is representative of intrinsic Galactic Be-BH binary population, it may indicate that standard evolutionary theory needs to be revised. This would pose another evolutionary problem in understanding black hole (BH) binaries, with BH X-ray novae formation issue being the prime example. Future evolution of MWC 656 with an ∼5 M⊙ BH and with an ∼13 M⊙ main-sequence companion on an ∼60 d orbit may lead to the formation of a coalescing BH–NS (neutron star) system. The estimated Advanced LIGO/Virgo detection rate of such systems is up to ∼0.2 yr−1. This empirical estimate is a lower limit as it is obtained with only one particular evolutionary scenario, the MWC 656 binary. This is only a third such estimate available (after Cyg X-1 and Cyg X-3), and it lends additional support to the existence of so far undetected BH–NS binaries.

Age constraints on the tectonic evolution of the Itremo region in Central Madagascar

The Itremo region in Central Madagascar comprises a deformed metasedimentary sequence (Itremo Group) that has undergone greenschist to lower amphibolite facies metamorphism. During a first phase of deformation (D1) Itremo Group sediments were deformed into a fold-and-thrust belt and transported toward the E to NE on top of migmatitic gneisses rocks of Anatananarivo block. A second phase of deformation (D2) affected both the fold-and-thrust belt and structurally underlying units, and formed large-scale N-S trending folds with steeply dipping axial planes. A Late Neoproterozoic Th–U–Pb XRF monazite age (565±17 Ma) dates the emplacement of a granite that truncates first-phase structures in the Itremo Group, and indicates that the fold-and-thrust belt formed prior to ≈565 Ma. Th–U–Pb electron microprobe dating was applied to elongated monazites that lie within the first-phase foliation of Itremo Group metapelites. The detrital cores of zoned monazites reveal two distinct age populations at ∼2000 and 1700 Ma, the latter age giving a maximum depositional age for the Itremo Group. Statistical analysis of ages determined from the rims of zoned monazites and from unzoned monazites indicates three Late Proterozoic–Early Paleozoic monazite growth events at about 565–540, 500 and 430 Ma. The oldest age population is contemporaneous within error, with the intrusion of the dated granite. The two younger age populations are found both in the Th–U–Pb and Ar–Ar data; together with the perturbation of the Rb–Sr system we interpret both ages as due to alteration related to fluid circulation events, possibly connected to the emplacement of pegmatite fields in Central Madagascar. Syn-D1 tectonic growth of contact metamorphism minerals such as andalusite has been observed locally in metapelites along the margin of Middle Neoproterozoic (≈800 Ma) granites, suggesting that D1 in the Itremo Group is contemporaneous with the intrusion of granites at ≈800 Ma. The N-S trending D2 folds are associated with ≈E-W shortening during the final assembly of Gondwana in Late Neoproterozoic–Early Cambrian times.

The evolution of microRNA in primates

Thesis (Ph.D.)--University of Washington, 2016-06

Piecing together evolution of the vertebrate endocrine system

One of the great challenges in biology is to understand how particular complex morphological and physiological characters originated in specific evolutionary lineages. In this article, we address the origin of the vertebrate hypothalamic-pituitary-peripheral gland (H-P-PG) endocrine system, a complex network of specialized tissues, ligands and receptors. Analysis of metazoan nucleotide and protein sequences reveals a patchwork pattern of H-P-PG gene conservation between vertebrates and closely related invertebrates (ascidians). This is consistent with a model of how the vertebrate H-P-PG endocrine system could have emerged in relatively few steps by gene family expansion and by regulatory and structural modifications to genes that are present in a chordate ancestor. Some of these changes might have resulted in new connections between metabolic or signaling pathways, such as the bridging of 'synthesis islands' to form an efficient system for steroid hormone synthesis.

Evolution of polydnaviruses as insect immune suppressors

Polydnaviruses (PDVs) are endogenous particles that are used by some endoparasitic hymenoptera to disrupt host immunity and development. Recent analyses of encapsidated PDV genes have increased the number of known PDV gene families, which are often closely related to insect genes. Several PDV proteins inactivate host haemocytes by damaging their actin cytoskeleton. These proteins share no significant sequence homology and occur in polyphyletic PDV genera, possibly indicating that convergent evolution has produced functionally similar immune-suppressive molecules causing a haemocyte phenotype characterised by damaged cytoskeleton and inactivation. These phenomena provide further insights into the immune-suppressive activity of PDVs and raise interesting questions about PDV evolution, a topic that has puzzled researchers ever since the discovery of PDVs.

Evolution of the regular zone of histone H1 in Fabaceae plants

An analysis of the historic H1 subtype, H1-1, in eight legumes belonging to four genera of the tribe Vicieae (Pisum, Lathyrus, Lens, and Vicia), revealed an extended region consisting of the tandemly repeated AKPAAK motifs. We named this region the Regular zone (RZ). The AKPAAK motifs are organized into two blocks separated by a short (two or six amino acids) intervening sequence (IS). The distal block contains six AKPAAK motifs, while the number of repeats in the proximal block varies from six in V. faba to seven in the other species. In V. hirsuta, the first two repeated units of the proximal block are octapeptides AKAKPAAK. The apparent rate of synonymous substitutions in the blocks of RZ is much higher than in the rest of the gene. This can be explained by repeat shuffling within each block. In the C-domain of the orthologous H1 subtype froth Medicago truncatula (tribe Trifolieae), a region corresponding to the RZ of Vicieae species was found. It also consists of two blocks of AKPAAK motifs (four and three repeats in the proximal and distal blocks, respectively). These blocks are separated by a 20-amino acid IS. The first 20 amino acids of the Medicago RZ are not part of AKPAAK repeats. We hypothesise that the RZ has most probably evolved as a result of an expansion of AKPAAK repeats from two separate sites in the C-domain. This process started tens of millions of years ago and was most likely directed by positive selection.

Evolution of resistance to sulfadoxine-pyrimethamine in Plasmodium falciparum

The development of resistance to sulfadoxine-pyrimethamine by Plasmodium parasites is a major problem for the effective treatment of malaria, especially P. falciparum malaria. Although the molecular basis for parasite resistance is known, the factors promoting the development and transmission of these resistant parasites are less clear. This paper reports the results of a quantitative comparison of factors previously hypothesized as important for the development of drug resistance, drug dosage, time of treatment, and drug elimination half-life, with an in-host dynamics model of P. falciparum malaria in a malaria-naive host. The results indicate that the development of drug resistance can be categorized into three stages. The first is the selection of existing parasites with genetic mutations in the dihydrofolate reductase or dihydropteroate synthetase gene. This selection is driven by the long half-life of the sulfadoxine-pyrimethamine combination. The second stage involves the selection of parasites with allelic types of higher resistance within the host during an infection. The timing of treatment relative to initiation of a specific anti-P. falciparum EMP1 immune response is an important factor during this stage, as is the treatment dosage. During the third stage, clinical treatment failure becomes prevalent as the parasites develop sufficient resistance mutations to survive therapeutic doses of the drug combination. Therefore, the model output reaffirms the importance of correct treatment of confirmed malaria cases in slowing the development of parasite resistance to sulfadoxine-pyrimethamine.

Sequence, circulating levels, and expression of C-type natriuretic peptide in a euryhaline elasmobranch, Carcharhinus leucas

The present study has examined expression and circulating levels of C-type natriuretic peptide (CNP) in the euryhaline bull shark, Carcharhinus leucas. Complementary DNA and deduced amino acid sequence for CNP in C leucas were determined by RACE methods. Homology of CNP amino acid sequence in C. leucas was high both for proCNP and for mature CNP when compared with previously identified elasmobranch CNPs. Mature CNP sequence in C. leucas was identical to that in Triakis seyllia and Seyliorhinus canicula. Levels of expression of CNP mRNA were significantly decreased in the atrium but did not change in either the brain or ventricle following acclimation to a SW environment. However, circulating levels of CNP significantly increased from 86.0 +/- 7.9 fmol ml(-1) in FW to 144.9 +/- 19.5 fmol ml(-1) in SW. The results presented demonstrate that changes in environmental salinity influences both synthesis of CNP from the heart and also circulating levels in C. leucas. Potential stimulus for release and modes of action are discussed. (c) 2005 Elsevier Inc. All rights reserved.

SeqDoC: rapid SNP and mutation detection by direct comparison of DNA sequence chromatograms

Background: This paper describes SeqDoC, a simple, web-based tool to carry out direct comparison of ABI sequence chromatograms. This allows the rapid identification of single nucleotide polymorphisms (SNPs) and point mutations without the need to install or learn more complicated analysis software. Results: SeqDoC produces a subtracted trace showing differences between a reference and test chromatogram, and is optimised to emphasise those characteristic of single base changes. It automatically aligns sequences, and produces straightforward graphical output. The use of direct comparison of the sequence chromatograms means that artefacts introduced by automatic base-calling software are avoided. Homozygous and heterozygous substitutions and insertion/deletion events are all readily identified. SeqDoC successfully highlights nucleotide changes missed by the Staden package 'tracediff' program. Conclusion: SeqDoC is ideal for small-scale SNP identification, for identification of changes in random mutagenesis screens, and for verification of PCR amplification fidelity. Differences are highlighted, not interpreted, allowing the investigator to make the ultimate decision on the nature of the change.

The prion protein gene: Identifying regulatory signals using marsupial sequence

The function of the prion protein gene (PRNP) and its normal product PrPC is elusive. We used comparative genomics as a strategy to understand the normal function of PRNP. As the reliability of comparisons increases with the number of species and increased evolutionary distance, we isolated and sequenced a 66.5 kb BAC containing the PRNP gene from a distantly related mammal, the model Australian marsupial Macropus eugenii (tammar wallaby). Marsupials are separated from eutherians such as human and mouse by roughly 180 million years of independent evolution. We found that tammar PRNP, like human PRNP, has two exons. Prion proteins encoded by the tammar wallaby and a distantly related marsupial, Monodelphis domestica (Brazilian opossum) PRNP contain proximal PrP repeats with a distinct, marsupial-specific composition and a variable number. Comparisons of tammar wallaby PRNP with PRNPs from human, mouse, bovine and ovine allowed us to identify non-coding gene regions conserved across the marsupial-eutherian evolutionary distance, which are candidates for regulatory regions. In the PRNP 3' UTR we found a conserved signal for nuclear-specific polyadenylation and the putative cytoplasmic polyadenylation element (CPE), indicating that post-transcriptional control of PRNP mRNA activity is important. Phylogenetic footprinting revealed conserved potential binding sites for the MZF-1 transcription factor in both upstream promoter and intron/intron 1, and for the MEF2, MyTI, Oct-1 and NFAT transcription factors in the intron(s). The presence of a conserved NFAT-binding site and CPE indicates involvement of PrPC in signal transduction and synaptic plasticity. (c) 2004 Elsevier B.V. All rights reserved.