Polymorphism and divergence at three duplicate genes in Brassica nigra

The CONSTANS-like gene family has been shown to evolve exceptionally fast in Brassicaceae. In the present study we analyzed sequence polymorphism and divergence of three genes from this family: COL1 (CONSTANS-LIKE 1) and two copies of CO (CONSTANS), COa and COb, in B. nigra. There was a significant fourfold difference in overall nucleotide diversity among the three genes, with BniCOb having twice as much variation as BniCOL1, which in turn was twice as variable as BniCOa. The ratio of nonsynonymous-to-synonymous substitutions (dN/dS) was high for all three genes, confirming previous studies. While we did not detect evidence of selection at BniCOa and BniCOb, there was a significant excess of polymorphic synonymous mutations in a McDonald-Kreitman test comparing COL1 in B. nigra and A. thaliana. This is apparently the result of an increase in selective constraint on COL1 in B. nigra combined with a decrease in A. thaliana. In conclusion, a complex scenario involving both demography and selection seems to have shaped the pattern of polymorphism at the three genes.     

Recent degeneration of an old duplicated flowering time gene in Brassica nigra

Gene and genome duplications play a major role in the evolution of plant species. The Brassica nigra genome is highly replicated as a result of ancient polyploidization events. Two copies of the flowering time gene CONSTANS (COa and COb) have been identified in B. nigra, and previous studies showed that COa is functional. In the present study, the polymorphism of 92 COb alleles sampled in seven populations was analyzed. Both polymorphism and recombination levels were elevated and varied strongly among populations and 8% of COb alleles exhibit apparently disabling mutations. Sequence data, however, do not provide unambiguous support for the presence of relaxed selective constraint on COb as compared to known functional CO genes. On the one hand, some of the disabling mutations reached high-frequency arguing for a loss of function but, on the other hand, the ratio of nonsynonymous to synonymous nucleotide polymorphism and diversity is low and similar to that observed in other B. nigra CO and CO-like genes, supporting the conservation of some function. We also showed that COb is still transcribed. Finally, the flowering time of Arabidopsis thaliana co mutant plants transformed with COb alleles with and without apparent disabling mutations was similar. We propose that COb was retained for a long period after duplication, but a recent fixation of a detrimental mutation, possibly as an effect of a bottleneck, resulted in its nonfunctionalization. We also speculate as to the presence of subsequent selection for rapid degeneration of the gene.     

Sequence divergence at the putative flowering time locus COL1 in Brassicaceae

An insertion/deletion polymorphism (Ind2) in the Brassica nigra CONSTANS LIKE 1 (Bni COL1) gene was previously found to be associated with variation in flowering time. In the present study we examine the inter-specific divergence of COL1 in the family Brassicaceae. Analysis of codon substitution models did not reveal evidence of positive Darwinian selection, but comparisons of the COL1 gene in different species revealed a surprising number of indels. A total of 24 indels were found in the 650 bp of the middle variable region of the gene. This high number of indels could reflect a lack of constraint on length of this region of the protein, or the effect of positive selection. The number of indels was close to that expected in non-coding DNA, but the indels were longer in COL1 than those observed in non-coding regions. Reconstruction of indel evolution indicated that most indels resulted from deletions rather than insertions. The Ind2 indel that has shown association with flowering time in Brassica nigra exhibited a remarkable distribution in the Brassicaceae family, indicating that the polymorphism may have persisted more than ten million years. Considering presumed historic populations sizes of Brassicaceae species, such a long persistence time seems unlikely for a neutral polymorphism.     

Comparative mapping in Arabidopsis and Brassica, fine scale genome collinearity and congruence of genes controlling flowering time

The model dicotyledonous plant, Arabidopsis thaliana , is closely related to Brassica crop species. It is intended that information concerning the genetic control of basic biological processes in Arabidopsis will be transferable to other species. Genome collinearity and its potential to facilitate the identification of candidate genes in Arabidopsis homologous to genes controlling important agronomic traits in Brassica was investigated. Genetic mapping in B. nigra identified two loci influencing flowering time (FT), with loci on linkage groups 2 and 8 explaining 53% and 12% of the total variation in FT, respectively. The CO gene exerts an important control over FT in A. thaliana , and B. nigra homologues of CO probably also play an important role in regulating FT. B. nigra homologues of CO were identified on linkage groups 2 and 8, the homologue on group 2 was coincident with the major locus controlling FT while the homologue on group 8 was within the 90% confidence interval of the weaker FT gene. The CO homologue on group 2 exhibits abundant allelic variation suggesting that it naturally controls a wide range of flowering times. Fine-scale A. thaliana/B. nigra comparative mapping demonstrated short-range collinearity between the genomes of Arabidopsis and Brassica . Eleven DNA fragments spaced over a 1.5 Mb contig in A. thaliana were used as RFLP probes in B. nigra . Three collinear representations of the A. thaliana contig were identified in B. nigra , with one interrupted by a large chromosomal inversion. Collinearity over this range will allow the resources generated by the Arabidopsis genome project to facilitate map-based cloning in Brassica crops.     

Conserved structure and function of the Arabidopsis flowering time gene CONSTANS in Brassica napus

The Arabidopsis thaliana CONSTANS (CO) gene which promotes flowering in long days was recently isolated by chromosome walking. The mapping of QTLs controlling flowering time in Brassica species has identified genomic regions that contain homologues of the CO gene. Four genes homologous to the Arabidopsis CO gene were isolated from a pair of homoeologous loci in each of two doubled-haploid Brassica napus lines displaying different flowering times, N-o-1 and N-o-9. The four genes, BnCOa1, BnCOa9, BnCOb1 and BnCOb9, are located on linkage groups N10 and N19, and are highly similar to each other and to the Arabidopsis CO gene. Two regions of the proteins are particularly well conserved, a N-terminal region with two putative zinc fingers and a C-terminal region which may contain a nuclear localization signal. All four genes appear to be expressed in B. napus. The BnCOa1 allele was shown to complement the co-2 mutation in Arabidopsis in a dosage-dependent manner causing earlier flowering than in wild type under both long- and short-day conditions.     

Structural organization of the estalpha3(1) gene in a Colombian strain of Culex quinquefasciatus differs from that in Cuba

In Culex mosquitoes (Diptera: Culicidae), the most common mechanism for resistance to organophosphorus (OP) insecticides involves amplification of one or more esterases. Two esterase loci are often involved, with different allelic forms co-amplified. Estalpha3(1) is co-amplified with estbeta1 in a Colombian (COL) strain of Culex quinquefasciatus Say. These two alleles co-migrate on acrylamide gels, often leading to misscoring of the phenotype as elevation of a single estbeta enzyme. By sequencing COL genomic DNA, we determined the estalpha3(1) gene length is 1623 nucleotides. The open reading frame of estalpha3(1) encodes a 540 amino acid protein, as for estalpha2(1) in strain Pel RR from Sri Lanka. The intron/exon boundaries of estalpha3(1) are identical to those of estalpha2(1), suggesting that they are alleles of the same locus. The COL estalpha3(1) gene differs from estaalpha3(2) in strain MRES from Cuba, although they have equivalent electrophoretic mobility, showing that these two strains contain distinct resistance-associated amplicons. Twenty nucleotide differences were scored between the MRES partial 495 bp sequence and that in the COL strain, with two amino acid changes, demonstrating distinct estalpha enzymes. Our sequencing data show 95% identity between the three estalpha genes (each has six introns and seven exons) in OP-resistant Cx. quinquefasciatus. Amplified estalpha3(1) and estbeta1 are at least 10kb apart in temephos-selected COL and 2.7kb apart in Pel RR, whereas these non-amplified genes are only 1.7kb apart in the nonselected parental COL stock, as in Pel SS (susceptible Sri Lankan strain), demonstrating that this region of the genome is susceptible to expansion and contraction.     

RIRE2, a novel gypsy-type retrotransposon from rice

The 441-bp DNA segment in a PCR-amplified fragment from Oryza sativa cv. IR36 was found to have a sequence with features characteristic of LTRs of retroelements, which was named RIRE2 (Rice retroelement #2) and further analyzed. Cloning and sequencing analyses of the DNA segments connected to LTR-like sequence showed that RIRE2 has a long internal region almost 10 kb long that is flanked by LTR-like sequences. This internal region carries a primer binding site (PBS) and polypurine tract (PPT) which are necessary for cDNA synthesis of retroelements. The PBS sequence is complementary to the 3' end region of tRNA(Arg). The internal region has an rt gene homologous to that of gypsy-type retrotransposons, evidence that RIRE2 is indeed a retrotransposon related to gypsy from Drosophila. RIRE2 has an extra sequence more than 4 kb long in the region downstream of gag-pol. Phylogenetic analysis of the putative amino-acid sequences of the rt gene as well as the int gene showed that RIRE2 is related to a group of gypsy-type retrotransposons of a large size that include Grande1-4 of teosinte, Tat4-1 and Athila1-1 of Arabidopsis thaliana, and Cyclops-2 of pea, but distantly related to any other group of gypsy-type retrotransposons, including RIRE3 and RIRE8 of rice. RIRE2 and Grande1-4 had the highest homology in the gag-pol region, but the nucleotide sequences of the LTR regions differed. Both elements had significant homology in the middle area of the extra regions downstream of gag-pol, in which they had an open reading frame encoding a protein with no known function on the opposite strand from that coding for gag-pol.     

DNA sequence variation in BpMADS2 gene in two populations of Betula pendula

The PISTILLATA (PI) homologue, BpMADS2, was isolated from silver birch (Betula pendula Roth) and used to study nucleotide polymorphism. Two regions (together about 2450 bp) comprising mainly untranslated sequences were sequenced from 10 individuals from each of two populations in Finland. The nucleotide polymorphism was low in the BpMADS2 locus, especially in the coding region. The synonymous site overall nucleotide diversity (pis) was 0.0043 and the nonsynonymous nucleotide diversity (pia) was only 0.000052. For the whole region, the pi values for the two populations were 0.0039 and 0.0045, and for the coding regions, the pi values were only 0 and 0.00066 (for the corresponding coding regions of Arabidopsis thaliana PI world-wide pi was 0.0021). Estimates of pi or theta did not differ significantly between the two populations, and the two populations were not diverged from each other. Two classes of BpMADS2 alleles were present in both populations, suggesting that this gene exhibits allelic dimorphism. In addition to the nucleotide site variation, two microsatellites were also associated within the haplotypes. This allelic dimorphism might be the result of postglacial re-colonization partly from northwestern, partly from southeastern/eastern refugia. The sequence comparison detected five recombination events in the regions studied. The large number of microsatellites in all of the three introns studied suggests that BpMADS2 is a hotspot for microsatellite formation.     

Differences in the nuclease sensitivity between the two alleles of the immunoglobulin kappa light chain genes in mouse liver and myeloma nuclei.

In mouse myeloma T the productive kappa light chain gene differs from its aberrantly rearranged allele in the patterns of DNAase I hypersensitive sites. In the region of the alleles where they are identical in sequence they have one site in common which lies 0.8 kb downstream of the coding region; but two sites upstream of and within the C gene segment (2) are found only on the non-productive allele. Within the region of different sequences both alleles have analogously located DNAase I hypersensitive sites; they lie 0.15 kb upstream of the respective leader segments and cover putative promoter sequences. Only one of the six DNAase I hypersensitive sites is also very sensitive towards micrococcal nuclease due to its particular DNA sequence. The non-rearranged gene studied in liver nuclei has no DNAase I hypersensitive sites but is preferentially cleaved in A/T rich regions.     

The CCAAT binding factor can mediate interactions between CONSTANS-like proteins and DNA

CONSTANS-Like (COL) proteins are plant-specific nuclear regulators of gene expression but do not contain a known DNA-binding motif. We tested whether a common DNA-binding protein can deliver these proteins to specific cis-acting elements. We screened for proteins that interact with two members of a subgroup of COL proteins. These COL proteins were Tomato COL1 (TCOL1), which does not seem to be involved in the control of flowering time, and the Arabidopsis thaliana CONSTANS (AtCO) protein which mediates photoperiodic induction of flowering. We show that the C-terminal plant-specific CCT (CO, CO-like, TIMING OF CAB EXPRESSION 1) domain of both proteins binds the trimeric CCAAT binding factor (CBF) via its HAP5/NF-YC component. Chromatin immunoprecipitation demonstrated that TCOL is recruited to the CCAAT motifs of the yeast CYC1 and HEM1 promoters by HAP5. In Arabidopsis, each of the three CBF components is encoded by several different genes that are highly transcribed. Under warm long days, high levels of expression of a tomato HAP5 (THAP5a) gene can reduce the flowering time of Arabidopsis. A mutation in the CCT domain of TCOL1 disrupts the interaction with THAP5 and the analogous mutation in AtCO impairs its function and delays flowering. CBFs are therefore likely to recruit COL proteins to their DNA target motifs in planta.     

Sequence variation and haplotype structure at the putative flowering-time locus COL1 of Brassica nigra

Motivated by a previous study indicating that polymorphism at an indel, Ind2, within the Brassica nigra COL1 gene is significantly associated with flowering time, we searched for evidence of selection in a sample of 41 complete sequences of B. nigra COL1. The within-gene population recombination rate is moderate, and all neutrality tests used in the present study failed to detect departure from the standard neutral model or evidence of selection. The haplotype structure of the 5'-half of the gene is primarily associated with the demographic history of the species and more specifically with the split between European and Ethiopian populations, whereas the structure of the 3'-half reflects the polymorphism at Ind2. This could be the result of selection or a combination of recombination and migration during the history of the sample of sequences. Without additional information on polymorphism in flanking areas, these two alternatives are difficult to tell apart. If selection acted on the gene, we suggest that if the indel itself is not the target of selection, among the polymorphic sites cosegregating with the polymorphism at Ind2, replacement polymorphisms around sites 890 and 1260 are the most likely quantitative trait nucleotides within the gene.     

Sequence variation and haplotype structure surrounding the flowering time locus FRI in Arabidopsis thaliana

Linkage disequilibrium in highly selfing organisms is expected to extend well beyond the scale of individual genes. The pattern of polymorphism in such species must thus be studied over a larger scale. We sequenced 14 short (0.5-1 kb) fragments from a 400-kb region surrounding the flowering time locus FRI in a sample of 20 accessions of Arabidopsis thaliana. The distribution of allele frequencies, as quantified by Tajima's D, varies considerably over the region and is incompatible with a standard neutral model. The region is characterized by extensive haplotype structure, with linkage disequilibrium decaying over 250 kb. In particular, recombination is evident within 35 kb of FRI in a haplotype associated with a functionally important allele. This suggests that A. thaliana may be highly suitable for linkage disequilibrium mapping.     

Photoperiod-regulated expression of the PpCOL1 gene encoding a homolog of CO/COL proteins in the moss Physcomitrella patens

The CONSTANS (CO) protein is a critical regulator of the photoperiodic control of flowering in Arabidopsis thaliana and Oryza sativa. We isolated a cDNA PpCOL1 encoding a homolog of the CO/CO-LIKE (COL) family proteins from a cryptogam Physcomitrella patens. The predicted PpCOL1 protein has N-terminal zinc finger and C-terminal CCT domains, which are conserved in the angiosperm CO/COL proteins. Structurally, PpCOL1 is the most closely related to the Group Ia or Ic proteins, which include AtCO and AtCOL1/2, among diverged members of the family. A transient expression assay using GFP showed that the CCT domain of PpCOL1 contains a nuclear-localizing signal. Northern blotting analyses revealed that the PpCOL1 expression is controlled by the circadian clock, and moreover, it is photoperiodically regulated at a gametophore stage when the rate of sporophyte formation is affected by day length. These observations indicate a possible involvement of PpCOL1 as a nuclear factor in the photoperiodic regulation of reproduction of Physcomitrella.     

Extensive variation in the ddl gene of penicillin-resistant Streptococcus pneumoniae results from a hitchhiking effect driven by the penicillin-binding protein 2b gene

An internal fragment of the ddl gene, encoding the cytoplasmic enzyme D-alanyl-D-alanine ligase, was sequenced from 566 isolates of Streptococcus pneumoniae and single isolates of Streptococcus mitis and Streptococcus oralis. The 52 alleles found among the S. pneumoniae isolates fell into two groups. Group A alleles were very uniform in sequence and were present in both penicillin-susceptible and penicillin-resistant pneumococci. Group B alleles were much more diverse and were found only in penicillin-resistant isolates. The Streptococcus oralis and Streptococcus mitis alleles were less diverged from group A alleles than some of the group B pneumococcal alleles, suggesting that the latter alleles contain interspecies recombinational replacements. The ddl gene was located 783 bp downstream of the penicillin-binding protein 2b gene (pbp2b). Sequencing of the pbp2b-recR-ddl-murF region of three penicillin-resistant pneumococci that had diverged ddl alleles showed that the whole region from pbp2b to ddl (or beyond) was highly diverged (about 8%) compared with the sequences from three penicillin-susceptible isolates. The high levels of diversity in the group B ddl alleles from penicillin-resistant isolates were ascribed to a hitchhiking effect whereby interspecies recombinational exchanges at pbp2b, selected by penicillin usage, often extend into, or through, the ddl gene. The data allow the average size of the interspecies recombinational replacements to be estimated at about 6 kb.