Positive selection on the mitochondrial ATP synthase 6 and the NADH dehydrogenase 2 genes across 22 hare species (genus Lepus)

Amino acid changes in mitochondrial (mt) oxidative phosphorylation (OXPHOS) genes have been suggested as a key adaptation to environmental variation. Here, we analyzed 416 sequences of ATPase synthase 6 (MT‐ATP6) and NADH dehydrogenase 2 (MT‐ND2) in 22 different hare (Lepus) species from across a wide range of habitats and climates. We used site‐ and branch‐based methods to test for positive selection on specific codons and lineages. We found four codons in MT‐ATP6 and five in MT‐ND2 under positive selection, affecting several species lineages. We investigated the association of protein variants at each locus with climate zone, using multinomial generalized linear models (glm), including species, regions, historical introgression events, and the co‐occurring protein variant at the other locus as additional explanatory variables. A significant climate effect as based on the “Köppen climate classification” was observed for MT‐ND2 protein variants as translated from our nucleotide sequences. Moreover, MT‐ND2 protein variants were significantly affected by the co‐occurring MT‐ATP6 protein variant in the same mtDNA molecule. Contrary to the expectation for non‐recombining mitochondrial DNA molecules, the presence of an evolutionarily relatively ancestral protein variant at one locus was associated with a relatively derived protein at the other locus in the same mitochondrial molecule, respectively. The relative evolutionary status of a protein variant was evaluated according to its positions relative to the respective out‐group protein variant in a network analysis of nucleotide sequences. All our results suggest a complex effect of various climatic parameters acting on multiple mtOXPHOS genes in a co‐adaptive way, favoring combinations of ancestral and derived variants.     

Revisiting the Molecular Evolutionary History of <Emphasis Type="Italic">Shigella</Emphasis> spp.

The theory that Shigella is derived from multiple independent origins of Escherichia coli (Pupo et al. 2000) has been challenged by recent findings that the virulence plasmids (VPs) and the chromosomes share a similar evolutionary history (Escobar-Paramo et al. 2003), which suggests that an ancestral VP entered an E. coli strain only once, which gave rise to Shigella spp. In an attempt to resolve these conflicting theories, we constructed three phylogenetic trees in this study: a robust chromosomal tree using 23 housekeeping genes from 46 strains of Shigella and enteroinvasive E. coli (EIEC), a chromosomal tree using 4 housekeeping genes from 19 EcoR strains and 46 Shigella/EIEC strains, and a VP tree using 5 genes outside of the VP cell-entry region from 38 Shigella/EIEC strains. Both chromosomal trees group Shigella into three main clusters and five outliers, and strongly suggest that Shigella has multiple origins within E. coli. Most strikingly, the VP tree shows that the VPs from two main Shigella clusters, C1 and C2, are more closely related, which contradicts the chromosomal trees that place C2 and C3 next to each other but C1 at a distance. Additionally, we have identified a complete tra operon of the F-plasmid in the genome sequence of an EIEC strain and found that two other EIEC strains are also likely to possess a complete tra operon. All lines of evidence support an alternative multiorigin theory that transferable diverse ancestral VPs entered diverse origins of E. coli multiple times during a prolonged period of time, resulting in Shigella species with diverse genomes but similar pathogenic properties. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Martin Kreitman] Jian Yang and Huan Nie contributed equally to this work.     

Expression of ASK1-like genes in arrested stamens of female Silene latifolia plants

Silene latifolia is a dioecious plant in which sex is determined by heteromorphic sex chromosomes. In female plants, stamen development is arrested before microspore mother cells are formed. In this study, we isolated four cDNAs (SlSKP1-1 to 4) encoding ASK1-like protein as expression markers to reveal when expression levels are reduced in arrested stamens of female flowers. Expression patterns of the SlSKP1 genes were analyzed by in-situ hybridization. We use the flower development classification of Grant et al. (in Plant J 6:471–480, 1994). SlSKP1 genes were highly expressed in primary parietal cells and primary sporogenous cells that develop into microspore mother cells in male flowers. Expression levels started to be reduced in the external stamens of the female flowers when stamen development was arrested at stage 7. Although microspore mother cells could not be developed in female flowers and SlSKP1 expression may be unnecessary in arrested stamens, SlSKP1 genes were still expressed in sporogenous cells of degenerated stamens at stage 8. Parietal cells stopped differentiating earlier than sporogenous cells in arrested stamens. These results suggest that not all types of cell are arrested simultaneously at a particular stage of stamen development during stamen suppression of S. latifolia.     

Hox gene expression in larval development of the polychaetes Nereis virens and Platynereis dumerilii (Annelida, Lophotrochozoa)

The bilaterian animals are divided into three great branches: the Deuterostomia, Ecdysozoa, and Lophotrochozoa. The evolution of developmental mechanisms is less studied in the Lophotrochozoa than in the other two clades. We have studied the expression of Hox genes during larval development of two lophotrochozoans, the polychaete annelids Nereis virens and Platynereis dumerilii. As reported previously, the Hox cluster of N. virens consists of at least 11 genes (de Rosa R, Grenier JK, Andreeva T, Cook CE, Adoutte A, Akam M, Carroll SB, Balavoine G, Nature, 399:772–776, 1999; Andreeva TF, Cook C, Korchagina NM, Akam M, Dondua AK, Ontogenez 32:225–233, 2001); we have also cloned nine Hox genes of P. dumerilii. Hox genes are mainly expressed in the descendants of the 2d blastomere, which form the integument of segments, ventral neural ganglia, pre-pygidial growth zone, and the pygidial lobe. Patterns of expression are similar for orthologous genes of both nereids. In Nereis, Hox2, and Hox3 are activated before the blastopore closure, while Hox1 and Hox4 are activated just after this. Hox5 and Post2 are first active during the metatrochophore stage, and Hox7, Lox4, and Lox2 at the late nectochaete stage only. During larval stages, Hox genes are expressed in staggered domains in the developing segments and pygidial lobe. The pattern of expression of Hox cluster genes suggests their involvement in the vectorial regionalization of the larval body along the antero-posterior axis. Hox gene expression in nereids conforms to the canonical patterns postulated for the two other evolutionary branches of the Bilateria, the Ecdysozoa and the Deuterostomia, thus supporting the evolutionary conservatism of the function of Hox genes in development. Milana Kulakova, Nadezhda Bakalenko and Elena Novikova contributed equally to this work.     

Camelid Ig V genes reveal significant human homology not seen in therapeutic target genes,providing for a powerful therapeutic antibody platform

Camelid immunoglobulin variable (IGV) regions were found homologous to their human counterparts; however, the germline V repertoires of camelid heavy and light chains are still incomplete and their therapeutic potential is only beginning to be appreciated. We therefore leveraged the publicly available HTG and WGS databases of Lama pacos and Camelus ferus to retrieve the germline repertoire of V genes using human IGV genes as reference. In addition, we amplified IGKV and IGLV genes to uncover the V germline repertoire of Lama glama and sequenced BAC clones covering part of the Lama pacos IGK and IGL loci. Our in silico analysis showed that camelid counterparts of all human IGKV and IGLV families and most IGHV families could be identified, based on canonical structure and sequence homology. Interestingly, this sequence homology seemed largely restricted to the Ig V genes and was far less apparent in other genes: 6 therapeutically relevant target genes differed significantly from their human orthologs. This contributed to efficient immunization of llamas with the human proteins CD70, MET, interleukin (IL)-1β and IL-6, resulting in large panels of functional antibodies. The in silico predicted human-homologous canonical folds of camelid-derived antibodies were confirmed by X-ray crystallography solving the structure of 2 selected camelid anti-CD70 and anti-MET antibodies. These antibodies showed identical fold combinations as found in the corresponding human germline V families, yielding binding site structures closely similar to those occurring in human antibodies. In conclusion, our results indicate that active immunization of camelids can be a powerful therapeutic antibody platform.     

A Voucher Specimen for <Emphasis Type="Italic">Macaca munzala</Emphasis>: Interspecific Affinities,Evolution, and Conservation of a Newly Discovered Primate

Sinha, A., Datta, A., Madhusudan, M. D., & Mishra, C. (2005. Macaca munzala: A new species from western Arunachal Pradesh, northeastern India. International Journal of Primatology, 26, 977–989) discovered Arunachal macaques (Macaca munzala), a species new to science, in the eastern Himalaya of Arunachal Pradesh in northeastern India. They depicted the holotype and paratypes of the species in photographs, and a specimen of the species had been unavailable for preservation and examination. In March 2005, we obtained an entire specimen of an adult male Macaca munzala, which we propose as a voucher specimen for the species. We provide detailed morphological and anatomical measurements of the specimen and examine its affinities with other macaques. Macaca munzala appears to be unique among macaques in craniodental size and structure, baculum, and aspects of caudal structure, while exhibiting affinities with the other members of the sinica-group to which it belongs. We summarize our insights on the origins and phylogeny of Macaca munzala. Finally, we review the current conservation status of the macaques, which are threatened by extensive hunting in the only 2 districts of Arunachal Pradesh where they are documented to occur.     

Distribution Patterns of Isomorphic Cold-Water Dinoflagellates (Scrippsiella/Woloszynskia Complex) Causing ‘red tides’ in the Baltic Sea

During the latest years medium-sized (15–30 μm), single-celled dinoflagellates have been reported to form blooms in the northern Baltic Proper and the Gulf of Finland in winter and spring. Recent studies (Kremp et al., 2003. Proceedings of the 7th International conference of Modern and Fossil Dinoflagellates, September 21–25, Nagasaki, Japan, 66 pp.) indicate that those blooms are caused by two isomorphic species – Scrippsiella hangoei (Schiller) Larsen, and a new species, tentatively belonging to the genus Woloszynskia. Until now there has been no report on how widely distributed these phytoplankton species are in the Baltic Sea. In this study, the occurrence of Scrippsiella/Woloszynskia complex in the entire Baltic Sea was investigated, by using monitoring data from 1997 to 2003. The species occurred in a salinity range from 2 to 8 PSU. Highest concentrations were observed at salinity 4.5–6.5 PSU. Maximum cell densities of Scrippsiella/Woloszynskia complex in the water column were mainly obtained in April or in the beginning of May by the water temperature <3 °C prior to stratification was formed. In the central Gulf of Finland, the second maximum was found in 1999 and 2002 by the temperature >6 °C. Bloom formations in the Baltic Proper and in the Gulf of Finland may not only be explained by optimum temperature and salinity, but also with other factors e.g. high nutrient concentrations and good seeding conditions from the sediments.     

Nematodes from saline and freshwater lakes of the Vestfold Hills,East Antarctica,including the description of <Emphasis Type="Italic">Hypodontolaimus antarcticus</Emphasis> sp. n.

The invertebrate fauna of many Antarctic ice-free areas, even those close to permanent research stations, can be poorly known. Here we describe some nematodes from freshwater and saline, marine-derived lakes of the Vestfold Hills, East Antarctica. The freshwater lakes contained the widespread East Antarctic endemic species, Plectus frigophilus Kirjanova, 1958. The saline lakes were inhabited by two recently described species, Halomonhystera halophila Andrássy, 2006 and Halomonhystera continentalis Andrássy, 2006, and by a new species described in this report, Hypodontolaimus antarcticus sp. n. Originally marine but now brackish Highway Lake contained a nematode fauna with both freshwater and marine-derived components. The nematode fauna of Antarctica now consists of 54 named species, 22 of which are found in East Antarctica.     

The rat Ruby (<Emphasis Type="Italic">R</Emphasis>) locus is <Emphasis Type="Italic">Rab38</Emphasis>: identical mutations in Fawn-hooded and Tester-Moriyama rats derived from an ancestral Long Evans rat sub-strain

Hermansky-Pudlak syndrome (HPS) is a group of rare, recessive disorders in which oculocutaneous albinism, progressive pulmonary fibrosis, bleeding diathesis, and other abnormalities result from defective biogenesis of multiple cytoplasmic organelles. Seven different HPS genes are known in humans; in mouse, at least 16 loci are associated with HPS-like mutant phenotypes. In the rat, only two HPS models are known, Fawn-hooded (FH) and Tester Moriyama (TM), non-complementing strains in which HPS-like hypopigmentation and platelet storage pool deficiency result from a mutation of the Ruby (red eyed dilution; R) locus on Chromosome (Chr) 1. We have identified the R locus as the Rab38 gene, establishing that rat R is homologous to mouse chocolate (cht). Further, we show that FH and TM rats have identical Rab38 Met1Ile mutations, occurring on an identical Chr 1 marker allele haplotype, indicating that these two strains derive from a common ancestor. This ancestor appears to have been a sub-strain of the outbred Long Evans (LE) strain, and several modern LE sub-strains carry the Rab38 Met1Ile R mutation on the same Chr 1 marker haplotype. These findings have significant implications for the many past and ongoing studies that involve the FH and LE-derivative rat strains. Hermansky-Pudlak syndrome (HPS; MIM 203300) is a group of autosomal recessive diseases in which oculocutaneous albinism (OCA), progressive and fatal pulmonary fibrosis, and bleeding diathesis due to platelet storage pool deficiency result from defects in the biogenesis of specific cytoplasmic organelles and granules: melanosomes, lysosomes, and platelet dense granules (reviewed in Spritz 1999, 2000; Spritz et al. 2003). In humans, seven different HPS genes are known (Oh et al. 1996; DellAngelica et al. 1999; Anikster et al. 2001; Suzuki et al. 2002; Li et al. 2003; Zhang et al. 2003). In the mouse, at least 16 loci associated with HPS-like mutant phenotypes are known, seven of which are homologous to the human HPS loci (Swank et al. 1998; Bennett and Lamoreux 2003). The nucleotide sequence data reported in this paper have been submitted to GenBank and have been assigned the accession number AY425759. (Naoki Oiso) Present address: Department of Dermatology, Saiseikai Tondabayashi Hospital, Tondabayashi, Osaka 584-0082, Japan.     

Analysis of a plastid gene cluster reveals a close relationship betweenCyanidioschyzon andCyanidium

Cyanidioschyzon merolae andCyanidium caldarium are representative species among of the most primitive algae, although the two species are distinctly different in various morphological traits. We determined the nucleotide sequence of therbcL gene and a flanking 8-kb region in the plastid genome of each of these algae. In both algae, 12 genes were identified in this region, in an identical order. This gene order is not conserved in the plastid genomes of other species of the kingdom Plantae that have been sequenced to data. An additional unidentified open reading frame was also found in the two algae that we analyzed, which has not been described in any other species of algae includingPorphyra purpurea. Comparison of the amino acid sequences of selected genes also supported the conclusion thatCyanidioschyzon merolae andCyanidium caldarium are closely related and that they are distinct from other rhodophytes. The nucleotide sequence data reported in this paper will appear in the DDBJ, EMBL and GenBank Nucleotide Sequence Databases under the accession numbers D63675 and D63676.     

Induction of a non-specific permeability transition in mitochondria from Yarrowia lipolytica and Dipodascus (Endomyces) magnusii yeasts

In this study we used tightly-coupled mitochondria from Yarrowia lipolytica and Dipodascus (Endomyces) magnusii yeasts, possessing a respiratory chain with the usual three points of energy conservation. High-amplitude swelling and collapse of the membrane potential were used as parameters for demonstrating induction of the mitochondrial permeability transition due to opening of a pore (mPTP). Mitochondria from Y. lipolytica, lacking a natural mitochondrial Ca²⁺ uptake pathway, and from D. magnusii, harboring a high-capacitive, regulated mitochondrial Ca²⁺ transport system (Bazhenova et al. J Biol Chem 273:4372–4377, 1998a; Bazhenova et al. Biochim Biophys Acta 1371:96–100, 1998b; Deryabina and Zvyagilskaya Biochemistry (Moscow) 65:1352–1356, 2000; Deryabina et al. J Biol Chem 276:47801–47806, 2001) were very resistant to Ca²⁺ overload. However, exposure of yeast mitochondria to 50–100 μM Ca²⁺ in the presence of the Ca²⁺ ionophore ETH129 induced collapse of the membrane potential, possibly due to activation of the fatty acid-dependent Ca²⁺/nH⁺-antiporter, with no classical mPTP induction. The absence of response in yeast mitochondria was not simply due to structural limitations, since large-amplitude swelling occurred in the presence of alamethicin, a hydrophobic, helical peptide, forming voltage-sensitive ion channels in lipid membranes. Ca²⁺- ETH129-induced activation of the Ca²⁺/H⁺-antiport system was inhibited and prevented by bovine serum albumin, and partially by inorganic phosphate and ATP. We subjected yeast mitochondria to other conditions known to induce the permeability transition in animal mitochondria, i.e., Ca²⁺ overload (in the presence of ETH129) combined with palmitic acid (Mironova et al. J Bioenerg Biomembr 33:319–331, 2001; Sultan and Sokolove Arch Biochem Biophys 386:37–51, 2001), SH-reagents, carboxyatractyloside (an inhibitor of the ADP/ATP translocator), depletion of intramitochondrial adenine nucleotide pools, deenergization of mitochondria, and shifting to acidic pH values in the presence of high phosphate concentrations. None of the above-mentioned substances or conditions induced a mPTP-like pore. It is thus evident that the permeability transition in yeast mitochondria is not coupled with Ca²⁺ uptake and is differently regulated compared to the mPTP of animal mitochondria.     

Rapid evolution of cytochrome c oxidase subunit II in camelids (Tylopoda, Camelidae)

Within cetartiodactyl species, both New and Old World camelids are uniquely adapted to the extremely hot and dry climates of African-Asian territories and to the high altitude cold and hypoxic environment of the whole Andean area. In order to investigate the potential association between these particular adaptations and mitochondrial aerobic energy production, we examined the camelid genes of cytochrome c oxidase subunits I, II, and III and the replacement of amino acids inferred. We found that all subunits had undergone a number of replacements in sites otherwise conserved in other cetartiodactyls. Changes of COXI and COXIII were mainly located in the transmembrane helices of proteins. For COXII, although most of the changes did not occur in sites directly involved in electron transfer, a shift of D by T at 115 position of Old World camelid might modify electrostatic interactions with cytochrome c. COXII also showed an increased relative evolutionary rate respect to other cetartiodactyls compared.     

Cloning and characterization of five cDNAs for genes differentially expressed during fruit development of kiwifruit (Actinidia deliciosa var.deliciosa)

Five cDNAs for genes differentially expressed during fruit development of kiwifruit (Actinidia deliciosa var.deliciosa cv. Hayward) were isolated from a library made from young fruit, 8–10 days after anthesis. One gene (pKIWI503) has low levels of expression in young fruit but is induced late in fruit development and during fruit ripening, and has some homology to plant metallothionein-like proteins. The other four genes are highly expressed in young fruit with reduced expression in the later stages of fruit development. pKIWI504 has strong homology to plant metallothionein-like proteins and pKIWI505 exhibits homology to the -subunit of the mitochondrial ATP synthase gene. The two other genes (pKIWI501 and 502) encode proteins with no significant homology to other known sequences.     

Organization and nucleotide sequence of the atp genes encoding the ATP synthase from alkaliphilic Bacillus firmus OF4

Summary The atp operon from the extreme alkaliphile Bacillus firmus OF4 was cloned and sequenced, and shown to contain genes for the eight structural subunits of the ATP synthase, preceded by a ninth gene predicted to encode a 14 kDa hydrophobic protein. The arrangement of genes is identical to that of the atp operons from Escherichia coli, Bacillus megaterium, and thermophilic Bacillus PS3. The deduced amino acid sequences of the subunits of the enzyme are also similar to their homologs in other ATP synthases, except for several unusual substitutions, particularly in the a and c subunits. These substitutions are in domains that have been implicated in the mechanism of proton translocation through F₀-ATPase, and therefore could contribute to the gating properties of the alkaliphile ATP synthase or its capacity for proton capture.     

Deep Origin of Plastid/Parasite ATP/ADP Translocases

Abstract Membrane proteins that transport ATP and ADP have been identified in mitochondria, plastids, and obligate intracellular parasites. The mitochondrial ATP/ADP transporters are derived from a broad-specificity transport family of eukaryotic origin, whereas the origin of the plastid/parasite ATP/ADP translocase is more elusive. Here we present the sequences of five genes coding for ATP/ADP translocases from four species of Rickettsia. The results are consistent with an early duplication and divergence of the five ATP/ADP translocases within the rickettsial lineage. A comparison of the phylogenetic depths of the mitochondrial and the plastid/parasite ATP/ADP translocases indicates a deep origin for both transporters. The results provide no evidence for a recent acquisition of the ATP/ADP transporters in Rickettsia via horizontal gene transfer, as previously suggested. A possible function of the two types of ATP/ADP translocases was to allow switches between glycolysis and aerobic respiration in the early eukaryotic cell and its endosymbiont.     

Foot differentiation and genomic plasticity in Hydra: lessons from the PPOD gene family

In Hydra, developmental processes are permanently active to maintain a simple body plan consisting of a two-layered, radially symmetrical tube with two differentiated structures, head and foot. Foot formation is a dynamic process and includes terminal differentiation of gastric epithelial cells into mucous secreting basal disc cells. A well-established marker for this highly specialized cell type is a locally expressed peroxidase (Hoffmeister et al. 1985). Based on the foot-specific peroxidase activity, the gene PPOD1 has been identified (Hoffmeister-Ullerich et al. 2002). Unexpectedly, this approach led to the identification of a second gene, PPOD2, with high sequence similarity to PPOD1 but a strikingly different expression pattern. Here, we characterize PPOD2 in more detail and show that both genes, PPOD1 and PPOD2, are members of a gene family with differential complexity and expression patterns in different Hydra species. At the genomic level, differences in gene number and structure within the PPOD gene family, even among closely related species, support a recently proposed phylogeny of the genus Hydra and point to unexpected genomic plasticity within closely related species of this ancient metazoan taxon. Electronic supplementary material Supplementary material is available in the online version of this article at