Complementing model species with model clades.

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant Tree of Life continues to improve. The intersection of these two research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade". These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis thaliana and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

Sex-linked gene expression and the emergence of hermaphrodites in Carica papaya.

PREMISE: One evolutionary path from hermaphroditism to dioecy is via a gynodioecious intermediate. The evolution of dioecy may also coincide with the formation of sex chromosomes that possess sex-determining loci that are physically linked in a region of suppressed recombination. Dioecious papaya (Carica papaya) has an XY chromosome system, where the presence of a Y chromosome determines maleness. However, in cultivation, papaya is gynodioecious, due to the conversion of the male Y chromosome to a hermaphroditic Yh chromosome during its domestication. METHODS: We investigated gene expression linked to the X, Y, and Yh chromosomes at different floral developmental stages to identify differentially expressed genes that may be involved in the sexual transition of males to hermaphrodites. RESULTS: We identified 309 sex-biased genes found on the sex chromosomes, most of which are found in the pseudoautosomal regions. Female (XX) expression in the sex-determining region was almost double that of X-linked expression in males (XY) and hermaphrodites (XYh ), which rules out dosage compensation for most sex-linked genes; although, an analysis of hemizygous X-linked loci found evidence of partial dosage compensation. Furthermore, we identified a candidate gene associated with sex determination and the transition to hermaphroditism, a homolog of the MADS-box protein SHORT VEGETATIVE PHASE. CONCLUSIONS: We identified a pattern of partial dosage compensation for hemizygous genes located in the papaya sex-determining region. Furthermore, we propose that loss-of-expression of the Y-linked SHORT VEGETATIVE PHASE homolog facilitated the transition from males to hermaphrodites in papaya.

Origin and domestication of papaya Yh chromosome.

Sex in papaya is controlled by a pair of nascent sex chromosomes. Females are XX, and two slightly different Y chromosomes distinguish males (XY) and hermaphrodites (XY(h)). The hermaphrodite-specific region of the Y(h) chromosome (HSY) and its X chromosome counterpart were sequenced and analyzed previously. We now report the sequence of the entire male-specific region of the Y (MSY). We used a BAC-by-BAC approach to sequence the MSY and resequence the Y regions of 24 wild males and the Y(h) regions of 12 cultivated hermaphrodites. The MSY and HSY regions have highly similar gene content and structure, and only 0.4% sequence divergence. The MSY sequences from wild males include three distinct haplotypes, associated with the populations' geographic locations, but gene flow is detected for other genomic regions. The Y(h) sequence is highly similar to one Y haplotype (MSY3) found only in wild dioecious populations from the north Pacific region of Costa Rica. The low MSY3-Y(h) divergence supports the hypothesis that hermaphrodite papaya is a product of human domestication. We estimate that Y(h) arose only ∼ 4000 yr ago, well after crop plant domestication in Mesoamerica >6200 yr ago but coinciding with the rise of the Maya civilization. The Y(h) chromosome has lower nucleotide diversity than the Y, or the genome regions that are not fully sex-linked, consistent with a domestication bottleneck. The identification of the ancestral MSY3 haplotype will expedite investigation of the mutation leading to the domestication of the hermaphrodite Y(h) chromosome. In turn, this mutation should identify the gene that was affected by the carpel-suppressing mutation that was involved in the evolution of males.

A wild origin of the loss-of-function lycopene beta cyclase (CYC-b) allele in cultivated, red-fleshed papaya (Carica papaya).

PREMISE OF THE STUDY: The red flesh of some papaya cultivars is caused by a recessive loss-of-function mutation in the coding region of the chromoplast-specific lycopene beta cyclase gene (CYC-b). We performed an evolutionary genetic analysis of the CYC-b locus in wild and cultivated papaya to uncover the origin of this loss-of-function allele in cultivated papaya. METHODS: We analyzed the levels and patterns of genetic diversity at the CYC-b locus and six loci in a 100-kb region flanking CYC-b and compared these to genetic diversity levels at neutral autosomal loci. The evolutionary relationships of CYC-b haplotypes were assessed using haplotype network analysis of the CYC-b locus and the 100-kb CYC-b region. KEY RESULTS: Genetic diversity at the recessive CYC-b allele (y) was much lower relative to the dominant Y allele found in yellow-fleshed wild and cultivated papaya due to a strong selective sweep. Haplotype network analyses suggest the y allele most likely arose in the wild and was introduced into domesticated varieties after the first papaya domestication event. The shared haplotype structure between some wild, feral, and cultivated haplotypes around the y allele supports subsequent escape of this allele from red cultivars back into wild populations through feral intermediates. CONCLUSIONS: Our study supports a protracted domestication process of papaya through the introgression of wild-derived traits and gene flow from cultivars to wild populations. Evidence of gene flow from cultivars to wild populations through feral intermediates has implications for the introduction of transgenic papaya into Central American countries.

The Evolutionary Tempo of Sex Chromosome Degradation in Carica papaya.

Genes on non-recombining heterogametic sex chromosomes may degrade over time through the irreversible accumulation of deleterious mutations. In papaya, the non-recombining male-specific region of the Y (MSY) consists of two evolutionary strata corresponding to chromosomal inversions occurring approximately 7.0 and 1.9 MYA. The step-wise recombination suppression between the papaya X and Y allows for a temporal examination of the degeneration progress of the young Y chromosome. Comparative evolutionary analyses of 55 X/Y gene pairs showed that Y-linked genes have more unfavorable substitutions than X-linked genes. However, this asymmetric evolutionary pattern is confined to the oldest stratum, and is only observed when recently evolved pseudogenes are included in the analysis, indicating a slow degeneration tempo of the papaya Y chromosome. Population genetic analyses of coding sequence variation of six Y-linked focal loci in the oldest evolutionary stratum detected an excess of nonsynonymous polymorphism and reduced codon bias relative to autosomal loci. However, this pattern was also observed for corresponding X-linked loci. Both the MSY and its corresponding X-specific region are pericentromeric where recombination has been shown to be greatly reduced. Like the MSY region, overall selective efficacy on the X-specific region may be reduced due to the interference of selective forces between highly linked loci, or the Hill-Robertson effect, that is accentuated in regions of low or suppressed recombination. Thus, a pattern of gene decay on the X-specific region may be explained by relaxed purifying selection and widespread genetic hitchhiking due to its pericentromeric location.

A polymorphic pseudoautosomal boundary in the Carica papaya sex chromosomes.

Sex chromosomes are defined by a non-recombining sex-determining region (SDR) flanked by one or two pseudoautosomal regions (PARs). The genetic composition and evolutionary dynamics of the PAR is also influenced by its linkage to the differentiated non-recombining SDR; however, understanding the effects of this linkage requires a precise definition of the PAR boundary. Here, we took a molecular population genetic approach to further refine the location of the PAR boundary of the evolutionary young sex chromosomes of the tropical plant, Carica papaya. We were able to map the position of the papaya PAR boundary A to a 100-kb region between two genetic loci approximately 2 Mb upstream of the previously genetically identified PAR boundary. Furthermore, this boundary is polymorphic within natural populations of papaya, with an approximately 100-130 kb expansion of the non-recombining SDR found in 16 % of individuals surveyed. The expansion of the PAR boundary in one Y haplotype includes at least one additional gene. Homologs of this gene are involved in male gametophyte and pollen development in other plant species.

Multiple developmental processes underlie sex differentiation in angiosperms.

The production of unisexual flowers has evolved numerous times in dioecious and monoecious plant taxa. Based on repeated evolutionary origins, a great variety of developmental and genetic mechanisms underlying unisexual flower development is predicted. Here, we comprehensively review the modes of development of unisexual flowers, test potential correlations with sexual system, and end with a synthesis of the genetics and hormonal regulation of plant sex determination. We find that the stage of organ abortion in male and female flowers is temporally correlated within species and also confirm that the arrest of development does not tend to occur preferentially at a particular stage, or via a common process.

Contrasting patterns of X/Y polymorphism distinguish Carica papaya from other sex chromosome systems.

The sex chromosomes of the tropical crop papaya (Carica papaya) are evolutionarily young and consequently allow for the examination of evolutionary mechanisms that drive early sex chromosome divergence. We conducted a molecular population genetic analysis of four X/Y gene pairs from a collection of 45 wild papaya accessions. These population genetic analyses reveal striking differences in the patterns of polymorphism between the X and Y chromosomes that distinguish them from other sex chromosome systems. In most sex chromosome systems, the Y chromosome displays significantly reduced polymorphism levels, whereas the X chromosome maintains a level of polymorphism that is comparable to autosomal loci. However, the four papaya sex-linked loci that we examined display diversity patterns that are opposite this trend: the papaya X alleles exhibit significantly reduced polymorphism levels, whereas the papaya Y alleles maintain greater than expected levels of diversity. Our analyses suggest that selective sweeps in the regions of the X have contributed to this pattern while also revealing geographically restricted haplogroups on the Y. We discuss the possible role sexual selection and/or genomic conflict have played in shaping the contrasting patterns of polymorphism found for the papaya X and Y chromosomes.

Asymmetric functional divergence of young, dispersed gene duplicates in Arabidopsis thaliana.

One prediction of the classic Ohno model of gene duplication predicts that new genes form from the asymmetric functional divergence of a newly arisen, redundant duplicate locus. In order to understand the mechanisms which give rise to functional divergence of newly formed dispersed duplicates, we assessed the expression and molecular evolutionary divergence of a suite of 19 highly similar dispersed duplicates in Arabidopsis thaliana. These duplicates have a K sil equal to or less than 5 % and are specific to the A. thaliana lineage; thus, they predictably represent some of the youngest duplicates in the A. thaliana genome. We found that the majority of young duplicate loci exhibit asymmetric expression patterns, with the daughter locus exhibiting reduced expression across all tissues analyzed relative to the progenitor locus or simply not expressed. Furthermore, daughter loci, on the whole, have significantly more nonsynonymous substitutions than the progenitor loci. We also identified four pairs of loci which exhibit significant (P < 0.05) evolutionary rate asymmetry, three of which exhibit elevated dN/dS in the duplicate copy. We suggest, based on these data, that functional diversification initially takes the form of asymmetric regulatory divergence that can be a direct consequence of the mode of duplication. The reduced and/or absence of expression in the daughter copy relaxes functional constraint on its protein coding sequence leading to the asymmetric accumulation of nonsynonymous mutations. Thus, our data both affirm Ohno's prediction while explaining the mechanism by which functional divergence initially occurs following duplication for dispersed gene duplicates.

The sex-specific region of sex chromosomes in animals and plants.

Our understanding of the evolution of sex chromosomes has increased greatly in recent years due to a number of molecular evolutionary investigations in divergent sex chromosome systems, and these findings are reshaping theories of sex chromosome evolution. In particular, the dynamics of the sex-determining region (SDR) have been demonstrated by recent findings in ancient and incipient sex chromosomes. Radical changes in genomic structure and gene content in the male specific region of the Y chromosome between human and chimpanzee indicated rapid evolution in the past 6 million years, defying the notion that the pace of evolution in the SDR was fast at early stages but slowed down overtime. The chicken Z and the human X chromosomes appeared to have acquired testis-expressed genes and expanded in intergenic regions. Transposable elements greatly contributed to SDR expansion and aided the trafficking of genes in the SDR and its X or Z counterpart through retrotransposition. Dosage compensation is not a destined consequence of sex chromosomes as once thought. Most X-linked microRNA genes escape silencing and are expressed in testis. Collectively, these findings are challenging many of our preconceived ideas of the evolutionary trajectory and fates of sex chromosomes.

High Prevalence of Clonal Reproduction and Low Genetic Diversity in Scutellaria floridana, a Federally Threatened Florida-Endemic Mint.

The threatened mint Florida skullcap (Scutellaria floridana) is endemic to four counties in the Florida panhandle. Because development and habitat modification extirpated several historical occurrences, only 19 remain to date. To inform conservation management and delisting decisions, a comprehensive investigation of the genetic diversity and relatedness, population structure, and clonal diversity was conducted using SNP data generated by ddRAD. Compared with other Lamiaceae, we detected low genetic diversity (HE = 0.125-0.145), low to moderate evidence of inbreeding (FIS = -0.02-0.555), and moderate divergence (FST = 0.05-0.15). We identified eight populations with most of the genetic diversity, which should be protected in situ, and four populations with low genetic diversity and high clonality. Clonal reproduction in our circular plots and in 92% of the sites examined was substantial, with average clonal richness of 0.07 and 0.59, respectively. Scutellaria floridana appears to have experienced a continued decline in the number of extant populations since its listing under the Endangered Species Act; still, the combination of sexual and asexual reproduction may be advantageous for maintaining the viability of extant populations. However, the species will likely require ongoing monitoring, management, and increased public awareness to ensure its survival and effectively conserve its genetic diversity.

The origin of the non-recombining region of sex chromosomes in Carica and Vasconcellea.

Carica and Vasconcellea are two closely related sister genera in the family Caricaceae, and were once classified as two sections under Carica. Sex chromosomes have been found in papaya and originated approximately 2-3 million years ago. The objectives of this study were to determine whether sex chromosomes have evolved in Vasconcellea. Six X/Y gene pairs were cloned, sequenced and analyzed from three dioecious, one trioecious and one monoecious species of Vasconcellea. The isolation of distinctive X and Y alleles in dioecious and trioecious species of Vasconcellea demonstrated that sex chromosomes have evolved in this genus. Phylogenetic analyses indicated a monophyletic relationship between the X/Y alleles of Carica and those of Vasconcellea. Distinctive clusters of X/Y alleles were documented in V. parviflora and V. pulchra for all available gene sequences, and in V. goudatinana and V. cardinamarcensis for some X/Y alleles. The X and Y alleles within each species shared most single nucleotide polymorphism haplotypes that differed from other species. Limited evidence of gene conversion was documented among the X/Y alleles of some species, but was not sufficient to cause the evolutionary patterns reported herein. The Carica and Vasconcellea sex chromosomes may have originated from the same autosomes bearing the X allelic form that still exist in the monoecious species V. monoica, and have evolved independently after the speciation event that separated Carica from Vasconcellea. Within Vasconcellea, sex chromosomes have evolved at the species level, at least for some species.

Population genomics of the Arabidopsis thaliana flowering time gene network.

The time to flowering is a key component of the life-history strategy of the model plant Arabidopsis thaliana that varies quantitatively among genotypes. A significant problem for evolutionary and ecological genetics is to understand how natural selection may operate on this ecologically significant trait. Here, we conduct a population genomic study of resequencing data from 52 genes in the flowering time network. McDonald-Kreitman tests of neutrality suggested a strong excess of amino acid polymorphism when pooling across loci. This excess of replacement polymorphism across the flowering time network and a skewed derived frequency spectrum toward rare alleles for both replacement and noncoding polymorphisms relative to synonymous changes is consistent with a large class of deleterious polymorphisms segregating in these genes. Assuming selective neutrality of synonymous changes, we estimate that approximately 30% of amino acid polymorphisms are deleterious. Evidence of adaptive substitution is less prominent in our analysis. The photoperiod regulatory gene, CO, and a gibberellic acid transcription factor, AtMYB33, show evidence of adaptive fixation of amino acid mutations. A test for extended haplotypes revealed no examples of flowering time alleles with haplotypes comparable in length to those associated with the null fri(Col) allele reported previously. This suggests that the FRI gene likely has a uniquely intense or recent history of selection among the flowering time genes considered here. Although there is some evidence for adaptive evolution in these life-history genes, it appears that slightly deleterious polymorphisms are a major component of natural molecular variation in the flowering time network of A. thaliana.

The evolutionary dynamics of plant duplicate genes.

Given the prevalence of duplicate genes and genomes in plant species, the study of their evolutionary dynamics has been a focus of study in plant evolutionary genetics over the past two decades. The past few years have been a particularly exciting time because recent theoretical and experimental investigations have led to a rethinking of the classic paradigm of duplicate gene evolution. By combining recent advances in genomic analysis with a new conceptual framework, researchers are determining the contributions of single-gene and whole-genome duplications to the diversification of plant species. This research provides insights into the roles that gene and genome duplications play in plant evolution.

Extremely low nucleotide diversity in the X-linked region of papaya caused by a strong selective sweep.

BACKGROUND: The papaya Y-linked region showed clear population structure, resulting in the detection of the ancestral male population that domesticated hermaphrodite papayas were selected from. The same populations were used to study nucleotide diversity and population structure in the X-linked region. RESULTS: Diversity is very low for all genes in the X-linked region in the wild dioecious population, with nucleotide diversity π syn = 0.00017, tenfold lower than the autosomal region (π syn = 0.0017) and 12-fold lower than the Y-linked region (π syn = 0.0021). Analysis of the X-linked sequences shows an undivided population, suggesting a geographically wide diversity-reducing event, whereas two subpopulations were observed in the autosomes separating gynodioecy and dioecy and three subpopulations in the Y-linked region separating three male populations. The extremely low diversity in the papaya X-linked region was probably caused by a recent, strong selective sweep before domestication, involving either the spread of a recessive mutation in an X-linked gene that is beneficial to males or a partially dominant mutation that benefitted females or both sexes. Nucleotide diversity in the domesticated X samples is about half that in the wild Xs, probably due to the bottleneck when hermaphrodites were selected during domestication. CONCLUSIONS: The extreme low nucleotide diversity in the papaya X-linked region is much greater than observed in humans, great apes, and the neo-X chromosome of Drosophila miranda, which show the expected pattern of Y-linked genes < X-linked genes < autosomal genes; papaya shows an unprecedented pattern of X-linked genes < autosomal genes < Y-linked genes.

The sex chromosomes of Silene latifolia revisited and revised.

Classical studies have established that, during meiosis, the X and Y chromosomes of the model dioecious plant Silene latifolia pair over a region at the ends of their q arms. We used fluorescence in situ hybridization of two molecular markers to demonstrate that this widely accepted model is incorrect. From these data we conclude that the homologous arm of the X chromosome is the p arm and that of the Y chromosome is the q arm. The establishment of the proper orientation of the pseudoautosomal region is essential for mapping and evolutionary studies.

Genetic and functional analysis of DD44, a sex-linked gene from the dioecious plant Silene latifolia, provides clues to early events in sex chromosome evolution.

Silene latifolia is a dioecious plant with heteromorphic sex chromosomes. The sex chromosomes of S. latifolia provide an opportunity to study the early events in sex chromosome evolution because of their relatively recent emergence. In this article, we present the genetic and physical mapping, expression analysis, and molecular evolutionary analysis of a sex-linked gene from S. latifolia, DD44 (Differential Display 44). DD44 is homologous to the oligomycin sensitivity-conferring protein, an essential component of the mitochondrial ATP synthase, and is ubiquitously expressed in both sexes. We have been able to genetically map DD44 to a region of the Y chromosome that is genetically linked to the carpel-suppressing locus. Although we have physically mapped DD44 to the distal end of the long arm of the X chromosome using fluorescence in situ hybridization (FISH), DD44 maps to the opposite arm of the Y chromosome as determined by our genetic map. These data suggest that chromosomal rearrangements have occurred on the Y chromosome, which may have contributed to the genetic isolation of the Y chromosome. We discuss the implications of these results with respect to the structural and functional evolution of the S. latifolia Y chromosome.

The genomic organization of Ty3/gypsy-like retrotransposons in Helianthus (Asteraceae) homoploid hybrid species.

The origin of new diploid, or homoploid, hybrid species is associated with rapid genomic restructuring in the hybrid neospecies. This mode of speciation has been best characterized in wild sunflower species in the genus Helianthus, where three homoploid hybrid species (H. anomalus, H. deserticola, and H. paradoxus) have independently arisen via ancient hybridization events between the same two parental species (H. annuus and H. petiolaris). Most previous work examining genomic restructuring in these sunflower hybrid species has focused on chromosomal rearrangements. However, the origin of all three homoploid hybrid sunflower species also is associated with massive proliferation events of Ty3/gypsy-like retrotransposons in the hybrid species' genomes. We compared the genomic organization of these elements in the parent species and two of the homoploid hybrid species using fluorescence in situ hybridization (FISH). We found a significant expansion of Ty3/gypsy-like retrotransposons confined to the pericentromeric regions of two hybrid sunflower species, H. deserticola and H. paradoxus. In contrast, we detected no significant increase in the frequency or extent of dispersed retrotransposon populations in the hybrid species within the resolution limits of our assay. We discuss the potential role that transposable element proliferation and localization plays in the evolution of homoploid hybrid species.

Local patterns of nucleotide polymorphism are highly variable in the selfing species Arabidopsis thaliana.

Neighboring genes predictably share similar evolutionary histories to an extent delineated by recombination. This correlation should extend across multiple linked genes in a selfing species such as Arabidopsis thaliana due to its low effective recombination rate. To test this prediction, we performed a molecular population genetics analysis of nucleotide polymorphism and divergence in chromosomal regions surrounding four low-diversity loci. Three of these loci, At1g67140, At3g03700, and TERMINAL FLOWER1 (TFL1), have been previously implicated as targets of selection and we would predict stronger correlations in polymorphism between neighboring loci due to genetic hitchhiking around these loci. The remaining locus, At1g04300, was identified in a study of linkage disequilibrium surrounding the CRYPTOCHROME2 (CRY2) locus. Although we found broad valleys of reduced nucleotide variation around two of our focal genes, At1g67140 and At3g03700, all chromosomal regions exhibited extreme variation in the patterns of polymorphism and evolution between neighboring loci. Although three of our four regions contained potential targets of selection, application of the composite-likelihood-ratio test of selection in conjunction with a goodness-of-fit test supports the selection hypothesis only for the region containing At3g03700. The degree of discordance in evolutionary histories between linked loci within each region generally correlated with estimates of recombination and linkage disequilibrium for that region, with the exception of the region containing At1g04300. We discuss the implications of these data for future population genetics analyses and genomics studies in A. thaliana.