The effect of mating system differences on nucleotide diversity at the phosphoglucose isomerase locus in the plant genus Leavenworthia.

To test the theoretical prediction that highly inbreeding populations should have low neutral genetic diversity relative to closely related outcrossing populations, we sequenced portions of the cytosolic phosphoglucose isomerase (PgiC) gene in the plant genus Leavenworthia, which includes both self-incompatible and inbreeding taxa. On the basis of sequences of intron 12 of this gene, the expected low diversity was seen in both populations of the selfers Leavenworthia uniflora and L. torulosa and in three highly inbreeding populations of L. crassa, while high diversity was found in self-incompatible L. stylosa, and moderate diversity in L. crassa populations with partial or complete self-incompatibility. In L. stylosa, the nucleotide diversity was strongly structured into three haplotypic classes, differing by several insertion/deletion sequences, with linkage disequilibrium between sequences of the three types in intron 12, but not in the adjacent regions. Differences between the three kinds of haplotypes are larger than between sequences of this gene region from different species. The haplotype divergence suggests the presence of a balanced polymorphism at this locus, possibly predating the split between L. stylosa and its two inbreeding sister taxa, L. uniflora and L. torulosa. It is therefore difficult to distinguish between different potential causes of the much lower sequence diversity at this locus in inbreeding than outcrossing populations. Selective sweeps during the evolution of these populations are possible, or background selection, or merely loss of a balanced polymorphism maintained by overdominance in the populations that evolved high selfing rates.     

DNA polymorphism in active gene and pseudogene of the cytosolic phosphoglucose isomerase (PgiC) loci in Arabidopsis halleri ssp. gemmifera

DNA variations in two PgiC loci were investigated in 15 strains of Arabidopsis halleri ssp. gemmifera. In a 5.5-kb region of the PgiC1 locus, 127 nucleotide substitutions and 33 length variations were observed. In a 6.0-kb region of the PgiC2 locus, 138 nucleotide substitutions and 33 length variations were observed. Frame shift, novel stop codons, and large length variations were observed in the PgiC2 coding region. These findings suggested that PgiC2 may be a pseudogene. The nucleotide diversities (pi) for the entire regions of both PgiC loci were approximately 0.0033. Tajima's test of both PgiC loci yielded significantly negative results. In the coding regions, the high proportions of replacement substitutions caused significant deviations from neutrality in McDonald and Kreitman's test. An excess of singletons and a high proportion of replacement polymorphic sites have been observed in the Adh and ChiA regions of A. halleri ssp. gemmifera. Thus, the A. halleri ssp. gemmifera population may not have reached equilibrium, and thus nonneutral patterns of DNA polymorphism were observed.     

DNA polymorphism at the cytosolic phosphoglucose isomerase (PgiC) locus of the wild plant Arabidopsis thaliana

DNA variation in a 4.7-kb region of the cytosolic phosphoglucose isomerase (PgiC) locus was investigated for 21 ecotypes of Arabidopsis thaliana. The estimated nucleotide diversity was 0.0038, which was one-third of those in previously investigated loci. Since most of the nucleotide variations (93%) were singleton and doubleton, Tajima's test statistic was significantly negative. About 50% of nucleotide polymorphisms in exons were replacement, which caused significance in McDonald and Kreitman's test when compared with Arabis gemmifera and Cardaminopsis petraea. These results indicated that DNA polymorphism at the PgiC locus was not under neutrality. There were two divergent sequence types in the PgiC region, which were associated with allozyme variation. The Fast allozyme was shown to have originated from the Slow allozyme, since two outgroup species had the Slow form. A phylogenetic tree of ecotypes with the Fast allozyme had the shape of a star phylogeny. Mismatch distribution of the Fast allozyme ecotypes resembled that expected under an expanding population model. These results suggest positive selection for the Fast allozyme of the PGIC in A. thaliana.     

A statistical test for the difference in the amounts of DNA variation between two populations

A statistical test for the difference in the amounts of DNA variation between two populations is developed. The test statistic involves the covariance of the amount of variation between two populations, which is given by a function of their divergence time, T0. Accordingly, the power (rejection probability) of the test depends on T0. In this article, T0 is treated as unknown because it is very difficult to estimate. The test is most conservative when T0 = infinity is assumed because the covariance is zero. If T0 = 0 is assumed, the largest value of the rejection probability is obtained. Thus, the test provides a range of rejection probability unless we have a reliable estimate of T0. The test is applied to the PgiC region in three mustard species: Leavenworthia stylosa, L. crassa and L. uniflora. The results of our test show that the level of variation in L. stylosa is significantly higher than those in the other species.     

The introgression of a functional nuclear gene from Poa to Festuca ovina

In sheep's fescue, Festuca ovina, genes coding for the cytosolic enzyme phosphoglucose isomerase, PGIC, are not only found at the standard locus, PgiC1, but also at a segregating second locus, PgiC2. We have used PCR-based sequencing to characterize the molecular structure and evolution of five PgiC1 and three PgiC2 alleles in F. ovina. The three PgiC2 alleles were complex in that they carried two gene copies: either two active genes or one active and one pseudogene. All the PgiC2 sequences were very similar to each other but highly diverged from the five PgiC1 sequences. We also sequenced PgiC genes from several other grass species. Phylogenetic analysis of these sequences indicates that PgiC2 has introgressed into F. ovina from the distant genus Poa. Such an introgression may, for example, follow from a non-standard fertilization with more than one pollen grain, or a direct horizontal gene transfer mediated by a plant virus.     

Phylogeny of Leavenworthia S-alleles suggests unidirectional mating system evolution and enhanced positive selection following an ancient population bottleneck

The adoption of self-fertilization from an ancestral outcrossing state is one of the most common evolutionary transitions in the flowering plants. In the mustard family, outcrossing is typically enforced by sporophytic self-incompatibility (SI), but there are also many self-compatible species. The genus Leavenworthia contains taxa that either possess or lack SI. Here, we present data showing that SI is associated with strict outcrossing and that there is widespread trans-specific sequence polymorphism at the locus involved in the recognition of self-pollen (the S-locus). This ancestral polymorphism is consistent with the presence of an outcrossing mating system in the common ancestor of Leavenworthia species, and suggests that there have been several independent losses of SI in the group. When compared with other mustard species, the bulk of Leavenworthia S-allele sequences are highly diverged from those found in other Brassicaceae and show relatively low levels of nucleotide diversity, a pattern that suggests the common ancestor of the genus likely underwent a strong population bottleneck. The hypothesis of postbottleneck S-locus rediversification is supported by tests showing stronger positive selection acting on S-alleles from Leavenworthia than those found in other Brassicaceae.     

Determining the physical limits of the Brassica S locus by recombinational analysis

A genetic analysis was performed to study the frequency of recombination for intervals across the Brassica S locus region. No recombination was observed between the S locus glycoprotein gene and the S receptor kinase gene in the segregating populations that we analyzed. However, a number of recombination breakpoints in regions flanking these genes were identified, allowing the construction of an integrated genetic and physical map of the genomic region encompassing one S haplotype. We identified, based on the pollination phenotype of plants homozygous for recombinant S haplotypes, a 50-kb region that encompasses all specificity functions in the S haplotype that we analyzed. Mechanisms that might operate to preserve the tight linkage of self-incompatibility specificity genes within the S locus complex are discussed in light of the relatively uniform recombination frequencies that we observed across the S locus region and of the structural heteromorphisms that characterize different S haplotypes.     

A recently silenced, duplicate PgiC locus in Clarkia

Previous electrophoretic analysis showed that 17 diploid species of the wildflower Clarkia (Onagraceae) have two cytosolic isozymes of phosphoglucose isomerase (PGIC; EC 5.3.1.9), whereas 15 other diploid species have a single PGIC. Molecular studies revealed that the two isozymes in the former species are encoded by duplicate genes, PgiC1 and PgiC2, whereas the single isozyme in the latter is always encoded by PgiC1. Phylogenetic analysis of the nucleotide sequences implied that PgiC2 was silenced four times independently in the genus. Here we describe a psi PgiC2 from C. mildrediae, a species in which only PgiC1 is expressed. The discovery of the psi PgiC2 is significant because it confirms a formal prediction of the phylogenetic analysis. The psi PgiC2 includes 5,039 nucleotides corresponding to 18 of the 23 exons of PgiC, as well as the intervening introns and 3' nontranslated region. The absence of an increase of nucleotide substitutions in its "exons" suggests that the gene was silenced recently. The present study appears to be the first to establish that a specific duplicate gene locus regularly expressed in a group of related plant species has been silenced in one of them. The multiple independent silencings of PgiC2 suggest that it remained functional but inessential in ancestral lineages. We discuss the possibility that PgiC2 may have been preserved in these lineages by selection against mutants causing defective PGIC1- PGIC2 heterodimers.      

Characterization of duplicated two cytosolic phosphoglucose isomerase (PgiC) loci in Arabidopsis halleri ssp. gemmifera

Arabidopsis halleri ssp. gemmifera has two cytosolic phosphoglucose isomerase (PgiC) loci. A 48-bp deletion was observed in the junction of exon 17 and intron 17 for a locus (PgiC2). PCR-RFLP analysis using cDNA template did not detect the PgiC2 locus. Another locus (PgiC1) has common structure with A. thaliana and expressed normally. A phylogenetic tree of PgiC sequences revealed that duplication of the two loci in A. gemmifera occurred after species splitting of A. thaliana and A. gemmifera. More than 12 kb region encompassing PgiC was sequenced for both loci. In both PgiC1 and PgiC2, sequence homologous to A. thaliana PgiC 5' upstream region was not detected. A gene located on chromosome 4 of A. thaliana was detected in the 5' upstream of PgiC2. This result suggested that the microsyntheny around the PgiC region between A. thaliana and A. gemmifera is not established.     

Structural and transcriptional analysis of the self-incompatibility locus of almond: identification of a pollen-expressed F-box gene with haplotype-specific polymorphism

Gametophytic self-incompatibility in Rosaceae, Solanaceae, and Scrophulariaceae is controlled by the S locus, which consists of an S-RNase gene and an unidentified "pollen S" gene. An approximately 70-kb segment of the S locus of the rosaceous species almond, the S haplotype-specific region containing the S-RNase gene, was sequenced completely. This region was found to contain two pollen-expressed F-box genes that are likely candidates for pollen S genes. One of them, named SFB (S haplotype-specific F-box protein), was expressed specifically in pollen and showed a high level of S haplotype-specific sequence polymorphism, comparable to that of the S-RNases. The other is unlikely to determine the S specificity of pollen because it showed little allelic sequence polymorphism and was expressed also in pistil. Three other S haplotypes were cloned, and the pollen-expressed genes were physically mapped. In all four cases, SFBs were linked physically to the S-RNase genes and were located at the S haplotype-specific region, where recombination is believed to be suppressed, suggesting that the two genes are inherited as a unit. These features are consistent with the hypothesis that SFB is the pollen S gene. This hypothesis predicts the involvement of the ubiquitin/26S proteasome proteolytic pathway in the RNase-based gametophytic self-incompatibility system.     

In Southern Africa, brown oculocutaneous albinism (BOCA) maps to the OCA2 locus on chromosome 15q: P-gene mutations identified

In southern Africa, brown oculocutaneous albinism (BOCA) is a distinct pigmentation phenotype. In at least two cases, it has occurred in the same families as tyrosinase-positive oculocutaneous albinism (OCA2), suggesting that it may be allelic, despite the fact that this phenotype was attributed to mutations in the TYRP1 gene in an American individual of mixed ancestry. Linkage analysis in five families mapped the BOCA locus to the same region as the OCA2 locus (maximum LOD 3.07; theta=0 using a six-marker haplotype). Mutation analysis of the human homologue of the mouse pink-eyed dilution gene (P), in 10 unrelated individuals with BOCA revealed that 9 had one copy of the 2.7-kb deletion. No other mutations were identified. Additional haplotype studies, based on closely linked markers (telomere to centromere: D15S1048, D15S1019, D15S1533, P-gene 2.7-kb deletion, D15S219, and D15S156) revealed several BOCA-associated P haplotypes. These could be divided into two core haplotypes, suggesting that a limited number of P-gene mutations give rise to this phenotype.     

DNA polymorphism haplotypes of the human apolipoprotein APOA1-APOC3-APOA4 gene cluster

Summary The genes coding for apolipoproteins A1, C3, and A4 (APOA1, APOC3, APOA4) are closely linked and tandemly organized within a 15-kilobase (kb) DNA segment on the long arm of human chromosome 11. The nucleotide variability of a 61-kb DNA segment containing these genes and their flanking sequences was studied by restriction analysis of a sample of 18 unrelated Northern Europeans using seven different genomic DNA probes. Eleven restriction site polymorphisms located within this DNA segment were used for haplotype analysis of 129 Mediterranean and 67 American black chromosomes. Estimation of the extent of nonrandom association between these polymorphisms indicated considerable linkage disequilibrium within the APOA1-APOC3-APOA4 gene cluster. Several haplotypes arose by recombination, and the rate of recombination within this gene cluster was estimated to be at least 4 times greater than that expected based on uniform recombination. The polymorphism information content of each of these polymorphisms, taken individually, ranges between 0.053 and 0.375, while that of their haplotypes ranges between 0.858 and 0.862. Therefore, DNA polymorphism haplotypes in the APOA1-APOC3-APOA4 gene cluster constitute a highly informative genetic marker on the long arm of human chromosome 11.     

Marked genetic polymorphism of the swine steroid 21-hydroxylase gene, and its location between the SLA class I and class II regions

Restriction fragment length polymorphism (RFLP) analysis of the swine 21-hydroxylase (CYP21) region was conducted on 31 unrelated SLA class I typed pigs, mainly Large Whites, including 15 haplotypes. Ten haplotypes were from SLA genotypic homozygotes and five were from SLA class I phenotypic homozygotes. DNA digestion with Hin dIII, TaqI and PstI, and hybridization to a 4.5-kb swine CYP21 genomic probe yielded respectively two, four and three RFLP patterns. Six patterns were identified with combined RFLP. In addition, analysis of the CYP21 region in families comprising several SLA recombinants demonstrated that the CYP21 gene lies in the DNA segment between the SLA class I and class II regions. These overall results reinforce our previous conclusion about the existence in the pig of a single 21-hydroxylase gene. The characterization of at least six CYP21 allelic patterns provides a new tool for studying the associations between the SLA region and zootechnical traits.     

Distinguishing the Forces Controlling Genetic Variation at the Xdh Locus in Drosophila Pseudoobscura

Fifty-eight isochromosomal lines sampled from two natural populations of Drosophila pseudoobscura in California and one from Bogota, Colombia, were examined using four-cutter restriction mapping. A 4.6-kb region of the xanthine dehydrogenase locus was probed and 66 of 135 restriction sites scored were polymorphic. This predicts that on average every 12th bp would be polymorphic in this region for the genes surveyed if polymorphism occurred randomly along the coding region. In addition, there were 12 insertion/deletion polymorphisms. Forty-nine distinct haplotypes were recognized in the 58 lines examined. The most common haplotype obtained a frequency of only 5%. Measures of base pair heterozygosity (0.0097) and linkage disequilibrium lead to a predicted population size in the range of 1.2-2.4 X 10(6) for the species. High levels of recombination (including gene conversion) can be inferred from the presence of all four gametic types in the data set.     

The extent of nucleotide polymorphism is highly variable across a 3-kb region on Plasmodium falciparum chromosome 2

Genomic nucleotide polymorphism in the virulent human malarial parasite Plasmodium falciparum was surveyed by sequencing a 3-kb region of chromosome 2 from 21 isolates, including the MSP4 and MSP5 genes. Extensive sequence polymorphism was observed in the coding regions of these genes and in the region downstream to MSP5, and the average pairwise divergence time of haplotypes in this region was estimated to be at least about 200,000 years. But nucleotide polymorphism was not found in the introns and was much reduced in the intergenic region. Over the entire region, nucleotide diversity was negatively correlated with a nucleotide content skewed toward thymine. Together with the previous evidence of limited nucleotide polymorphism in introns of P. falciparum, these data suggest the existence of a mechanism suppressing single-nucleotide polymorphism in regions of the P. falciparum genome with highly skewed nucleotide content.     

Molecular cloning of the t complex responder genetic locus

Although mouse t haplotypes carry recessive mutations causing male sterility and embryonic lethality, they persist in wild mouse populations via male transmission ratio distortion (TRD). Genetic evidence suggests that at least five t-haplotype-encoded loci combine to cause TRD. One of these loci, called the t complex responder (Tcr), is absolutely required for any deviation from Mendelian segregation to occur. A candidate for the Tcr gene has previously been identified. Evidence that this gene represents Tcr is its localization to the appropriate genomic subregion and testis-specific expression pattern. Here, we report the molecular cloning of the region between recombinant chromosome breakpoints defining the Tcr locus. These results circumscribe Tcr to a 150- to 220-kb region of DNA, including the 22-kb candidate responder gene. This gene and two other homologs were created by large genomic duplications, each involving segments of DNA 10-fold larger than the individual genes.     

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.     

Chromosomal elements regulate gene activity and chromatin structure of the human serpin gene cluster at 14q32.1

The human serine protease inhibitor (serpin) gene cluster at 14q32.1 contains a number of genes that are specifically expressed in hepatic cells. Cell-specific enhancers have been identified in several of these genes, but elements involved in locus-wide gene and chromatin control have yet to be defined. To identify regulatory elements in this region, we prepared a series of mutant chromosomal alleles by homologous recombination and transferred the specifically modified human chromosomes to hepatic cells for functional tests. We report that deletion of an 8-kb DNA segment upstream of the human alpha1-antitrypsin gene yields a mutant serpin allele that fails to be activated in hepatic cells. Within this region, a 2.3-kb DNA segment between kb -8.1 and -5.8 contains a previously unrecognized control region that is required not only for serpin gene activation but also for chromatin remodeling of the entire locus.