Search

Additional file 1.

Additional file 2.

Additional file 4.

Additional file 3.

The relative rate of evolution for sex‐biased genes has often been used as a measure of the strength of sex‐specific selection. In contrast to studies in a wide variety of animals, far less is known about the molecular evolution of sex‐biased genes in plants, particularly in dioecious angiosperms. Here, we investigate the gene expression patterns and evolution of sex‐biased genes in the dioecious plant Salix viminalis. We observe lower rates of sequence evolution for male‐biased genes expressed in the reproductive tissue compared to unbiased and female‐biased genes. These results could be partially explained by the lower codon usage bias for male‐biased genes leading to elevated rates of synonymous substitutions compared to unbiased genes. However, the stronger haploid selection in the reproductive tissue of plants, together with pollen competition, would also lead to higher levels of purifying selection acting to remove deleterious variation. Future work should focus on the differential evolution of haploid‐ and diploid‐specific genes to understand the selective dynamics acting on these loci.

Background Biased sex ratios are common among dioecious plant species despite the theoretical prediction of selective advantage of even sex ratios. Albeit the high prevalence of deviations from even sex ratios, the genetic causes to sex biases are rarely known outside of a few model species. Here we present a mechanism underlying the female biased sex ratio in the dioecious willow species Salix viminalis. Results We compared the segregation pattern of genome-wide single nucleotide polymorphism markers in two contrasting bi-parental pedigree populations, the S3 with even sex ratio and the S5 with a female biased sex ratio. With the segregation analysis and comparison between the two populations, we were able to demonstrate that sex determination and sex ratio distortion are controlled by different genetic mechanisms. We furthermore located the sex ratio distorter locus to a Z/W-gametologous region on chromosome 15, which was in close linkage with the sex determination locus. Interestingly, all males in the population with biased sex ratio have in this sex ratio distorter locus the same genotype, meaning that males with the Z1/Z3-genotype were missing from the population, thereby creating the 2:1 female biased sex ratio. Conclusions We attribute the absence of Z1/Z3 males to an allelic incompatibility between maternally and paternally inherited alleles in this sex ratio distorter locus. Due to the tight linkage with the sex determination locus only male individuals are purged from the population at an early age, presumably before or during seed development. We showed that such allelic incompatibility could be stably maintained over evolutionary times through a system of overdominant or pseudooverdominant alleles. Thus, it is possible that the same mechanism generates the female biased sex ratio in natural willow populations. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3634-5) contains supplementary material, which is available to authorized users.

Supplementary genetic model description. (PDF 46 kb)

Genetic map and syntheny analysis. S. viminalis whole genome genetic map and analysis of syntheny with the P. trichocarpa genome sequence. (PDF 200Â kb)

Sex chromosomes can evolve when recombination is halted between a pair of chromosomes, and this can lead to degeneration of the sex-limited chromosome. In the early stages of differentiation sex chromosomes are homomorphic, and even though homomorphic sex chromosomes are very common throughout animals and plants, we know little about the evolutionary forces shaping these types of sex chromosomes. We used DNA- and RNA-Seq data from females and males to explore the sex chromosomes in the female heterogametic willow, Salix viminalis, a species with ancient dioecy but with homomorphic sex chromosomes. We detected no major sex differences in read coverage in the sex determination (SD) region, indicating that the W region has not significantly degenerated. However, single nucleotide polymorphism densities in the SD region are higher in females compared with males, indicating very recent recombination suppression, followed by the accumulation of sex-specific single nucleotide polymorphisms. Interestingly, we identified two female-specific scaffolds that likely represent W-chromosome-specific sequence. We show that genes located in the SD region display a mild excess of male-biased expression in sex-specific tissue, and we use allele-specific gene expression analysis to show that this is the result of masculinization of expression on the Z chromosome rather than degeneration of female-expression on the W chromosome. Together, our results demonstrate that insertion of small DNA fragments and accumulation of sex-biased gene expression can occur before the detectable decay of the sex-limited chromosome.

Background Drought is a major environmental stress that can have severe impacts on plant productivity and survival. Understanding molecular mechanisms of drought responses is crucial in order to breed for drought adapted plant cultivars. The aim of the present study was to investigate phenotypic and transcriptional drought responses in two willow genotypes (520 and 592) originating from an experimental cross between S. viminalis × (S. viminalis × S. schwerinii). Willows are woody perennials in the Salicaceae plant family that are grown as bioenergy crops worldwide. Methods An experiment was conducted where plants were exposed to drought and different eco-physiological parameters were assessed. RNA-seq data was furthermore generated with the Illumina technology from root tips and leaves from plants grown in drought and well-watered (WW) conditions. The RNA-seq data was assembled de novo with the Trinity assembler to create a reference gene set to which the reads were mapped in order to obtain differentially expressed genes (DEGs) between the drought and WW conditions. To investigate molecular mechanisms involved in the drought response, GO enrichment analyses were conducted. Candidate genes with a putative function in the drought response were also identified. Results A total of 52,599 gene models were obtained and after filtering on gene expression (FPKM ≥ 1), 35,733 gene models remained, of which 24,421 contained open reading frames. A total of 5,112 unique DEGs were identified between drought and WW conditions, of which the majority were found in the root tips. Phenotypically, genotype 592 displayed less growth reduction in response to drought compared to genotype 520. At the transcriptional level, genotype 520 displayed a greater response in the leaves as more DEGs were found in genotype 520 compared to genotype 592. In contrast, the transcriptional responses in the root tips were rather similar between the two genotypes. A core set of candidate genes encoding proteins with a putative function in drought response was identified, for example MYBs and bZIPs as well as chlorophyll a/b binding proteins. Discussion We found substantial differences in drought responses between the genotypes, both at the phenotypic and transcriptional levels. In addition to the genotypic variation in several traits, we also found indications for genotypic variation in trait plasticity, which could play a role in drought adaptation. Furthermore, the two genotypes displayed overall similar transcriptional responses in the root tips, but more variation in the leaves. It is thus possible that the observed phenotypic differences could be a result of transcriptional differences mostly at the leaf level. Conclusions This study has contributed to a better general understanding of drought responses in woody plants, specifically in willows, and has implications for breeding research towards more drought adapted plants. Electronic supplementary material The online version of this article (doi:10.1186/s12870-015-0630-2) contains supplementary material, which is available to authorized users.

lnPD and deltaK from the Bayesian cluster analyses with the STRUCTURE software. (PDF 672 kb)

Background: In woody plants from temperate regions, adaptation to the local climate results in annual cycles of growth and dormancy, and optimal regulation of these cycles are critical for growth, long-term survival, and competitive success. In this study we have investigated the genetic background to growth phenology in a Salix pedigree by assessing genetic and phenotypic variation in growth cessation, leaf senescence and bud burst in different years and environments. A previously constructed linkage map using the same pedigree and anchored to the annotated genome of P. trichocarpa was improved in target regions and used for QTL analysis of the traits. The major aims in this study were to map QTLs for phenology traits in Salix, and to identify candidate genes in QTL hot spots through comparative mapping with the closely related Populus trichocarpa. Results: All traits varied significantly among genotypes and the broad-sense heritabilities ranged between 0.5 and 0.9, with the highest for leaf senescence. In total across experiment and years, 80 QTLs were detected. For individual traits, the QTLs explained together from 21.5 to 56.5% of the variation. Generally each individual QTL explained a low amount of the variation but three QTLs explained above 15% of the variation with one QTL for leaf senescence explaining 34% of the variation. The majority of the QTLs were recurrently identified across traits, years and environments. Two hotspots were identified on linkage group (LG) II and X where narrow QTLs for all traits co-localized. Conclusions: This study provides the most detailed analysis of QTL detection for phenology in Salix conducted so far. Several hotspot regions were found where QTLs for different traits and QTLs for the same trait but identified during different years co-localised. Many QTLs co-localised with QTLs found in poplar for similar traits that could indicate common pathways for these traits in Salicaceae. This study is an important first step in identifying QTLs and candidate genes for phenology traits in Salix.

We present here a draft genome sequence of the red jungle fowl, Gallus gallus. Because the chicken is a modern descendant of the dinosaurs and the first non-mammalian amniote to have its genome sequenced, the draft sequence of its genome--composed of approximately one billion base pairs of sequence and an estimated 20,000-23,000 genes--provides a new perspective on vertebrate genome evolution, while also improving the annotation of mammalian genomes. For example, the evolutionary distance between chicken and human provides high specificity in detecting functional elements, both non-coding and coding. Notably, many conserved non-coding sequences are far from genes and cannot be assigned to defined functional classes. In coding regions the evolutionary dynamics of protein domains and orthologous groups illustrate processes that distinguish the lineages leading to birds and mammals. The distinctive properties of avian microchromosomes, together with the inferred patterns of conserved synteny, provide additional insights into vertebrate chromosome architecture.

By comparing sequence diversity and divergence on sex chromosomes one can study how the rate of evolution in affected by mutation and/or selection. The rate of mutation in male biased, meaning that relatively more mutations are created in the male germ line than in the female. Since the male mutation bias (αm) most likely is a consequence of the difference in cell divisions between male and female germ lines, life history characters that affect this difference should covary with αm. Indeed, we found a positive correlation between estimates of αm and increased generation times and increased intensity of sperm competition. We have also found that estimates of αm varied significantly between gametologous introns located on the sex chromosomes. This could be a consequence of the variation in substitution rates between loci. Population genetics theory predicts that both positive and negative selection reduce genetic diversity around a selected locus at a distance determined by the rate of recombination. Consequently, a non-recombining chromosome, like the female specific W chromosome in birds, selection is expected to have a large impact on sequence diversity. Indeed, in a large sequence screening we found only one segregating site among 7643 base pairs sequenced in 47 chicken females. Furthermore, we also found that deleterious substitutions are fixed in a higher rate for W- than Z-linked sequences, which is in agreement with the lack of recombination and strong genetic drift due to the low effective population size. Rarely non-synonymous mutations are beneficial for an individual, but when it happens, the mutation is positively selected and rapidly reaches fixation in a population. We have found that positive selection has been acting on the female reproductive protein, zona pellucida c in birds. This rapid evolution is likely a mechanism to prevent hybridisation.