miR-430 microRNA Family in Fishes: Molecular Characterization and Evolution

The miR-430 microRNA family has been described in multiple fish species as one of the first microRNAs expressed by the zygote. It has been suggested that this family is implicated in maternal mRNA elimination, but may also play a role in steroidogenesis, sexual differentiation, and flatfish metamorphosis. The miR-430 sequences have been found in multiple-copy tandem clusters but evidence of their conservation outside of teleost fishes is scarce. In the present study, we have characterized the tandem repeats organization of these microRNAs in different fish species, both model and of interest in aquaculture. A phylogenetic analysis of this family has allowed us to identify that the miR-430 duplication, which took place before the Chondrostei and Neopterygii groups' divergence, has resulted in three variants ("a", "b", and "c"). According to our data, variant "b" is the most closely related to the ancestral sequence. Furthermore, we have detected isolated instances of the miR-430 repeat subunit in some species, which suggests that this microRNA family may be affected by DNA rearrangements. This study provides new data about the abundance, variability, and organization of the miR-430 family in fishes.

Pistachio genomes provide insights into nut tree domestication and ZW sex chromosome evolution

Pistachio is a nut crop domesticated in the Fertile Crescent and a dioecious species with ZW sex chromosomes. We sequenced the genomes of Pistacia vera cultivar (cv.) Siirt, the female parent, and P. vera cv. Bagyolu, the male parent. Two chromosome-level reference genomes of pistachio were generated, and Z and W chromosomes were assembled. The ZW chromosomes originated from an autosome following the first inversion, which occurred approximately 8.18 Mya. Three inversion events in the W chromosome led to the formation of a 12.7-Mb (22.8% of the W chromosome) non-recombining region. These W-specific sequences contain several genes of interest that may have played a pivotal role in sex determination and contributed to the initiation and evolution of a ZW sex chromosome system in pistachio. The W-specific genes, including defA, defA-like, DYT1, two PTEN1, and two tandem duplications of six VPS13A paralogs, are strong candidates for sex determination or differentiation. Demographic history analysis of resequenced genomes suggest that cultivated pistachio underwent severe domestication bottlenecks approximately 7640 years ago, dating the domestication event close to the archeological record of pistachio domestication in Iran. We identified 390, 211, and 290 potential selective sweeps in 3 cultivar subgroups that underlie agronomic traits such as nut development and quality, grafting success, flowering time shift, and drought tolerance. These findings have improved our understanding of the genomic basis of sex determination/differentiation and horticulturally important traits and will accelerate the improvement of pistachio cultivars and rootstocks.

A chromosome-level genome assembly enables the identification of the follicule stimulating hormone receptor as the master sex-determining gene in the flatfish Solea senegalensis

Sex determination (SD) shows huge variation among fish and a high evolutionary rate, as illustrated by the Pleuronectiformes (flatfishes). This order is characterized by its adaptation to demersal life, compact genomes and diversity of SD mechanisms. Here, we assembled the Solea senegalensis genome, a flatfish of great commercial value, into 82 contigs (614 Mb) combining long- and short-read sequencing, which were next scaffolded using a highly dense genetic map (28,838 markers, 21 linkage groups), representing 98.9% of the assembly. Further, we established the correspondence between the assembly and the 21 chromosomes by using BAC-FISH. Whole genome resequencing of six males and six females enabled the identification of 41 single nucleotide polymorphism variants in the follicle stimulating hormone receptor (fshr) consistent with an XX/XY SD system. The observed sex association was validated in a broader independent sample, providing a novel molecular sexing tool. The fshr gene displayed differential expression between male and female gonads from 86 days post-fertilization, when the gonad is still an undifferentiated primordium, concomitant with the activation of amh and cyp19a1a, testis and ovary marker genes, respectively, in males and females. The Y-linked fshr allele, which included 24 nonsynonymous variants and showed a highly divergent 3D protein structure, was overexpressed in males compared to the X-linked allele at all stages of gonadal differentiation. We hypothesize a mechanism hampering the action of the follicle stimulating hormone driving the undifferentiated gonad toward testis.

Detection and variability analyses of CRISPR-like loci in the H. pylori genome

Helicobacter pylori is a human pathogenic bacterium with a high genomic plasticity. Although the functional CRISPR-Cas system has not been found in its genome, CRISPR-like loci have been recently identified. In this work, 53 genomes from different geographical areas are analyzed for the search and analysis of variability of this type of structure. We confirm the presence of a locus that was previously described in the VlpC gene in al lgenomes, and we characterize new CRISPR-like loci in other genomic locations. By studying the variability and gene location of these loci, the evolution and the possible roles of these sequences are discussed. Additionally, the usefulness of this type of sequences as a phylogenetic marker has been demonstrated, associating the different strains by geographical area.

Centromeric Satellite DNA in Flatfish (Order Pleuronectiformes) and Its Relation to Speciation Processes

Two new centromeric satellite DNAs in flatfish (Order Pleuronectiformes) have been characterized. The SacI-family from Hippoglossus hippoglossus, restricted to this species, had a monomeric size of 334 base pair (bp) and was located in most of the centromeres of its karyotype. The PvuII-family, with a monomeric size of 177 bp, was initially isolated from the genome of Solea senegalensis, and fluorescent in situ hybridization (FISH) localized the repeat to centromeres of most of the chromosomes. This family could only be amplified in 2 other species of the genus Solea (Solea solea and Solea lascaris). Molecular features and chromosomal location indicated a possible structural and/or functional role of these sequence repeats. The presence of species-specific satellite-DNA families in the centromeres and their possible role in the speciation processes in this group of fishes is discussed.

The Molecular Cytogenetic Characterization of Pistachio (Pistacia vera L.) Suggests the Arrest of Recombination in the Largest Heteropycnotic Pair HC1

This paper represents the first molecular cytogenetic characterization of the strictly dioecious pistachio tree (Pistacia vera L.). The karyotype was characterized by fluorescent in situ hybridization (FISH) with probes for 5S and 45S rDNAs, and the pistachio specific satellite DNAs PIVE-40, and PIVE-180, together with DAPI-staining. PIVE-180 has a monomeric unit of 176–178 bp and high sequence homology between family members; PIVE-40 has a 43 bp consensus monomeric unit, and is most likely arranged in higher order repeats (HORs) of two units. The P. vera genome is highly heterochromatic, and prominent DAPI positive blocks are detected in most chromosomes. Despite the difficulty in classifying chromosomes according to morphology, 10 out of 15 pairs (2n = 30) could be distinguished by their unique banding patterns using a combination of FISH probes. Significantly, the largest pair, designated HC1, is strongly heteropycnotic, shows differential condensation, and has massive enrichment in PIVE-40 repeats. There are two types of HC1 chromosomes (type-I and type-II) with differing PIVE-40 hybridization signal. Only type-I/II heterozygotes and type-I homozygotes individuals were found. We speculate that the differentiation between the two HC1 chromosomes is due to suppression of homologous recombination at meiosis, reinforced by the presence of PIVE-40 HORs and differences in PIVE-40 abundance. This would be compatible with a ZW sex-determination system in the pistachio tree.

First haploid genetic map based on microsatellite markers in Senegalese sole (Solea senegalensis, Kaup 1858)

The Senegalese sole (Solea senegalensis, Kaup 1858) is a flatfish species of great value for aquaculture. In this study, we develop the first linkage map in this species based on microsatellite markers characterized from genomic DNA libraries and EST databases of Senegalese sole and from other flatfish species. Three reference gynogenetic families were obtained by chromosome-manipulation techniques: two haploid gynogenetics, used to assign and order microsatellites to linkage groups and another diploid gynogenetic family, used for estimating marker-centromere distances. The consensus map consists of 129 microsatellites distributed in 27 linkage groups (LG), with an average density of 4.7 markers per LG and comprising 1,004 centimorgans (cM). Additionally, 15 markers remained unlinked. Through half-tetrad analysis, we were able to estimate the centromere distance for 81 markers belonging to 24 LG, representing an average of 3 markers per LG. Comparative mapping was performed between flatfish species LG and model fish species chromosomes (stickleback, Tetraodon, medaka, fugu and zebrafish). The usefulness of microsatellite markers and the genetic map as tools for comparative mapping and evolution studies is discussed.

Chromosomal manipulation in Senegalese sole (Solea senegalensis Kaup, 1858): induction of triploidy and gynogenesis

In this study we have developed protocols for induced triploidy and gynogenesis of Senegalese sole (Solea senegalensis), a promising flatfish species for marine aquaculture, in order to: 1) identify the sex-determination mechanism; and 2) to improve its production by generating a) sterile fish, avoiding problems related with sexual maturation, and b) all-female stocks, of higher growth rate. Triploidy was induced by means of a cold shock. Gynogenesis was induced by activating eggs with UV-irradiated sperm, and to prompt diploid gynogenesis, a cold-shock step was also used. Ploidy of putative triploid larvae and gynogenetic embryos were determined by means of karyotyping and microsatellite analysis. Haploid gynogenetic embryos showed the typical "haploid syndrome". As expected, triploid and gynogenetic groups showed lower fertilization, hatching, and survival rates than in the diploid control group. Survival rate, calculated 49 days after hatching, for haploid and diploid gynogenetic groups was similar to those observed in other fish species (0% and 62.5%, respectively), whereas triploids showed worse values (45%). Sex was determined macroscopically and by histological procedures, revealing that all the diploid gynogenetic individuals were females. In conclusion, we have successfully applied chromosomal-manipulation techniques in the flatfish species Senegalese sole in order to produce triploid, haploid, and diploid gynogenetic progenies.

Sequencing papaya X and Yh chromosomes reveals molecular basis of incipient sex chromosome evolution

Sex determination in papaya is controlled by a recently evolved XY chromosome pair, with two slightly different Y chromosomes controlling the development of males (Y) and hermaphrodites (Y(h)). To study the events of early sex chromosome evolution, we sequenced the hermaphrodite-specific region of the Y(h) chromosome (HSY) and its X counterpart, yielding an 8.1-megabase (Mb) HSY pseudomolecule, and a 3.5-Mb sequence for the corresponding X region. The HSY is larger than the X region, mostly due to retrotransposon insertions. The papaya HSY differs from the X region by two large-scale inversions, the first of which likely caused the recombination suppression between the X and Y(h) chromosomes, followed by numerous additional chromosomal rearrangements. Altogether, including the X and/or HSY regions, 124 transcription units were annotated, including 50 functional pairs present in both the X and HSY. Ten HSY genes had functional homologs elsewhere in the papaya autosomal regions, suggesting movement of genes onto the HSY, whereas the X region had none. Sequence divergence between 70 transcripts shared by the X and HSY revealed two evolutionary strata in the X chromosome, corresponding to the two inversions on the HSY, the older of which evolved about 7.0 million years ago. Gene content differences between the HSY and X are greatest in the older stratum, whereas the gene content and order of the collinear regions are identical. Our findings support theoretical models of early sex chromosome evolution.

An Expressed Sequence Tag (EST)-enriched genetic map of turbot (Scophthalmus maximus): a useful framework for comparative genomics across model and farmed teleosts

Background

The turbot (Scophthalmus maximus) is a relevant species in European aquaculture. The small turbot genome provides a source for genomics strategies to use in order to understand the genetic basis of productive traits, particularly those related to sex, growth and pathogen resistance. Genetic maps represent essential genomic screening tools allowing to localize quantitative trait loci (QTL) and to identify candidate genes through comparative mapping. This information is the backbone to develop marker-assisted selection (MAS) programs in aquaculture. Expressed sequenced tag (EST) resources have largely increased in turbot, thus supplying numerous type I markers suitable for extending the previous linkage map, which was mostly based on anonymous loci. The aim of this study was to construct a higher-resolution turbot genetic map using EST-linked markers, which will turn out to be useful for comparative mapping studies.

Results

A consensus gene-enriched genetic map of the turbot was constructed using 463 SNP and microsatellite markers in nine reference families. This map contains 438 markers, 180 EST-linked, clustered at 24 linkage groups. Linkage and comparative genomics evidences suggested additional linkage group fusions toward the consolidation of turbot map according to karyotype information. The linkage map showed a total length of 1402.7 cM with low average intermarker distance (3.7 cM; ~2 Mb). A global 1.6:1 female-to-male recombination frequency (RF) ratio was observed, although largely variable among linkage groups and chromosome regions. Comparative sequence analysis revealed large macrosyntenic patterns against model teleost genomes, significant hits decreasing from stickleback (54%) to zebrafish (20%). Comparative mapping supported particular chromosome rearrangements within Acanthopterygii and aided to assign unallocated markers to specific turbot linkage groups.

Conclusions

The new gene-enriched high-resolution turbot map represents a useful genomic tool for QTL identification, positional cloning strategies, and future genome assembling. This map showed large synteny conservation against model teleost genomes. Comparative genomics and data mining from landmarks will provide straightforward access to candidate genes, which will be the basis for genetic breeding programs and evolutionary studies in this species.

Validation and comparison of microsatellite markers derived from Senegalese sole (Solea senegalensis, Kaup) genomic and expressed sequence tags libraries

In this work, we tested 100 potential new microsatellites (SSRs) equally derived from expressed sequence tag (EST) and enriched genomic-DNA libraries from Senegalese sole (Solea senegalensis, Kaup), a valuable cultured flatfish species. A final set of 69 new polymorphic microsatellites were validated after a population analysis, 37 of which corresponded to the first EST library constructed for Senegalese sole (EST-SSR). Although differences were not significant, EST sequences provided a higher proportion of quality markers (74%) than anonymous ones (64%). Most of the rejected anonymous SSRs (17 loci) were discarded because they did not generate PCR products; only one was monomorphic. On the contrary, all EST-SSRs gave PCR products, although monomorphism was more frequent (26%). Altogether, the number of alleles per locus was fairly similar in both SSR types, ranging from 2 to 19. The observed and expected heterozygosities varied from 0.105 to 1 and from 0.108 to 0.937, respectively. The main difference between the two sets was the percentage of annotated loci, being higher in EST-SSRs, as expected. Within the EST-SSRs, 46% of them showed flanking regions that significantly matched with EST sequences from other three flatfish species; however, the microsatellite itself was present only on half of these cases. These two new SSR sets constitute a suitable tool for fingerprinting, gene flow, genetic diversity, genome mapping studies and molecular-assisted breeding in this species.

Construction of physical maps for the sex-specific regions of papaya sex chromosomes

Background

Papaya is a major fruit crop in tropical and subtropical regions worldwide. It is trioecious with three sex forms: male, female, and hermaphrodite. Sex determination is controlled by a pair of nascent sex chromosomes with two slightly different Y chromosomes, Y for male and Y^h for hermaphrodite. The sex chromosome genotypes are XY (male), XY^h (hermaphrodite), and XX (female). The papaya hermaphrodite-specific Y^h chromosome region (HSY) is pericentromeric and heterochromatic. Physical mapping of HSY and its X counterpart is essential for sequencing these regions and uncovering the early events of sex chromosome evolution and to identify the sex determination genes for crop improvement.

Results

A reiterate chromosome walking strategy was applied to construct the two physical maps with three bacterial artificial chromosome (BAC) libraries. The HSY physical map consists of 68 overlapped BACs on the minimum tiling path, and covers all four HSY-specific Knobs. One gap remained in the region of Knob 1, the only knob structure shared between HSY and X, due to the lack of HSY-specific sequences. This gap was filled on the physical map of the HSY corresponding region in the X chromosome. The X physical map consists of 44 BACs on the minimum tiling path with one gap remaining in the middle, due to the nature of highly repetitive sequences. This gap was filled on the HSY physical map. The borders of the non-recombining HSY were defined genetically by fine mapping using 1460 F₂ individuals. The genetically defined HSY spanned approximately 8.5 Mb, whereas its X counterpart extended about 5.4 Mb including a 900 Kb region containing the Knob 1 shared by the HSY and X. The 8.5 Mb HSY corresponds to 4.5 Mb of its X counterpart, showing 4 Mb (89%) DNA sequence expansion.

Conclusion

The 89% increase of DNA sequence in HSY indicates rapid expansion of the Y^h chromosome after genetic recombination was suppressed 2–3 million years ago. The genetically defined borders coincide with the common BACs on the minimum tiling paths of HSY and X. The minimum tiling paths of HSY and its X counterpart are being used for sequencing these X and Y^h-specific regions.

Exploitation of a turbot (Scophthalmus maximus L.) immune-related expressed sequence tag (EST) database for microsatellite screening and validation

In this study, we identified and characterized 160 microsatellite loci from an expressed sequence tag (EST) database generated from immune-related organs of turbot (Scophthalmus maximus). A final set of 83 new polymorphic microsatellites were validated after the analysis of 40 individuals of Atlantic origin including both wild and farmed individuals. The allele number and the expected heterozygosity ranged from 2 to 18 and from 0.021 to 0.951, respectively. Evidences of null alleles at moderate-high frequencies were detected at six loci using population data. None of the analysed loci showed deviations from Mendelian segregation after the analysis of five full-sib families including approximately 92 individuals/family. The markers are used to consolidate the turbot genetic map, and because they are mostly EST-derived, they will be very useful for comparative genomic studies within flatfishes and with model fish species. Using an in silico approach, we detected significant homologies of microsatellite sequences with the EST databases of the flatfish species with highest genomic resources (Senegalese sole, Atlantic halibut, bastard halibut) in 31% of these turbot markers. The conservation of these microsatellites within Pleuronectiformes will pave the way for anchoring genetic maps of different species and identifying genomic regions related to productive traits.

A satellite DNA evolutionary analysis in the North American endemic dioecious plant Rumex hastatulus (Polygonaceae)

We studied the evolution of RAE180 satellite DNA family in the North American endemic dioecious plant Rumex hastatulus. In this species, the Texas race is characterized by a single XX/XY sex chromosome system, whereas the North Carolina race has evolved a derived complex XX/XY(1)Y(2) sex chromosome system. RAE180 repeats were autosomic and poorly represented (2 × 10(-4)% of the genome) with no differences between individuals of different genders or different races of R. hastatulus. In fact, the sex chromosomes of the North Carolina race are still euchromatic, and they have not accumulated satellite DNA sequences, which contrasts with that occurring in the rest of dioecious XX/XY(1)Y(2) Rumex species. In R. hastatulus, we detected the existence of three RAE180 subfamilies. Notwithstanding, while in the Texas race the TX1/NC1 subfamily is the most frequent, the TX2/NC2 subfamily is the most abundant in the North Carolina race. Additionally, the third, less represented subfamily (TX3/NC3) appears currently as relict sequences in both genomes. A common feature of RAE180 satellite is the sudden replacement of one sequence variant by another in different species (or populations as in R. hastatulus races). Thus, the phylogenetic analysis of RAE180 repeats from six dioecious Rumex species supports the "library" hypothesis. According to this hypothesis, we assume that a set of divergent RAE180 variants were present in the ancestral genome of dioecious Rumex species, from which novel tandem arrays originated by the amplification of different variants in different lineages. Differential levels of RAE180 satellite DNA amplification in each lineage, at different evolutionary times, and in different chromosomal positions gave rise to differential patterns of sequence evolution.

Micro-collinearity and genome evolution in the vicinity of an ethylene receptor gene of cultivated diploid and allotetraploid coffee species (Coffea)

Arabica coffee (Coffea arabica L.) is a self-compatible perennial allotetraploid species (2n=4x=44), whereas Robusta coffee (C. canephora L.) is a self-incompatible perennial diploid species (2n=2x=22). C. arabica (C(a) C(a) E(a) E(a) ) is derived from a spontaneous hybridization between two closely related diploid coffee species, C. canephora (CC) and C. eugenioides (EE). To investigate the patterns and degree of DNA sequence divergence between the Arabica and Robusta coffee genomes, we identified orthologous bacterial artificial chromosomes (BACs) from C. arabica and C. canephora, and compared their sequences to trace their evolutionary history. Although a high level of sequence similarity was found between BACs from C. arabica and C. canephora, numerous chromosomal rearrangements were detected, including inversions, deletions and insertions. DNA sequence identity between C. arabica and C. canephora orthologous BACs ranged from 93.4% (between E(a) and C(a) ) to 94.6% (between C(a) and C). Analysis of eight orthologous gene pairs resulted in estimated ages of divergence between 0.046 and 0.665 million years, indicating a recent origin of the allotetraploid species C. arabica. Analysis of transposable elements revealed differential insertion events that contributed to the size increase in the C(a) sub-genome compared to its diploid relative. In particular, we showed that insertion of a Ty1-copia LTR retrotransposon occurred specifically in C. arabica, probably shortly after allopolyploid formation. The two sub-genomes of C. arabica, C(a) and E(a) , showed sufficient sequence differences, and a whole-genome shotgun approach could be suitable for sequencing the allotetraploid genome of C. arabica.

Molecular characterization and evolution of an interspersed repetitive DNA family of oysters

When genomic DNA from the European flat oyster Ostrea edulis L. was digested by BclI enzyme, a band of about 150 bp was observed in agarose gel. After cloning and sequencing this band and analysing their molecular characteristics and genomic organization by means of Southern blot, in situ hybridisation, and polymerase chain reaction (PCR) protocols, we concluded that this band is an interspersed highly repeated DNA element, which is related in sequence to the flanking regions of (CT)-microsatellite loci of the species O. edulis and Crassostrea gigas. Furthermore, we determined that this element forms part of a longer repetitive unit of 268 bp in length that, at least in some loci, is present in more than one copy. By Southern blot hybridisation and PCR amplifications-using primers designed for conserved regions of the 150-bp BclI clones of O. edulis-we determined that this repetitive DNA family is conserved in five other oyster species (O. stentina, C. angulata, C. gigas, C. ariakensis, and C. sikamea) while it is apparently absent in C. gasar. Finally, based on the analysis of the repetitive units in these oyster species, we discuss the slow degree of concerted evolution in this interspersed repetitive DNA family and its use for phylogenetic analysis.

Effect of location, organization, and repeat-copy number in satellite-DNA evolution

Here, we analyze the evolutionary dynamics of a satellite-DNA family in an attempt to understand the effect of factors such as location, organization, and repeat-copy number in the molecular drive process leading to the concerted-evolution pattern found in this type of repetitive sequences. The presence of RAE180 satellite-DNA in the dioecious species of the plant genus Rumex is a noteworthy feature at this respect, as RAE180 satellite repeats have accumulated differentially, showing a distinct distribution pattern in different species. The evolution of dioecious Rumex gave rise to two phylogenetic clades: one clade composed of species with an ancestral XX/XY sex chromosome system and a second, derived clade of species with a multiple sex-chromosome system XX/XY(1)Y(2). While in the XX/XY dioecious species, the RAE180 satellite-DNA is located only in a small autosomal locus, the RAE180 repeats are present also in a small autosomal locus and additionally have been massively amplified in the Y chromosomes of XX/XY(1)Y(2) species. Here, we have found that the RAE180 repeats of the autosomal locus of XX/XY species are characterized by intra-specific sequence homogeneity and inter-specific divergence and that the comparison of individual nucleotide positions between related species shows a general pattern of concerted evolution. On the contrary, both in the autosomal and the Y-linked loci of XX/XY(1)Y(2) species, ancestral variability has remained with reduced rates of sequence homogenization and of evolution. Thus, this study demonstrates that molecular mechanisms of non-reciprocal exchange are key factors in the molecular drive process; the satellite DNAs in the non-recombining Y chromosomes show low rates of concerted evolution and intra-specific variability increase with no inter-specific divergence. By contrast, freely recombining loci undergo concerted evolution with genetic differentiation between species as occurred in the autosomal locus of XX/XY species. However, evolutionary periods of rapid sequence change might alternate with evolutionary periods of stasis with variability remaining by the reduced action of molecular mechanisms of non-reciprocal exchange as occurred in XX/XY(1)Y(2) species, which could depend on repeat-copy number and the processes involved in their amplification.

Characterization of RUSI, a telomere-associated satellite DNA, in the genus Rumex (Polygonaceae)

A satellite-DNA family (RUSI) has been isolated and characterized in Rumexinduratus Boiss and Reuter (Polygonaceae), an Iberian endemic polygamous sorrel. The RUSI repeats are 170 bp in length and approximately 68% AT-rich containing different variants of degenerate telomere motifs--(TT)(n)AN(GG)(n) -, a typical feature of subtelomeric DNA repeats adjacent to telomeres, which have been referred to as telomere-associated sequences or TASs. In fact, fluorescent in situhybridization showed that this satellite DNA is located in subtelomeric positions of most of the chromosomes of R. induratus, with some centromeric loci. PCR and Southern-blot hybridization assays for sequence conservation in the genus Rumex, indicated that the RUSI sequences are restricted to the genomes of R. induratus and R. scutatus, both species of the section Scutati, suggesting that they are recently evolved. Sequence variation within the two species is high (mean value of sequence differences between repeats of 15% for R. induratus and 7.5% for R. scutatus) and the degree of sequence differentiation between species is low with no species-specific variants, postulated to be due to slowed rates of spreading of sequence variants by molecular homogenizing mechanisms. Characteristics of RUSI sequences are discussed in the light of their chromosomal location and analyzed for their evolutionary and phylogenetic implications.

Patterns of tandem repetition in plant whole genome assemblies

Tandem repeats often confound large genome assemblies. A survey of tandemly arrayed repetitive sequences was carried out in whole genome sequences of the green alga Chlamydomonas reinhardtii, the moss Physcomitrella patens, the monocots rice and sorghum, and the dicots Arabidopsis thaliana, poplar, grapevine, and papaya, in order to test how these assemblies deal with this fraction of DNA. Our results suggest that plant genome assemblies preferentially include tandem repeats composed of shorter monomeric units (especially dinucleotide and 9-30-bp repeats), while higher repetitive units pose more difficulties to assemble. Nevertheless, notwithstanding that currently available sequencing technologies struggle with higher arrays of repeated DNA, major well-known repetitive elements including centromeric and telomeric repeats as well as high copy-number genes, were found to be reasonably well represented. A database including all tandem repeat sequences characterized here was created to benefit future comparative genomic analyses.

The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus)

Papaya, a fruit crop cultivated in tropical and subtropical regions, is known for its nutritional benefits and medicinal applications. Here we report a 3x draft genome sequence of 'SunUp' papaya, the first commercial virus-resistant transgenic fruit tree to be sequenced. The papaya genome is three times the size of the Arabidopsis genome, but contains fewer genes, including significantly fewer disease-resistance gene analogues. Comparison of the five sequenced genomes suggests a minimal angiosperm gene set of 13,311. A lack of recent genome duplication, atypical of other angiosperm genomes sequenced so far, may account for the smaller papaya gene number in most functional groups. Nonetheless, striking amplifications in gene number within particular functional groups suggest roles in the evolution of tree-like habit, deposition and remobilization of starch reserves, attraction of seed dispersal agents, and adaptation to tropical daylengths. Transgenesis at three locations is closely associated with chloroplast insertions into the nuclear genome, and with topoisomerase I recognition sites. Papaya offers numerous advantages as a system for fruit-tree functional genomics, and this draft genome sequence provides the foundation for revealing the basis of Carica's distinguishing morpho-physiological, medicinal and nutritional properties.

The evolution of sex chromosomes in the genus Rumex (Polygonaceae): Identification of a new species with heteromorphic sex chromosomes

The structural features and evolutionary state of the sex chromosomes of the XX/XY species of Rumex are unknown. Here, we report a study of the meiotic behaviour of the XY bivalent in Rumex acetosella and R. suffruticosus, a new species which we describe cytogenetically for the first time in this paper, and also that of the XY(1)Y(2) trivalent of R. acetosa by both conventional cytogenetic techniques and analysis of synaptonemal complex formation. Fluorescent in situ hybridization with satellite DNA and rDNA sequences as probes was used to analyse the degree of cytogenetic differentiation between the X and Y chromosomes in order to depict their evolutionary stage in the three species. Contrasting with the advanced state of genetic differentiation between the X and the Y chromosomes in R. acetosa, we have found that R. acetosella and R. suffruticosus represent an early stage of genetic differentiation between sex chromosomes. Our findings further demonstrate the usefulness of the genus Rumex as a model for analysing the evolution of sex chromosomes in plants, since within this genus it is now possible to study the different levels of genetic differentiation between the sex chromosomes and to analyse their evolutionary history from their origin.

SatDNA Analyzer: a computing tool for satellite-DNA evolutionary analysis

satDNA Analyzer is a program, implemented in C++, for the analysis of the patterns of variation at each nucleotide position considered independently amongst all units of a given satellite-DNA family when comparing it between a pair of species. The program classifies each site accordingly as monomorphic or polymorphic, discriminates shared from non-shared polymorphisms and classifies each non-shared polymorphism according to the model proposed by Strachan et al. in six different stages of transition during the spread of a variant repeat unit toward its fixation. Furthermore, this program implements several other utilities for satellite-DNA analysis evolution such as the design of the average consensus sequences, the average base pair contents, the distribution of variant sites, the transition to transversion ratio and different estimates of intra-specific variation and inter-specific variation. Aprioristic hypotheses on factors influencing the molecular drive process and the rates and biases of concerted evolution can be tested with this program. Additionally, satDNA Analyzer generates an output file containing a sequence alignment without shared polymorphisms to be used for further evolutionary analysis by using different phylogenetic softwares.

AVAILABILITY

satDNA Analyzer is freely available at http://satdna.sourceforge.net/. SatDNA Analyzer has been designed to operate on Windows, Linux and Mac OS X.

Cloning and characterization of dispersed repetitive DNA derived from microdissected sex chromosomes of Rumex acetosa

Rumex acetosa is characterized by a multiple chromosome system (2n = 12 + XX for females, and 2n = 12 + XY1Y2 for males), in which sex is determined by the ratio between the number of X chromosomes and autosome sets. For a better understanding of the molecular structure and evolution of plant sex chromosomes, we have generated a sex chromosome specific library of R. acetosa by microdissection. The screening of this library has allowed us to identify 5 repetitive DNA families that have been characterized in detail. One of these families, DOP-20, has shown no homology with other sequences in databases. Nevertheless, the putative proteins encoded by the other 4 families, DOP-8, DOP-47, DOP-60, and DOP-61, show homology with proteins from different plant retroelements, including poly proteins from Ty3-gypsy- and Ty1-copia-like long terminal repeat (LTR) retroelements, and reverse transcriptase from non-LTR retro elements. Results indicate that sequences from these 5 families are dispersed throughout the genome of both males and females, but no appreciable accumulation or differentiation of these types of sequences have been found in the Y chromosomes. These repetitive DNA sequences are more conserved in the genome of other dioecious species such as Rumex papillaris, Rumex intermedius, Rumex thyrsoides, Rumex hastatulus, and Rumex suffruticosus, than in the polygamous, gynodioecious, or hermaphrodite species Rumex induratus, Rumex lunaria, Rumex con glom er atus, Rumex crispus, and Rumex bucephalo phorus, which supports a single origin of dioecious species in this genus. The implication of these transposable elements in the origin and evolution of the heteromorphic sex chromosomes of R. acetosa is discussed.

The origin and evolution of the variability in a Y-specific satellite-DNA of Rumex acetosa and its relatives

In this paper, we analyze a satellite-DNA family, the RAYSI family, which is specific of the Y chromosomes of Rumex acetosa, a dioecious plant species with a multiple sex-chromosome system in which the females are XX and the males are XY(1)Y(2). Here, we demonstrate that this satellite DNA is common to other relatives of R. acetosa, including Rumex papillaris, Rumex intermedius, Rumex thyrsoides and Rumex tuberosus that are also dioecious species with a multiple system of sex chromosomes. This satellite-DNA family is absent from the genomes of other dioecious Rumex species having an XX/XY sex-chromosome system. Our data confirm recent molecular phylogenies that support a unique origin for all dioecious species of Rumex and two separate lineages for species with single or complex sex-chromosome systems. Our data also support an accelerated degeneration of Y-chromosome in XX/XY(1)Y(2) species by the accumulation of satellite-DNA sequences. On the other hand, the particular non-recombining nature of the Y chromosomes of R. acetosa and their closest relatives lead to a particular mode of evolution of RAYSI sequences. Thus, mechanisms leading to the suppression of recombination between the Y chromosomes reduced the rate of concerted evolution and gave rise to the apparition of different RAYSI subfamilies. Thus, R. acetosa and R. intermedius have two subfamilies (the RAYSI-S and RAYSI-J subfamilies and the INT-A and INT-B subfamilies, respectively), while R. papillaris only has one, the RAYSI-J subfamily. The RAYSI-S and RAYSI-J subfamilies of R. acetosa differ in 83 fixed diagnostic sites and several diagnostic deletions while the INT-A and the INT-B of R. intermedius differ in 27 fixed diagnostic sites. Pairwise comparisons between RAYSI-S and RAYSI-J sequences or between INT-A and INT-B sequences revealed these sites to be shared mutations detectable in repeats of the same variant in same positions. Evolutionary comparisons suggest that the subfamily RAYSI-J has appeared in the common ancestor of R. acetosa and R. papillaris, in which RAYSI-J has replaced totally (R. papillaris) or almost totally the ancestral sequence (R. acetosa). This scenario assumes that RAYSI-S sequences should be considered ancestral sequences and that a secondary event of subfamily subdivision should be occurring in R. intermedius, with their RAYSI subfamilies more closely related to one another than with other RAYSI sequences. Our analysis suggests that the different subfamilies diverged by a gradual and cohesive way probably mediated by sister-chromatid interchanges while their expansion or contraction in number might be explained by alternating cycles of sudden mechanisms of amplification or elimination.

The controversial telomeres of lily plants

The molecular structure of the exceptional telomeres of six plant species belonging to the order Asparagales and two species of the order Liliales was analyzed using Southern blot and fluorescence in situ hybridization. Three different situations were found, namely: i) In the two Liliales species, Tulipa australis (Liliaceae) and Merendera montana (Colchicaceae), the chromosome ends display hybridization signals with oligonucleotides resembling telomere repeats of both plants (TTTAGGG)n and vertebrates (TTAGGG)n. ii) Asparagales species such as Phormium tenax (Hemerocallidaceae), Muscari comosum (Hyacinthaceae), Narcissus jonquilla (Amaryllidaceae) and Allium sativum (Alliaceae) lack both the plant telomere repeats and the vertebrate telomere repeats. iii) Two other Asparagales species, Aloe vera (Asphodelaceae) and an Iris hybrid (Iridaceae), display positive hybridization with the vertebrate telomere repeats but not with the plant telomere repeats. Southern blot hybridization revealed concurring results. On this basis, the composition of the telomere structure in this plant group is discussed.

Reduced rates of sequence evolution of Y-linked satellite DNA in Rumex (Polygonaceae)

One characteristic of sex chromosomes is the accumulation of a set of different types of repetitive DNA sequences in the Y chromosomes. However, little is known about how this occurs or about how the absence of recombination affects the subsequent evolutionary fate of the repetitive sequences in the Y chromosome. Here we compare the evolutionary pathways leading to the appearance of three different families of satellite-DNA sequences within the genomes of Rumex acetosa and R. papillaris, two dioecious plant species with a complex XX/XY(1)Y(2) sex-chromosome system. We have found that two of these families, one autosomic (the RAE730 family) and one Y-linked (the RAYSI family), arose independently from the ancestral duplication of the same 120-bp repeat unit. Conversely, a comparative analysis of the three satellite-DNA families reveals no evolutionary relationships between these two and the third, RAE180, also located in the Y chromosomes. However, we have demonstrated that, regardless of the mechanisms that gave rise to these families, satellite-DNA sequences have different evolutionary fates according to their location in different types of chromosomes. Specifically, those in the Y chromosomes have evolved at half the rate of those in the autosomes, our results supporting the hypothesis that satellite DNAs in nonrecombining Y chromosomes undergo lower rates of sequence evolution and homogenization than do satellite DNAs in autosomes.

The evolution of reproductive systems and sex-determining mechanisms within rumex (polygonaceae) inferred from nuclear and chloroplastidial sequence data

The genus Rumex includes hermaphroditic, polygamous, gynodioecious, monoecious, and dioecious species, with the dioecious species being represented by different sex-determining mechanisms and sex-chromosome systems. Therefore, this genus represents an exceptional case study to test several hypotheses concerning the evolution of both mating systems and the genetic control of sex determination in plants. Here, we compare nuclear intergenic transcribed spacers and chloroplast intergenic sequences of 31 species of Rumex. Our phylogenetic analysis supports a systematic classification of the genus, which differs from that currently accepted. In contrast to the current view, this new phylogeny suggests a common origin for all Eurasian and American dioecious species of Rumex, with gynodioecy as an intermediate state on the way to dioecy. Our results support the contention that sex determination based on the balance between the number of X chromosomes and the number of autosomes (X/A balance) has evolved secondarily from male-determining Y mechanisms and that multiple sex-chromosome systems, XX/XY1Y2, were derived twice from an XX/XY system. The resulting phylogeny is consistent with a classification of Rumex species according to their basic chromosome number, implying that the evolution of Rumex species might have followed a process of chromosomal reduction from x = 10 toward x = 7 through intermediate stages (x = 9 and x = 8).