mtDNA Diversity and genetic lineages of eighteen cattle breeds from <Emphasis Type="Italic">Bos taurus </Emphasis>and<Emphasis Type="Italic"> Bos indicus</Emphasis> in China

In order to clarify the origin and genetic diversity of indigenous cattle breeds in China, we carried out phylogenetic analysis of representatives of those breeds by employing mitochondrial gene polymorphism. Complete cyt b gene sequences, 1140 bp in length, were determined for a total of 136 individuals from 18 different breeds and these sequences were clustered into two distinct genetic lineages: taurine (Bos taurus) and zebu (Bos indicus). In analysis of the cyt b gene diversity, Chinese cattle showed higher nucleotide (0.00923) and haplotype diversity (0.848) than the reports from other studies, and the animals from the taurine lineage indicated higher nucleotide diversity (0.00330) and haplotype diversity (0.746) than the ones from the zebu lineage (0.00136; 0.661). The zebu mtDNA dominated in the southern breeds (63.3–100%), while the taurine dominated in the northern breeds (81.8–100%). Six cattle breeds from the central area of China exhibited intermediate frequencies of zebu mtDNA (25–71.4%). This polymorphism revealed a declining south-to-north gradient of female zebu introgression and a geographical hybrid zone of Bos taurus and Bos indicus in China.     

Phylogeny of members of the <Emphasis Type="Italic">Frankia</Emphasis> genus based on <Emphasis Type="Italic">gyr</Emphasis>B, <Emphasis Type="Italic">nif</Emphasis>H and <Emphasis Type="Italic">gln</Emphasis>II sequences

To construct an evolutionary hypothesis for the genus Frankia, gyrB (encoding gyrase B), nifH (encoding nitrogenase reductase) and glnII (encoding glutamine synthetase II) gene sequences were considered for 38 strains. The overall clustering pattern among Frankia strains based on the three analyzed sequences varied among themselves and with the previously established 16S rRNA gene phylogeny and they did not reliably reflect clear evolution of the four discerned Frankia clusters (1, 2, 3 and 4). Based on concatenated gyrB, nifH and glnII, robust phylogenetic trees were observed with the three treeing methods (Maximum Likelihood, Parsimony and Neighbor-Joining) and supported by strong bootstrap and posterior probability values (>75%) for overall branching. Cluster 4 (non-infective and/or non-nitrogen-fixing Frankia) was positioned at a deeper branch followed by cluster 3 (Rhamnaceae and Elaeagnaceae infective Frankia), while cluster 2 represents uncultured Frankia microsymbionts of the Coriariaceae, Datiscaceae, Rosaceae and of Ceanothus sp. (Rhamnaceae); Cluster 1 (Betulaceae, Myricaceae and Casuarinaceae infective Frankia) appears to have diverged more recently. The present study demonstrates the utility of phylogenetic analyses based upon concatenated gyrB, nifH and glnII sequences to help resolve previously unresolved or poorly resolved nodes and will aid in describing species among the genus Frankia.     

Phylogenetic Studies of <Emphasis Type="Italic">Gordonia</Emphasis> Species Based on <Emphasis Type="Italic">gyrB</Emphasis> and <Emphasis Type="Italic">secA1</Emphasis> Gene Analyses

Phylogenetic relations within the genus Gordonia were analyzed using partial gyrB and secA1 gene sequences of 23 type species in comparison with those of 16S rRNA gene. The gyrB and secA1 phylogenies showed agreement with that constructed using 16S rRNA gene sequences. The degrees of divergence of the gyrB and secA1 genes were approximately 3.4 and 1.7 times greater, respectively, than that of 16S rRNA gene. The gyrB gene showed more discriminatory power than either the secA1 or 16S rRNA gene, facilitating clear differentiation of any two Gordonia species using gyrB gene analysis. Our data indicate that gyrB and secA1 gene sequences are useful as markers for phylogenetic study and identification at the species level of the genus Gordonia.     

Application of <Emphasis Type="Italic">gyrB</Emphasis> and <Emphasis Type="Italic">parE</Emphasis> sequence similarity analyses for differentiation of species within the genus <Emphasis Type="Italic">Geobacillus</Emphasis>

The primary structures of the genes encoding the β-subunits of a type II topoisomerase (gyrase, gyrB) and a type IV topoisomerase (parE) were determined for 15 strains of thermophilic bacteria of the genus Geobacillus. The obtained sequences were used for analysis of the phylogenetic similarity between members of this genus. Comparison of the phylogenetic trees of geobacilli constructed on the basis of the 16S rRNA, gyrB, and parE gene sequences demonstrated that the level of genetic distance between the sequences of the genes encoding the β-subunits of type II topoisomerases significantly exceeded the values obtained by comparative analysis of the 16S rRNA gene sequences of Geobacillus strains. It was shown that, unlike the 16S rRNA gene analysis, comparative analysis of the gyrB and parE gene sequences provided a more precise determination of the phylogenetic position of bacteria at the species level. The data obtained suggest the possibility of using the genes encoding the β-subunits of type II topoisomerases as phylogenetic markers for determination of the species structure of geobacilli.     

Polymorphism analysis of <Emphasis Type="Italic">csd</Emphasis> gene in six <Emphasis Type="Italic">Apis mellifera</Emphasis> subspecies

The complementary sex determination (csd) gene is the primary gene determining the gender of honey bees (Apis spp). In this study we analyzed the polymorphism of csd gene in six Apis mellifera subspecies. The genomic region 3 of csd gene in these six A. mellifera was cloned, and identified. A total of 79 haplotypes were obtained from these six subspecies. Analysis showed that region 3 of csd gene has a high level of polymorphism in all the six A. mellifera subspecies. The A. m. anatolica subspecies has a slightly higher nucleotide diversity (π) than other subspecies, while the π values showed no significant difference among the other five subspecies. The phylogenetic tree showed that all the csd haplotypes from different A. mellifera subspecies are scattered throughout the tree, without forming six different clades. Population differentiation analysis showed that there are significant genetic differentiations among some of the subspecies. The NJ phylogenetic tree showed that the A. m. caucasica and A. m. carnica have the closest relationship, followed by A. m. ssp, A. m. ligustica, A. m. carpatica and A. m. anatolica that were gathered in the tree in turn.     

Phylogenetic study of Oryzoideae species and related taxa of the Poaceae based on <Emphasis Type="Italic">atpB</Emphasis>-<Emphasis Type="Italic">rbcL</Emphasis> and <Emphasis Type="Italic">ndhF</Emphasis> DNA sequences

Oryzoideae (Poaceae) plants have economic and ecological value. However, the phylogenetic position of some plants is not clear, such as Hygroryza aristata (Retz.) Nees. and Porteresia coarctata (Roxb.) Tateoka (syn. Oryza coarctata). Comprehensive molecular phylogenetic studies have been carried out on many genera in the Poaceae. The different DNA sequences, including nuclear and chloroplast sequences, had been extensively employed to determine relationships at both higher and lower taxonomic levels in the Poaceae. Chloroplast DNA ndhF gene and atpB-rbcL spacer were used to construct phylogenetic trees and estimate the divergence time of Oryzoideae, Bambusoideae, Panicoideae, Pooideae and so on. Complete sequences of atpB-rbcL and ndhF were generated for 17 species representing six species of the Oryzoideae and related subfamilies. Nicotiana tabacum L. was the outgroup species. The two DNA datasets were analyzed, using Maximum Parsimony and Bayesian analysis methods. The molecular phylogeny revealed that H. aristata (Retz.) Nees was the sister to Chikusichloa aquatica Koidz. Moreover, P. coarctata (Roxb.) Tateoka was in the genus Oryza. Furthermore, the result of evolution analysis, which based on the ndhF marker, indicated that the time of origin of Oryzoideae might be 31 million years ago.     

Association between expression levels and growth trait-related SNPs located in promoters of the <Emphasis Type="Italic">MC4R</Emphasis> and <Emphasis Type="Italic">MSTN</Emphasis> genes in <Emphasis Type="Italic">Spinibarbus hollandi</Emphasis>

Melanocortin 4 receptor: (MC4R) and Myostatin (MSTN) are two important growth trait-related genes in animals. In this study, we showed that two SNPs, MC4R-719A>G and MSTN-519C>T, found in the promoters of the MC4R and MSTN genes, respectively, are both associated with growth traits in Spinibarbus hollandi. Furthermore, we observed that there were significant associations between the expression levels of the MC4R and MSTN genes and these two growth trait-related SNPs. The expression level of MC4R gene in brain was lower in GG genotype fish with extremely high growth performance than that in AA genotype fish with extremely low growth performance. Expression level of the MSTN gene in muscle was lower in TT genotype fish with extremely high growth performance than that in CC and CT genotype fish with lower growth performance. The results indicated that these SNPs located in the promoters of MC4R and MSTN are associated with growth-related traits through modification of gene expression levels. The MSTN and MC4R SNPs may have useful application in effective marker-assisted selection aimed to increase output in S. hollandi.     

Extreme conservation of non-repetitive non-coding regions near <Emphasis Type="Italic">HoxD</Emphasis>complex of vertebrates

Homeotic gene complexes determine the anterior-posterior body axis in animals. The expression pattern and function of hox genes along this axis is colinear with the order in which they are organized in the complex. This 'chromosomal organization and functional correspondence' is conserved in all bilaterians investigated. Although the molecular basis of this 'colinearity' in not yet understood, it is possible that there are control elements within or in the proximity of these complexes that establish and maintain the expression patterns of hox genes in a coordinated fashion. We report here an unprecedented conservation of non-coding DNA sequences adjacent to the HoxD complex of vertebrates. Stretches of hundreds of base pairs in a 7 kb region, upstream of HoxD complex, show 100% conservation from fish to human. Using primers designed from these sequences of human HoxD complex, we amplified the corresponding regions from different vertebrates, including mammals, aves, reptiles, amphibians and pisces. Such a high degree of conservation, where no variation was allowed during ~500 million years of evolution, suggests critical function for these sequences in the regulation of the HoxD complex. Furthermore, these sequences provide a molecular handle to gain insight into the mechanism of regulation of this complex.     

Arbuscular mycorrhizal fungi in roots of non-photosynthetic plants, <Emphasis Type="Italic">Sciaphila japonica</Emphasis> and <Emphasis Type="Italic">Sciaphila tosaensis</Emphasis> (Triuridaceae)

The mycorrhizal fungi in the roots of achlorophyllous Sciaphila japonica and S. tosaensis (Triuridaceae) were identified by molecular methods. The habitats of S. japonica were in a tree plantation of Japanese cypress, Chamaecyparis obtusa, and bamboo forests, and those of S. tosaensis were in a camellia forest and a bamboo forest. In the root cortical cells of both plants, aseptate hyphal coils were observed, which suggested the Paris-type arbuscular mycorrhiza (AM). A phylogenetic analysis based on a partial sequence of an AM fungal nuclear small subunit ribosomal RNA gene showed that the fungal DNA sequences of S. japonica were separated into three closely related clades. Those of S. tosaensis were separated into two clades, which were also closely related to each other. The AM fungi of S. japonica and S. tosaensis were completely separated in the phylogenetic tree even among those found in the same habitat, which suggests the high specificities in the plant-fungal partnerships. All the detected AM fungi in these plants belonged to Glomus-group A. Even though the habitats are in quite common environments, both plant species are known as endangered species in Japan. Such a definite specificity in AM symbioses seems to restrict the distribution of the myco-heterotrophic plants.     

Gene Sequence Phylogenies of the Family <Emphasis Type="Italic">Microbacteriaceae</Emphasis>

The type strains of 32 species of 13 genera of the family Microbacteriaceae were analysed with respect to gene-coding phylogeny for DNA gyrase subunit B (gyrB), RNA-polymerase subunit B (rpoB), recombinase A (recA), and polyphosphate kinase (ppk). The resulting gene trees were compared with the 16S rRNA gene phylogeny of the same strains. The topology of neighbour-joining and maximum parsimony phylogenetic trees, based on nucleic-acid sequences and protein sequences of housekeeping genes, differed from one another, and no gene tree was identical to that of the 16S rRNA gene tree. Most genera analysed containing >1 strain formed phylogenetically coherent taxa. The three pathovars of Curtobacterium flaccumfaciens clustered together to the exclusion of the type strains of other Curtobacterium species in all DNA - and protein-based analyses. In no tree did the distribution of a major taxonomic marker, i.e., diaminobutyric acid versus lysine and/or ornithine in the peptidoglycan, or acyl type of peptidoglycan, correlate with the phylogenetic position of the organisms. The changing phylogenetic position of Agrococcus jenensis was unexpected: This strain defined individual lineages in the trees based on 16S rRNA and gyrB and showed identity with Microbacterium saperdae in the other three gene trees.     

<Emphasis Type="Italic">Wolbachia</Emphasis> Infections of the Whitefly <Emphasis Type="Italic">Bemisia tabaci</Emphasis>

We report the first systematic survey for the presence of Wolbachia endosymbionts in aphids and whiteflies, particularly different populations and biotypes of Bemisia tabaci. Additional agriculturally important species included were predator species, leafhoppers, and lepidopterans. We used a polymerase chain reaction (PCR)-based detection assay with ribosomal 16S rDNA and Wolbachia cell surface protein (wsp) gene primers. Wolbachia were detected in a number of whitefly populations and species, whitefly predators, and one leafhopper species; however, none of the aphid species tested were found infected. Single, double, and triple infections were detected in some of the B. tabaci populations. PCR and phylogenetic analysis of wsp gene sequences indicated that all Wolbachia strains found belong to group B. Topologies of the optimal tree derived by maximum likelihood (ML) and a ML tree in which Wolbachia sequences from B. tabaci are constrained to be monophyletic are significantly different. Our results indicate that there have been at least four independent Wolbachia infection events in B. tabaci. The importance of the presence of Wolbachia infections in B. tabaci is discussed. RID= ID= <E5>Correspondence to: </E5>K. Bourtzis; <E5>email:</E5> kbourtz&commat;cc.uoi.gr Received: 9 September 2002 / Accepted: 25 September 2002     

Characterization of the horse (<Emphasis Type="Italic">Equus caballus</Emphasis>)<Emphasis Type="Italic"> IGHA</Emphasis> gene

Nucleotide sequences of the immunoglobulin constant heavy chain genes of the horse have been described for IGHM, IGHG and IGHE genes, but not for IGHA. Here, we provide the nucleotide sequence of the genomic IGHA gene of the horse (Equus caballus), including its secretion region and the transmembrane exon. The equine IGHA gene shows the typical structure of a mammalian IGHA gene, with only three exons, separated by two introns of similar size. The hinge exon is located at the 5 end of the CH2 exon and encodes a hinge region of 11 amino acids, which contains five proline residues. The coding nucleotide sequence of the secreted form of the equine IGHA gene shares around 72% identity with the human IGHA1 and IGHA2 genes, as well as the bovine, ovine, porcine and canine IGHA genes, without distinct preference for any of these species. The same species also cluster together in a phylogenetic tree of the IGHA coding regions of various mammals, whereas rodent, rabbit, marsupial and monotreme IGHA genes each build a separate cluster.The nucleotide sequences reported in this paper have been assigned the EMBL/GenBank accession numbers AY247966 and AY351982     

Phylogeny of Saururaceae based on mitochondrial <Emphasis Type="Italic">matR</Emphasis> gene sequence data

DNA sequences of matR gene from three species of Saururaceae and the selected outgroups, Chloranthus holostegius and Zippelia begoniaefolia, are reported. All DNA sequences of six species in four genera of Saururaceae and the two outgroups are analyzed on PAUP 4.0 8b to reconstruct the phylogeny. A single matR gene tree is generated from parsimony, distance, and likelihood analyses, respectively. The three trees with the same topology are slightly different in bootstrapping support for some clades. The result indicates that Saururaceae is monophyletic. Anemopsis is sister to Houttuynia, and the two genera form the first diverging lineage of the family. The sister group relationship between Saururus and Gymnotheca is also supported by a relatively high bootstrap value. The result is different from all the former phylogenetic opinions on Saururaceae based on morphology, but it is supported by the evolution of flower-bract stalk in Saururaceae. In addition, some characteristics of the matR gene are analyzed. The MatR gene is a relatively better tool to reconstruct the molecular clock because the base substitution bias greatly decreases in the gene. Received: January 19, 2001 / Accepted: October 22, 2001     

Polymorphism and selection in the major histocompatibility complex <Emphasis Type="Italic">DRA</Emphasis> and <Emphasis Type="Italic">DQA</Emphasis> genes in the family Equidae

The major histocompatibility complex genes coding for antigen binding and presenting molecules are the most polymorphic genes in the vertebrate genome. We studied the DRA and DQA gene polymorphism of the family Equidae. In addition to 11 previously reported DRA and 24 DQA alleles, six new DRA sequences and 13 new DQA alleles were identified in the genus Equus. Phylogenetic analysis of both DRA and DQA sequences provided evidence for trans-species polymorphism in the family Equidae. The phylogenetic trees differed from species relationships defined by standard taxonomy of Equidae and from trees based on mitochondrial or neutral gene sequence data. Analysis of selection showed differences between the less variable DRA and more variable DQA genes. DRA alleles were more often shared by more species. The DQA sequences analysed showed strong amongst-species positive selection; the selected amino acid positions mostly corresponded to selected positions in rodent and human DQA genes.     

Isolation and sequence analysis of <Emphasis Type="Italic">Sox</Emphasis> genes from lizard <Emphasis Type="Italic">Eremias</Emphasis> multiocellata

The Sox (SRY-related high-mobility-group box) family of genes shares a conserved HMG box and is involved in a diverse range of developmental processes and sex determination in vertebrates. Twenty Sox genes are present in the genomes of humans and mice, but far less is known about the Sox gene family in reptiles. Using two pairs of highly degenerate primers designed from a multiple alignment of Sox amino acid sequences in several species, different positive clones were obtained from male and female Eremias multiocellata, a viviparous lizard which is subject to TSD (temperature-dependent sex determination). These clones were sequenced and identified. They are members of the SoxB (Sox2, Sox14), SoxC (Sox11, Sox12) and SoxE (Sox9a, Sox9b, Sox10) groups. No sex-specific differences were observed. Based on the amino acid sequence similarities, the phylogenetic analysis was carried out and these genes clustered with their orthologues. In addition, we found the gene duplication in E. multiocellata, it may be a mechanism to produce new functional genes.     

Genetic population structure of <Emphasis Type="Italic">Lefua echigonia</Emphasis> inferred from allozymic and mitochondrial cytochrome <Emphasis Type="Italic">b</Emphasis> variations

We examined the genetic population structure of Lefua echigonia (Japanese name, hotoke-dojo) using polymerase chain reaction restriction-fragment length polymorphisms (PCR-RFLPs) of the mitochondrial cytochrome b gene and two allozymic loci. The phylogenetic relationships of L. echigonia and those among L. echigonia, Lefua sp. (nagare hotoke-dojo), and L. nikkonis (ezo hotoke-dojo) were also investigated based on the nucleotide sequences of the cytochrome b gene. PCR-RFLP analysis revealed 18 mitotypes in L. echigonia, 2 in Lefua sp., and 1 in L. nikkonis. Phylogenetic trees based on the cytochrome b sequences indicated that the 18 mitotypes in L. echigonia were divided into five distinct groups (South-Kanto, North-Kanto, Tohoku, Echigo, and Tokai-Kinki clades) that differed by 8.5–15.3%, reflecting region-specific geographic distributions. The distributions of alleles in two allozymic loci roughly corresponded to those of the mitotype groups. The divergence times of the five groups were estimated to be about 3.4–7.7 million years ago by applying a general rate for mitochondrial DNA, suggesting that the divergence among them might have occurred in the late Tertiary. It can be inferred that the regional differentiation of each group was mainly due to geographic isolation and that this has been maintained, because the boundaries among the groups corresponded to geological features. The trees also supported the existence of three taxa, L. echigonia, Lefua sp., and L. nikkonis. We concluded that Lefua sp. was distinguished from other species in Lefua by morphological and ecological characters and also by genetic divergences of the cytochrome b gene. Our study also demonstrated the superior efficacy and simplicity of PCR-RFLP analysis as a method for detecting genetic variation in L. echigonia.     

Phylogenetic studies in <Emphasis Type="Italic">Ravenelia esculenta</Emphasis> and related rust fungi

Ravenelia esculenta Naras. and Thium. is a rust fungus, which infects mostly thorns, inflorescences, flowers and fruits of Acacia eburnea Willd. Aecial stages of the rust produce hypertrophy in infected parts. DNA of the rust fungus was isolated from aeciospores by ‘freeze thaw’ method. 18S rDNA was amplified and sequenced by automated DNA sequencer. BLAST of the sequence at NCBI retrieved 96 sequences producing significant alignments. Multiple sequence alignment of these sequences was done by ClustalW. Phylogenetic analysis was done by using MEGA 3.1. UPGMA Minimum Evolution tree with bootstrap value of 1000 replicates was constructed using these sequences. From phylogenetic tree it is observed that Ravenelia esculenta and the genus Gymnosporangium share a common ancestry, though Ravenelia esculenta is autoecious on angiosperm and the genus Gymnosporangium is heteroecious with pycnia, aecia on angiosperm and uredia, telia on gymnosperm. Two major clades are recognized which are based on the nature of aecial host (gymnosperm or angiosperm). These clades were also showing shift from pteridophytes to angiosperms as telial hosts. The tree can be interpreted in the other way also where there is separation of 14 families of Uredinales depending upon nature of teliospores, nature of aeciospores and structure of pycnia. These studies determine the phylogenetic position of Ravenelia esculenta among other rust fungi besides broad separation of Uredinales into two clades. These studies also show that there is phylogenetic correlation between molecular and morphological data. This is first report of DNA sequencing and phylogenetic positioning in genus Ravenelia from India.     

Molecular phylogeny of the blind cavefish <Emphasis Type="Italic">Phreatichthys andruzzii</Emphasis> and <Emphasis Type="Italic">Garra barreimiae</Emphasis> within the family Cyprinidae

The phylogenetic relationships of two cavefish, Phreatichthys andruzzii and Garra barreimiae, belonging to the family Cyprinidae, were investigated by sequencing the mitochondrial cytochrome b gene. These cavefish species are native to Somalia (eastern Africa) and Oman (southeastern Arabian peninsula), respectively, and so far no molecular support to their taxonomy and phylogenetic position was ever provided. The analysis of cytochrome b sequences showed that the species are monophyletic taxa, closely related to each other and to other species of the genus Garra. Molecular clock calculations allowed to date the origin of these hypogaean species back to the Plio-Pleistocene and support the hypothesis that African cyprinids originated from Miocenic immigrations of Asian ancestors.     

Molecular evolution of a family of resistance gene analogs of nucleotide-binding site sequences in <Emphasis Type="Italic">Solanum lycopersicum</Emphasis>

Nucleotide-binding site-leucine-rich repeats (NBS-LRR) gene families are one of the major plant resistance genes. Genomic NBS evolution was studied in many plant species for diverse arrays of NBS gene families. In this study, we focused on one family of NBS sequences in an attempt to understand how closely related NBS sequences evolved in the light of selection in domesticated plant species. A phylogenetic analysis revealed five major clades (A–E) and five subclades (A1–A5) within clade A of cloned NBS sequences. Positive selection was only detected in newly evolved NBS lineages in subclades of clade A. Positively selected codon sites were found among NBS sequences of clade A. A sliding-window analysis revealed that regions with Ka/Ks ratios of >1 were in the inter-motifs when paired clades were compared, but regions with Ka/Ks ratios of >1 were found across NBS sequences when subclades of clade A were compared. Our results based on a family of closely related NBS sequences showed that positive selection was first exerted on specific lineages across all NBS sequences after selective constraints. Subsequently, sequences with mutations in commonly conserved motifs were scrutinized by purifying selection. In the long term, conserved high frequency alleles in commonly conserved motifs and changes in inter-motifs were maintained in the investigated family of NBS sequences. Moreover, codons identified to be under positive selection in the inter-motifs were mainly located in regions involved in functions of ATP binding or hydrolysis.