BACKGROUND
sex chromosomes represent a unique “chromosomal elements” appearing in most gonochoristic animals and some plant dioecious species  <cit> .
since most living organisms are still far from being sequenced, mining of data about structure and function of individual sex chromosomes remains a largely unaddressed issue.
various molecular methods have been employed to identify sex specific markers and genes.
in these approaches, dna from individual sex chromosome is either mechanically separated by microdissection or flow-sorting of sex chromosomes  <cit> , followed by construction of libraries and sequencing  <cit> , or alternatively, sex-linked markers have been isolated by various molecular methods such as restriction fragment length polymorphism , amplified fragment length polymorphism  or representational difference analysis   <cit> .
next generation sequencing has become a unique platform for large-scale analysis of genome structure and dynamics.
the use of massive sequencing data has broadly been adopted in genome and transcriptome comparative studies, for characterization of species and/or tissue specific expression profiles, identification of disease associated genes, and comparative genomics analysis  <cit> .
in principle, rna-seq data from each individual of a cross can be processed separately or, alternatively, progeny data can be pooled  to minimize sequencing costs and extent of computation.
here, we present a fully automated pipeline for detection of sex-linked genes that combines classical genetic segregation analysis with rna-seq data for male and female parent generation and pooled male and female progeny.
this pipeline has been used for the identification of x- and y-linked genes in rumex acetosa , a model dioecious plant species with xx/xy1y <dig> sex chromosome system.
the generated data were verified using pcr with newly identified sex-specific primers.
biological material and patterns of segregation used by our algorithm
possible patterns of sex gene segregation are shown in figure  <dig>  in pattern a, the mother is homozygous x1x <dig> , while the father has two different alleles .
thus, all sons possess the same allele variant as the mother - x <dig>  and y allele inherited from the father.
all daughters inherit one allele from each parent, resulting in the x1x <dig> genotype.
in pattern b, the mother is heterozygous x1x <dig> while the father has a different x linked variant .
according to segregation rules, the sons will be x1y or x3y and daughters x1x <dig> or x3x <dig> .
in pattern c, the mother is heterozygous and shares one allele  with the father.figure  <dig>
segregation patterns.
all possible combinations of sex-linked alleles in parent generation are listed.
alleles on the x chromosome exhibit criss-cross inheritance, while alleles on the y chromosome can only be transmitted from father to son .
segregation patterns are used as the starting point for the design of filtering rules.
implementation
workflow of the pipeline used in the present study is visualized in figure  <dig>  the feasibility of our approach was tested using the cross of male and female individuals of r. acetosa.
the total rna was isolated from young leaves  <cit>  and oligo dt primers were used in order to enrich for polya rna.
two very distant lineages were chosen for the study - mother lineage being r. acetosa, almería, spain while the father lineage was r. acetosa, brno-reckovice, czech republic.
<dig> male and  <dig> female individuals from f <dig> generation were chosen for separate isolation, test of integrity and concentration measurements.
finally, male and female rna was pooled in the same proportions  and sequenced along with mother and father samples on illumina hiseq with resulting paired-end reads displaying a length peak of 87 bp for both forward and reverse reads.
with an estimate of 2% of the r. acetosa genome being transcribed, this leads us to an approximation of 62x coverage for the assembly of reference x-linked genes and 44x coverage for the assembly of reference y-linked genes  <cit> .figure  <dig>
simplified workflow.
mother’s reads are assembled with trinity assembler and newly created contigs are then used as a reference to which all reads are mapped.
this yields into  <dig> bam files storing alignments.
after applying filters, x-linked contigs are identified.
for detection of y-linked genes, reference contigs are assembled from male reads
the pipeline therefore requires four datasets with preprocessed rna-seq reads  - mother reads, father reads, reads of pooled female progeny and reads of pooled male progeny.
both, fasta and fastq files are accepted and can be easily interchangeable in galaxy workflow.
for maximal efficiency of our pipeline, we recommend the following guidelines:reads from the mother should be the most abundant since they are essential for deriving x-linked contigs,
attention should be paid such that even numbers of individuals are pooled as it increases probability that the observed allele ratio will be concordantly to segregation rules 50%: 50%,
if possible, a cross of distant populations is recommended  as it increases the probability of allele variability,
users are encouraged to use long pair-end reads for sequencing in order to lower the number of reads that cannot be mapped uniquely,
if possible, strand-specific rna-seq protocol should be used preferably
for this project, we designed filtering rules for each segregation pattern and implemented rules for pattern a .
according to these rules, only contigs representing high quality sex-linked markers are preserved.table  <dig>
list of filtering rules used for identification of sex-linked contigs; filtering rules for identification of x-linked genes following segregation pattern a
mother must be homozygous .
sons must be same as mother .
father must differ from mother.
daughters must have mother and father alleles .
daughters and sons cannot have any intersection.
sufficient coverage.
the pipeline is accessible under the list of published galaxy workflows on the webpage: usegalaxy.org/workflow/list_published.
it incorporates custom tools, such as linkyx_x, linkyx_y and linkyx_mpileup for identification of x-linked genes, identification of y-linked genes and custom mpileup, respectively.
all these tools can be installed from galaxy tool shed to any galaxy instance  <cit> .
in our experience, a typical dataset with a total of  <dig> million reads requires a computation time in the range of 18– <dig> hours, with most of the time dedicated to de-novo assembly.
the project home page at github.com/biomonika/linkyx/ describes how to install the tools and built-in wrappers of the tools should install all of the dependencies automatically.
identification of x-linked genes
in order to obtain transcripts of putative x-linked genes, all mother reads were assembled into contigs with trinity assembler  <cit> .
these contigs therefore represented the reference.
this is a crucial step of the procedure since in previous works  <cit> , the reference was derived from all individuals.
including all sequences from all individuals for building reference contigs necessarily leads to the creation of an artificial individual.
all polymorphisms then need to be compared to this artificial individual and subsequently between each other .
also, some contigs will come from male-specific sequences only and the method will fail when diverged x and y homologs are not assembled together into one contig.
therefore, building reference contigs exclusively from the mother’s dataset significantly aids the gene detection process.
in the next step, reads from all family members are mapped back onto the contigs with bwa   <cit> .
four sam/bam files containing alignment of reads are created with samtools and sorted  <cit> .
duplicate reads, which may result from artificial pcr amplification step during sequencing procedure, are removed.
three software solutions have been tested for this procedure - picard’s markduplicates function , samtools’ rmdup and fastx collapse .
since fastx collapse cannot produce fastq files and markduplicates repeatedly reported biggest number of properly paired reads and lowest number of singletons for our data, it was chosen for further use in the pipeline.
subsequently, variants were called.
samtools mpileup was set to -uf and bcftools view was used with following parameters: “-p  <dig>  -cgv -”.
resulting vcf files contain snps, indels and short structural rearrangements.
in principle, they contain all information with regard to how a sample differs from the reference sequence.
samtools parameters have been chosen by rule of thumb after extensive testing.
it must be taken into account that reporting many variants that will contain sequencing errors or misalignments can lead to excluding a true x-linked gene for technical reasons, reporting fewer variants may lead to reporting autosomal genes that break one of our filtering criteria with these fraudulences not being recognized.
the selected criteria should thus represent a balance between both the above-mentioned problems – the first of which raises the level of false negatives and the second which raises the level of false positives.
beside variants generated by samtools, the program itself checks for variants that could not be reported by samtools and still should be forbidden at certain positions.
for example, if there are  <dig> reads covering a certain position,  <dig> of them reporting adenin and  <dig> reporting cytosin, samtools may not report the  <dig> cytosins.
however, if according to filtering rules cytosins are forbidden at that position, the program excludes the whole contig from further computation .
in this way, the initial dataset  is filtered to contain only candidate transcriptionally active sex-linked markers.
custom bash and perl scripts were written for the filtering step with rules described in tables  <dig>   <dig> and  <dig>  genes exhibiting pattern a are filtered based on  <dig> rules described in table  <dig> and the principle of filtering is visualized in figure  <dig> table  <dig>
list of filtering rules used for identification of sex-linked contigs; filtering rules for identification of y-linked genes
minimal and maximal number of reads.
fraction of female and male reads.
male reads must be more abundant than female reads.
1
depth in number of reads:
tolerated depth of alternative allele  in number of reads:
0-4
5-24
> = 25schematic representation of identification of x-linked genes exhibiting the a pattern.
mother is homozygote , father has a different variant .
approximately half of the daughter’s reads follow the father’s variant  while sons inherit mother’s variant.
variants specific only for males suggesting the y origin are neglected.
for every set of genes, the following information is provided: contig name, sequence, isoforms listed .
the list of genes is sorted.
for x-linked genes, their list is sorted by the number of uniquely mapped reads for each family member, following the expectation that contigs with the least number of such reads are the easiest to be confirmed as truly x-linked .
identification of y-linked genes
in order to obtain transcripts of putative y-linked genes, all male reads were assembled into contigs with trinity assembler.
all reads are then mapped back, sorted and deduplicated in a manner that is similar to the previous case.
if the contigs come from the y chromosome or are expressed specifically in males, then the female reads do not map to those contigs, while male reads do .
four bam files  are therefore filtered according to designed rules described in table  <dig>  before applying these rules, read numbers are normalized to the number of mother reads.figure  <dig>
schematic representation of identification of y-linked genes.
reference contigs are assembled from male reads only.
then, female reads do not map to contigs of y-chromosome origin while male reads do.
short fragments in red color visualized on the figure represent limited tolerance of mapping/assemble errors.
to test that proposed filters sufficiently separated y-linked contigs from contigs present in both sexes, all reads were mapped again on the male reference contigs with highly relaxed conditions represented by bwa parameter aln -n  <dig> .
with these conditions,  <dig> % of reads were mapped instead of  <dig> %, which confirmed reliability of chosen parameters and showed that identified contigs were either very divergent y copies or standalone genes.
again, putative genes are sorted from the most relevant ones.
the key for sorting is coverage, with the first contig covered by the biggest number of reads.
RESULTS
pcr confirmation of sex chromosome-linkage of identified genes
<dig> candidate x-linked genes and  <dig> candidate y-linked genes have been obtained by the pipeline presented here.
from these,  <dig> x-linked genes and  <dig> y-linked genes were randomly chosen for wet-lab confirmation.
first, pcr was performed for the  <dig> candidate x-linked genes, pcr products were cloned and sequenced.
manual inspection of these sequences  confirmed their x-linkage .
in order to verify the y-linkage of predicted genes, pcr was performed using dna from male and female plants of parental generation.
a majority of genes  showed a prominent male pcr product suggesting their y origin .
in the case of one contig, pcr products could be seen in both male and female samples suggesting that this contig might be present in both sexes, while manifesting a sex-dependent expression.figure  <dig>
experimental laboratory verification of sex-linked contigs.
a) pcr products of candidate x-linked contigs were sequenced and clustered.
father and daughter variants cluster together, which confirms x-linkage of a selected gene.
note that snp marked in blue color is shared only among sequences of father and daughters.
another snp in yellow represents sequencing error.
b) pcr products of candidate y-linked contigs.
for every contig/gene genomic dna of  <dig> individuals  was used as a template: father , brother of father , two sisters of father , mother , two brothers of mother .
product is present only in male individuals.
in order to exclude the possibility that the patterns of 7/ <dig> contigs were caused by difference between the spanish  and the czech  populations, new individuals were included from both populations.
this step confirmed that  <dig> out of  <dig> contigs were truly of y chromosome origin.
in one case, pcr products were missing in the spanish population which suggests that this population has different variants on primer-binding sites.
taken together, 6/ <dig> genes were reliably confirmed and none were certainly refuted.
for x-linked verification, we used the dna from father, mother and pooled male and female progeny.
we sequenced 20– <dig> sequences for each contig after amplification.
after alignment, contigs where father allele was passed only to daughters  were considered x-linked.
we were able to reliably confirm that  <dig> contigs were truly x-linked.
in  <dig> contigs, a conclusion could not be made due to the low coverage and the possibility of not sampling all alleles.
one contig represented false positive result and was probably of autosomal origin.
discussion
recently, chibalina et al.
<cit>  used an rna-seq based approach for identification of sex-linked genes in silene latifolia, another dioecious plant with xy sex chromosomes.
they used two different crosses and subsequent sequencing of all members of the family  for identification of marker chromosome linkage.
the same approach was used by hough et al.
<cit>  to identify sex-linked genes in the dioecious plant r. hastatulus with a neo-y sex chromosome system.
they used transcriptomes from parents and f <dig> progeny from two within-population crosses, one from a population with xy males and one from a population with xy1y <dig> males.
similar strategy was used by muyle et al.
<cit>  who analysed a 10-generations inbred line of silene latifolia for determination of sex linkage of genes using rna-seq.
three males and three females were pooled and massively sequenced.
subsequently, contigs representing consensus sequence for all individuals were assembled together.
in such an inbred line, a very low heterozygosity was expected and thus the contigs were searched for a pattern where all males showed the heterozygous xy pattern, and all females presented the homozygous xx pattern.
it was expected that in males, reads of the x chromosome origin would be the same as those in females, whereas reads from the y chromosome would exhibit a new variant.
all three reports predicted hundreds of new x- or y-linked genes.
however, we assume that some predicted genes could be false positives, since segregation of any snp according to expected pattern was sufficient for considering contigs as x/y-chromosome linked.
in contrast, in our approach all snps in analyzed expressed markers exhibit expected segregation.
therefore, the number of predicted genes is lower but reliability is higher.
for this reason, genes predicted by our tool are more suitable as an initial dataset for isolation of y- and x-linked bac clones and physical mapping of sex chromosomes.
it should be noted that our approach would not identify sex-linked genes that are not recombinationally trapped  or are recent additions to non-recombining regions.
we conclude that this approach together with experimental laboratory verification represents a reliable tool for obtaining sex-linked markers for non-model species.
CONCLUSIONS
in this work, we developed a bioinformatics approach for detection of sex-linked genes using rna-seq data.
such an approach could be used for identification of both y  chromosome-linked or x  chromosome-linked genes in animal and plant species.
no bioinformatics skills are needed for users of our pipeline since fastq files  are the only input for data processing.
the tool is freely accessible to the scientific community via github.com/biomonika/linkyx/.
as a proof of concept, rna-seq data were used in combination with a cross of parents and f <dig> male and female offspring in r. acetosa.
after de novo assembly of illumina based data we identified x- or y- linked patterns of segregation based on nucleotide polymorphism analysis.
for validation of our prediction, we selected y-specific est for pcr experiments that proved the reliability of our approach for high throughput identification of sex-linked markers.
moreover, we showed for the first time that r. acetosa heterochromatic y chromosomes contain transcriptionally active genes.
prior to the present work, only repetitive dna were shown to be gathered on the both y chromosomes in this species  <cit> .
availability and requirements
project name: fully automated pipeline for detection of sex linked genes using rna-seq data.
project home page:github.com/biomonika/linkyx/.
operating system: unix-based, primarily tested on debian.
programming language: bash, perl, python.
other requirements: java  <dig> .
license: academic free license  <dig> .
any restrictions to use by non-academics: no.
additional file
additional file 1: table s <dig>  x-linked and y-linked primers and sequences used for verification.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
ek and bv designed experiment and critically revised manuscript.
mm and rh wrote manuscript, mm implemented pipeline and carried out verification of y-linked genes and zk calibrated the parameters of the pipeline and carried out verification of x-linked genes.
all authors read and approved the final manuscript.

 sex chromosomes represent a unique “chromosomal elements” appearing in most gonochoristic animals and some plant dioecious species. since most living organisms are still far from being sequenced, mining of data about structure and function of individual sex chromosomes remains a largely unaddressed issue. various molecular methods have been employed to identify sex specific markers and genes. in these approaches, dna from individual sex chromosome is either mechanically separated by microdissection or flow-sorting of sex chromosomes, followed by construction of libraries and sequencing, or alternatively, sex-linked markers have been isolated by various molecular methods such as restriction fragment length polymorphism, amplified fragment length polymorphism or representational difference analysis. next generation sequencing has become a unique platform for large-scale analysis of genome structure and dynamics. the use of massive sequencing data has broadly been adopted in genome and transcriptome comparative studies, for characterization of species and/or tissue specific expression profiles, identification of disease associated genes, and comparative genomics analysis. in principle, rna-seq data from each individual of a cross can be processed separately or, alternatively, progeny data can be pooled to minimize sequencing costs and extent of computation. here, we present a fully automated pipeline for detection of sex-linked genes that combines classical genetic segregation analysis with rna-seq data for male and female parent generation and pooled male and female progeny. this pipeline has been used for the identification of x- and y-linked genes in rumex acetosa, a model dioecious plant species with xx/xy1y sex chromosome system. the generated data were verified using pcr with newly identified sex-specific primers. biological material and patterns of segregation used by our algorithm possible patterns of sex gene segregation are shown in figure  in pattern a, the mother is homozygous x1x, while the father has two different alleles. thus, all sons possess the same allele variant as the mother - x and y allele inherited from the father. all daughters inherit one allele from each parent, resulting in the x1x genotype. in pattern b, the mother is heterozygous x1x while the father has a different x linked variant. according to segregation rules, the sons will be x1y or x3y and daughters x1x or x3x. in pattern c, the mother is heterozygous and shares one allele with the father.figure segregation patterns. all possible combinations of sex-linked alleles in parent generation are listed. alleles on the x chromosome exhibit criss-cross inheritance, while alleles on the y chromosome can only be transmitted from father to son. segregation patterns are used as the starting point for the design of filtering rules. implementation workflow of the pipeline used in the present study is visualized in figure  the feasibility of our approach was tested using the cross of male and female individuals of r. acetosa. the total rna was isolated from young leaves and oligo dt primers were used in order to enrich for polya rna. two very distant lineages were chosen for the study - mother lineage being r. acetosa, almería, spain while the father lineage was r. acetosa, brno-reckovice, czech republic. male and female individuals from f generation were chosen for separate isolation, test of integrity and concentration measurements. finally, male and female rna was pooled in the same proportions and sequenced along with mother and father samples on illumina hiseq with resulting paired-end reads displaying a length peak of 87 bp for both forward and reverse reads. with an estimate of 2% of the r. acetosa genome being transcribed, this leads us to an approximation of 62x coverage for the assembly of reference x-linked genes and 44x coverage for the assembly of reference y-linked genes.figure simplified workflow. mother’s reads are assembled with trinity assembler and newly created contigs are then used as a reference to which all reads are mapped. this yields into bam files storing alignments. after applying filters, x-linked contigs are identified. for detection of y-linked genes, reference contigs are assembled from male reads the pipeline therefore requires four datasets with preprocessed rna-seq reads - mother reads, father reads, reads of pooled female progeny and reads of pooled male progeny. both, fasta and fastq files are accepted and can be easily interchangeable in galaxy workflow. for maximal efficiency of our pipeline, we recommend the following guidelines:reads from the mother should be the most abundant since they are essential for deriving x-linked contigs, attention should be paid such that even numbers of individuals are pooled as it increases probability that the observed allele ratio will be concordantly to segregation rules 50%: 50%, if possible, a cross of distant populations is recommended as it increases the probability of allele variability, users are encouraged to use long pair-end reads for sequencing in order to lower the number of reads that cannot be mapped uniquely, if possible, strand-specific rna-seq protocol should be used preferably for this project, we designed filtering rules for each segregation pattern and implemented rules for pattern a. according to these rules, only contigs representing high quality sex-linked markers are preserved.table list of filtering rules used for identification of sex-linked contigs; filtering rules for identification of x-linked genes following segregation pattern a mother must be homozygous. sons must be same as mother. father must differ from mother. daughters must have mother and father alleles. daughters and sons cannot have any intersection. sufficient coverage. the pipeline is accessible under the list of published galaxy workflows on the webpage: usegalaxy.org/workflow/list_published. it incorporates custom tools, such as linkyx_x, linkyx_y and linkyx_mpileup for identification of x-linked genes, identification of y-linked genes and custom mpileup, respectively. all these tools can be installed from galaxy tool shed to any galaxy instance. in our experience, a typical dataset with a total of million reads requires a computation time in the range of 18– hours, with most of the time dedicated to de-novo assembly. the project home page at github.com/biomonika/linkyx/ describes how to install the tools and built-in wrappers of the tools should install all of the dependencies automatically. identification of x-linked genes in order to obtain transcripts of putative x-linked genes, all mother reads were assembled into contigs with trinity assembler. these contigs therefore represented the reference. this is a crucial step of the procedure since in previous works, the reference was derived from all individuals. including all sequences from all individuals for building reference contigs necessarily leads to the creation of an artificial individual. all polymorphisms then need to be compared to this artificial individual and subsequently between each other. also, some contigs will come from male-specific sequences only and the method will fail when diverged x and y homologs are not assembled together into one contig. therefore, building reference contigs exclusively from the mother’s dataset significantly aids the gene detection process. in the next step, reads from all family members are mapped back onto the contigs with bwa. four sam/bam files containing alignment of reads are created with samtools and sorted. duplicate reads, which may result from artificial pcr amplification step during sequencing procedure, are removed. three software solutions have been tested for this procedure - picard’s markduplicates function, samtools’ rmdup and fastx collapse. since fastx collapse cannot produce fastq files and markduplicates repeatedly reported biggest number of properly paired reads and lowest number of singletons for our data, it was chosen for further use in the pipeline. subsequently, variants were called. samtools mpileup was set to -uf and bcftools view was used with following parameters: “-p -cgv -”. resulting vcf files contain snps, indels and short structural rearrangements. in principle, they contain all information with regard to how a sample differs from the reference sequence. samtools parameters have been chosen by rule of thumb after extensive testing. it must be taken into account that reporting many variants that will contain sequencing errors or misalignments can lead to excluding a true x-linked gene for technical reasons, reporting fewer variants may lead to reporting autosomal genes that break one of our filtering criteria with these fraudulences not being recognized. the selected criteria should thus represent a balance between both the above-mentioned problems – the first of which raises the level of false negatives and the second which raises the level of false positives. beside variants generated by samtools, the program itself checks for variants that could not be reported by samtools and still should be forbidden at certain positions. for example, if there are reads covering a certain position, of them reporting adenin and reporting cytosin, samtools may not report the cytosins. however, if according to filtering rules cytosins are forbidden at that position, the program excludes the whole contig from further computation. in this way, the initial dataset is filtered to contain only candidate transcriptionally active sex-linked markers. custom bash and perl scripts were written for the filtering step with rules described in tables  and genes exhibiting pattern a are filtered based on rules described in table  and the principle of filtering is visualized in figure  table list of filtering rules used for identification of sex-linked contigs; filtering rules for identification of y-linked genes minimal and maximal number of reads. fraction of female and male reads. male reads must be more abundant than female reads. 1 depth in number of reads: tolerated depth of alternative allele in number of reads: 0-4 5-24 > = 25schematic representation of identification of x-linked genes exhibiting the a pattern. mother is homozygote, father has a different variant. approximately half of the daughter’s reads follow the father’s variant while sons inherit mother’s variant. variants specific only for males suggesting the y origin are neglected. for every set of genes, the following information is provided: contig name, sequence, isoforms listed. the list of genes is sorted. for x-linked genes, their list is sorted by the number of uniquely mapped reads for each family member, following the expectation that contigs with the least number of such reads are the easiest to be confirmed as truly x-linked. identification of y-linked genes in order to obtain transcripts of putative y-linked genes, all male reads were assembled into contigs with trinity assembler. all reads are then mapped back, sorted and deduplicated in a manner that is similar to the previous case. if the contigs come from the y chromosome or are expressed specifically in males, then the female reads do not map to those contigs, while male reads do. four bam files are therefore filtered according to designed rules described in table  before applying these rules, read numbers are normalized to the number of mother reads.figure schematic representation of identification of y-linked genes. reference contigs are assembled from male reads only. then, female reads do not map to contigs of y-chromosome origin while male reads do. short fragments in red color visualized on the figure represent limited tolerance of mapping/assemble errors. to test that proposed filters sufficiently separated y-linked contigs from contigs present in both sexes, all reads were mapped again on the male reference contigs with highly relaxed conditions represented by bwa parameter aln -n. with these conditions, % of reads were mapped instead of %, which confirmed reliability of chosen parameters and showed that identified contigs were either very divergent y copies or standalone genes. again, putative genes are sorted from the most relevant ones. the key for sorting is coverage, with the first contig covered by the biggest number of reads.  pcr confirmation of sex chromosome-linkage of identified genes candidate x-linked genes and candidate y-linked genes have been obtained by the pipeline presented here. from these, x-linked genes and y-linked genes were randomly chosen for wet-lab confirmation. first, pcr was performed for the candidate x-linked genes, pcr products were cloned and sequenced. manual inspection of these sequences confirmed their x-linkage. in order to verify the y-linkage of predicted genes, pcr was performed using dna from male and female plants of parental generation. a majority of genes showed a prominent male pcr product suggesting their y origin. in the case of one contig, pcr products could be seen in both male and female samples suggesting that this contig might be present in both sexes, while manifesting a sex-dependent expression.figure experimental laboratory verification of sex-linked contigs. a) pcr products of candidate x-linked contigs were sequenced and clustered. father and daughter variants cluster together, which confirms x-linkage of a selected gene. note that snp marked in blue color is shared only among sequences of father and daughters. another snp in yellow represents sequencing error. b) pcr products of candidate y-linked contigs. for every contig/gene genomic dna of individuals was used as a template: father, brother of father, two sisters of father, mother, two brothers of mother. product is present only in male individuals. in order to exclude the possibility that the patterns of 7/ contigs were caused by difference between the spanish and the czech populations, new individuals were included from both populations. this step confirmed that out of contigs were truly of y chromosome origin. in one case, pcr products were missing in the spanish population which suggests that this population has different variants on primer-binding sites. taken together, 6/ genes were reliably confirmed and none were certainly refuted. for x-linked verification, we used the dna from father, mother and pooled male and female progeny. we sequenced 20– sequences for each contig after amplification. after alignment, contigs where father allele was passed only to daughters were considered x-linked. we were able to reliably confirm that contigs were truly x-linked. in contigs, a conclusion could not be made due to the low coverage and the possibility of not sampling all alleles. one contig represented false positive result and was probably of autosomal origin. discussion recently, chibalina et al. used an rna-seq based approach for identification of sex-linked genes in silene latifolia, another dioecious plant with xy sex chromosomes. they used two different crosses and subsequent sequencing of all members of the family for identification of marker chromosome linkage. the same approach was used by hough et al. to identify sex-linked genes in the dioecious plant r. hastatulus with a neo-y sex chromosome system. they used transcriptomes from parents and f progeny from two within-population crosses, one from a population with xy males and one from a population with xy1y males. similar strategy was used by muyle et al. who analysed a 10-generations inbred line of silene latifolia for determination of sex linkage of genes using rna-seq. three males and three females were pooled and massively sequenced. subsequently, contigs representing consensus sequence for all individuals were assembled together. in such an inbred line, a very low heterozygosity was expected and thus the contigs were searched for a pattern where all males showed the heterozygous xy pattern, and all females presented the homozygous xx pattern. it was expected that in males, reads of the x chromosome origin would be the same as those in females, whereas reads from the y chromosome would exhibit a new variant. all three reports predicted hundreds of new x- or y-linked genes. however, we assume that some predicted genes could be false positives, since segregation of any snp according to expected pattern was sufficient for considering contigs as x/y-chromosome linked. in contrast, in our approach all snps in analyzed expressed markers exhibit expected segregation. therefore, the number of predicted genes is lower but reliability is higher. for this reason, genes predicted by our tool are more suitable as an initial dataset for isolation of y- and x-linked bac clones and physical mapping of sex chromosomes. it should be noted that our approach would not identify sex-linked genes that are not recombinationally trapped or are recent additions to non-recombining regions. we conclude that this approach together with experimental laboratory verification represents a reliable tool for obtaining sex-linked markers for non-model species.  in this work, we developed a bioinformatics approach for detection of sex-linked genes using rna-seq data. such an approach could be used for identification of both y chromosome-linked or x chromosome-linked genes in animal and plant species. no bioinformatics skills are needed for users of our pipeline since fastq files are the only input for data processing. the tool is freely accessible to the scientific community via github.com/biomonika/linkyx/. as a proof of concept, rna-seq data were used in combination with a cross of parents and f male and female offspring in r. acetosa. after de novo assembly of illumina based data we identified x- or y- linked patterns of segregation based on nucleotide polymorphism analysis. for validation of our prediction, we selected y-specific est for pcr experiments that proved the reliability of our approach for high throughput identification of sex-linked markers. moreover, we showed for the first time that r. acetosa heterochromatic y chromosomes contain transcriptionally active genes. prior to the present work, only repetitive dna were shown to be gathered on the both y chromosomes in this species. availability and requirements project name: fully automated pipeline for detection of sex linked genes using rna-seq data. project home page:github.com/biomonika/linkyx/. operating system: unix-based, primarily tested on debian. programming language: bash, perl, python. other requirements: java. license: academic free license. any restrictions to use by non-academics: no. additional file additional file 1: table s x-linked and y-linked primers and sequences used for verification. competing interests the authors declare that they have no competing interests. authors’ contributions ek and bv designed experiment and critically revised manuscript. mm and rh wrote manuscript, mm implemented pipeline and carried out verification of y-linked genes and zk calibrated the parameters of the pipeline and carried out verification of x-linked genes. all authors read and approved the final manuscript. 
for this reason, only datasets are needed – reads from high-throughput sequencing platforms for parent generation and f generation. 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. sex chromosomesrna-seqsegregation analysisrumex acetosaissue-copyright-statement© the author 2015 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. our tools are present in the galaxy tool shed, from which they can be installed to any local galaxy instance. we present a pipeline for detection of sex linked genes based on nucleotide polymorphism analysis. thus, polymorphisms following segregation pattern can be easily visualized, which significantly enhances primer design and subsequent steps of wet-lab verification. 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. in our approach, tracking of nucleotide polymorphisms is carried out using a cross of preferably distant populations. 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.
 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 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.  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 datasets are needed – reads from high-throughput sequencing platforms for parent generation and f generation. 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.  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. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords sex chromosomesrna-seqsegregation analysisrumex acetosaissue-copyright-statement© the author 2015 