Identification of male-specific amh duplication,sexually differentially expressed genes and microRNAs at early embryonic development of Nile tilapia (Oreochromis niloticus)

Background

The probable influence of genes and the environment on sex determination in Nile tilapia suggests that it should be regarded as a complex trait. Detection of sex determination genes in tilapia has both scientific and commercial importance. The main objective was to detect genes and microRNAs that were differentially expressed by gender in early embryonic development.

Results

Artificial fertilization of Oreochromis niloticus XX females with either sex-reversed ΔXX males or genetically-modified YY ‘supermales’ resulted in all-female and all-male embryos, respectively. RNA of pools of all-female and all-male embryos at 2, 5 and 9 dpf were used as template for a custom Agilent eArray hybridization and next generation sequencing. Fifty-nine genes differentially expressed between genders were identified by a false discovery rate of p < 0.05. The most overexpressed genes were amh and tspan8 in males, and cr/20β-hsd, gpa33, rtn4ipl and zp3 in females (p < 1 × 10⁻⁹). Validation of gene expression using qPCR in embryos and gonads indicated copy number variation in tspan8, gpa33, cr/20β-hsd and amh. Sequencing of amh identified a male-specific duplication of this gene, denoted amhy, differing from the sequence of amh by a 233 bp deletion on exonVII, hence lacking the capability to encode the protein motif that binds to the transforming growth factor beta receptor (TGF-β domain). amh and amhy segregated in the mapping family in full concordance with SD-linked marker on LG23 signifying the QTL for SD. We discovered 831 microRNAs in tilapia embryos of which nine had sexually dimorphic expression patterns by a false discovery rate of p < 0.05. An up-regulated microRNA in males, pma-mir-4585, was characterized with all six predicted target genes including cr/20β-hsd, down-regulated in males.

Conclusions

This study reports the first discovery of sexually differentially expressed genes and microRNAs at a very early stage of tilapia embryonic development, i.e. from 2 dpf. Genes with sexually differential expression patterns are enriched for copy number variation. A novel male-specific duplication of amh, denoted amhy, lacking the TGF-β domain was identified and mapped to the QTL region on LG23 for SD, thus indicating its potential role in SD.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-774) contains supplementary material, which is available to authorized users.     

Sex-dependent genetic effects on immune responses to a parasitic nematode

Background

Many disease aetiologies have sex specific effects, which have important implications for disease management. It is now becoming increasingly evident that such effects are the result of the differential expression of autosomal genes rather than sex-specific genes. Such sex-specific variation in the response to Trichuris muris, a murine parasitic nematode infection and model for the human parasitic nematode T. trichiura, has been well documented, however, the underlying genetic causes of these differences have been largely neglected. We used the BXD mouse set of recombinant inbred strains to identify sex-specific loci that contribute to immune phenotypes in T. muris infection.

Results

Response phenotypes to T. muris infection were found to be highly variable between different lines of BXD mice. A significant QTL on chromosome 5 (TM5) associated with IFN-γ production was found in male mice but not in female mice. This QTL was in the same location as a suggestive QTL for TNF-α and IL-6 production in male mice suggesting a common control of these pro-inflammatory cytokines. A second QTL was identified on chromosome 4 (TM4) affecting worm burden in both male and female cohorts. We have identified several genes as potential candidates for modifying responses to T. muris infection.

Conclusions

We have used the largest mammalian genetic model system, the BXD mouse population, to identify candidate genes with sex-specific effects in immune responses to T. muris infection. Some of these genes may be differentially expressed in male and female mice leading to the difference in immune response between the sexes reported in previous studies. Our study further highlights the importance of considering sex as an important factor in investigations of immune response at the genome-wide level, in particular the bias that can be introduced when generalizing results obtained from only one sex or a mixed sex population. Rather, analyses of interaction effects between sex and genotype should be part of future studies.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-193) contains supplementary material, which is available to authorized users.     

The genetics of domestication of rice bean, Vigna umbellata

Background and Aims

The Asian genus Vigna, to which four cultivated species (rice bean, azuki bean, mung bean and black gram) belong, is suitable for comparative genomics. The aims were to construct a genetic linkage map of rice bean, to identify the genomic regions associated with domestication in rice bean, and to compare these regions with those in azuki bean.

Methods

A genetic linkage map was constructed by using simple sequence repeat and amplified fragment length polymorphism markers in the BC₁F₁ population derived from a cross between cultivated and wild rice bean. Using this map, 31 domestication-related traits were dissected into quantitative trait loci (QTLs). The genetic linkage map and QTLs of rice bean were compared with those of azuki bean.

Key Results

A total of 326 markers converged into 11 linkage groups (LGs), corresponding to the haploid number of rice bean chromosomes. The domestication-related traits in rice bean associated with a few major QTLs distributed as clusters on LGs 2, 4 and 7. A high level of co-linearity in marker order between the rice bean and azuki bean linkage maps was observed. Major QTLs in rice bean were found on LG4, whereas major QTLs in azuki bean were found on LG9.

Conclusions

This is the first report of a genetic linkage map and QTLs for domestication-related traits in rice bean. The inheritance of domestication-related traits was so simple that a few major QTLs explained the phenotypic variation between cultivated and wild rice bean. The high level of genomic synteny between rice bean and azuki bean facilitates QTL comparison between species. These results provide a genetic foundation for improvement of rice bean; interchange of major QTLs between rice bean and azuki bean might be useful for broadening the genetic variation of both species.     

QTL analysis of heterostyly in Primula sieboldii and its application for morph identification in wild populations

Background and Aims

Primula sieboldii is a perennial clonal herb that is distributed around the Sea of Japan and is endangered in Japan. Its breeding system is characterized by heteromorphic self-incompatibility, and the morph ratio within a population is very important for reproductive success. The aims of this study were to construct a linkage map, map the S locus as a qualitative trait and quantitative trait loci (QTLs) for floral morphological traits related to heterostyly, and predict the morph type in wild populations by using molecular markers for devising a conservation strategy.

Methods

A linkage map was constructed with 126 markers. The QTLs for four floral traits and the S locus were mapped. Using the genotypes of loci that were located near both the S locus and the QTLs with large effects, morphs of 59 wild genets were predicted.

Key Results

The linkage map consisted of 14 linkage groups (LGs). The S locus was mapped to LG 7. Major QTLs for stigma and anther heights were detected in the same region as the S locus. These QTLs exhibited high logarithm of the odds scores and explained a high percentage of the phenotypic variance (>85 %). By analysing these two traits within each morph, additional QTLs for each trait were detected. Using the four loci linked to the S locus, the morphs of 43 genets in three wild populations could be predicted.

Conclusions

This is the first report of a linkage map and QTL analysis for floral morphology related to heterostyly in P. sieboldii. Floral morphologies related to heterostyly are controlled by the S locus in LG 7 and by several QTLs in other LGs. Additionally, this study showed that molecular markers are effective tools for investigating morph ratios in a population containing the non-flowering individuals or during the non-flowering seasons.     

Germination variation in Arabidopsis thaliana accessions under moderate osmotic and salt stresses

Background

Water and salt stresses are two important environmental factors that limit the germination of seeds in most ecological environments. Most studies conducted so far to address the genetic basis of the above phenomenon have used stress conditions that are much more extreme than those found in natural environments. Furthermore, although an excess of ions and water restrictions have similar osmotic effects on germination, the common and divergent signalling components mediating the effects of both factors remain unknown.

Methods

The germination of seeds was compared under solutions of NaCl (50 mm) and polyethylene glycol (PEG, −0·6 MPa), that establish mild stress conditions, in 28 Arabidopsis thaliana accessions. Because Bayreuth (Bay) and Shadara (Sha) accessions showed contrasting sensitivity responses to both stresses, a quantitative trait locus (QTL) analysis was carried out using Bay × Sha recombinant inbred lines (RILs) to identify loci involved in the control of germination under mild salt and osmotic stresses.

Key Results

Two loci associated with the salt sensitivity response, named SSR1 and SSR2 QTLs, and four loci for the osmotic sensitivity response, named OSR1–OSR4 QTLs, were mapped. The effects of the SSR1 QTL on toxic salt sensitivity, and the osmotic contribution of OSR1, were confirmed by heterogeneous inbred families (HIFs). Whilst the SSR1 QTL had a significant effect under a wide range of NaCl concentrations, the OSR1 QTL was confirmed only under moderate drought stress. Interestingly the OSR1 QTL also showed pleiotropic effects on biomass accumulation in response to water deficit.

Conclusions

The regulation of germination under moderate salt and osmotic stresses involves the action of independent major loci, revealing the existence of loci specifically associated with the toxic component of salt and not just its osmotic effect. Furthermore, this work demonstrates that novel loci control germination under osmotic stress conditions simulating more realistic ecological environments as found by populations of seeds in nature.     

Maternal sex effects and inbreeding depression under varied environmental conditions in gynodioecious Fragaria vesca subsp. bracteata

Background and Aims

Gynodioecy (coexistence of females and hermaphrodites) is a sexual system that occurs in numerous flowering plant lineages. Thus, understanding the features that affect its maintenance has wide importance. Models predict that females must have a seed fitness advantage over hermaphrodites, and this may be achieved via seed quality or quantity. Females in a population of Fragaria vesca subsp. bracteata, a long-lived gynodioecious perennial, do not demonstrate a seed quantity advantage, so this study explored whether females produced better quality seed via maternal sex effects or avoidance of inbreeding depression (IBD).

Methods

Families of selfed and outcrossed seed were created using hermaphrodite mothers and families of outcrossed seed were created using female mothers. The effects of these pollination treatments were assessed under benign conditions early in life and under varied conditions later in life. To test for an effect of maternal sex, fitness components and traits associated with acclimation to variable environments of progeny of outbred hermaphrodites and females were compared. To test for expression of IBD, fitness parameters between inbred and outbred progeny of hermaphrodites were compared.

Key Results

Offspring of females were more likely to germinate in benign conditions and survive in harsh resource environments than outbred progeny of hermaphrodites. IBD was low across most life stages, and both the effect of maternal sex on progeny quality and the expression of IBD depended on both maternal family and resource condition of the progeny.

Conclusions

The effect of maternal sex and IBD on progeny quality depended on resource conditions, maternal lineage and progeny life stage. In conjunction with known lack of differences in seed quantity, the quality advantages and IBD observed here are still unlikely to be sufficient for maintenance of gynodioecy under nuclear inheritance of male sterility.     

Sexual dimorphism in a dioecious population of the wind-pollinated herb Mercurialis annua: the interactive effects of resource availability and competition

Background and Aims

Male-biased sex allocation commonly occurs in wind-pollinated hermaphroditic plants, and is often positively associated with size, notably in terms of height. Currently, it is not well established whether a corresponding pattern holds for dioecious plants: do males of wind-pollinated species exhibit greater reproductive allocation than females? Here, sexual dimorphism is investigated in terms of life history trade-offs in a dioecious population of the wind-pollinated ruderal herb Mercurialis annua.

Methods

The allocation strategies of males and females grown under different soil nutrient availability and competitive (i.e. no, male or female competitor) regimes were compared.

Key Results

Male reproductive allocation increased disproportionately with biomass, and was greater than that of females when grown in rich soils. Sexual morphs differentially adjusted their reproductive allocation in response to local environmental conditions. In particular, males reduced their reproductive allocation in poor soils, whereas females increased theirs, especially when competing with another female rather than growing alone. Finally, males displayed smaller above-ground vegetative sizes than females, but neither nutrient availability nor competition had a strong independent effect on relative size disparities between the sexes.

Conclusions

Selection appears to favour plasticity in reproductive allocation in dioecious M. annua, thereby maintaining a relatively constant size hierarchy between sexual morphs. In common with other dioecious species, there seems to be little divergence in the niches occupied by males and females of M. annua.     

A new SNP-based vision of the genetics of sex determination in European sea bass (Dicentrarchus labrax)

Background

European sea bass (Dicentrarchus labrax) is one of the most important farmed species in Mediterranean aquaculture. The observed sexual growth and maturity dimorphism in favour of females adds value towards deciphering the sex determination system of this species. Current knowledge indicates the existence of a polygenic sex determining determination system that interacts with temperature. This was explored by restriction-site associated DNA (RAD) marker analysis in a test panel of 175 offspring that originated from a factorial cross between two dams and four sires from a single full-sib family.

Results

The first high-density single nucleotide polymorphism (SNP) based linkage map for sea bass was constructed, consisting of 6706 SNPs on 24 linkage groups. Indications for putative sex-determining QTL (quantitative trait loci) that were significant at the genome-wide threshold were detected on linkage groups 6, 11 and 18 to 21, although a genome-wide association study (GWAS) did not identify individual significant SNPs at a genome-wide threshold. A preliminary genomic prediction approach that tested the efficiency of SNP-based selection for female sea bass showed a slight advantage compared to traditional pedigree-based selection. However, when the same models were tested on the same animals for selection for greater length, a clear advantage of the SNP-based selection was observed.

Conclusions

Overall, the results of this study provide additional support to the polygenic sex determination hypothesis in sea bass. In addition, identification of sex-ratio QTL may provide new opportunities for sex-ratio control in sea bass.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0148-y) contains supplementary material, which is available to authorized users.     

High-throughput root phenotyping screens identify genetic loci associated with root architectural traits in Brassica napus under contrasting phosphate availabilities

Background and Aims

Phosphate (Pi) deficiency in soils is a major limiting factor for crop growth worldwide. Plant growth under low Pi conditions correlates with root architectural traits and it may therefore be possible to select these traits for crop improvement. The aim of this study was to characterize root architectural traits, and to test quantitative trait loci (QTL) associated with these traits, under low Pi (LP) and high Pi (HP) availability in Brassica napus.

Methods

Root architectural traits were characterized in seedlings of a double haploid (DH) mapping population (n = 190) of B. napus [‘Tapidor’ × ‘Ningyou 7’ (TNDH)] using high-throughput phenotyping methods. Primary root length (PRL), lateral root length (LRL), lateral root number (LRN), lateral root density (LRD) and biomass traits were measured 12 d post-germination in agar at LP and HP.

Key Results

In general, root and biomass traits were highly correlated under LP and HP conditions. ‘Ningyou 7’ had greater LRL, LRN and LRD than ‘Tapidor’, at both LP and HP availability, but smaller PRL. A cluster of highly significant QTL for LRN, LRD and biomass traits at LP availability were identified on chromosome A03; QTL for PRL were identified on chromosomes A07 and C06.

Conclusions

High-throughput phenotyping of Brassica can be used to identify root architectural traits which correlate with shoot biomass. It is feasible that these traits could be used in crop improvement strategies. The identification of QTL linked to root traits under LP and HP conditions provides further insights on the genetic basis of plant tolerance to P deficiency, and these QTL warrant further dissection.     

On the roles of colour and scent in a specialized floral mimicry system

Background and Aims

Sexually deceptive orchids achieve cross-pollination by mimicking the mating signals of female insects, generally hymenopterans. This pollination mechanism is often highly specific as it is based primarily on the mimicry of mating signals, especially the female sex pheromones of the targeted pollinator. Like many deceptive orchids, the Mediterranean species Ophrys arachnitiformis shows high levels of floral trait variation, especially in the colour of the perianth, which is either green or white/pinkinsh within populations. The adaptive significance of perianth colour polymorphism and its influence on pollinator visitation rates in sexually deceptive orchids remain obscure.

Methods

The relative importance of floral scent versus perianth colour in pollinator attraction in this orchid pollinator mimicry system was evaluated by performing floral scent analyses by gas chromatography-mass spectrometry (GC-MS) and behavioural bioassays with the pollinators under natural conditions were performed.

Key Results

The relative and absolute amounts of behaviourally active compounds are identical in the two colour morphs of O. arachnitiformis. Neither presence/absence nor the colour of the perianth (green versus white) influence attractiveness of the flowers to Colletes cunicularius males, the main pollinator of O. arachnitiformis.

Conclusion

Chemical signals alone can mediate the interactions in highly specialized mimicry systems. Floral colour polymorphism in O. arachnitiformis is not subjected to selection imposed by C. cunicularius males, and an interplay between different non-adaptive processes may be responsible for the maintenance of floral colour polymorphism both within and among populations.     

Limits on the reproducibility of marker associations with southern leaf blight resistance in the maize nested association mapping population

Background

A previous study reported a comprehensive quantitative trait locus (QTL) and genome wide association study (GWAS) of southern leaf blight (SLB) resistance in the maize Nested Association Mapping (NAM) panel. Since that time, the genomic resources available for such analyses have improved substantially. An updated NAM genetic linkage map has a nearly six-fold greater marker density than the previous map and the combined SNPs and read-depth variants (RDVs) from maize HapMaps 1 and 2 provided 28.5 M genomic variants for association analysis, 17 fold more than HapMap 1. In addition, phenotypic values of the NAM RILs were re-estimated to account for environment-specific flowering time covariates and a small proportion of lines were dropped due to genotypic data quality problems. Comparisons of original and updated QTL and GWAS results confound the effects of linkage map density, GWAS marker density, population sample size, and phenotype estimates. Therefore, we evaluated the effects of changing each of these parameters individually and in combination to determine their relative impact on marker-trait associations in original and updated analyses.

Results

Of the four parameters varied, map density caused the largest changes in QTL and GWAS results. The updated QTL model had better cross-validation prediction accuracy than the previous model. Whereas joint linkage QTL positions were relatively stable to input changes, the residual values derived from those QTL models (used as inputs to GWAS) were more sensitive, resulting in substantial differences between GWAS results. The updated NAM GWAS identified several candidate genes consistent with previous QTL fine-mapping results.

Conclusions

The highly polygenic nature of resistance to SLB complicates the identification of causal genes. Joint linkage QTL are relatively stable to perturbations of data inputs, but their resolution is generally on the order of tens or more Mbp. GWAS associations have higher resolution, but lower power due to stringent thresholds designed to minimize false positive associations, resulting in variability of detection across studies. The updated higher density linkage map improves QTL estimation and, along with a much denser SNP HapMap, greatly increases the likelihood of detecting SNPs in linkage with causal variants. We recommend use of the updated genetic resources and results but emphasize the limited repeatability of small-effect associations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1068) contains supplementary material, which is available to authorized users.     

Genetic architecture of resistance to Septoria tritici blotch in European wheat

Background

Septoria tritici blotch is an important leaf disease of European winter wheat. In our survey, we analyzed Septoria tritici blotch resistance in field trials with a large population of 1,055 elite hybrids and their 87 parental lines. Entries were fingerprinted with the 9 k SNP array. The accuracy of prediction of Septoria tritici blotch resistance achieved with different genome-wide mapping approaches was evaluated based on robust cross validation scenarios.

Results

Septoria tritici blotch disease severities were normally distributed, with genotypic variation being significantly (P < 0.01) larger than zero. The cross validation study revealed an absence of large effect QTL for additive and dominance effects. Application of genomic selection approaches particularly designed to tackle complex agronomic traits allowed to double the accuracy of prediction of Septoria tritici blotch resistance compared to calculation methods suited to detect QTL with large effects.

Conclusions

Our study revealed that Septoria tritici blotch resistance in European winter wheat is controlled by multiple loci with small effect size. This suggests that the currently achieved level of resistance in this collection is likely to be durable, as involvement of a high number of genes in a resistance trait reduces the risk of the resistance to be overcome by specific pathogen isolates or races.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-858) contains supplementary material, which is available to authorized users.     

Sexual differences in physiological integration in the dioecious shrub Lindera triloba: a field experiment using girdling manipulation

Background and Aims

It is important to consider the modular level when verifying sexual dimorphism in dioecious plants. Nevertheless, between-sex differences in resource translocation among modules (i.e. physiological integration) have not been tested at the whole-plant level. In this study, sexual differences in physiological integration were examined among ramets, within a genet in the dioecious sprouting shrub Lindera triloba, by a field experiment with girdling manipulation.

Methods

Female and male genets were randomly assigned to girdled or intact groups. Girdling of the main ramets was conducted in May 2009 by removing a ring of bark and cambium approx. 1 cm wide at a height of 80–100 cm. The effects of treatment and sex on ramet dynamics (mortality, recruitment and diameter growth) and inflorescence production during 1 year after girdling were examined.

Key Results

The diameter growth rate of main ramets of both sexes was lower at ground level (D₀) but higher at breast height (dbh) in girdled than in intact groups. In sprouted ramets with a dbh of 0–2 cm, males in girdled groups had lower growth rates at D₀ than those of intact groups, whereas no girdling effect was found for females. The main ramets in girdled groups produced more inflorescences than intact groups, irrespective of sex, but male ramets showed a greater response to the treatment than females.

Conclusions

In L. triloba, physiological integration exists at the whole-plant level, and sprouted ramets are dependent on assimilates translocated from main ramets, but this dependence weakens as sprouted ramets get larger. Female sprouted ramets can grow in a physiologically independent manner from the main ramet earlier than those of males. This study highlights the importance of considering modular structures and physiological integration when evaluating sexual differences in demographic patterns of clonal plants.     

Mapping of a milk production quantitative trait locus to a 1.056 Mb region on bovine chromosome 5 in the Fleckvieh dual purpose cattle breed

Background

In a previous study in the Fleckvieh dual purpose cattle breed, we mapped a quantitative trait locus (QTL) affecting milk yield (MY1), milk protein yield (PY1) and milk fat yield (FY1) during first lactation to the distal part of bovine chromosome 5 (BTA5), but the confidence interval was too large for positional cloning of the causal gene. Our objective here was to refine the position of this QTL and to define the candidate region for high-throughput sequencing.

Methods

In addition to those previously studied, new Fleckvieh families were genotyped, in order to increase the number of recombination events. Twelve new microsatellites and 240 SNP markers covering the most likely QTL region on BTA5 were analysed. Based on haplotype analysis performed in this complex pedigree, families segregating for the low frequency allele of this QTL (minor allele) were selected. Single- and multiple-QTL analyses using combined linkage and linkage disequilibrium methods were performed.

Results

Single nucleotide polymorphism haplotype analyses on representative family sires and their ancestors revealed that the haplotype carrying the minor QTL allele is rare and most probably originates from a unique ancestor in the mapping population. Analyses of different subsets of families, created according to the results of haplotype analysis and availability of SNP and microsatellite data, refined the previously detected QTL affecting MY1 and PY1 to a region ranging from 117.962 Mb to 119.018 Mb (1.056 Mb) on BTA5. However, the possibility of a second QTL affecting only PY1 at 122.115 Mb was not ruled out.

Conclusion

This study demonstrates that targeting families segregating for a less frequent QTL allele is a useful method. It improves the mapping resolution of the QTL, which is due to the division of the mapping population based on the results of the haplotype analysis and to the increased frequency of the minor allele in the families. Consequently, we succeeded in refining the region containing the previously detected QTL to 1 Mb on BTA5. This candidate region contains 27 genes with unknown or partially known function(s) and is small enough for high-throughput sequencing, which will allow future detailed analyses of candidate genes.     

Does the silver moss Bryum argenteum exhibit sex-specific patterns in vegetative growth rate,asexual fitness or prezygotic reproductive investment?

Background and Aims

Expected life history trade-offs associated with sex differences in reproductive investment are often undetected in seed plants, with the difficulty arising from logistical issues of conducting controlled experiments. By controlling genotype, age and resource status of individuals, a bryophyte was assessed for sex-specific and location-specific patterns of vegetative, asexual and sexual growth/reproduction across a regional scale.

Methods

Twelve genotypes (six male, six female) of the dioecious bryophyte Bryum argenteum were subcultured to remove environmental effects, regenerated asexually to replicate each genotype 16 times, and grown over a period of 92 d. Plants were assessed for growth rates, asexual and sexual reproductive traits, and allocation to above- and below-ground regenerative biomass.

Key Results

The degree of sexual versus asexual reproductive investment appears to be under genetic control, with three distinct ecotypes found in this study. Protonemal growth rate was positively correlated with asexual reproduction and sexual reproduction, whereas asexual reproduction was negatively correlated (appeared to trade-off) with vegetative growth (shoot production). No sex-specific trade-offs were detected. Female sex-expressing shoots were longer than males, but the sexes did not differ in growth traits, asexual traits, sexual induction times, or above- and below-ground biomass. Males, however, had much higher rates of inflorescence production than females, which translated into a significantly higher (24x) prezygotic investment for males relative to females.

Conclusions

Evidence for three distinct ecotypes is presented for a bryophyte based on regeneration traits. Prior to zygote production, the sexes of this bryophyte did not differ in vegetative growth traits but significantly differed in reproductive investment, with the latter differences potentially implicated in the strongly biased female sex ratio. The disparity between males and females for prezygotic reproductive investment is the highest known for bryophytes.     

Gender plasticity and sexual system stability in Wurmbea

Background and aims

Sexually dimorphic populations are often located in drier habitats than cosexual populations. Gender plasticity (GP), whereby hermaphrodites alter female and male functions depending on resources, and sex-differential plasticity (SDP) between hermaphrodites and unisexuals are predicted to affect sexual system stability. Here, GP and SDP are evaluated in cosexual and gynodioecious Wurmbea biglandulosa and sub-dioecious and dioecious W. dioica.

Methods

GP was evaluated under two resource conditions, compared among sexual systems and assessed as to whether (1) males produced perfect flowers and (2) hermaphrodites altered investment in perfect (female function) and total (male function) flowers. SDP was assessed within sexual systems as differences between sex functions of hermaphrodites vs. unisexuals. Males and hermaphrodites were compared to assess whether size thresholds for female function differed among sexual systems. Plasticity costs were evaluated using correlations between female function and male traits in hermaphrodites, and in W. dioica by comparing hermaphrodite and male regressions between plant size and pollen production.

Key Results

In dioecious W. dioica no males exhibited GP, whereas 100 % did in gynodioecious and cosexual W. biglandulosa. In sub-dioecious W. dioica, resources affected GP (high, 66 %; low, 42 %). Hermaphrodites in all sexual systems reduced perfect but not total flowers under low resources. Unisexuals were unaffected, demonstrating SDP for female function only. Thresholds for female function were greater in sub-dioecious W. dioica than in W. biglandulosa. Plasticity costs were detected only in sub-dioecious W. dioica.

Conclusions

SDP for female function could assist female establishment in cosexual populations and maintain females in gynodioecious and sub-dioecious populations. Although the absence of male SDP should stabilize sub-dioecy, plasticity costs would render sub-dioecy unstable, favouring canalized males over hermaphrodites. This study highlights the importance of interactions between environmental conditions and hermaphrodite sex expression for the stability of dimorphic sexual systems.