Sexy ways: the methodical approaches to study plant sex chromosomes.

Sex chromosomes have evolved in many plant species with separate sexes. Current plant research is shifting from examining the structure of sex chromosomes to exploring their functional aspects. New studies are progressively unveiling the specific genetic and epigenetic mechanisms responsible for shaping distinct sexes in plants. While the fundamental methods of molecular biology and genomics are generally employed for the analysis of sex chromosomes, it is often necessary to modify classical procedures not only to simplify and expedite analyses but sometimes to make them possible at all. In this review, we demonstrate how, at the level of structural and functional genetics, cytogenetics, and bioinformatics, it is essential to adapt established procedures for sex chromosome analysis.

HiC-TE: a computational pipeline for Hi-C data analysis to study the role of repeat family interactions in the genome 3D organization.

MOTIVATION: The role of repetitive DNA in the 3D organization of the interphase nucleus is a subject of intensive study. In studies of 3D nucleus organization, mutual contacts of various loci can be identified by Hi-C sequencing. Typical analyses use binning of read pairs by location to reduce noise. We use binning by repeat families instead to make similar conclusions about repeat regions. RESULTS: To achieve this, we combined Hi-C data, reference genome data and tools for repeat analysis into a Nextflow pipeline identifying and quantifying the contacts of specific repeat families. As an output, our pipeline produces heatmaps showing contact frequency and circular diagrams visualizing repeat contact localization. Using our pipeline with tomato data, we revealed the preferential homotypic interactions of ribosomal DNA, centromeric satellites and some LTR retrotransposon families and, as expected, little contact between organellar and nuclear DNA elements. While the pipeline can be applied to any eukaryotic genome, results in plants provide better coverage, since the built-in TE-greedy-nester software only detects tandems and LTR retrotransposons. Other repeats can be fed via GFF3 files. This pipeline represents a novel and reproducible way to analyze the role of repetitive elements in the 3D organization of genomes. AVAILABILITY AND IMPLEMENTATION: https://gitlab.fi.muni.cz/lexa/hic-te/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Adaptive changes of the autosomal part of the genome in a dioecious clade of Silene.

The genus Silene brings many opportunities for the study of various processes involved in the evolution of dioecy and young sex chromosomes. Here we focus on a dioecious clade in Silene subgenus Silene and closely related species. This study provides improved support for monophyly of this clade (based on inclusion of further dioecious species) and a new estimate of its age (ca 2.3 million years). We observed a rise in adaptive evolution in the autosomal and pseudoautosomal parts of the genome on the branch where dioecy originated. This increase is not a result of the accumulation of sexually antagonistic genes in the pseudoautosomal region. It is also not caused by the coevolution of genes acting in mitochondria (despite the possibility that dioecy along this branch could have evolved from a nucleo-cytoplasmic male sterility-based system). After considering other possibilities, the most parsimonious explanation for the increase seen in the number of positively selected codons is the adaptive evolution of genes involved in the adaptation of the autosomal part of the genome to dioecy, as described in Charnov's sex-allocation theory. As the observed coincidence cannot prove causality, studies in other dioecious clades are necessary to allow the formation of general conclusions. This article is part of the theme issue 'Sex determination and sex chromosome evolution in land plants'.

Fundamentally different repetitive element composition of sex chromosomes in Rumex acetosa.

BACKGROUND AND AIMS: Dioecious species with well-established sex chromosomes are rare in the plant kingdom. Most sex chromosomes increase in size but no comprehensive analysis of the kind of sequences that drive this expansion has been presented. Here we analyse sex chromosome structure in common sorrel (Rumex acetosa), a dioecious plant with XY1Y2 sex determination, and we provide the first chromosome-specific repeatome analysis for a plant species possessing sex chromosomes. METHODS: We flow-sorted and separately sequenced sex chromosomes and autosomes in R. acetosa using the two-dimensional fluorescence in situ hybridization in suspension (FISHIS) method and Illumina sequencing. We identified and quantified individual repeats using RepeatExplorer, Tandem Repeat Finder and the Tandem Repeats Analysis Program. We employed fluorescence in situ hybridization (FISH) to analyse the chromosomal localization of satellites and transposons. KEY RESULTS: We identified a number of novel satellites, which have, in a fashion similar to previously known satellites, significantly expanded on the Y chromosome but not as much on the X or on autosomes. Additionally, the size increase of Y chromosomes is caused by non-long terminal repeat (LTR) and LTR retrotransposons, while only the latter contribute to the enlargement of the X chromosome. However, the X chromosome is populated by different LTR retrotransposon lineages than those on Y chromosomes. CONCLUSIONS: The X and Y chromosomes have significantly diverged in terms of repeat composition. The lack of recombination probably contributed to the expansion of diverse satellites and microsatellites and faster fixation of newly inserted transposable elements (TEs) on the Y chromosomes. In addition, the X and Y chromosomes, despite similar total counts of TEs, differ significantly in the representation of individual TE lineages, which indicates that transposons proliferate preferentially in either the paternal or the maternal lineage.

Sex and the flower - developmental aspects of sex chromosome evolution.

Background: The evolution of dioecious plants is occasionally accompanied by the establishment of sex chromosomes: both XY and ZW systems have been found in plants. Structural studies of sex chromosomes are now being followed up by functional studies that are gradually shedding light on the specific genetic and epigenetic processes that shape the development of separate sexes in plants. Scope: This review describes sex determination diversity in plants and the genetic background of dioecy, summarizes recent progress in the investigation of both classical and emerging model dioecious plants and discusses novel findings. The advantages of interspecies hybrids in studies focused on sex determination and the role of epigenetic processes in sexual development are also overviewed. Conclusions: We integrate the genic, genomic and epigenetic levels of sex determination and stress the impact of sex chromosome evolution on structural and functional aspects of plant sexual development. We also discuss the impact of dioecy and sex chromosomes on genome structure and expression.

The slowdown of Y chromosome expansion in dioecious Silene latifolia due to DNA loss and male-specific silencing of retrotransposons.

BACKGROUND: The rise and fall of the Y chromosome was demonstrated in animals but plants often possess the large evolutionarily young Y chromosome that is thought has expanded recently. Break-even points dividing expansion and shrinkage phase of plant Y chromosome evolution are still to be determined. To assess the size dynamics of the Y chromosome, we studied intraspecific genome size variation and genome composition of male and female individuals in a dioecious plant Silene latifolia, a well-established model for sex-chromosomes evolution. RESULTS: Our genome size data are the first to demonstrate that regardless of intraspecific genome size variation, Y chromosome has retained its size in S. latifolia. Bioinformatics study of genome composition showed that constancy of Y chromosome size was caused by Y chromosome DNA loss and the female-specific proliferation of recently active dominant retrotransposons. We show that several families of retrotransposons have contributed to genome size variation but not to Y chromosome size change. CONCLUSIONS: Our results suggest that the large Y chromosome of S. latifolia has slowed down or stopped its expansion. Female-specific proliferation of retrotransposons, enlarging the genome with exception of the Y chromosome, was probably caused by silencing of highly active retrotransposons in males and represents an adaptive mechanism to suppress degenerative processes in the haploid stage. Sex specific silencing of transposons might be widespread in plants but hidden in traditional hermaphroditic model plants.

Impact of Repetitive Elements on the Y Chromosome Formation in Plants.

In contrast to animals, separate sexes and sex chromosomes in plants are very rare. Although the evolution of sex chromosomes has been the subject of numerous studies, the impact of repetitive sequences on sex chromosome architecture is not fully understood. New genomic approaches shed light on the role of satellites and transposable elements in the process of Y chromosome evolution. We discuss the impact of repetitive sequences on the structure and dynamics of sex chromosomes with specific focus on Rumex acetosa and Silene latifolia. Recent papers showed that both the expansion and shrinkage of the Y chromosome is influenced by sex-specific regulation of repetitive DNA spread. We present a view that the dynamics of Y chromosome formation is an interplay of genetic and epigenetic processes.

Satellite DNA and Transposable Elements in Seabuckthorn (Hippophae rhamnoides), a Dioecious Plant with Small Y and Large X Chromosomes.

Seabuckthorn (Hippophae rhamnoides) is a dioecious shrub commonly used in the pharmaceutical, cosmetic, and environmental industry as a source of oil, minerals and vitamins. In this study, we analyzed the transposable elements and satellites in its genome. We carried out Illumina DNA sequencing and reconstructed the main repetitive DNA sequences. For data analysis, we developed a new bioinformatics approach for advanced satellite DNA analysis and showed that about 25% of the genome consists of satellite DNA and about 24% is formed of transposable elements, dominated by Ty3/Gypsy and Ty1/Copia LTR retrotransposons. FISH mapping revealed X chromosome-accumulated, Y chromosome-specific or both sex chromosomes-accumulated satellites but most satellites were found on autosomes. Transposable elements were located mostly in the subtelomeres of all chromosomes. The 5S rDNA and 45S rDNA were localized on one autosomal locus each. Although we demonstrated the small size of the Y chromosome of the seabuckthorn and accumulated satellite DNA there, we were unable to estimate the age and extent of the Y chromosome degeneration. Analysis of dioecious relatives such as Shepherdia would shed more light on the evolution of these sex chromosomes.

Overexpression of a flower-specific aerolysin-like protein from the dioecious plant Rumex acetosa alters flower development and induces male sterility in transgenic tobacco.

Sex determination in Rumex acetosa, a dioecious plant with a complex XY1 Y2 sex chromosome system (females are XX and males are XY1 Y2 ), is not controlled by an active Y chromosome but depends on the ratio between the number of X chromosomes and autosomes. To gain insight into the molecular mechanisms of sex determination, we generated a subtracted cDNA library enriched in genes specifically or predominantly expressed in female floral buds in early stages of development, when sex determination mechanisms come into play. In the present paper, we report the molecular and functional characterization of FEM32, a gene encoding a protein that shares a common architecture with proteins in different plants, animals, bacteria and fungi of the aerolysin superfamily; many of these function as ß pore-forming toxins. The expression analysis, assessed by northern blot, RT-PCR and in situ hybridization, demonstrates that this gene is specifically expressed in flowers in both early and late stages of development, although its transcripts accumulate much more in female flowers than in male flowers. The ectopic expression of FEM32 under both the constitutive promoter 35S and the flower-specific promoter AP3 in transgenic tobacco showed no obvious alteration in vegetative development but was able to alter floral organ growth and pollen fertility. The 35S::FEM32 and AP3::FEM32 transgenic lines showed a reduction in stamen development and pollen viability, as well as a diminution in fruit set, fruit development and seed production. Compared with other floral organs, pistil development was, however, enhanced in plants overexpressing FEM32. According to these effects, it is likely that FEM32 functions in Rumex by arresting stamen and pollen development during female flower development. The aerolysin-like pore-forming proteins of eukaryotes are mainly involved in defence mechanisms against bacteria, fungi and insects and are also involved in apoptosis and programmed cell death (PCD), a mechanism that could explain the role of FEM32 in Rumex sex determination.

Impact of repetitive DNA on sex chromosome evolution in plants.

Structurally and functionally diverged sex chromosomes have evolved in many animals as well as in some plants. Sex chromosomes represent a specific genomic region(s) with locally suppressed recombination. As a consequence, repetitive sequences involving transposable elements, tandem repeats (satellites and microsatellites), and organellar DNA accumulate on the Y (W) chromosomes. In this paper, we review the main types of repetitive elements, their gathering on the Y chromosome, and discuss new findings showing that not only accumulation of various repeats in non-recombining regions but also opposite processes form Y chromosome. The aim of this review is also to discuss the mechanisms of repetitive DNA spread involving (retro) transposition, DNA polymerase slippage or unequal crossing-over, as well as modes of repeat removal by ectopic recombination. The intensity of these processes differs in non-recombining region(s) of sex chromosomes when compared to the recombining parts of genome. We also speculate about the relationship between heterochromatinization and the formation of heteromorphic sex chromosomes.

Fully automated pipeline for detection of sex linked genes using RNA-Seq data.

BACKGROUND: Sex chromosomes present a genomic region which to some extent, differs between the genders of a single species. Reliable high-throughput methods for detection of sex chromosomes specific markers are needed, especially in species where genome information is limited. Next generation sequencing (NGS) opens the door for identification of unique sequences or searching for nucleotide polymorphisms between datasets. A combination of classical genetic segregation analysis along with RNA-Seq data can present an ideal tool to map and identify sex chromosome-specific expressed markers. To address this challenge, we established genetic cross of dioecious plant Rumex acetosa and generated RNA-Seq data from both parental generation and male and female offspring. RESULTS: We present a pipeline for detection of sex linked genes based on nucleotide polymorphism analysis. In our approach, tracking of nucleotide polymorphisms is carried out using a cross of preferably distant populations. For this reason, only 4 datasets are needed - reads from high-throughput sequencing platforms for parent generation (mother and father) and F1 generation (male and female progeny). Our pipeline uses custom scripts together with external assembly, mapping and variant calling software. Given the resource-intensive nature of the computation, servers with high capacity are a requirement. Therefore, in order to keep this pipeline easily accessible and reproducible, we implemented it in Galaxy - an open, web-based platform for data-intensive biomedical research. Our tools are present in the Galaxy Tool Shed, from which they can be installed to any local Galaxy instance. As an output of the pipeline, user gets a FASTA file with candidate transcriptionally active sex-linked genes, sorted by their relevance. At the same time, a BAM file with identified genes and alignment of reads is also provided. Thus, polymorphisms following segregation pattern can be easily visualized, which significantly enhances primer design and subsequent steps of wet-lab verification. CONCLUSIONS: Our pipeline presents a simple and freely accessible software tool for identification of sex chromosome linked genes in species without an existing reference genome. Based on combination of genetic crosses and RNA-Seq data, we have designed a high-throughput, cost-effective approach for a broad community of scientists focused on sex chromosome structure and evolution.

Identification of a novel retrotransposon with sex chromosome-specific distribution in Silene latifolia.

Silene latifolia is a dioecious plant species with chromosomal sex determination. Although the evolution of sex chromosomes in S. latifolia has been the subject of numerous studies, a global view of X chromosome structure in this species is still missing. Here, we combine X chromosome microdissection and BAC library screening to isolate new X chromosome-linked sequences. Out of 8 identified BAC clones, only BAC 86M14 showed an X-preferential signal after FISH experiments. Further analysis revealed the existence of the Athila retroelement which is enriched in the X chromosome and nearly absent in the Y chromosome. Based on previous data, the Athila retroelement belongs to the CL3 group of most repetitive sequences in the S. latifolia genome. Structural, transcriptomics and phylogenetic analyses revealed that Athila CL3 represents an old clade in the Athila lineage. We propose a mechanism responsible for Athila CL3 distribution in the S. latifolia genome.

Possible mechanisms responsible for absence of a retrotransposon family on a plant Y chromosome.

Some transposable elements (TEs) show extraordinary variance in abundance along sex chromosomes but the mechanisms responsible for this variance are unknown. Here, we studied Ogre long terminal repeat (LTR) retrotransposons in Silene latifolia, a dioecious plant with evolutionarily young heteromorphic sex chromosomes. Ogre elements are ubiquitous in the S. latifolia genome but surprisingly absent on the Y chromosome. Bacterial artificial chromosome (BAC) library analysis and fluorescence in situ hybridization (FISH) were used to determine Ogre structure and chromosomal localization. Next generation sequencing (NGS) data were analysed to assess the transcription level and abundance of small RNAs. Methylation of Ogres was determined by bisulphite sequencing. Phylogenetic analysis was used to determine mobilization time and selection forces acting on Ogre elements. We characterized three Ogre families ubiquitous in the S. latifolia genome. One family is nearly absent on the Y chromosome despite all the families having similar structures and spreading mechanisms. We showed that Ogre retrotransposons evolved before sex chromosomes appeared but were mobilized after formation of the Y chromosome. Our data suggest that the absence of one Ogre family on the Y chromosome may be caused by 24-nucleotide (24-nt) small RNA-mediated silencing leading to female-specific spreading. Our findings highlight epigenetic silencing mechanisms as potentially crucial factors in sex-specific spreading of some TEs, but other possible mechanisms are also discussed.

Expansion of microsatellites on evolutionary young Y chromosome.

Sex chromosomes are an ideal system to study processes connected with suppressed recombination. We found evidence of microsatellite expansion, on the relatively young Y chromosome of the dioecious plant sorrel (Rumex acetosa, XY1Y2 system), but no such expansion on the more ancient Y chromosomes of liverwort (Marchantia polymorpha) and human. The most expanding motifs were AC and AAC, which also showed periodicity of array length, indicating the importance of beginnings and ends of arrays. Our data indicate that abundance of microsatellites in genomes depends on the inherent expansion potential of specific motifs, which could be related to their stability and ability to adopt unusual DNA conformations. We also found that the abundance of microsatellites is higher in the neighborhood of transposable elements (TEs) suggesting that microsatellites are probably targets for TE insertions. This evidence suggests that microsatellite expansion is an early event shaping the Y chromosome where this process is not opposed by recombination, while accumulation of TEs and chromosome shrinkage predominate later.

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 Yh for hermaphrodite. The sex chromosome genotypes are XY (male), XYh (hermaphrodite), and XX (female). The papaya hermaphrodite-specific Yh 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 F2 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 Yh 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 Yh-specific regions.

Genomic diversity in two related plant species with and without sex chromosomes--Silene latifolia and S. vulgaris.

BACKGROUND: Genome size evolution is a complex process influenced by polyploidization, satellite DNA accumulation, and expansion of retroelements. How this process could be affected by different reproductive strategies is still poorly understood. METHODOLOGY/PRINCIPAL FINDINGS: We analyzed differences in the number and distribution of major repetitive DNA elements in two closely related species, Silene latifolia and S. vulgaris. Both species are diploid and possess the same chromosome number (2n = 24), but differ in their genome size and mode of reproduction. The dioecious S. latifolia (1C = 2.70 pg DNA) possesses sex chromosomes and its genome is 2.5× larger than that of the gynodioecious S. vulgaris (1C = 1.13 pg DNA), which does not possess sex chromosomes. We discovered that the genome of S. latifolia is larger mainly due to the expansion of Ogre retrotransposons. Surprisingly, the centromeric STAR-C and TR1 tandem repeats were found to be more abundant in S. vulgaris, the species with the smaller genome. We further examined the distribution of major repetitive sequences in related species in the Caryophyllaceae family. The results of FISH (fluorescence in situ hybridization) on mitotic chromosomes with the Retand element indicate that large rearrangements occurred during the evolution of the Caryophyllaceae family. CONCLUSIONS/SIGNIFICANCE: Our data demonstrate that the evolution of genome size in the genus Silene is accompanied by the expansion of different repetitive elements with specific patterns in the dioecious species possessing the sex chromosomes.

Survey of repetitive sequences in Silene latifolia with respect to their distribution on sex chromosomes.

We carried out a global survey of all major types of transposable elements in Silene latifolia, a model species with sex chromosomes that are in the early stages of their evolution. A shotgun genomic library was screened with genomic DNA to isolate and characterize the most abundant elements. We found that the most common types of elements were the subtelomeric tandem repeat X-43.1 and Gypsy retrotransposons, followed by Copia retrotransposons and LINE non-LTR elements. SINE elements and DNA transposons were less abundant. We also amplified transposable elements with degenerate primers and used them to screen the library. The localization of elements by FISH revealed that most of the Copia elements were accumulated on the Y chromosome. Surprisingly, one type of Gypsy element, which was similar to Ogre elements known from legumes, was almost absent on the Y chromosome but otherwise uniformly distributed on all chromosomes. Other types of elements were ubiquitous on all chromosomes. Moreover, we isolated and characterized two new tandem repeats. One of them, STAR-C, was localized at the centromeres of all chromosomes except the Y chromosome, where it was present on the p-arm. Its variant, STAR-Y, carrying a small deletion, was specifically localized on the q-arm of the Y chromosome. The second tandem repeat, TR1, co-localized with the 45S rDNA cluster in the subtelomeres of five pairs of autosomes. FISH analysis of other Silene species revealed that some elements (e.g., Ogre-like elements) are confined to the section Elisanthe while others (e.g. Copia or Athila-like elements) are present also in more distant species. Similarly, the centromeric satellite STAR-C was conserved in the genus Silene whereas the subtelomeric satellite X-43.1 was specific for Elisanthe section. Altogether, our data provide an overview of the repetitive sequences in Silene latifolia and revealed that genomic distribution and evolutionary dynamics differ among various repetitive elements. The unique pattern of repeat distribution is found on the Y chromosome, where some elements are accumulated while other elements are conspicuously absent, which probably reflects different forces shaping the Y chromosome.

Microsatellite accumulation on the Y chromosome in Silene latifolia.

The dioecious plant Silene latifolia possesses evolutionarily young sex chromosomes, and so serves as a model system to study the early stages of sex chromosome evolution. Sex chromosomes often differ distinctly from autosomes in both their structure and their patterns of evolution. The S. latifolia Y chromosome is particularly unique owing to its large size, which contrasts with the size of smaller, degenerate mammalian Y chromosomes. It is thought that the suppression of recombination on the S. latifolia Y chromosome could have resulted in the accumulation of repetitive sequences that account for its large size. Here we used fluorescence in situ hybridization (FISH) to study the chromosomal distribution of various microsatellites in S. latifolia including all possible mono-, di-, and tri-nucleotides. Our results demonstrate that a majority of microsatellites are accumulated on the q arm of the Y chromosome, which stopped recombining relatively recently and has had less time to accumulate repetitive DNA sequences compared with the p arm. Based on these results we can speculate that microsatellites have accumulated in regions that predate the genome expansion, supporting the view that the accumulation of repetitive DNA sequences occurred prior to, not because of, the degeneration of genes.

High intrachromosomal similarity of retrotransposon long terminal repeats: evidence for homogenization by gene conversion on plant sex chromosomes?

Retrotransposons are ubiquitous in the plant genomes and are responsible for their plasticity. Recently, we described a novel family of gypsy-like retrotransposons, named Retand, in the dioecious plant Silene latifolia possessing evolutionary young sex chromosomes of the mammalian type (XY). Here we have analyzed long terminal repeats (LTRs) of Retand that were amplified from laser microdissected X and Y sex chromosomes and autosomes of S. latifolia. A majority of X and Y-derived LTRs formed a few separate clades in phylogenetic analysis reflecting their high intrachromosomal similarity. Moreover, the LTRs localized on the Y chromosome were less divergent than the X chromosome-derived or autosomal LTRs. These data can be explained by a homogenization process, such as gene conversion, working more intensively on the Y chromosome.

Quantitative structure-activity relationships for toxicity and genotoxicity of halogenated aliphatic compounds: wing spot test of Drosophila melanogaster.

Halogenated aliphatic compounds were evaluated for toxic and genotoxic effects in the somatic mutation and recombination test employing Drosophila melanogaster. The tested chemicals included chlorinated, brominated and iodinated; mono-, di- and tri-substituted; saturated and unsaturated alkanes: 1,2-dibromoethane, 1-bromo-2-chloroethane, 1-iodopropane, 2,3-dichloropropene, 3-bromo-1-propene, epibromohydrin, 2-iodobutane, 3-chloro-2-methylpropene, 1,2,3-trichloropropane, 1,2-dichloroethane, 1,2-dichlorobutane, 1-chloro-2-methylpropane, 1,3-dichloropropane, 1,2-dichloropropane, 2-chloroethymethylether, 1-bromo-2-methylpropane and 1-chloropentane. N-methyl-N-nitrosourea served as the positive and distilled water as the negative control. The set of chemicals for the toxicological testing was selected by the use of statistical experiment design. Group of unsaturated aliphatic hydrocarbons were generally more toxic than saturated analogues. The genotoxic effect was observed with 14 compounds in the wing spot test, while 3 substances did not show any genotoxicity by using the wing spot test at 50% lethal concentration. The highest number of wing spots was observed in genotoxicity assay with 1-bromo-2-chloroethane, 1,2-dichloroethane, 1,2-dibromoethane and 1-iodopropane. Nucleophilic superdelocalizability calculated by quantum mechanics appears to be a good parameter for prediction of both toxicity and genotoxicity effects of halogenated aliphatic compounds.