BACKGROUND
phytochromes are specialized photoreceptors that perceive and interpret light signals from the environment to regulate virtually all aspects of plant development, including seed germination, chloroplast development, tropisms, shade avoidance responses, floral initiation, circadian rhythms, pigmentation, and senescence  <cit> .
the phytochromes have a primary role in sensing red  and far-red  light, and also play a role in the perception of blue  and ultraviolet  light  <cit> .
the active phytochrome molecule consists of a large  apoprotein bound to a phycobilin chromophore  <cit> .
the phytochrome apoproteins are encoded by a small gene family in all plant taxonomic divisions, including parasitic plants, mosses, cryptogams, and green algae  <cit> .
in angiosperms, the phytochrome apoprotein genes have been classified into four or five gene sub-families based on sequence similarity to the five phytochrome genes of arabidopsis: phya, phyb, phyc, phyd, and phye  <cit> .
all five arabidopsis phytochromes share an amino acid sequence similarity of 46-56%, with the exception phyb and phyd--which are the result of recent gene duplication and share ~80% amino acid identity  <cit> .
thus, the five arabidopsis genes are often assigned to four subfamilies: phya, phyb/d, phyc, and phye  <cit> .
the arabidopsis phyb/d subfamily is more closely related to phye gene  than to the phya and phyc genes , which together form a separate ancient evolutionary clade  <cit> .
having presumably arisen by gene duplication and subsequent subfunctionalization and/or neofunctionalization, the phytochrome gene family in toto performs a complex network of redundant, partially redundant, non-overlapping, and in some cases antagonistic regulatory functions throughout plant development  <cit> .
for example, all arabidopsis phytochromes play diverse and interacting roles in photoperiodic regulation of floral initiation.
phya, phyb, phyd and phye act partially redundantly in the light-dependent entrainment of the circadian clock  <cit> , which in turn regulates transcription of the floral inducer constans  in a circadian manner  <cit> .
in arabidopsis, phya, in conjunction with blue-light dependent cryptochrome photoreceptors cry <dig> and cry <dig>  promotes flowering by inhibiting the degradation of co protein, while phyb acts antagonistically to stimulate co degradation  <cit> .
in addition, phyb, phyd and phye act partially redundantly as repressors of flowering that are dependent on r/fr ratio  <cit> .
in this role, phyb also acts downstream of co as a negative regulator of transcription of the 'florigen' molecule ft  in a tissue specific manner  <cit> .
mutant analyses indicate that phyc also plays a role in photoperiodic flowering  <cit> .
further, genetic variation at the phyc locus underlies some of the natural phenotypic variation in flowering time in arabidopsis  <cit> .
in angiosperms, the composition of phytochrome gene family varies significantly among taxonomic lineages.
although a single phya gene is present in most flowering plants, some plant families, such as carnation  and legumes , have two distinct phya genes  <cit> .
similarly, several plant lineages have gained multiple phyb-like genes through independent gene duplications of phyb  <cit> .
for example, tomato has two phyb genes  that are not directly orthologous to arabidopsis phyb and phyd, respectively  <cit> .
while most angiosperms have a single phyc gene, species in some families such as fabaceae and salicaceae appear to have lost phyc during evolution  <cit> .
although a single phye-like gene is present in most flowering plants, phye is completely absent in poplar , in the piperales, and some monocots such as maize  <cit> .
finally, the novel phyf subfamily, which groups with phya/c clade, has been identified in tomato  <cit> .
little is known about the composition of the phytochrome gene family in cultivated cottons or their wild relatives  in the malvaceae family.
this is despite the fact that physiological experiments suggest that phytochromes regulate economically important aspects of cotton development, including drought resistance, seed germination, plant architecture, photoperiodic flowering, and fiber elongation  <cit> .
for example, r/fr photon ratio influences the length and diameter of developing seed fiber; fibers exposed to a high r/fr photon ratio during development were longer than those that received lower r/fr ratio, implicating the involvement of a phytochrome  <cit> .
while modern domesticated varieties of the major cultivated cottons g. hirsutum l. and g. barbadense l. exhibit photoperiod independent flowering, wild and 'primitive' accessions of g. hirsutum and g. barbadense flower under short-day photoperiodic control  <cit> .
an understanding of the molecular-genetic basis of differences in photoperiodic flowering in cottons will accelerate strategies for improvement of cultivated varieties through the introgression of valuable genetic traits from wild germplasm  <cit> .
in this regard, it is important to note that mutational changes in phytochrome function have been implicated in the loss of photoperiod sensitivity in several major crops including sorghum, barley, rice, and soy  <cit> .
a thorough characterization of the phytochrome gene family in cotton species is necessary for understanding the molecular basis of photoperiodic flowering, the influences of light quality on cotton fiber elongation, and other aspects of cotton development.
any inventory of phytochrome genes of cottons is complicated by the fact that the major cultivated species, g. hirsutum and g. barbadense are allotetraploids.
diploid species in the genus gossypium are categorized into eight genome groups  based on cytogenetic and phylogenetic criteria  <cit> .
the old-world a genome group and the new world d genome group diverged from each other on the order of 1- <dig> mya  <cit> , then underwent hybridization and polyploidization creating an ad allopolyploid lineage ancestral to g. hirsutum  and g. barbadense  on the order of  <dig> mya  <cit> .
in this study, we utilized a pcr-based approach with low-degeneracy primers to obtain gene fragments, or 'genome sequence tags'  that yield an initial description of the composition and evolution of the phytochrome gene family in the new world allotetraploid cottons gossypium hirsutum and g. barbadense, and in the old-world diploids g herbaceum l. and g. raimondii ulbr., which are considered to be extant relatives of the a- and d-genome diploid ancestors  of the allotetraploid lineage.
this study provides a necessary foundation for studies of the specific biological functions of each of the phytochrome genes in cotton species, and helps to illuminate the evolutionary patterns of duplicated genes in complex genomes, as well as the evolutionary history of the world's most important fiber crop species.
RESULTS
because our results were derived from pcr, our inventory of the phytochrome gene family in gossypium spp.
is provisional.
all sequences have been submitted to genbank .
phytochrome hinge amplification using 'universal' primers
between n-terminal 'photoperception domain' and c-terminal 'signaling domain' of the phytochrome apoprotein is a short 'hinge region'  that shows relatively high sequence variation, and has proven useful for characterization of the phytochrome gene complement in a variety of plant species, and for robust phylogenetic analyses  <cit> .
to amplify the hinge region of all cotton phytochromes, we used an alignment of eudicot phytochrome sequences to design a 768-fold degenerate pcr primer  based on the conserved hypatdip peptide in the n-terminal domain, and a  <dig> -fold degenerate pcr primer , based on the conserved pfplryac peptide in the c-terminal domain .
i corresponds to inosine.
r, y, m, k, s, w correspond to the iupac-iub ambiguity set.
amplification across the hinge region using taq dna polymerase yielded pcr products from all taxa.
we cloned the amplification products from each taxon into an e. coli vector, then sequenced ~ <dig> clones for each taxon.
for all taxa, a majority  of clones showed the highest similarity in blast searches to arabidopsis phye .
for each taxon, only a minority of clones showed high-scoring similarity to arabidopsis phya or phyb.
this apparently skewed distribution of amplification products -- observed across all taxa -- suggested an amplification bias in favor of phye amplicons.
no clones were obtained from any taxon that had high-scoring similarity to arabidopsis phyc or phyd.
no new phytochrome sub-families were observed.
amplification of the phya gene sub-family
because of possible biased amplification, we designed new less-degenerate hinge-region primer sets for the phya, phyb/d, and phyc sub-families  using available phytochrome sequences from species in the rosid clade, which includes both cotton and arabidopsis  <cit> .
the hinge regions of phya genes were amplified using phyabnondeg-f and phyadeg-r , yielding a ~ <dig> bp amplification product from all accessions.
in blast database searches, all clones had a high-scoring pair relationship with arabidopsis phya .
sequences from a total of more than  <dig> clones across all taxa yielded two distinct consensus contigs from each of the diploids g. herbaceum and g. raimondii, and four distinct contigs from the allotetraploids g. barbadense and g. hirsutum.
when aligned across all taxa, these contigs yielded a  <dig> bp consensus alignment that had an average pairwise sequence similarity of  <dig> %, with  <dig> sites  identical across all taxa, and no stop codons or indels in any taxa.
distance analysis  showed two well-separated gene sub-clades .
these sub-clades were designated phya <dig> and phya <dig>  the level of hinge-region differentiation between these two sub-clades was far greater than that seen in other cotton phytochrome gene sub-families , with an uncorrected "p" distance of  <dig> , corresponding to  <dig> nt changes  based on parsimony.
these data indicated that a single phya gene underwent duplication after the divergence of the cotton and arabidopsis lineages, but prior to the divergence of a-genome and d-genome lineages, leaving each of the modern diploids in our study  with a complement of two phya paralogs .
indeed, four distinct contigs were observed in both the inbred g. hirsutum cultivar tm- <dig> and in the doubled-haploid line g. barbadense 3- <dig>  for each allotetraploid taxon, two contigs fell into each of the phya- <dig> and phya- <dig> clades .
a conservative inventory of available est sequences indicated that at least two distinct phya loci are expressed in g. hirsutum .
within each of the phya <dig> and phya <dig> clades, the level of nucleotide diversity was very low, with at most four parsimonious nucleotide changes separating each contig.
however, within the phya <dig> clade, the contigs resolved into two subclades  that each contained a single contig from one of the diploid taxa and one contig from each of the allotetraploids.
for example, g. raimondii  phya <dig> grouped in a single contig from each of g. hirsutum and g. barbadense.
based on this grouping, the latter contigs were assigned the provisional designation of phya <dig> d. similarly, g. herbaceum  grouped with g. hirsutum phya <dig> a and g. barbadense phya <dig> a. based on similar criteria, the phya <dig> clade was also divided into phya <dig> a and phya <dig> d subclades .
the phylogenetic resolution of a- and d-genome subclades supported the hypothesis that each of the a- and d-genome diploids contributed both phya <dig> and phya <dig> to the allotetraploid lineage.
thus, although hinge-region nucleotide diversity within each of the phya <dig> and phya <dig> clades was low, it was sufficient to resolve a tentative phya gene complement for each taxon, as well as the pattern of gene inheritance through the allopolyploidization event.
amplification of the phyb/d gene sub-family
a ~ <dig> bp fragment from the phyb/d hinge region was obtained by amplification using primers phyabnondeg-f and phybdeg-r .
sequences from a total of  <dig> clones yielded a single consensus contig from each of the diploid cottons g. herbaceum and g. raimondii, and from the allotetraploid g. hirsutum.
two distinct contigs were assembled from clones derived from the allotetraploid g. barbadense.
these clone sequences shared ~85% nucleotide identity with the arabidopsis phyb gene and ~78% nt identity with arabidopsis phyd.
all clones had a high-scoring pair relationship with the arabidopsis phyb gene  as well as significant similarity to the arabidopsis phyd gene .
consensus sequences were aligned across all taxa, yielding a  <dig> bp alignment with an average pairwise sequence similarity of  <dig> %, with  <dig> sites  identical across all taxa, no stop codons and no indels.
although these data indicated the presence of at least one phyb gene in each of the a- and d-genome diploid plants and in g. hirsutum, and at least two genes phyb genes in the g. barbadense, the low level of nucleotide differentiation observed within the hinge region yielded insufficient phylogenetic information to characterize the phyb gene complement in any of the study taxa.
to obtain better resolution of the phyb gene complement, additional low degeneracy primers 1010-f, 1910-f, 1910-r, and 2848-r  were used along with primer phyabnondeg-f to create a  <dig>  kb long series of overlapping amplicons corresponding to approximately  <dig>  kb of the arabidopsis phyb gene and extending from the hinge, through the first intron and into the second exon .
after amplification, cloning and sequencing, the amplicons were assembled for each taxon.
in all gossypium taxa examined, the first intron was ~ <dig> bp longer than the first intron of phyb from arabidopsis.
unlike the other phytochrome amplicons, we detected a high frequency of putative pcr-mediated recombination events  <cit>  within the phyb <dig>  kb fragment from amplifications using g. barbadense as template.
the recombination detection algorithm rdp <dig>  <cit>  identified a number of clones resulting from apparent recombination between the a-genome and d-genome derived homeologous sequences, with predicted breakpoints  between nucleotides  <dig> and  <dig> of the alignment.
after omission of these recombinant clones, composite amplicon sequences from each taxon were aligned, creating a consensus alignment of  <dig>  bp with  <dig> % average pairwise similarity and  <dig>  identical sites .
overall, the cotton phyb genes shared 65% nucleotide identity with the arabidopsis phyb ortholog.
no stop codons or indels were detected in exon sequences.
a  <dig> bp putative deletion was observed in one contig  from g. hirsutum.
in addition, a  <dig> bp indel was polymorphic between the phyb.a and phyb.d clades.
finally, phyb of g. raimondii had an additional  <dig> bp insertion.
all indel polymorphisms were located within first introns.
detailed phylogenetic analyses of the  <dig>  bp contigs from a-, d-, and ad-genome cottons  indicated the presence of least one phyb locus in the two diploid cottons, g. herbaceum and g. raimondii, and at least two phyb loci in both allotetraploid cottons.
the g. hirsutum and g. barbadense sequence contigs each grouped into two sub-clades .
the single phyb contig from g. herbaceum was used to define the phyb.a cluster , while the single phyb contig from g. raimondii anchored the phyb.d cluster.
from these results, we concluded that phyb.a and phyb.d, which shared ~98% nucleotide sequence identity, arose as orthologs at the time of divergence of the a- and d-genome diploid lineages.
we surmised that phyb.a was contributed to the allotetraploids via the a-genome ancestor and phyb.d was contributed via the d-genome ancestor.
available est sequences indicated that at least one phyb locus is expressed in g. hirsutum .
amplification from the phyc gene sub-family
several sets of degenerate primer pairs that were designed on the basis of the conserved hypatdip and pfplryac regions -- including several designed from rosid phyc nucleotide sequences -- failed to produce detectable pcr amplification products from the gossypium species tested .
however, the identification of a small est clone  with similarity to arabidopsis phyc  in a library from g. raimondii floral tissue  <cit> , allowed us to design the primer phyc_1r_dfci within the c-terminal domain .
when used in combination with phydeg-f, this primer amplified a ~ <dig> kb fragment composed entirely of coding sequence from the first exon of phyc, including the hinge .
all clones obtained using this primer pair had a high-scoring similarity to arabidopsis phyc .
from these clones, we assembled a single consensus contig from each of the diploid species g. herbaceum and g. raimondii, and two distinct consensus contigs from each of the allotetraploids g. hirsutum and g. barbadense.
consensus sequences for each of the putative phyc contigs were aligned across all taxa, yielding a  <dig>  bp alignment with an average pairwise sequence similarity of  <dig> %,  <dig>  sites  identical across all taxa, with no indels or stop codons in any taxa.
in phylogenetic analyses , the phyc consensus sequences grouped into two major clades .
one of these clades contained the g. herbaceum contig and one contig from each of g. hirsutum and g. barbadense.
this clade was designated phyc.a.
the other clade, designated phyc.d, included the g. raimondii contig along with the other of the two contigs from each of g. hirsutum and g. barbadense.
these data indicated that both the a- and d-genome ancestors had one phyc gene, and that upon hybridization and polyploidization, this gene was contributed from each diploid to the allotetraploid ancestor of g. hirsutum and g. barbadense.
for comparison with the other phytochromes, we also analyzed a portion of the phyc alignment corresponding to the hinge region only.
this alignment was  <dig> nucleotide pairs in length, with pairwise sequence similarity of  <dig> %,  <dig> sites  identical across all taxa, with no indels.
although it encompassed fewer variable nucleotides, nj analysis of the hinge region alone could be used to differentiate the phyc.a and phyc.d clades  and to infer the composition and evolutionary inheritance of the phyc gene family in cottons .
our failure to obtain phyc hinge amplification with several sets of both universal  and rosid specific primers was entirely due to substantial nucleotide differentiation in phyc, particularly within the hinge region.
for example, the  <dig> nt long phydeg-r primer had six nucleotide mismatches with the cotton phyc genes, including three transitions and three transversions.
five of the six mismatches occurred at what are considered to be invariant  nucleotide positions.
it should be noted that these divergent nucleotides in the conserved primer-binding site did not alter the amino acid sequence .
the phye gene sub-family
phye hinge region consensus contigs from our study taxa formed a  <dig> bp alignment with an average pairwise similarity of  <dig> %, with  <dig>  invariant sites, no indels, and no stop codons in any taxa.
the consensus of the aligned phye sequences had 80% nucleotide similarity to the corresponding fragment of the arabidopsis phye gene.
based on maximum parsimony, nucleotide diversity in the cotton phye hinge sequences could be explained by a minimum of six nucleotide changes, all of which were synonymous.
nj analysis of the cotton phye hinge region showed two distinct clades  corresponding to the a- and d-genome derived orthologs , a finding consistent with a hypothesis in which each diploid ancestor contributed a single phye ortholog to the allotetraploid lineage .
interestingly, while two distinct phye contigs were obtained from g. hirsutum, only a single contig, which grouped with the d-genome clade, was obtained from g. barbadense.
available est sequences indicated that at least one phye locus is expressed in g. hirsutum .
a global hinge-based alignment of arabidopsis and cotton phytochromes
phya, phyb, phyc and phye hinge regions from arabidopsis and gossypium spp.
were aligned to create a global phytochrome alignment  <dig> nucleotides in length, with an average pairwise similarity of  <dig> % and  <dig> identical sites .
the gene phylogeny generated from this alignment  reflected divergence of phya, phyb, phyc and phye as a result of speciation  and gene duplication .
the level of nucleotide divergence of each of the gene sub-families after nodes 1a, 1b, 1c and 1e  was similar, with a mean of  <dig>  ±  <dig>  nucleotide substitutions per site.
however, the synonymous  and non-synonymous  substitution rates were both significantly more variable among the various gene sub-families defined by nodes 1a, 1b, 1c and 1d than were simple nucleotide distances .
despite this variation, all sub-families showed a ka/ks ratio < <dig> , implying that each remains under purifying selection for function.
further, excessively long branch-lengths, which are often found in pseudogenes, were not observed.
in the phyb, phyc and phye clades, the branch lengths leading to the arabidopsis orthologs, which have known biological functions, were longer than the branches leading to their respective cotton orthologs.
considered together, these lines of evidence indicate that each of the phytochrome sub-families retains some biological function in gossypium, as they do in arabidopsis  <cit> .
further, our topology supports the conclusion that phyd is the result of a relatively recent gene duplication that may be exclusive to the brassicaceae family  <cit> .
nodes refer to the nj tree in figure  <dig>  k-2p indicates the mean kimura 2-parameter distances between arabidopsis and cotton gene sequences.
discussion
resolution of the phytochrome gene family
in three out of four cases, we were able to successfully resolve the inventory and evolutionary relationships of the phytochrome genes in diploid and allotetraploid cottons using the hinge region only.
this finding supports the general utility of employing the hinge region for identifying gsts for phytochromes.
in only one case  was additional gene sequence required for sufficient phylogenetic resolution.
in another case , nucleotide divergence at a commonly used primer-binding site prevented the characterization of the hinge region by the typical strategy of using primers based on conserved flanking peptides hypatdip and pfplryac.
however, nucleotide diversity within the phyc hinge region itself was sufficiently informative to resolve the pattern of evolutionary inheritance through allotetraploidization event.
the sequencing of phytochrome gene fragments from a- and d-genome diploids, as well as from ad allotetraploid taxa, provides an essential foundation for all subsequent analysis of phytochrome function and evolution in gossypium.
the sequenced fragments provide sufficient information  to unequivocally identify or 'tag' various orthologs, homeologs and paralogs, as well as monitor their patterns of nucleotide divergence, and trace their evolutionary inheritance through the allopolyploidization event.
this information will serve as a foundation for further sequence assembly and annotation, and will be used to design locus-specific primer sets for quantitative rt-pcr assays that will measure transcript levels for each gene family member.
in some cases  levels of sequence divergence are high enough to support studies of gene function using rnai or amirna approaches to create gene-specific knockouts  <cit> .
the use of well characterized 'candidate genes' of agronomic interest is becoming an integral component of marker-assisted selection efforts in plants  <cit> .
several snp-based molecular markers  <cit>  are now being developed using the diagnostic nucleotide characters identified in this study, and are being mapped in experimental cotton populations that show segregation of phytochrome-controlled traits such as fiber length and flowering time.
the ancestral phytochrome gene complement of the malvales and brassicales
our study indicated that the diploid ancestors to the world's major fiber crops  had a complement of phytochrome apoprotein genes that was very similar to that of the model plant arabidopsis thaliana.
this was not entirely unexpected given the relatively close phylogenetic relationship of the two lineages  <cit> .
the most-simple evolutionary scenario is that the last common ancestor of arabidopsis and cotton, possibly an arborescent species in the late cretaceous period  <cit> , had a phytochrome gene complement consisting of one functional gene in each of the phya, phyb/d, phyc and phye subfamilies.
phya duplication in gossypium
after the divergence of the malvales and brassicales, the ancestral phya gene underwent duplication resulting in the observed phya- <dig> and phya- <dig> paralogs of modern gossypium spp.
as the a- and d-genome diploids have both paralogs, the duplication event occurred prior to the divergence of the a- and d-genome lineages.
using  <dig> mya  as a rough estimate of the time of divergence of the malvales and brassicales  <cit> , along with our observed ks of  <dig>  in the phya hinge region in this time interval, we can derive a crude estimate of  <dig>  substitutions/synonymous-site/million years, and an estimate of the time of phya duplication of ~ <dig> mya.
this estimate places the duplication well within the crown group of malvales and the malvaceae family  <cit> .
given our time estimate, the phya duplication may be exclusive to the genus gossypium, but would have occurred prior to the estimated time of divergence of the a and d genome groups  <cit> .
as neither we nor others  <cit>  have observed evidence of additional nuclear gene duplications or chromosomal duplications in this time period, the phya event was likely a tandem or segmental duplication, rather than whole genome duplication.
after a gene duplication event, one of the two newly duplicated genes is theoretically unconstrained by selection for function, and is thus free to accumulate mutations leading to a pseudogene fate, subfunctionalization, or neofunctionalization  <cit> .
although we did not obtain definitive evidence of pseudogenic sequences in any of the phytochromes or taxa studied , we did observe significant variation in ka/ks ratios in pairwise interspecific comparisons , leaving open the possibility of pseudogene outcomes.
alternatively, one of the duplicated genes may undergo positive selection to gain a novel function .
further, duplicated gene-pairs may subdivide the function of ancestral gene .
perhaps the most intriguing fate, which has been observed empirically, but not yet explained in theory, is the situation in which both gene copies may be retained for a lengthy period under what appears to be purifying or negative selection  <cit> .
one approach to understanding the evolutionary fates of duplicated genes is through an analysis of the signature of natural selection on amino acid encoding sequences.
although the hinge regions of phytochromes display relatively high levels of nucleotide diversity  <cit> , they do not evolve under neutrality.
the hinge region participates in inter-domain communication in phytochrome molecules  <cit> .
for example, phosphorylation of a serine residue in the phya hinge plays a likely role in regulating protein-protein interactions between phytochrome and downstream signal-transducing molecules  <cit> .
compared to wild-type, a mutation in the hinge region of arabidopsis phyb is deficient in localization into distinct nuclear bodies  <cit> .
further, a single nucleotide polymorphism  in the hinge of one of two phyb genes in aspen  was associated with natural geographic variation in the timing of bud-set  <cit> .
in comparisons between cotton and arabidopsis , the ka/ks ratio for the phya hinge region was  <dig>  -- a value that is typical for genes under purifying selection  <cit> .
in contrast, the ka/ks ratio for phya after gene duplication  was  <dig> , or ~ <dig> -fold higher.
this value is also ~ <dig> -fold greater than the mean ka/ks ratio of all phytochrome hinge regions  of approximately  <dig>  ±  <dig> .
this significantly elevated ka/ks ratio after the phya duplication could be attributed to a relaxation of stabilizing selection and/or subfunctionalization of the nascent phya paralogs .
the possible functional divergence of phya <dig> and phya <dig> may be more pronounced after the separation of the a- and d-genome lineages .
a comparison of phya <dig> in the two diploids yields a ka/ks ratio of ~ <dig> , primarily due to amino acid substitutions in phya <dig> d, while phya <dig> has a ka/ks ratio of  <dig>  in the same taxonomic comparisons.
although this difference is suggestive of possible differential rates of functional evolution in the paralogs, it is not statistically significant in fisher's exact test .
it will be of interest to determine whether the cotton phya paralogs have distinct functions.
experiments are underway to determine the respective biological functions of each phya- <dig> and phya- <dig> in g. hirsutum and g. barbadense using paralog-specific rt-pcr, rnai gene knockout, and tests for genetic associations between phytochrome-controlled phenotypic traits and phya- <dig> and phya- <dig> specific molecular markers.
a 'candidate gene' approach has recently been used in soy  to uncover a genetic linkage between the photoperiod insensitivity locus e <dig> and one of the two the phya genes, designated gmphya <dig> and gmphya <dig>  <cit> .
loss of photoperiodic flowering is associated with a ty1/copia-like retrotransposon insertion into exon  <dig> of gmphya <dig>  the authors argue that gene duplication and partial redundancy of the phya genes may have facilitated the loss of photoperiod sensitivity by allowing the gmphya <dig>  mutant to avoid the major deleterious phenotypic effects that would have been caused by complete deficiency of phya gene function.
d-d indicates means of comparisons within extant diploids, d-t indicates means of comparisons of extant diploids with tetraploids, t-t indicates means of comparisons within tetraploids.
s dif, synonymous differences; s pos, synonymous positions; ns dif, non-synonymous differences; ns pos, non-synonymous positions.
p indicates significance as determined by fisher's exact test.
persistence and loss of phytochrome paralogs after allopolyploidization
all phytochromes underwent gene duplication by polyploidization at the time of formation of the ad allotetraploids, on the order of  <dig> - <dig>  mya  <cit> .
for example, in g. hirsutum, we detected a minimum set of ten distinct phytochrome genes, including four phya genes.
in order to assess the evolutionary trajectory of these recently duplicated genes, we examined the synonymous and non-synonymous divergence rates of a- and d-genome phytochrome orthologs and homeologs  in pairwise comparisons of 1) diploids with diploids , 2) diploids with tetraploids , and 3) tetraploids with tetraploids .
given that the allotetraploid cottons had both a- and d-genome derived copies of each gene on the order of hundreds of thousands of years, we hypothesized that there may be a relaxation of selection in the allotetraploids, as one of the two copies should no longer be evolutionarily constrained.
however, in comparisons of a- vs. d-genome derived orthologs or homeologs for six gsts , we did not observe dramatic differences in ka/ks between diploid and allotetraploids in any gst except the hinge region of phya <dig> .
because of low levels of nucleotide divergence, we employed fisher's exact test  <cit>  and found no significant differences in the patterns of nucleotide evolution in allotetraploids vs. diploids.
thus, there was no broad evidence of relaxation of natural selection on gene function after gene duplication by allotetraploidization.
further, the generally low ka/ks ratios across all genes and taxa support a model in which that the phytochrome homeologs are largely evolving independently by a birth-and-death model rather than concerted evolution  <cit> .
the coding sequences of the phyb  <dig>  kb fragment also appeared be evolving under stabilizing selection in both the diploids  and allotetraploids  reflecting continued selective constraint on coding sequence evolution after polyploidization.
however, there was a significant excess of non-synonymous substitutions in both diploids and allotetraploids  indicating a partial relaxation of negative selection and/or functional divergence of the phyb homeologs.
in the allotetraploid cottons, both phyc.a and phyc.d are also evolving in a pattern consistent with purifying selection .
however, it should be noted that the phyc.d clade appears to be evolving at distinctly faster rate  than the phyc.a clade .
this suggests either a relaxation of purifying selection in, or functional divergence of phyc.d.
in a similar study of phytochromes in cultivated sorghum  and its wild congeneric relatives  <cit> , phyc was undergoing faster amino acid evolution than phya or phyb.
in the both the phyb and phyc gene subfamilies of cotton, the sequences of the c-terminal signaling domain had higher ka/ks ratios than the corresponding hinge region alone.
this may reflect the co-evolution of protein-protein interactions with downstream signaling partners, which are mediated by the c-terminal 'signal transduction' domain  <cit> .
while phye-related contigs had low ka/ks values , indicating purifying selection, no contig corresponding to an expected g. barbadense phye.a ortholog was observed.
this may have been due to under-sampling of g. barbadense clones for sequencing, or due to nucleotide divergence in primer sites .
of the  <dig> phye-like clone sequences obtained from g. barbadense, all were in the d-genome derived clade, which would be an unlikely result  assuming equal amplification efficiencies for phye.a and phye.d.
alternatively, the apparent lack of a phye.a ortholog in g. barbadense could be explained by concerted evolution, gene conversion, or by pcr-mediated recombination  <cit> .
overall, the phye genes, like the other cotton phytochromes, had more synonymous than non-synonymous nucleotide substitutions, favoring a birth-and-death model of gene evolution.
CONCLUSIONS
our preliminary efforts to obtain an inventory of the cotton phytochrome gene family  indicated that diploid a- and d-genome diploid cottons have two paralogous phya genes , and one each of phyb, phyc, and phye gene sub-families.
coding sequence evolution in phya <dig> was significantly elevated, suggesting loss of selection for function, or incipient subfunctionalization.
other than this duplication and the lack of a separate phyd gene, the phytochrome complement of diploid cottons was very similar to that observed in the closely related model plant arabidopsis thaliana, which greatly facilitates cross-species comparisons.
whole genome duplication via allopolyploidization  resulted in additive amalgamation of phytochrome genes within a single nucleus in the allotetraploid, retaining complete gene complements of at least four phya genes, two genes of each phyb, phyc and phye in ad-genome g. hirsutum.
g. barbadense may lack the phye gene contributed by the a-genome ancestor.
strong purifying selection on nearly all of the phytochrome genes suggests some level of conservation of function of each of the genes after polyploidization.
with the possible exception of one of the phye.a homeologs in g. barbadense, we did not see evidence of gene loss.
we did not observe any convincing evidence of concerted evolution by gene conversion.
rather, the genes duplicated by allopolyploidy appear to be largely retained, and evolving independently as observed in  <dig> other nuclear genes in allotetraploid cottons  <cit> .
these results further our understanding of the evolutionary fates of duplicate genes following allopolyploidization.
information on key evolutionary events , as well as rates and patterns of evolutionary change, are an important component of the functional annotation of genes and genomes  <cit> .
these data provide the foundation for more comprehensive studies of the biological functions of each of the cotton paralogs and homeologs.
the development of phytochrome 'candidate gene' markers based on the gsts identified here may prove useful in the mobilization of valuable genes from photoperiodic wild and primitive cottons into elite cotton varieties, in order to improve stress tolerance, disease resistance, fiber quality, and other traits.

 phytochromes are specialized photoreceptors that perceive and interpret light signals from the environment to regulate virtually all aspects of plant development, including seed germination, chloroplast development, tropisms, shade avoidance responses, floral initiation, circadian rhythms, pigmentation, and senescence. the phytochromes have a primary role in sensing red and far-red light, and also play a role in the perception of blue and ultraviolet light. the active phytochrome molecule consists of a large apoprotein bound to a phycobilin chromophore. the phytochrome apoproteins are encoded by a small gene family in all plant taxonomic divisions, including parasitic plants, mosses, cryptogams, and green algae. in angiosperms, the phytochrome apoprotein genes have been classified into four or five gene sub-families based on sequence similarity to the five phytochrome genes of arabidopsis: phya, phyb, phyc, phyd, and phye. all five arabidopsis phytochromes share an amino acid sequence similarity of 46-56%, with the exception phyb and phyd--which are the result of recent gene duplication and share ~80% amino acid identity. thus, the five arabidopsis genes are often assigned to four subfamilies: phya, phyb/d, phyc, and phye. the arabidopsis phyb/d subfamily is more closely related to phye gene than to the phya and phyc genes, which together form a separate ancient evolutionary clade. having presumably arisen by gene duplication and subsequent subfunctionalization and/or neofunctionalization, the phytochrome gene family in toto performs a complex network of redundant, partially redundant, non-overlapping, and in some cases antagonistic regulatory functions throughout plant development. for example, all arabidopsis phytochromes play diverse and interacting roles in photoperiodic regulation of floral initiation. phya, phyb, phyd and phye act partially redundantly in the light-dependent entrainment of the circadian clock, which in turn regulates transcription of the floral inducer constans in a circadian manner. in arabidopsis, phya, in conjunction with blue-light dependent cryptochrome photoreceptors cry and cry promotes flowering by inhibiting the degradation of co protein, while phyb acts antagonistically to stimulate co degradation. in addition, phyb, phyd and phye act partially redundantly as repressors of flowering that are dependent on r/fr ratio. in this role, phyb also acts downstream of co as a negative regulator of transcription of the'florigen' molecule ft in a tissue specific manner. mutant analyses indicate that phyc also plays a role in photoperiodic flowering. further, genetic variation at the phyc locus underlies some of the natural phenotypic variation in flowering time in arabidopsis. in angiosperms, the composition of phytochrome gene family varies significantly among taxonomic lineages. although a single phya gene is present in most flowering plants, some plant families, such as carnation and legumes, have two distinct phya genes. similarly, several plant lineages have gained multiple phyb-like genes through independent gene duplications of phyb. for example, tomato has two phyb genes that are not directly orthologous to arabidopsis phyb and phyd, respectively. while most angiosperms have a single phyc gene, species in some families such as fabaceae and salicaceae appear to have lost phyc during evolution. although a single phye-like gene is present in most flowering plants, phye is completely absent in poplar, in the piperales, and some monocots such as maize. finally, the novel phyf subfamily, which groups with phya/c clade, has been identified in tomato. little is known about the composition of the phytochrome gene family in cultivated cottons or their wild relatives in the malvaceae family. this is despite the fact that physiological experiments suggest that phytochromes regulate economically important aspects of cotton development, including drought resistance, seed germination, plant architecture, photoperiodic flowering, and fiber elongation. for example, r/fr photon ratio influences the length and diameter of developing seed fiber; fibers exposed to a high r/fr photon ratio during development were longer than those that received lower r/fr ratio, implicating the involvement of a phytochrome. while modern domesticated varieties of the major cultivated cottons g. hirsutum l. and g. barbadense l. exhibit photoperiod independent flowering, wild and'primitive' accessions of g. hirsutum and g. barbadense flower under short-day photoperiodic control. an understanding of the molecular-genetic basis of differences in photoperiodic flowering in cottons will accelerate strategies for improvement of cultivated varieties through the introgression of valuable genetic traits from wild germplasm. in this regard, it is important to note that mutational changes in phytochrome function have been implicated in the loss of photoperiod sensitivity in several major crops including sorghum, barley, rice, and soy. a thorough characterization of the phytochrome gene family in cotton species is necessary for understanding the molecular basis of photoperiodic flowering, the influences of light quality on cotton fiber elongation, and other aspects of cotton development. any inventory of phytochrome genes of cottons is complicated by the fact that the major cultivated species, g. hirsutum and g. barbadense are allotetraploids. diploid species in the genus gossypium are categorized into eight genome groups based on cytogenetic and phylogenetic criteria. the old-world a genome group and the new world d genome group diverged from each other on the order of 1- mya, then underwent hybridization and polyploidization creating an ad allopolyploid lineage ancestral to g. hirsutum and g. barbadense on the order of mya. in this study, we utilized a pcr-based approach with low-degeneracy primers to obtain gene fragments, or'genome sequence tags' that yield an initial description of the composition and evolution of the phytochrome gene family in the new world allotetraploid cottons gossypium hirsutum and g. barbadense, and in the old-world diploids g herbaceum l. and g. raimondii ulbr., which are considered to be extant relatives of the a- and d-genome diploid ancestors of the allotetraploid lineage. this study provides a necessary foundation for studies of the specific biological functions of each of the phytochrome genes in cotton species, and helps to illuminate the evolutionary patterns of duplicated genes in complex genomes, as well as the evolutionary history of the world's most important fiber crop species.  because our results were derived from pcr, our inventory of the phytochrome gene family in gossypium spp. is provisional. all sequences have been submitted to genbank. phytochrome hinge amplification using'universal' primers between n-terminal'photoperception domain' and c-terminal'signaling domain' of the phytochrome apoprotein is a short'hinge region' that shows relatively high sequence variation, and has proven useful for characterization of the phytochrome gene complement in a variety of plant species, and for robust phylogenetic analyses. to amplify the hinge region of all cotton phytochromes, we used an alignment of eudicot phytochrome sequences to design a 768-fold degenerate pcr primer based on the conserved hypatdip peptide in the n-terminal domain, and a -fold degenerate pcr primer, based on the conserved pfplryac peptide in the c-terminal domain. i corresponds to inosine. r, y, m, k, s, w correspond to the iupac-iub ambiguity set. amplification across the hinge region using taq dna polymerase yielded pcr products from all taxa. we cloned the amplification products from each taxon into an e. coli vector, then sequenced ~ clones for each taxon. for all taxa, a majority of clones showed the highest similarity in blast searches to arabidopsis phye. for each taxon, only a minority of clones showed high-scoring similarity to arabidopsis phya or phyb. this apparently skewed distribution of amplification products -- observed across all taxa -- suggested an amplification bias in favor of phye amplicons. no clones were obtained from any taxon that had high-scoring similarity to arabidopsis phyc or phyd. no new phytochrome sub-families were observed. amplification of the phya gene sub-family because of possible biased amplification, we designed new less-degenerate hinge-region primer sets for the phya, phyb/d, and phyc sub-families using available phytochrome sequences from species in the rosid clade, which includes both cotton and arabidopsis. the hinge regions of phya genes were amplified using phyabnondeg-f and phyadeg-r, yielding a ~ bp amplification product from all accessions. in blast database searches, all clones had a high-scoring pair relationship with arabidopsis phya. sequences from a total of more than clones across all taxa yielded two distinct consensus contigs from each of the diploids g. herbaceum and g. raimondii, and four distinct contigs from the allotetraploids g. barbadense and g. hirsutum. when aligned across all taxa, these contigs yielded a bp consensus alignment that had an average pairwise sequence similarity of %, with sites identical across all taxa, and no stop codons or indels in any taxa. distance analysis showed two well-separated gene sub-clades. these sub-clades were designated phya and phya the level of hinge-region differentiation between these two sub-clades was far greater than that seen in other cotton phytochrome gene sub-families, with an uncorrected "p" distance of, corresponding to nt changes based on parsimony. these data indicated that a single phya gene underwent duplication after the divergence of the cotton and arabidopsis lineages, but prior to the divergence of a-genome and d-genome lineages, leaving each of the modern diploids in our study with a complement of two phya paralogs. indeed, four distinct contigs were observed in both the inbred g. hirsutum cultivar tm- and in the doubled-haploid line g. barbadense 3- for each allotetraploid taxon, two contigs fell into each of the phya- and phya- clades. a conservative inventory of available est sequences indicated that at least two distinct phya loci are expressed in g. hirsutum. within each of the phya and phya clades, the level of nucleotide diversity was very low, with at most four parsimonious nucleotide changes separating each contig. however, within the phya clade, the contigs resolved into two subclades that each contained a single contig from one of the diploid taxa and one contig from each of the allotetraploids. for example, g. raimondii phya grouped in a single contig from each of g. hirsutum and g. barbadense. based on this grouping, the latter contigs were assigned the provisional designation of phya d. similarly, g. herbaceum grouped with g. hirsutum phya a and g. barbadense phya a. based on similar criteria, the phya clade was also divided into phya a and phya d subclades. the phylogenetic resolution of a- and d-genome subclades supported the hypothesis that each of the a- and d-genome diploids contributed both phya and phya to the allotetraploid lineage. thus, although hinge-region nucleotide diversity within each of the phya and phya clades was low, it was sufficient to resolve a tentative phya gene complement for each taxon, as well as the pattern of gene inheritance through the allopolyploidization event. amplification of the phyb/d gene sub-family a ~ bp fragment from the phyb/d hinge region was obtained by amplification using primers phyabnondeg-f and phybdeg-r. sequences from a total of clones yielded a single consensus contig from each of the diploid cottons g. herbaceum and g. raimondii, and from the allotetraploid g. hirsutum. two distinct contigs were assembled from clones derived from the allotetraploid g. barbadense. these clone sequences shared ~85% nucleotide identity with the arabidopsis phyb gene and ~78% nt identity with arabidopsis phyd. all clones had a high-scoring pair relationship with the arabidopsis phyb gene as well as significant similarity to the arabidopsis phyd gene. consensus sequences were aligned across all taxa, yielding a bp alignment with an average pairwise sequence similarity of %, with sites identical across all taxa, no stop codons and no indels. although these data indicated the presence of at least one phyb gene in each of the a- and d-genome diploid plants and in g. hirsutum, and at least two genes phyb genes in the g. barbadense, the low level of nucleotide differentiation observed within the hinge region yielded insufficient phylogenetic information to characterize the phyb gene complement in any of the study taxa. to obtain better resolution of the phyb gene complement, additional low degeneracy primers 1010-f, 1910-f, 1910-r, and 2848-r were used along with primer phyabnondeg-f to create a kb long series of overlapping amplicons corresponding to approximately kb of the arabidopsis phyb gene and extending from the hinge, through the first intron and into the second exon. after amplification, cloning and sequencing, the amplicons were assembled for each taxon. in all gossypium taxa examined, the first intron was ~ bp longer than the first intron of phyb from arabidopsis. unlike the other phytochrome amplicons, we detected a high frequency of putative pcr-mediated recombination events within the phyb kb fragment from amplifications using g. barbadense as template. the recombination detection algorithm rdp identified a number of clones resulting from apparent recombination between the a-genome and d-genome derived homeologous sequences, with predicted breakpoints between nucleotides and of the alignment. after omission of these recombinant clones, composite amplicon sequences from each taxon were aligned, creating a consensus alignment of bp with % average pairwise similarity and identical sites. overall, the cotton phyb genes shared 65% nucleotide identity with the arabidopsis phyb ortholog. no stop codons or indels were detected in exon sequences. a bp putative deletion was observed in one contig from g. hirsutum. in addition, a bp indel was polymorphic between the phyb.a and phyb.d clades. finally, phyb of g. raimondii had an additional bp insertion. all indel polymorphisms were located within first introns. detailed phylogenetic analyses of the bp contigs from a-, d-, and ad-genome cottons indicated the presence of least one phyb locus in the two diploid cottons, g. herbaceum and g. raimondii, and at least two phyb loci in both allotetraploid cottons. the g. hirsutum and g. barbadense sequence contigs each grouped into two sub-clades. the single phyb contig from g. herbaceum was used to define the phyb.a cluster, while the single phyb contig from g. raimondii anchored the phyb.d cluster. from these results, we concluded that phyb.a and phyb.d, which shared ~98% nucleotide sequence identity, arose as orthologs at the time of divergence of the a- and d-genome diploid lineages. we surmised that phyb.a was contributed to the allotetraploids via the a-genome ancestor and phyb.d was contributed via the d-genome ancestor. available est sequences indicated that at least one phyb locus is expressed in g. hirsutum. amplification from the phyc gene sub-family several sets of degenerate primer pairs that were designed on the basis of the conserved hypatdip and pfplryac regions -- including several designed from rosid phyc nucleotide sequences -- failed to produce detectable pcr amplification products from the gossypium species tested. however, the identification of a small est clone with similarity to arabidopsis phyc in a library from g. raimondii floral tissue, allowed us to design the primer phyc_1r_dfci within the c-terminal domain. when used in combination with phydeg-f, this primer amplified a ~ kb fragment composed entirely of coding sequence from the first exon of phyc, including the hinge. all clones obtained using this primer pair had a high-scoring similarity to arabidopsis phyc. from these clones, we assembled a single consensus contig from each of the diploid species g. herbaceum and g. raimondii, and two distinct consensus contigs from each of the allotetraploids g. hirsutum and g. barbadense. consensus sequences for each of the putative phyc contigs were aligned across all taxa, yielding a bp alignment with an average pairwise sequence similarity of %, sites identical across all taxa, with no indels or stop codons in any taxa. in phylogenetic analyses, the phyc consensus sequences grouped into two major clades. one of these clades contained the g. herbaceum contig and one contig from each of g. hirsutum and g. barbadense. this clade was designated phyc.a. the other clade, designated phyc.d, included the g. raimondii contig along with the other of the two contigs from each of g. hirsutum and g. barbadense. these data indicated that both the a- and d-genome ancestors had one phyc gene, and that upon hybridization and polyploidization, this gene was contributed from each diploid to the allotetraploid ancestor of g. hirsutum and g. barbadense. for comparison with the other phytochromes, we also analyzed a portion of the phyc alignment corresponding to the hinge region only. this alignment was nucleotide pairs in length, with pairwise sequence similarity of %, sites identical across all taxa, with no indels. although it encompassed fewer variable nucleotides, nj analysis of the hinge region alone could be used to differentiate the phyc.a and phyc.d clades and to infer the composition and evolutionary inheritance of the phyc gene family in cottons. our failure to obtain phyc hinge amplification with several sets of both universal and rosid specific primers was entirely due to substantial nucleotide differentiation in phyc, particularly within the hinge region. for example, the nt long phydeg-r primer had six nucleotide mismatches with the cotton phyc genes, including three transitions and three transversions. five of the six mismatches occurred at what are considered to be invariant nucleotide positions. it should be noted that these divergent nucleotides in the conserved primer-binding site did not alter the amino acid sequence. the phye gene sub-family phye hinge region consensus contigs from our study taxa formed a bp alignment with an average pairwise similarity of %, with invariant sites, no indels, and no stop codons in any taxa. the consensus of the aligned phye sequences had 80% nucleotide similarity to the corresponding fragment of the arabidopsis phye gene. based on maximum parsimony, nucleotide diversity in the cotton phye hinge sequences could be explained by a minimum of six nucleotide changes, all of which were synonymous. nj analysis of the cotton phye hinge region showed two distinct clades corresponding to the a- and d-genome derived orthologs, a finding consistent with a hypothesis in which each diploid ancestor contributed a single phye ortholog to the allotetraploid lineage. interestingly, while two distinct phye contigs were obtained from g. hirsutum, only a single contig, which grouped with the d-genome clade, was obtained from g. barbadense. available est sequences indicated that at least one phye locus is expressed in g. hirsutum. a global hinge-based alignment of arabidopsis and cotton phytochromes phya, phyb, phyc and phye hinge regions from arabidopsis and gossypium spp. were aligned to create a global phytochrome alignment nucleotides in length, with an average pairwise similarity of % and identical sites. the gene phylogeny generated from this alignment reflected divergence of phya, phyb, phyc and phye as a result of speciation and gene duplication. the level of nucleotide divergence of each of the gene sub-families after nodes 1a, 1b, 1c and 1e was similar, with a mean of ± nucleotide substitutions per site. however, the synonymous and non-synonymous substitution rates were both significantly more variable among the various gene sub-families defined by nodes 1a, 1b, 1c and 1d than were simple nucleotide distances. despite this variation, all sub-families showed a ka/ks ratio <, implying that each remains under purifying selection for function. further, excessively long branch-lengths, which are often found in pseudogenes, were not observed. in the phyb, phyc and phye clades, the branch lengths leading to the arabidopsis orthologs, which have known biological functions, were longer than the branches leading to their respective cotton orthologs. considered together, these lines of evidence indicate that each of the phytochrome sub-families retains some biological function in gossypium, as they do in arabidopsis. further, our topology supports the conclusion that phyd is the result of a relatively recent gene duplication that may be exclusive to the brassicaceae family. nodes refer to the nj tree in figure k-2p indicates the mean kimura 2-parameter distances between arabidopsis and cotton gene sequences. discussion resolution of the phytochrome gene family in three out of four cases, we were able to successfully resolve the inventory and evolutionary relationships of the phytochrome genes in diploid and allotetraploid cottons using the hinge region only. this finding supports the general utility of employing the hinge region for identifying gsts for phytochromes. in only one case was additional gene sequence required for sufficient phylogenetic resolution. in another case, nucleotide divergence at a commonly used primer-binding site prevented the characterization of the hinge region by the typical strategy of using primers based on conserved flanking peptides hypatdip and pfplryac. however, nucleotide diversity within the phyc hinge region itself was sufficiently informative to resolve the pattern of evolutionary inheritance through allotetraploidization event. the sequencing of phytochrome gene fragments from a- and d-genome diploids, as well as from ad allotetraploid taxa, provides an essential foundation for all subsequent analysis of phytochrome function and evolution in gossypium. the sequenced fragments provide sufficient information to unequivocally identify or'tag' various orthologs, homeologs and paralogs, as well as monitor their patterns of nucleotide divergence, and trace their evolutionary inheritance through the allopolyploidization event. this information will serve as a foundation for further sequence assembly and annotation, and will be used to design locus-specific primer sets for quantitative rt-pcr assays that will measure transcript levels for each gene family member. in some cases levels of sequence divergence are high enough to support studies of gene function using rnai or amirna approaches to create gene-specific knockouts. the use of well characterized'candidate genes' of agronomic interest is becoming an integral component of marker-assisted selection efforts in plants. several snp-based molecular markers are now being developed using the diagnostic nucleotide characters identified in this study, and are being mapped in experimental cotton populations that show segregation of phytochrome-controlled traits such as fiber length and flowering time. the ancestral phytochrome gene complement of the malvales and brassicales our study indicated that the diploid ancestors to the world's major fiber crops had a complement of phytochrome apoprotein genes that was very similar to that of the model plant arabidopsis thaliana. this was not entirely unexpected given the relatively close phylogenetic relationship of the two lineages. the most-simple evolutionary scenario is that the last common ancestor of arabidopsis and cotton, possibly an arborescent species in the late cretaceous period, had a phytochrome gene complement consisting of one functional gene in each of the phya, phyb/d, phyc and phye subfamilies. phya duplication in gossypium after the divergence of the malvales and brassicales, the ancestral phya gene underwent duplication resulting in the observed phya- and phya- paralogs of modern gossypium spp. as the a- and d-genome diploids have both paralogs, the duplication event occurred prior to the divergence of the a- and d-genome lineages. using mya as a rough estimate of the time of divergence of the malvales and brassicales, along with our observed ks of in the phya hinge region in this time interval, we can derive a crude estimate of substitutions/synonymous-site/million years, and an estimate of the time of phya duplication of ~ mya. this estimate places the duplication well within the crown group of malvales and the malvaceae family. given our time estimate, the phya duplication may be exclusive to the genus gossypium, but would have occurred prior to the estimated time of divergence of the a and d genome groups. as neither we nor others have observed evidence of additional nuclear gene duplications or chromosomal duplications in this time period, the phya event was likely a tandem or segmental duplication, rather than whole genome duplication. after a gene duplication event, one of the two newly duplicated genes is theoretically unconstrained by selection for function, and is thus free to accumulate mutations leading to a pseudogene fate, subfunctionalization, or neofunctionalization. although we did not obtain definitive evidence of pseudogenic sequences in any of the phytochromes or taxa studied, we did observe significant variation in ka/ks ratios in pairwise interspecific comparisons, leaving open the possibility of pseudogene outcomes. alternatively, one of the duplicated genes may undergo positive selection to gain a novel function. further, duplicated gene-pairs may subdivide the function of ancestral gene. perhaps the most intriguing fate, which has been observed empirically, but not yet explained in theory, is the situation in which both gene copies may be retained for a lengthy period under what appears to be purifying or negative selection. one approach to understanding the evolutionary fates of duplicated genes is through an analysis of the signature of natural selection on amino acid encoding sequences. although the hinge regions of phytochromes display relatively high levels of nucleotide diversity, they do not evolve under neutrality. the hinge region participates in inter-domain communication in phytochrome molecules. for example, phosphorylation of a serine residue in the phya hinge plays a likely role in regulating protein-protein interactions between phytochrome and downstream signal-transducing molecules. compared to wild-type, a mutation in the hinge region of arabidopsis phyb is deficient in localization into distinct nuclear bodies. further, a single nucleotide polymorphism in the hinge of one of two phyb genes in aspen was associated with natural geographic variation in the timing of bud-set. in comparisons between cotton and arabidopsis, the ka/ks ratio for the phya hinge region was -- a value that is typical for genes under purifying selection. in contrast, the ka/ks ratio for phya after gene duplication was, or ~ -fold higher. this value is also ~ -fold greater than the mean ka/ks ratio of all phytochrome hinge regions of approximately ±. this significantly elevated ka/ks ratio after the phya duplication could be attributed to a relaxation of stabilizing selection and/or subfunctionalization of the nascent phya paralogs. the possible functional divergence of phya and phya may be more pronounced after the separation of the a- and d-genome lineages. a comparison of phya in the two diploids yields a ka/ks ratio of ~, primarily due to amino acid substitutions in phya d, while phya has a ka/ks ratio of in the same taxonomic comparisons. although this difference is suggestive of possible differential rates of functional evolution in the paralogs, it is not statistically significant in fisher's exact test. it will be of interest to determine whether the cotton phya paralogs have distinct functions. experiments are underway to determine the respective biological functions of each phya- and phya- in g. hirsutum and g. barbadense using paralog-specific rt-pcr, rnai gene knockout, and tests for genetic associations between phytochrome-controlled phenotypic traits and phya- and phya- specific molecular markers. a'candidate gene' approach has recently been used in soy to uncover a genetic linkage between the photoperiod insensitivity locus e and one of the two the phya genes, designated gmphya and gmphya. loss of photoperiodic flowering is associated with a ty1/copia-like retrotransposon insertion into exon of gmphya the authors argue that gene duplication and partial redundancy of the phya genes may have facilitated the loss of photoperiod sensitivity by allowing the gmphya mutant to avoid the major deleterious phenotypic effects that would have been caused by complete deficiency of phya gene function. d-d indicates means of comparisons within extant diploids, d-t indicates means of comparisons of extant diploids with tetraploids, t-t indicates means of comparisons within tetraploids. s dif, synonymous differences; s pos, synonymous positions; ns dif, non-synonymous differences; ns pos, non-synonymous positions. p indicates significance as determined by fisher's exact test. persistence and loss of phytochrome paralogs after allopolyploidization all phytochromes underwent gene duplication by polyploidization at the time of formation of the ad allotetraploids, on the order of - mya. for example, in g. hirsutum, we detected a minimum set of ten distinct phytochrome genes, including four phya genes. in order to assess the evolutionary trajectory of these recently duplicated genes, we examined the synonymous and non-synonymous divergence rates of a- and d-genome phytochrome orthologs and homeologs in pairwise comparisons of 1) diploids with diploids, 2) diploids with tetraploids, and 3) tetraploids with tetraploids. given that the allotetraploid cottons had both a- and d-genome derived copies of each gene on the order of hundreds of thousands of years, we hypothesized that there may be a relaxation of selection in the allotetraploids, as one of the two copies should no longer be evolutionarily constrained. however, in comparisons of a- vs. d-genome derived orthologs or homeologs for six gsts, we did not observe dramatic differences in ka/ks between diploid and allotetraploids in any gst except the hinge region of phya. because of low levels of nucleotide divergence, we employed fisher's exact test and found no significant differences in the patterns of nucleotide evolution in allotetraploids vs. diploids. thus, there was no broad evidence of relaxation of natural selection on gene function after gene duplication by allotetraploidization. further, the generally low ka/ks ratios across all genes and taxa support a model in which that the phytochrome homeologs are largely evolving independently by a birth-and-death model rather than concerted evolution. the coding sequences of the phyb kb fragment also appeared be evolving under stabilizing selection in both the diploids and allotetraploids reflecting continued selective constraint on coding sequence evolution after polyploidization. however, there was a significant excess of non-synonymous substitutions in both diploids and allotetraploids indicating a partial relaxation of negative selection and/or functional divergence of the phyb homeologs. in the allotetraploid cottons, both phyc.a and phyc.d are also evolving in a pattern consistent with purifying selection. however, it should be noted that the phyc.d clade appears to be evolving at distinctly faster rate than the phyc.a clade. this suggests either a relaxation of purifying selection in, or functional divergence of phyc.d. in a similar study of phytochromes in cultivated sorghum and its wild congeneric relatives, phyc was undergoing faster amino acid evolution than phya or phyb. in the both the phyb and phyc gene subfamilies of cotton, the sequences of the c-terminal signaling domain had higher ka/ks ratios than the corresponding hinge region alone. this may reflect the co-evolution of protein-protein interactions with downstream signaling partners, which are mediated by the c-terminal'signal transduction' domain. while phye-related contigs had low ka/ks values, indicating purifying selection, no contig corresponding to an expected g. barbadense phye.a ortholog was observed. this may have been due to under-sampling of g. barbadense clones for sequencing, or due to nucleotide divergence in primer sites. of the phye-like clone sequences obtained from g. barbadense, all were in the d-genome derived clade, which would be an unlikely result assuming equal amplification efficiencies for phye.a and phye.d. alternatively, the apparent lack of a phye.a ortholog in g. barbadense could be explained by concerted evolution, gene conversion, or by pcr-mediated recombination. overall, the phye genes, like the other cotton phytochromes, had more synonymous than non-synonymous nucleotide substitutions, favoring a birth-and-death model of gene evolution.  our preliminary efforts to obtain an inventory of the cotton phytochrome gene family indicated that diploid a- and d-genome diploid cottons have two paralogous phya genes, and one each of phyb, phyc, and phye gene sub-families. coding sequence evolution in phya was significantly elevated, suggesting loss of selection for function, or incipient subfunctionalization. other than this duplication and the lack of a separate phyd gene, the phytochrome complement of diploid cottons was very similar to that observed in the closely related model plant arabidopsis thaliana, which greatly facilitates cross-species comparisons. whole genome duplication via allopolyploidization resulted in additive amalgamation of phytochrome genes within a single nucleus in the allotetraploid, retaining complete gene complements of at least four phya genes, two genes of each phyb, phyc and phye in ad-genome g. hirsutum. g. barbadense may lack the phye gene contributed by the a-genome ancestor. strong purifying selection on nearly all of the phytochrome genes suggests some level of conservation of function of each of the genes after polyploidization. with the possible exception of one of the phye.a homeologs in g. barbadense, we did not see evidence of gene loss. we did not observe any convincing evidence of concerted evolution by gene conversion. rather, the genes duplicated by allopolyploidy appear to be largely retained, and evolving independently as observed in other nuclear genes in allotetraploid cottons. these results further our understanding of the evolutionary fates of duplicate genes following allopolyploidization. information on key evolutionary events, as well as rates and patterns of evolutionary change, are an important component of the functional annotation of genes and genomes. these data provide the foundation for more comprehensive studies of the biological functions of each of the cotton paralogs and homeologs. the development of phytochrome'candidate gene' markers based on the gsts identified here may prove useful in the mobilization of valuable genes from photoperiodic wild and primitive cottons into elite cotton varieties, in order to improve stress tolerance, disease resistance, fiber quality, and other traits. 
we identified two paralogous phya genes in diploid cottons, the result of a malvaceae-specific phya gene duplication that occurred approximately million years ago, before the divergence of the a- and d-genome ancestors. detailed sequence analysis suggests that phytochrome genes retained after duplication by segmental duplication and allopolyploidy appear to be evolving independently under a birth-and-death-process with strong purifying selection. our study provides a preliminary phytochrome gene inventory that is necessary and sufficient for further characterization of the biological functions of each of the cotton phytochrome genes, and for the development of'candidate gene' markers that are potentially useful for cotton improvement via modern marker-assisted selection strategies. little is known about the composition and evolution of the phytochrome gene family in diploid or allotetraploid cotton species. the objective of this study was to obtain a preliminary inventory and molecular-evolutionary characterization of the phytochrome gene family in cotton. the allotetraploid genomes have largely retained the complete gene complements inherited from both of the diploid genome ancestors, with at least four phya genes and two genes encoding phyb, phyc and phye in the ad-genomes. we identified a single gene copy of phyb, phyc, and phye in diploid cottons. we used comparative sequence resources to design low-degeneracy pcr primers that amplify genomic sequence tags for members of the phya, phyb/d, phyc and phye gene sub-families from a- and d-genome diploid and ad-genome allotetraploid gossypium species. phytochromes are a family of red/far-red photoreceptors that regulate a number of important developmental traits in cotton, including plant architecture, fiber development, and photoperiodic flowering.
 phytochromes are a family of red/far-red photoreceptors that regulate a number of important developmental traits in cotton, including plant architecture, fiber development, and photoperiodic flowering. little is known about the composition and evolution of the phytochrome gene family in diploid or allotetraploid cotton species. the objective of this study was to obtain a preliminary inventory and molecular-evolutionary characterization of the phytochrome gene family in cotton.  we used comparative sequence resources to design low-degeneracy pcr primers that amplify genomic sequence tags for members of the phya, phyb/d, phyc and phye gene sub-families from a- and d-genome diploid and ad-genome allotetraploid gossypium species. we identified two paralogous phya genes in diploid cottons, the result of a malvaceae-specific phya gene duplication that occurred approximately million years ago, before the divergence of the a- and d-genome ancestors. we identified a single gene copy of phyb, phyc, and phye in diploid cottons. the allotetraploid genomes have largely retained the complete gene complements inherited from both of the diploid genome ancestors, with at least four phya genes and two genes encoding phyb, phyc and phye in the ad-genomes. we did not identify a phyd gene in any cotton genomes examined.  detailed sequence analysis suggests that phytochrome genes retained after duplication by segmental duplication and allopolyploidy appear to be evolving independently under a birth-and-death-process with strong purifying selection. our study provides a preliminary phytochrome gene inventory that is necessary and sufficient for further characterization of the biological functions of each of the cotton phytochrome genes, and for the development of'candidate gene' markers that are potentially useful for cotton improvement via modern marker-assisted selection strategies. 