Sequence Evolution and Expression Regulation of Stress-Responsive Genes in Natural Populations of Wild Tomato

The wild tomato species Solanum chilense and S. peruvianum are a valuable non-model system for studying plant adaptation since they grow in diverse environments facing many abiotic constraints. Here we investigate the sequence evolution of regulatory regions of drought and cold responsive genes and their expression regulation. The coding regions of these genes were previously shown to exhibit signatures of positive selection. Expression profiles and sequence evolution of regulatory regions of members of the Asr (ABA/water stress/ripening induced) gene family and the dehydrin gene pLC30-15 were analyzed in wild tomato populations from contrasting environments. For S. chilense, we found that Asr4 and pLC30-15 appear to respond much faster to drought conditions in accessions from very dry environments than accessions from more mesic locations. Sequence analysis suggests that the promoter of Asr2 and the downstream region of pLC30-15 are under positive selection in some local populations of S. chilense. By investigating gene expression differences at the population level we provide further support of our previous conclusions that Asr2, Asr4, and pLC30-15 are promising candidates for functional studies of adaptation. Our analysis also demonstrates the power of the candidate gene approach in evolutionary biology research and highlights the importance of wild Solanum species as a genetic resource for their cultivated relatives.     

Contrasting Patterns of Sequence Evolution at the Functionally Redundant <Emphasis Type="Italic">bric à brac</Emphasis> Paralogs in <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

Genes with overlapping expression and function may gradually diverge despite retaining some common functions. To test whether such genes show distinct patterns of molecular evolution within species, we examined sequence variation at the bric à brac (bab) locus of Drosophila melanogaster. This locus is composed of two anciently duplicated paralogs, bab1 and bab2, which are involved in patterning the adult abdomen, legs, and ovaries. We have sequenced the 148 kb genomic region spanning the bab1 and bab2 genes from 94 inbred lines of D. melanogaster sampled from a single location. Two non-coding regions, one in each paralog, appear to be under selection. The strongest evidence of directional selection is found in a region of bab2 that has no known functional role. The other region is located in the bab1 paralog and is known to contain a cis-regulatory element that controls sex-specific abdominal pigmentation. The coding region of bab1 appears to be under stronger functional constraint than the bab2 coding sequences. Thus, the two paralogs are evolving under different selective regimes in the same natural population, illuminating the different evolutionary trajectories of partially redundant duplicate genes.     

Asr genes belong to a gene family comprising at least three closely linked loci on chromosome 4 in tomato

Asr1, Asr2 andAsr3 are three homologous clones isolated from tomato whose expression is believed to be regulated by abscisic acid (ABA); the corresponding genes thus participate in physiological and developmental processes such as responses of leaf and root to water stress, and fruit ripening. In this report, results obtained with Near Isogenic Lines reveal thatAsr1, Asr2 andAsr3 represent three different loci. In addition, we map these genes on the restriction fragment length polymorphism (RFLP) map of the tomato genome by using an F₂ population derived from an interspecific hybrid crossL. esculentum × L. penelli. RFLP data allow us to map these genes on chromosome 4, suggesting that they belong to a gene family. The elucidation of the genomic organization of theAsr gene family may help in understanding the role of its members in the response to osmotic stress, as well as in fruit ripening, at the molecular level.     

Biochemical variation to determine phylogenetic relationships betweenHordeum chilense and other American species of the genusHordeum (Poaceae)

Slab gel electrophoresis techniques have been applied to the study of isozyme and kernel protein patterns in 20 accessions ofHordeum chilense and related species in order to elucidate their phylogenetic relationships. On the basis of our results we can conclude that: (1) Conventional classification based on morphological characters does not totally agree with biochemical data. (2) Sectt.Anisolepis andCritesion seem to be clearly differentiated. (3) The accessions classified asH. compressum present biochemical phenotypes quite different from the rest of the species. (4)H. stenostachys, H. muticum andH. chilense constitute a group of variable species with many biochemical similarities and close phylogenetic relationships. (5) The evolutionary pattern of these American species seems to follow a model of reticulate evolution.     

Evolutionary history of an MHC gene in two leporid species: characterisation of <Emphasis Type="Italic">Mhc-DQA</Emphasis> in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning–sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. The nucleotide sequence data reported in this paper have been submitted to Genbank and have been assigned the accession numbers FJ225335–FJ225346.     

Comparisons of the hybrids Hordeum chilensexH. vulgare,H. chilensexH. bulbosum,H. chilensexSecale cereale and the amphidiploid of H. chilensexH. vulgare

Summary Diploid hybrids between Hordeum chilense and three other species, namely H. vulgare, H. bulbosum and Secale cereale, are described together with the amphidiploid of H. chilensexH. vulgare. Both the diploid hybrid and the amphidiploid of H. chilensexH. vulgare were chromosomally unstable, H. chilensexH. bulbosum was less so, while H. chilensexS. cereale was stable. Differential amphiplasty was found in all combinations. No homoeologous pairing was found in the Hordeum hybrids but in H. chilensexS. cereale there was chromosome pairing both within the two genomes and between the genomes.     

First evidence of methylation in the genome of <Emphasis Type="Italic">Drosophila willistoni</Emphasis>

DNA methylation has been studied abundantly in vertebrates and recent evidence confirms that this phenomenon could be disseminated among some invertebrates groups, including Drosophila species. In this paper, we used the Methylation-Sensitive Restriction Endonuclease (MSRE) technique and Southern blot with specific probes, to detect methylation in the Drosophila willistoni species. We found differential cleavage patterns between males and females that cannot be explained by Mendelian inheritance, pointing to a DNA methylation phenomenon different from the Drosophila melanogaster one. The sequencing of some of these bands showed that these fragments were formed by different DNA elements, among which rDNA. We also characterized the D. willitoni dDnmt2 sequence, through a Mega Blast search against the D. willistoni Trace Archive Database using the D. melanogaster dDnmt2 nucleotide sequence as query. The complete analysis of D. willistoni dDnmt2 sequence showed that its promoter region is larger, its dDnmt2 nucleotide sequence is 33% divergent from the D. melanogaster one, Inverted Terminal Repeats (ITRs) are absent and only the B isoform of the enzyme is produced. In contrast, ORF2 is more conserved. Comparing the D. willistoni and D. melanogaster dDnmt2 protein sequences, we found higher conservation in motifs from the large domain, responsible for the catalysis of methyl transfer, and great variability in the region that carries out the recognition of specific DNA sequences (TRD). Globally, our results reveal that methylation of the D. willistoni genome could be involved in a singular process of species-specific dosage compensation and that the DNA methylation in the Drosophila genus can have diverse functions. This could be related to the evolutionary history of each species and also to the acquisition time of the dDnmt2 gene.     

Nucleotide substitutions in <Emphasis Type="Italic">rolC</Emphasis> and <Emphasis Type="Italic">nptII</Emphasis> gene sequences during long-term cultivation of <Emphasis Type="Italic">Panax ginseng</Emphasis> cell cultures

It has been shown previously that the rolC gene from Agrobacterium tumefaciens gene was stably and highly expressed in 15-year-old Panax ginseng transgenic cell cultures. In the present report, we analyze in detail the nucleotide composition of the rolC and nptII (neomycin phosphotransferase) genes, which is the selective marker used for transgenic cell cultures of P. ginseng. It has been established that the nucleotide sequences of the rolC and nptII genes underwent mutagenesis during cultivation. Particularly, 1–4 nucleotide substitutions were found per sequence in the 540 and 798 bp segments of the complete rolC and nptII genes, respectively. Approximately half of these nucleotide substitutions caused changes in the structure of the predicted gene product. In addition, we attempted to determine the rate of accumulation of these changes by comparison of DNA extracted from P. ginseng cell cultures from 1995 to 2007. It was observed that the frequency of nucleotide substitutions for the rolC and nptII genes in 1995 was 1.21 ± 0.02 per 1,000 nucleotides analyzed, while in 2007, the nucleotide substitutions significantly increased (1.37 ± 0.07 per 1,000 nucleotides analyzed). Analyzing the nucleotide substitutions, we found that substitution to G or to C nucleotides significantly increased (in 1.9 times) in the rolC and nptII genes compared with P. ginseng actin gene. Finally, the level of nucleotide substitutions in the rolC gene was 1.1-fold higher when compared with the nptII gene. Thus, for the first time, we have experimentally demonstrated the level of nucleotide substitutions in transferred genes in transgenic plant cell cultures.     

Nucleotide diversity patterns of local adaptation at drought‐related candidate genes in wild tomatoes

We surveyed nucleotide diversity at two candidate genes LeNCED1 and pLC30‐15, involved in an ABA (abscisic acid) signalling pathway, in two closely related tomato species Solanum peruvianum and Solanum chilense. Our six population samples (three for each species) cover a range of mesic to very dry habitats. The ABA pathway plays an important role in the plants’ response to drought stress. LeNCED1 is an upstream gene involved in ABA biosynthesis, and pLC30‐15 is a dehydrin gene positioned downstream in the pathway. The two genes show very different patterns of nucleotide variation. LeNCED1 exhibits very low nucleotide diversity relative to the eight neutral reference loci that were previously surveyed in these populations. This suggests that strong purifying selection has been acting on this gene. In contrast, pLC30‐15 exhibits higher levels of nucleotide diversity and, in particular in S. chilense, higher genetic differentiation between populations than the reference loci, which is indicative of local adaptation. In the more drought‐tolerant species S. chilense, one population (from Quicacha) shows a significant haplotype structure, which appears to be the result of positive (diversifying) selection.     

The sequence of the isoepoxydon dehydrogenase gene of the patulin biosynthetic pathway in <Emphasis Type="Italic">Penicillium</Emphasis> species

Interest in species of the genus Penicillium is related to their ability to produce the mycotoxin patulin and to cause spoilage of fruit products worldwide. The sequence of the isoepoxydon dehydrogenase (idh) gene, a gene in the patulin biosynthetic pathway, was determined for 28 strains representing 12 different Penicillium species known to produce the mycotoxin patulin. Isolates of Penicillium carneum, Penicillium clavigerum, Penicillium concentricum, Penicillium coprobium, Penicillium dipodomyicola, Penicillium expansum, Penicillium gladioli, Penicillium glandicola, Penicillium griseofulvum, Penicillium paneum, Penicillium sclerotigenum and Penicillium vulpinum were compared. Primer pairs for DNA amplification and sequencing were designed from the P. griseofulvum idh gene (GenBank AF006680). The two introns present were removed from the nucleotide sequences, which were translated to produce the IDH sequences of the 12 species for comparison. Phylogenetic relationships among the species were determined from rDNA (ITS1, 5.8 S, ITS2 and partial sequence of 28S rDNA) and from the idh nucleotide sequences minus the two introns. Maximum parsimony analysis showed trees based on rDNA and idh sequences to be congruent. It is anticipated that the genetic information obtained in the present study will aid in the design of probes, specific for patulin biosynthetic pathway genes, to identify the presence of these mycotoxigenic fungi. The U.S. Government's right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.     

Analysis of nucleotide sequence polymorphism of internal transcribed spacers of ribosomal genes in diploid <Emphasis Type="Italic">Aegilops</Emphasis> (L.) species

The nucleotide sequence of the fragment of the internal transcribed spacer (ITS) of rDNA comprising the full-length ITS1, the gene encoding 5.8S rRNA, and part of the ITS2 sequence was determined in 22 samples of five diploid Aegilops species. The full alignment length of compared sequences was 524 bp. Species-specific substitutions were found in the ITS nucleotide sequence of rDNA of different Aegilops species. Intraspecific differences in ITS structure in diploid Aegilops species were detected for the first time. Polymorphism of the ITS nucleotide sequence within the same sample was revealed, which might be due either to differences between the genomes of individual plants comprising the sample or to the presence of several types of ribosomal genes in the genome of one plant. In general, both interspecific and intraspecific variability of the ITS nucleotide sequences of rDNA is extremely low. In total, 26 variable sites, twelve of which were informative, were identified.__________Translated from Molekulyarnaya Biologiya, Vol. 39, No. 2, 2005, pp. 193–197.Original Russian Text Copyright © 2005 by Goryunova, Chikida, Gori, Kochieva.     

New species and a new combination in the <Emphasis Type="Italic">Hyphopichia</Emphasis> and <Emphasis Type="Italic">Yarrowia</Emphasis> yeast clades

Three new species of Candida and a new combination in the genus Hyphopichia are proposed from phylogenetic analysis of nucleotide divergence in domains D1/D2 of the large subunit (26S) rDNA. The new taxa and their type strains are the following: Candida bentonensis sp. nov. (NRRL YB-2364, CBS 9994), Candida hispaniensis sp. nov. (NRRL Y-5580, CBS 9996), Candida pseudorhagii sp. nov. (NRRL YB-2076, CBS 9998) and Hyphopichia heimii comb. nov. (NRRL Y-7502, CBS 6139), basionym Pichia heimii Pignal. Phylogenetic analysis placed C. pseudorhagii and H. heimii in the Hyphopichia clade whereas C. bentonensis and C. hispaniensis are members of the Yarrowia clade.     

<Emphasis Type="Italic">hedgehog</Emphasis> is a segment polarity gene in a crustacean and a chelicerate

The evolution of arthropod segmentation has been studied by comparing expression patterns of pair-rule and segment polarity genes in various species. In Drosophila, the formation and maintenance of the parasegmental boundaries depend on the interactions between the wingless (wg), engrailed (en) and hedgehog (hh) genes. Until now, the expression pattern of hh has not been analysed to such a great extent as en or wg. We report the cloning and expression analysis of hh genes from Euscorpius flavicaudis, a chelicerate, and Artemia franciscana, a branchiopod crustacean. Our data provide evidence that hh, being expressed in the posterior part of every segment, is a segment polarity gene in both organisms. Additional hh expression sites were observed in the rostrum and appendages of Euscorpius and in the gut of Artemia. From the available data on hh expression in various bilaterians, we review the various hypotheses on the evolution of hh function and we suggest an ancestral role of hh in proctodeum specification and gut formation.Edited by D. Tautz     

Phylogeography and population structure of the endangered Tehuantepec jackrabbit <Emphasis Type="Italic">Lepus flavigularis</Emphasis>: implications for conservation

The Tehuantepec jackrabbit (Lepus flavigularis) is an endangered species restricted to a small area in the Isthmus of Tehuantepec, Oaxaca, Mexico. To evaluate its phylogeographic structure, population genetics, and demographic history we sequenced the mitochondrial Control Region hypervariable domain (CR-1) for 42 individuals representing the entire species range. Phylogenetic patterns indicated that this species is subdivided into two highly divergent clades, with an average nucleotide genetic distance of 3.7% (TrN) between them. Clades A and B are geographically distributed in non-overlapping areas to the west and to the east of the Isthmus of Tehuantepec, respectively. Genetic diversity indices showed reduced genetic variability in L. flavigularis when compared to other species of Lepus within main clades and within populations. This low genetic diversity coupled with the restricted distribution to very small areas of occurrence and limited gene flow suggest that genetic drift has played an important role in the evolution of this species. Historical demographic analysis also pointed out that these two clades underwent a recent population expansion that started about 9,000 years ago for clade A and 3,200 years ago for clade B during the Holocene. Consequently, from the conservation perspective our results suggest that populations included in clades A and B should be regarded as distinct evolutionary lineages.     

Isolation and characterization of genome-specific DNA sequences in Triticeae species

Two contrasting genome-specific DNA sequences were isolated from Aegilops speltoides (wild goat grass) and Hordeum chilense (wild barley), each representing more than 1 % of the genomes. These repetitive DNA fragments were identified as being genome-specific before cloning by genomic Southern hybridization (using total genomic DNA as a probe), and hence extensive screening of clones was not required. For each fragment, up to six recombinant plasmid clones were screened and about half were genome-specific. Clone pAesKB52 from Ae. speltoides was a 763 by EcoRI fragment, physically organized in simple tandem repeats and shown to localize to sub-telomerec chromosome regions of species with the Triticeae S-genome by in situ hybridization to chromosomes. The sequence data showed an internal duplication of some 280 bp, which presumably occurred before sequence amplification and dispersion, perhaps by unequal crossing-over or reciprocal translocation. In situ hybridization showed that the sequence distribution varied between closely related (S-genome) species. Clone pHcKB6 was a 339 by DraI fragment from H. chilense, also tandemly repeated but more variable; loss of the DraI site resulting in a ladder pattern in Southern blots which had little background smear. In situ hybridization showed that the tandem repeats were present as small clusters dispersed along all chromosome arms except at a few discrete regions including the centromeres and telomeres. The clone hybridized essentially specifically to the H-genome of H. chilense and hence was able to identify the origin of chromosomes in a H. chilense x Secale africanum hybrid by in situ hybridization. It has a high A + T content (66%), small internal duplications, and a 50 by degenerate inverted repeat. We speculate that it has dispersed by retrotransposition in association with other sequences carrying coding domains. The organization and evolution of such sequences are important in understanding long-range genome organization and the types of change that can occur on evolutionary and plant breeding timescales. Genome-specific sequences are also useful as markers for alien chromatin in plant breeding.     

<Emphasis Type="Italic">RNase</Emphasis>-Based Gametophytic Self-Incompatibility Evolution: Questioning the Hypothesis of Multiple Independent Recruitments of the <Emphasis Type="Italic">S</Emphasis>-Pollen Gene

Multiple independent recruitments of the S-pollen component (always an F-box gene) during RNase-based gametophytic self-incompatibility evolution have recently been suggested. Therefore, different mechanisms could be used to achieve the rejection of incompatible pollen in different plant families. This hypothesis is, however, mainly based on the interpretation of phylogenetic analyses, using a small number of divergent nucleotide sequences. In this work we show, based on a large collection of F-box S-like sequences, that the inferred relationship of F-box S-pollen and F-box S-like sequences is dependent on the sequence alignment software and phylogenetic method used. Thus, at present, it is not possible to address the phylogenetic relationship of F-box S-pollen and S-like sequences from different plant families. In Petunia and Malus/Pyrus the putative S-pollen gene(s) show(s) variability patterns different than expected for an S-pollen gene, raising the question of false identification. Here we show that in Petunia, the unexpected features of the putative S-pollen gene are not incompatible with this gene’s being the S-pollen gene. On the other hand, it is very unlikely that the Pyrus SFBB-gamma gene is involved in specificity determination. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Phylogeny of <Emphasis Type="Italic">Banana Streak Virus</Emphasis> Reveals Recent and Repetitive Endogenization in the Genome of Its Banana Host (<Emphasis Type="Italic">Musa</Emphasis> sp.)

Banana streak virus (BSV) is a plant dsDNA pararetrovirus (family Caulimoviridae, genus badnavirus). Although integration is not an essential step in the BSV replication cycle, the nuclear genome of banana (Musa sp.) contains BSV endogenous pararetrovirus sequences (BSV EPRVs). Some BSV EPRVs are infectious by reconstituting a functional viral genome. Recent studies revealed a large molecular diversity of episomal BSV viruses (i.e., nonintegrated) while others focused on BSV EPRV sequences only. In this study, the evolutionary history of badnavirus integration in banana was inferred from phylogenetic relationships between BSV and BSV EPRVs. The relative evolution rates and selective pressures (d_N/d_S ratio) were also compared between endogenous and episomal viral sequences. At least 27 recent independent integration events occurred after the divergence of three banana species, indicating that viral integration is a recent and frequent phenomenon. Relaxation of selective pressure on badnaviral sequences that experienced neutral evolution after integration in the plant genome was recorded. Additionally, a significant decrease (35%) in the EPRV evolution rate was observed compared to BSV, reflecting the difference in the evolution rate between episomal dsDNA viruses and plant genome. The comparison of our results with the evolution rate of the Musa genome and other reverse-transcribing viruses suggests that EPRVs play an active role in episomal BSV diversity and evolution.