The optimization principle in phylogenetic analysis tends to give incorrect topologies when the number of nucleotides or amino acids used is small

In the maximum parsimony (MP) and minimum evolution (ME) methods of phylogenetic inference, evolutionary trees are constructed by searching for the topology that shows the minimum number of mutational changes required (M) and the smallest sum of branch lengths (S), respectively, whereas in the maximum likelihood (ML) method the topology showing the highest maximum likelihood (A) of observing a given data set is chosen. However, the theoretical basis of the optimization principle remains unclear. We therefore examined the relationships of M, S, and A for the MP, ME, and ML trees with those for the true tree by using computer simulation. The results show that M and S are generally greater for the true tree than for the MP and ME trees when the number of nucleotides examined (n) is relatively small, whereas A is generally lower for the true tree than for the ML tree. This finding indicates that the optimization principle tends to give incorrect topologies when n is small. To deal with this disturbing property of the optimization principle, we suggest that more attention should be given to testing the statistical reliability of an estimated tree rather than to finding the optimal tree with excessive efforts. When a reliability test is conducted, simplified MP, ME, and ML algorithms such as the neighbor-joining method generally give conclusions about phylogenetic inference very similar to those obtained by the more extensive tree search algorithms.

Atom density in protein structures

The residue environment in protein structures is studied with respect to the density of carbon (C), oxygen (O), and nitrogen (N) atoms within a certain distance (say 5 Å) of each residue. Two types of environments are evaluated: one based on side-chain atom contacts (abbreviated S-S) and the other based on all atom (side-chain + backbone) contacts (abbreviated A-A). Different atom counts are observed about nine-residue structural categories defined by three solvent accessibility levels and three secondary structure states. Among the structural categories, the S-S atom count ratios generally vary more than the A-A atom count ratios because of the fact that the backbone (O) and (N) atoms contribute equal counts. Secondary structure affects the (C) density for the A-A contacts whereas secondary structure has little influence on the (C) density for the S-S contacts. For S-S contacts, a greater density of (O) over (N) atom neighbors stands out in the environment of most amino acid types. By contrast, for A-A contacts, independent of the solvent accessibility levels, the ratio (O)/(N) is ≈1 in helical states, consistent with the geometry of α-helical residues whose side-chains tilt oppositely to the amino to carboxy α-helical axis. The highest ratio of neighbor (O)/(N) is achieved under solvent exposed conditions. This (O) vs. (N) prevalence is advantageous at the protein surface that generally exhibits an acid excess that helps to enhance protein solubility in the cell and to avoid nonspecific interactions with phosphate groups of DNA, RNA, and other plasma constituents.

Dual role of the nuclear factor of activated T cells insert region in DNA recognition and cooperative contacts to activator protein 1

The transcription factors nuclear factor of activated T cells (NFAT) and activator protein 1 (AP-1) coordinately regulate cytokine gene expression in activated T-cells by binding to closely juxtaposed sites in cytokine promoters. The structural basis for cooperative binding of NFAT and AP-1 to these sites, and indeed for the cooperative binding of transcription factors to composite regulatory elements in general, is not well understood. Mutagenesis studies have identified a segment of AP-1, which lies at the junction of its DNA-binding and dimerization domains (basic region and leucine zipper, respectively), as being essential for protein–protein interactions with NFAT in the ternary NFAT/AP-1/DNA complex. In a model of the ternary complex, the segment of NFAT nearest AP-1 is the Rel insert region (RIR), a feature that is notable for its hypervariability in size and in sequence amongst members of the Rel transcription factor family. Here we have used mutational analysis to study the role of the NFAT RIR in binding to DNA and AP-1. Parallel yeast one-hybrid screening assays in combination with alanine-scanning mutagenesis led to the identification of four amino acid residues in the RIR of NFAT2 (also known as NFATC1 or NFATc) that are essential for cooperativity with AP-1 (Ile-544, Glu-545, Thr-551, and Ile-553), and three residues that are involved in interactions with DNA (Lys-538, Arg-540, and Asn-541). These results were confirmed and extended through in vitro binding assays. We thus conclude that the NFAT RIR plays an essential dual role in DNA recognition and cooperative binding to AP-1 family transcription factors.

Electrostatic interaction between helical macromolecules in dense aggregates: An impetus for DNA poly- and meso-morphism

DNA exhibits a surprising multiplicity of structures when it is packed into dense aggregates. It undergoes various polymorphous transitions (e.g., from the B to A form) and mesomorphous transformations (from hexagonal to orthorhombic or monoclinic packing, changes in the mutual alignment of nearest neighbors, etc). In this report we show that such phenomena may have their origin in the specific helical symmetry of the charge distribution on DNA surface. Electrostatic interaction between neighboring DNA molecules exhibits strong dependence on the patterns of molecular surface groups and adsorbed counter-ions. As a result, it is affected by such structural parameters as the helical pitch, groove width, the number of base pairs per helical turn, etc. We derive expressions which relate the energy of electrostatic interaction with these parameters and with the packing variables characterizing the axial and azimuthal alignment between neighboring macromolecules. We show, in particular, that the structural changes upon the B-to-A transition reduce the electrostatic energy by ≈kcal/mol per base pair, at a random adsorption of counter ions. Ion binding into the narrow groove weakens or inverts this effect, stabilizing B-DNA, as it is presumably the case in Li+-DNA assemblies. The packing symmetry and molecular alignment in DNA aggregates are shown to be affected by the patterns of ion binding.

Did homeodomain proteins duplicate before the origin of angiosperms, fungi, and metazoa?

Homeodomain proteins are transcription factors that play a critical role in early development in eukaryotes. These proteins previously have been classified into numerous subgroups whose phylogenetic relationships are unclear. Our phylogenetic analysis of representative eukaryotic sequences suggests that there are two major groups of homeodomain proteins, each containing sequences from angiosperms, metazoa, and fungi. This result, based on parsimony and neighbor-joining analyses of primary amino acid sequences, was supported by two additional features of the proteins. The two protein groups are distinguished by an insertion/deletion in the homeodomain, between helices I and II. In addition, an amphipathic alpha-helical secondary structure in the region N terminal of the homeodomain is shared by angiosperm and metazoan sequences in one group. These results support the hypothesis that there was at least one duplication of homeobox genes before the origin of angiosperms, fungi, and metazoa. This duplication, in turn, suggests that these proteins had diverse functions early in the evolution of eukaryotes. The shared secondary structure in angiosperm and metazoan sequences points to an ancient conserved functional domain.

Accurate reconstruction of a known HIV-1 transmission history by phylogenetic tree analysis.

Phylogenetic analyses are increasingly used in attempts to clarify transmission patterns of human immunodeficiency virus type 1 (HIV-1), but there is a continuing discussion about their validity because convergent evolution and transmission of minor HIV variants may obscure epidemiological patterns. Here we have studied a unique HIV-1 transmission cluster consisting of nine infected individuals, for whom the time and direction of each virus transmission was exactly known. Most of the transmissions occurred between 1981 and 1983, and a total of 13 blood samples were obtained approximately 2-12 years later. The p17 gag and env V3 regions of the HIV-1 genome were directly sequenced from uncultured lymphocytes. A true phylogenetic tree was constructed based on the knowledge about when the transmissions had occurred and when the samples were obtained. This complex, known HIV-1 transmission history was compared with reconstructed molecular trees, which were calculated from the DNA sequences by several commonly used phylogenetic inference methods [Fitch-Margoliash, neighbor-joining, minimum-evolution, maximum-likelihood, maximum-parsimony, unweighted pair group method using arithmetic averages (UPGMA), and a Fitch-Margoliash method assuming a molecular clock (KITSCH)]. A majority of the reconstructed trees were good estimates of the true phylogeny; 12 of 13 taxa were correctly positioned in the most accurate trees. The choice of gene fragment was found to be more important than the choice of phylogenetic method and substitution model. However, methods that are sensitive to unequal rates of change performed more poorly (such as UPGMA and KITSCH, which assume a constant molecular clock). The rapidly evolving V3 fragment gave better reconstructions than p17, but a combined data set of both p17 and V3 performed best. The accuracy of the phylogenetic methods justifies their use in HIV-1 research and argues against convergent evolution and selective transmission of certain virus variants.

Evolutionary relationships of the coelacanth, lungfishes, and tetrapods based on the 28S ribosomal RNA gene.

The origin of land vertebrates was one of the major transitions in the history of vertebrates. Yet, despite many studies that are based on either morphology or molecules, the phylogenetic relationships among tetrapods and the other two living groups of lobe-finned fishes, the coelacanth and the lungfishes, are still unresolved and debated. Knowledge of the relationships among these lineages, which originated back in the Devonian, has profound implications for the reconstruction of the evolutionary scenario of the conquest of land. We collected the largest molecular data set on this issue so far, about 3,500 base pairs from seven species of the large 28S nuclear ribosomal gene. All phylogenetic analyses (maximum parsimony, neighbor-joining, and maximum likelihood) point toward the hypothesis that lungfishes and coelacanths form a monophyletic group and are equally closely related to land vertebrates. This evolutionary hypothesis complicates the identification of morphological or physiological preadaptations that might have permitted the common ancestor of tetrapods to colonize land. This is because the reconstruction of its ancestral conditions would be hindered by the difficulty to separate uniquely derived characters from shared derived characters in the coelacanth/lungfish and tetrapod lineages. This molecular phylogeny aids in the reconstruction of morphological evolutionary steps by providing a framework; however, only paleontological evidence can determine the sequence of morphological acquisitions that allowed lobe-finned fishes to colonize land.

Interresidue hydrogen bonding in a peptide nucleic acid.RNA heteroduplex.

Previous molecular mechanics calculations suggest that strands of peptide nucleic acids (PNAs) and complementary oligonucleotides form antiparallel duplexes stabilized by interresidue hydrogen bonds. In the computed structures, the amide carbonyl oxygen nearest the nucleobase (O7') forms an interresidue hydrogen bond with the backbone amide proton of the following residue, (n + 1)H1'. Of the 10 published two dimensional 1H NMR structures of a hexameric PNA.RNA heteroduplex. PNA(GAACTC).r(GAGUUC), 9 exhibit two to five potential interresidue hydrogen bonds. In our minimized average structure, created from the coordinates of these 10 NMR structures, three of the five possible interresidue hydrogen bond sites within the PNA backbone display the carbonyl oxygen (O7') and the amide proton (n + 1)H1' distances and N1'-H1'-(n - 1)O7' angles optimal for hydrogen bond formation. The finding of these interresidue hydrogen bonds supports the results of our previous molecular mechanics calculations.

Genetic heterogeneity in cystinuria: the SLC3A1 gene is linked to type I but not to type III cystinuria.

Cystinuria is an autosomal recessive amino-aciduria where three urinary phenotypes have been described (I, II, and III). An amino acid transporter gene, SLC3A1 (formerly rBAT), was found to be responsible for this disorder. To assess whether mutations in SLC3A1 are involved in different cystinuria phenotypes, linkage with this gene and its nearest marker (D2S119) was analyzed in 22 families with type I and/or type III cystinuria. Linkage with heterogeneity was proved (alpha = 0.45; P < 0.008). Type I/I families showed homogeneous linkage to SLC3A1 (Zmax > 3.0 at theta = 0.00; alpha = 1), whereas types I/III and III/III were not linked. Our data suggest that type I cystinuria is due to mutations in the SLC3A1 gene, whereas another locus is responsible for type III. This result establishes genetic heterogeneity for cystinuria, classically considered as a multiallelic monogenic disease.

Nucleotide sequence of the cDNA encoding the proenzyme of phenol oxidase A1 of Drosophila melanogaster.

Clones encoding pro-phenol oxidase [pro-PO; zymogen of phenol oxidase (monophenol, L-dopa:oxygen oxidoreductase, EC 1.14.18.1)] A1 were isolated from a lambda gt10 library that originated from Drosophila melanogaster strain Oregon-R male adults. The 2294 bp of the cDNA included a 13-bp 5'-noncoding region, a 2070-bp encoding open reading frame of 690 amino acids, and a 211-bp 3'-noncoding region. A hydrophobic NH2-terminal sequence for a signal peptide is absent in the protein. Furthermore, there are six potential N-glycosylation sites in the sequence, but no amino sugar was detected in the purified protein by amino acid analysis, indicating the lack of an N-linked sugar chain. The potential copper-binding sites, amino acids 200-248 and 359-414, are highly homologous to the corresponding sites of hemocyanin of the tarantula Eurypelma californicum, the horseshoe crab Limulus polyphemus, and the spiny lobster Panulirus interruptus. On the basis of the phylogenetic tree constructed by the neighbor-joining method, vertebrate tyrosinases and molluscan hemocyanins constitute one family, whereas pro-POs and arthropod hemocyanins group with another family. It seems, therefore, likely that pro-PO originates from a common ancestor with arthropod hemocyanins, independently to the vertebrate and microbial tyrosinases.

Cultural variation in Africa: role of mechanisms of transmission and adaptation.

Cultural inheritance can be considered as a mechanism of adaptation made possible by communication, which has reached its greatest development in humans and can allow long-term conservation or rapid change of culturally transmissible traits depending on circumstances and needs. Conservativeness/flexibility is largely modulated by mechanisms of sociocultural transmission. An analysis was carried out by testing the fit of three models to 47 cultural traits (classified in six groups) in 277 African societies. Model A (demic diffusion) is conservation over generations, as shown by correlations of cultural traits with language, used as a measure of historical connection. Model B (environmental adaptation) is measured by correlation to the natural environment. Model C (cultural diffusion) is the spread to neighbors by social contact in an epidemic-like fashion and was tested by measuring the tightness of geographic clustering of the traits. Most traits examined, in particular those affecting family structure and kinship, showed great conservation over generations, as shown by the fit of model A. They are most probably transmitted by family members. This is in agreement with the theoretical demonstration that cultural transmission in the family (vertical) is the most conservative one. Some traits show environmental effects, indicating the importance of adaptation to physical environment. Only a few of the 47 traits showed tight geographic clustering indicating that their spread to nearest neighbors follows model C, as is usually the case for transmission among unrelated people (called horizontal transmission).

NusA interferes with interactions between the nascent RNA and the C-terminal domain of the alpha subunit of RNA polymerase in Escherichia coli transcription complexes.

The effects of NusA on the RNA polymerase contacts made by nucleotides at internal positions in the nascent RNA in Escherichia coli transcription complexes were analyzed by using the photocrosslinking nucleotide analog 5-[(4-azidophenacyl) thio]-UMP. It was placed at nucleotides between +6 and +15 in RNA transcribed from the phage lambda PR' promoter. Crosslinks of analog in these positions in RNAs which contained either 15, 28, 29, or 49 nt were examined. Contacts between the nascent RNA and proteins in the transcription complex were analyzed as the RNA was elongated, by placing the crosslinker nearest the 5' end of the RNA 10, 23, 24, or 44 nt away from the 3' end. The beta or beta' subunit of polymerase, and NusA when added, were contacted by RNA from 15 to 49 nt long. When the upstream crosslinker was 24 nt from the 3" end of the RNA (29-nt RNA), alpha was also contacted in the absence of NusA. The addition of NusA prevented RNA crosslinking to alpha. When the crosslinker was 44 nt from the 3' end (49-nt RNA), alpha crosslinks were still observed, but crosslinks to beta or beta' and NusA were greatly diminished. RNA crosslinking to alpha, and loss of this crosslink when NusA was added, was observed in the presence of NusB, NusE, and NusG and when transcription was carried out in the presence of an E. coli S100 cell extract. Peptide mapping localized the RNA interactions to the C-terminal domain of alpha.

Trans-Pacific migrations of the loggerhead turtle (Caretta caretta) demonstrated with mitochondrial DNA markers.

Juvenile loggerhead turtles (Caretta caretta) have recently been documented in the vicinity of Baja California and thousands of these animals have been captured in oceanic fisheries of the North Pacific. The presence of loggerhead turtles in the central and eastern North Pacific is a prominent enigma in marine turtle distribution because the nearest documented nesting concentrations for this species are in Australia and Japan, over 10,000 km from Baja California. To determine the origin of the Baja California feeding aggregate and North Pacific fishery mortalities, samples from nesting areas and pelagic feeding aggregates were compared with genetic markers derived from mtDNA control region sequences. Overall, 57 of 60 pelagic samples (95%) match haplotypes seen only in Japanese nesting areas, implicating Japan as the primary source of turtles in the North Pacific Current and around Baja California. Australian nesting colonies may contribute the remaining 5% of these pelagic feeding aggregates. Juvenile loggerhead turtles apparently traverse the entire Pacific Ocean, approximately one-third of the planet, in the course of developmental migrations, but mortality in high-seas fisheries raises concern over the future of this migratory population.

Significados da reinserção escolar de crianças após queimaduras: um estudo etnográfico

As queimaduras na infância compreendem importante causa de atendimento hospitalar e internação e podem desencadear sequelas físicas e emocionais. Dependendo da gravidade e complexidade das lesões, a hospitalização da criança é prolongada para realização de tratamento adequado e ela é afastada tanto do convívio familiar como do social. O número de crianças em idade escolar que sofrem queimaduras, e consequentemente interrompem as atividades escolares, por períodos curtos ou até as abandonam, é significativo, tornando digno de preocupação o processo de reinserção escolar. Este estudo teve como objetivo interpretar os significados do processo de reinserção escolar de crianças sobreviventes de queimaduras. Realizou-se uma pesquisa com abordagem metodológica qualitativa, fundamentada no referencial teórico da Antropologia Interpretativa e no método etnográfico. Após a aprovação do estudo pelos Comitês de Ética em Pesquisa com Seres Humanos das instituições envolvidas, convidaram-se crianças sobreviventes de queimaduras em seguimento em um hospital-escola do interior paulista. Coletaram-se dados no hospital, no domicílio e na escola das crianças, de janeiro de 2012 a dezembro de 2013, por meio de entrevistas em profundidade audiogravadas e observação participante, complementadas pelo diário de campo. Participaram da pesquisa 14 crianças e os atores sociais envolvidos neste processo, como familiares, profissionais de saúde, professores, vizinha e amigo, totalizando 57 participantes. A coleta e análise dos dados ocorreram concomitantemente, e esta última seguiu os pressupostos da análise temática indutiva. Identificaram-se códigos, os quais, posteriormente, embasaram a construção das duas unidades de sentidos: \"Fatos e fatores prévios ao retorno escolar\" e \"A volta à escola\". A partir destas, identificou-se o núcleo temático, \"O olho puxa essas coisas. A gente olha mesmo\": enxergando o outro como diferente, o qual apresenta a explicação compreensiva do processo de reinserção escolar, fornecendo os significados da experiência. Os significados foram explicados por meio dos conceitos antropológicos de estigma, identidade e corporeidade. O culto ao corpo presente no contexto cultural brasileiro foi fator intensificador do estigma sofrido pelas crianças. Identificaram-se fatores facilitadores do processo, como: dar continuidade às atividades escolares e manter contato com professores e colegas durante a hospitalização; preparar a escola para receber a criança; abordar sobre queimadura na escola para evitar curiosidade; encarar a presença da discriminação e trabalhar as diferenças. Estas estratégias fazem com que todos os envolvidos se sintam ao menos mais confortáveis durante esse difícil processo que merece atenção e envolvimento dos familiares, profissionais de saúde e equipe escolar. A partir dos resultados desta pesquisa, esperamos empoderar as crianças e seus familiares acerca da melhor maneira para lidar com as queimaduras durante a reinserção escolar e contribuir para o desenvolvimento de ações e estratégias baseadas na cultura para que este processo ocorra da forma menos traumática possível

Propriedades magnéticas de trímeros de FexCo1-x depositados em Pt(111)

Quais propriedades magnéticas são modificadas quando se agrupam átomos de Fe/Co para formar estruturas quasi-2D, se comparadas aos nanofios (quasi-1D) de FexCo1-x? E como estas propriedades reagem com a variação da proporção de Fe/Co nos aglomerados? A fim de responder a estas questões, trímeros de FexCo1-x depositados em Pt(111) são investigados utilizando o método de primeiros princípios Real Space-Linear Muffin-Tin Orbital-Atomic Sphere Approximation (RS-LMTO-ASA) no âmbito da Teoria do Funcional da Densidade (DFT). Diferentes configurações de trímeros triangulares são consideradas, variando-se as posições e a concentração dos átomos de Fe/Co. Neste trabalho, demonstra-se a ocorrência de uma tendência não-linear estritamente decrescente dos momentos orbitais médios como função da concentração de Fe, distinta do encontrado tanto para os nanofios de FexCo1-x (dependência linear) quanto para a monocamada correspondente (dependência não-linear). Os resultados obtidos mostram ainda que os momentos orbitais variam com o ambiente local e com a direção de magnetização, especialmente quando associados aos átomos de Co, em concordância com publicações anteriores. A mudança de dimensionalidade quasi-1D (nanofios) para quasi-2D (trímeros compactos) não afeta o comportamento dos momentos de spin, que permanecem descritos por uma função linear com respeito à proporção de Fe/Co. Ambos o formato e a concentração de Fe nos sistemas apresentam um papel importante nos valores de energia de anisotropia magnética. Em adição, observou-se que o subtrato de Pt opera ativamente na definição das propriedades magnéticas dos aglomerados. Embora todas as configurações lineares e compactas dos aglomerados de FexCo1-x sejam estáveis e exibam interações fortemente ferromagnéticas entre os primeiros vizinhos, nem todas revelaram o ordenamento colinear como estado fundamental, apresentando uma interação de Dzyaloshinskii-Moriya não-desprezível induzida pelo acoplamento spin-órbita. Estes casos específicos são: o trímero triangular de Co puro e o trímero linear (nanofio) de Fe puro, para o qual foi verificado o acoplamento do tipo Ruderman-Kittel-Kasuya-Yosida entre os átomos de Fe constituintes. Os resultados obtidos contribuem para o entendimento de quais mecanismos definem o magnetismo nos trímeros de FexCo1-x/Pt(111), e discutem as questões presentes atualmente na literatura no contexto destes sistemas.