Patterns of Sequence Divergence and Evolution of the S1 Orthologous Regions between Asian and African Cultivated Rice Species

A strong postzygotic reproductive barrier separates the recently diverged Asian and African cultivated rice species, Oryza sativa and O. glaberrima. Recently a model of genetic incompatibilities between three adjacent loci: S₁A, S₁ and S₁B (called together the S₁ regions) interacting epistatically, was postulated to cause the allelic elimination of female gametes in interspecific hybrids. Two candidate factors for the S₁ locus (including a putative F-box gene) were proposed, but candidates for S₁A and S₁B remained undetermined. Here, to better understand the basis of the evolution of regions involved in reproductive isolation, we studied the genic and structural changes accumulated in the S₁ regions between orthologous sequences. First, we established an 813 kb genomic sequence in O. glaberrima, covering completely the S₁A, S₁ and the majority of the S₁B regions, and compared it with the orthologous regions of O. sativa. An overall strong structural conservation was observed, with the exception of three isolated regions of disturbed collinearity: (1) a local invasion of transposable elements around a putative F-box gene within S₁, (2) the multiple duplication and subsequent divergence of the same F-box gene within S₁A, (3) an interspecific chromosomal inversion in S₁B, which restricts recombination in our O. sativa×O. glaberrima crosses. Beside these few structural variations, a uniform conservative pattern of coding sequence divergence was found all along the S₁ regions. Hence, the S₁ regions have undergone no drastic variation in their recent divergence and evolution between O. sativa and O. glaberrima, suggesting that a small accumulation of genic changes, following a Bateson-Dobzhansky-Muller (BDM) model, might be involved in the establishment of the sterility barrier. In this context, genetic incompatibilities involving the duplicated F-box genes as putative candidates, and a possible strengthening step involving the chromosomal inversion might participate to the reproductive barrier between Asian and African rice species.     

QTL detection for rice grain quality traits using an interspecific backcross population derived from cultivated Asian (O. sativa L.) and African (O. glaberrima S.) rice.

An interspecific advanced backcross population derived from a cross between Oryza sativa "V20A" (a popular male-sterile line used in Chinese rice hybrids) and Oryza glaberrima (accession IRGC No. 103544 from Mali) was used to identify quantitative trait loci (QTL) associated with grain quality and grain morphology. A total of 308 BC3F1 hybrid families were evaluated for 16 grain-related traits under field conditions in Changsha, China, and the same families were evaluated for RFLP and SSR marker segregation at Cornell University (Ithaca, N.Y.). Eleven QTL associated with seven traits were detected in six chromosomal regions, with the favorable allele coming from O. glaberrima at eight loci. Favorable O. glaberrima alleles were associated with improvements in grain shape and appearance, resulting in an increase in kernel length, transgressive variation for thinner grains, and increased length to width ratio. Oryza glaberrima alleles at other loci were associated with potential improvements in crude protein content and brown rice yield. These results suggested that genes from O. glaberrima may be useful in improving specific grain quality characteristics in high-yielding O. sativa hybrid cultivars.     

Conservation and purifying selection of transcribed genes located in a rice centromere

Recombination is strongly suppressed in centromeric regions. In chromosomal regions with suppressed recombination, deleterious mutations can easily accumulate and cause degeneration of genes and genomes. Surprisingly, the centromere of chromosome8 (Cen8) of rice (Oryza sativa) contains several transcribed genes. However, it remains unclear as to what selective forces drive the evolution and existence of transcribed genes in Cen8. Sequencing of orthologous Cen8 regions from two additional Oryza species, Oryza glaberrima and Oryza brachyantha, which diverged from O. sativa 1 and 10 million years ago, respectively, revealed a set of seven transcribed Cen8 genes conserved across all three species. Chromatin immunoprecipitation analysis with the centromere-specific histone CENH3 confirmed that the sequenced orthologous regions are part of the functional centromere. All seven Cen8 genes have undergone purifying selection, representing a striking phenomenon of active gene survival within a recombination-free zone over a long evolutionary time. The coding sequences of the Cen8 genes showed sequence divergence and mutation rates that were significantly reduced from those of genes located on the chromosome arms. This suggests that Oryza has a mechanism to maintain the fidelity and functionality of Cen8 genes, even when embedded in a sea of repetitive sequences and transposable elements.     

Analysis of collinear regions of Oryza AA and CC genomes

Comparative analyses of genome structure and sequence of closely related species have yielded insights into the evolution and function of plant genomes. A total of 103,844 BAC end sequences delegated -73.8 Mb of O. officinalis that belongs to the CC genome type of the rice genus Oryza were obtained and compared with the genome sequences office cultivar, O. sativa ssp.japonica cv. Nipponbare. We found that more than 45% of O. officinalis genome consists of repeat sequences, which is higher than that of Nipponbare cultivar. To further investigate the evolutionary divergence of AA and CC genomes, two BAC-contigs of O. officinalis were compared with the collinear genomic regions of Nipponbare. Of 57 genes predicted in the AA genome orthologous regions, 39 had orthologs in the regions of the CC genome. Alignment of the orthologous regions indicated that the CC genome has undergone expansion in both genic and intergenic regions through primarily retroelement insertion. Particularly, the density of RNA transposable elements was 17.95% and 1.78% in O. officinalis and O. sativa, respectively. This explains why the orthologous region is about 100 kb longer in the CC genome in comparison to the AA genome.     

RAPD, RFLP and SSLP analyses of phylogenetic relationships between cultivated and wild species of rice.

RAPD, RFLP, nuclear SSLP and chloroplast SSLP analyses were carried out to clarify the phylogenetic relationships among A-genome species of rice. In total, 12 cultivars of Oryza sativa (4 Japonica, 3 Javanica and 5 Indica), one cultivar of O. glaberrima, and 17 wild accessions (12 O. rufipogon, 2 O. glumaepatula, 1 O. longistaminata, 1 O. meridionalis and 1 O. barthii) were used. Their banding patterns were scored and compared to evaluate the similarity between accessions. Genetic differentiation within and between taxa was examined based on the average similarity indices. Except for chloroplast SSLP analysis, the average similarities were higher within O. sativa than within O. rufipogon, and O. sativa Indica had greater intrasubspecific variation than Japonica and Javanica. Comparisons between cultivated and wild species showed that O. sativa was closely related to O. rufipogon, while O. glaberrima was closely related to O. barthii. This indicated that two cultivated species, O. sativa and O. glaberrima, originated from O. rufipogon and O. barthii, respectively. Domestication of O. sativa seemed to be diphyletic, since strong similarity was observed between O. sativa Japonica-Javanica and O. rufipogon from China and between O. sativa Indica and O. rufipogon from tropical Asia. In addition, dendrograms for RAPD, RFLP, and nuclear and chloroplast SSLP analyses were constructed to reveal the overall genetic relationships among A-genome species. In all analyses, O. sativa and O. glaberrima formed groups with O. rufipogon and O. barthii, respectively. However, their manners of clustering with other wild species were not the same. The results of RAPD and RFLP analyses indicate that O. glumaepatula was relatively close to the groups of O. sativa and O. glaberrima whereas O. longistaminata and O. meridionalis were highly differentiated from other A-genome species. On the other hand, clear interspecific relationships were not obtained by nuclear or chloroplast SSLP analyses.     

Sex-independent transmission ratio distortion system responsible for reproductive barriers between Asian and African rice species

* A sex-independent transmission ratio distortion (siTRD) system detected in the interspecific cross in rice was analyzed in order to understand its significance in reproductive barriers. The S(1) gene, derived from African rice Oryza glaberrima, induced preferential abortion of both male and female gametes possessing its allelic alternative (), from Asian rice O. sativa, only in the heterozygote. * The siTRD was characterized by resolving it into mTRD and fTRD occurring through male and female gametes, respectively, cytological analysis of gametophyte development, and mapping of the S(1) locus using near-isogenic lines. The allelic distribution of the S(1) locus in Asian and African rice species complexes was also analyzed. * The siTRD system involved at least two components affecting male and female gametogeneses, respectively, including a modifier(s) that enhances fTRD. The chromosomal location of the major component causing the mTRD was delimited within an approx. 40 kb region. The S(1) locus induced hybrid sterility in any pairwise combination between Asian and African rice species complexes. * The allelic state of the S(1) locus has diverged between Asian and African rice species complexes, suggesting that the TRD system has a significant role in the reproductive barriers in rice.     

Diversification of the rice Waxy gene by insertion of mobile DNA elements into introns.

The waxy (wx) gene of Oryza glaberrima was cloned, and its nucleotide sequence was determined. A waxy mutant of O. glaberrima showing a glutinous phenotype was found to contain a substitution mutation generating a termination codon in the coding region of the wx gene. The Wx sequence of O. glaberrima was different from that of Oryza sativa by substitutions and insertions/deletions, among which only a few substitutions occurred in several exons not to severely alter the amino acid sequence of the Wx protein. The most striking difference observed in introns was a 139-bp deletion (or insertion) in intron 10 of O. glaberrima (or O. sativa). In O. sativa, 125 bp of the 139-bp sequence was flanked by direct repeats of a 14-bp sequence. A sequence homologous to the 125-bp sequence was found in the region preceding exon 2; this sequence was also flanked by direct repeats of another 14-bp sequence. This result and the observation that the 125-bp sequence was interspersed in rice genomes indicate that they are SINEs (short interspersed elements) in the plant system. We also identified a DNA sequence with long terminal inverted repeats in intron 13 of both O. glaberrima and O. sativa. This sequence was present in multiple copies in rice genomes, suggesting that it is a transposable element. These results obtained suggest that mobile DNA elements have diversified the rice Waxy gene by inserting into introns, each of which may originally have a length of about 100 bp.     

Rate and polarity of gene fusion and fission in Oryza sativa and Arabidopsis thaliana

Eukaryotic gene fusion and fission events are mechanistically more complicated than in prokaryotes, and their quantitative contributions to genome evolution are still poorly understood. We have identified all differentially composite or split genes in 2 fully sequenced plant genomes, Oryza sativa and Arabidopsis thaliana. Out of 10,172 orthologous gene pairs, 60 (0.6% of the total) revealed a verified fusion or fission event in either lineage after the divergence of O. sativa and A. thaliana. Polarizing these events by outgroup comparison revealed differences in the rate of gene fission but not of gene fusion in the rice and Arabidopsis lineages. Gene fission occurred at a higher rate than gene fusion in the O. sativa lineage and was furthermore more common in rice than in Arabidopsis. Nucleotide insertion bias has promoted gene fission in the O. sativa lineage, consistent with its generally longer nucleotide sequences than A. thaliana in selectively neutral regions, and with the abundance of transposable elements in rice. The divergence time of monocots and dicots (140-200 Myr) indicates that gene fusion/fission events occur at an average rate of 1x10(-11) to 2x10(-11) events per gene per year, approximately 100-fold slower than the average per site nuclear nucleotide substitution rate in these lineages. Gene fusion and fission are thus rare and slow processes in higher plant genomes; they should be of utility to address deeper evolutionary relationships among plants--and the relationship of plants to other eukaryotic lineages--where sequence-based phylogenies provide equivocal or conflicting results.     

The Genic Nature of Gamete Eliminator in Rice

The two cultivated rice species, Oryza sativa and Oryza glaberrima, are morphologically alike but are reproductively isolated from each other by hybrid sterility. The hybrid is male sterile but partially female fertile. Backcross experiments were conducted to introduce an alien factor controlling hybrid sterility from O. glaberrima (W025) into O. sativa (T65wx) and examine the genetic basis. An extracted sterility factor, closely linked to the wx locus, induced gametic abortion due to allelic interaction and was tentatively designated as S(t). The segregation patterns for infertility was explained by assuming that W025 and T65wx carried S(t) and S(t)a, respectively, and gametes with S(t)a aborted only in the heterozygote (S(t)/S(t)a) although the elimination of female gametes was incomplete. Thus, S(t) seemed to be intermediate between a gamete eliminator and pollen killer. However, S(t) was proven to be likely the same as S1 which was formerly reported as gamete eliminator in a different genetic background of O. sativa. In addition, a chromosomal segment containing S1 (or S(t] caused a marked suppression of crossing over around it, suggesting the presence of an inversion. Further, female transmission of S1a increased as the segment containing S1 became small by recombination. After S1 was further purified by successive backcrosses up to the BC15 generation, it became pollen killer. The present results give evidence that a profound sterility gene such as gamete eliminator can be made from accumulation of pollen killer and its modifier(s) when pollen killer and modifier(s) are linked, they behave as a gene complex in the hybrid.     

Evolutionary relationships among rice species with AA genome based on SINE insertion analysis

Previous studies based on morphological and molecular markers indicated that there are two cultivated and five wild rice species within the Oryza genus with the AA genome. In the cultivated rice species, Oryza sativa, a retroposon named p-SINE1 has been identified. Some of the p-SINE1 members characterized previously showed interspecific insertion polymorphisms in the species with the AA genome. In this study, we identified new p-SINE1 members showing interspecific insertion polymorphisms from representative strains of four wild rice species with the AA genome: O. barthii, O. glumaepatula, O. longistaminata, and O. meridionalis. Some of these members were present only in strains of one species, whereas the others were present in strains of two or more species. The p-SINE1 insertion patterns in the strains of the Asian and African cultivated rice species O. sativa and O. glaberrima were very similar to those of the Asian and African wild rice species O. rufipogon and O. barthii, respectively. This is consistent with the previous hypothesis that O. sativa and O. glaberrima are derived from specific wild rice species. Phylogenetic analysis based on the p-SINE1 insertion patterns showed that the strains of each of the five wild rice species formed a cluster. The strains of O. longistaminata appear to be distantly related to those of O. meridionalis. The strains of these two species appear to be distantly related to those of three other species, O. rufipogon, O. barthii and O. glumaepatula. The latter three species are closely related to one another with O. barthii and O. glumaepatula being most closely related. A phylogenetic tree including a hypothetical ancestor with all loci empty for p-SINE1 insertion showed that the strains of O. longistaminata are related most closely to the hypothetical ancestor. This indicates that O. longistaminata and O. meridionalis diverged early on, whereas the other species diverged relatively recently, and suggests that the Oryza genus with AA genome might have originated in Africa, rather than in Asia.     

The population structure of African cultivated rice oryza glaberrima (Steud.): evidence for elevated levels of linkage disequilibrium caused by admixture with O. sativa and ecological adaptation

Genome-wide linkage disequilibrium (LD) was investigated for 198 accessions of Oryza glaberrima using 93 nuclear microsatellite markers. Significantly elevated levels of LD were detected, even among distantly located markers. Free recombination among loci at the population genetic level was shown (1) by a lack of decay in LD among markers on the same chromosome and (2) by a strictly increasing composite likelihood function for the recombination parameter. This suggested that the elevation in LD was due not to physical linkage but to other factors, such as population structure. A Bayesian clustering analysis confirmed this hypothesis, indicating that the sample of O. glaberrima in this study was subdivided into at least five cryptic subpopulations. Two of these subpopulations clustered with control samples of O. sativa, subspecies indica and japonica, indicating that some O. glaberrima accessions represent admixtures. The remaining three O. glaberrima subpopulations were significantly associated with specific combinations of phenotypic traits-possibly reflecting ecological adaptation to different growing environments.     

Long-range and targeted ectopic recombination between the two homeologous chromosomes 11 and 12 in Oryza species

Whole genome duplication (WGD) and subsequent evolution of gene pairs have been shown to have shaped the present day genomes of most, if not all, plants and to have played an essential role in the evolution of many eukaryotic genomes. Analysis of the rice (Oryza sativa ssp. japonica) genome sequence suggested an ancestral WGD ～50-70 Ma common to all cereals and a segmental duplication between chromosomes 11 and 12 as recently as 5 Ma. More recent studies based on coding sequences have demonstrated that gene conversion is responsible for the high sequence conservation which suggested such a recent duplication. We previously showed that gene conversion has been a recurrent process throughout the Oryza genus and in closely related species and that orthologous duplicated regions are also highly conserved in other cereal genomes. We have extended these studies to compare megabase regions of genomic (coding and noncoding) sequences between two cultivated (O. sativa, Oryza glaberrima) and one wild (Oryza brachyantha) rice species using a novel approach of topological incongruency. The high levels of intraspecies conservation of both gene and nongene sequences, particularly in O. brachyantha, indicate long-range conversion events less than 4 Ma in all three species. These observations demonstrate megabase-scale conversion initiated within a highly rearranged region located at ～2.1 Mb from the chromosome termini and emphasize the importance of gene conversion in cereal genome evolution.     

Evolutionary analysis of the Sub1 gene cluster that confers submergence tolerance to domesticated rice

BACKGROUND AND AIMS: Tolerance of complete submergence is recognized in a small number of accessions of domesticated Asian rice (Oryza sativa) and can be conferred by the Sub1A-1 gene of the polygenic Submergence-1 (Sub1) locus. In all O. sativa varieties, the Sub1 locus encodes the ethylene-responsive factor (ERF) genes Sub1B and Sub1C. A third paralogous ERF gene, Sub1A, is limited to a subset of indica accessions. It is thought that O. sativa was domesticated from the gene pools of the wild perennial species O. rufipogon Griff. and/or the annual species O. nivara Sharma et Shastry. The aim of this study was to evaluate the orthologues of the Sub1 locus in the closest relatives of O. sativa to provide insight into the origin of the gene and allelic variation of the Sub1 locus. METHODS: Orthologues of the Sub1 genes were isolated from O. rufipogon and O. nivara by use of oligonucleotide primers corresponding to the most highly conserved regions of the Sub1 genes of domesticated rice. The phylogenetic relatedness of Sub1 genes of O. sativa and its wild relatives was evaluated. KEY RESULTS AND CONCLUSIONS: Both O. rufipogon and O. nivara possess two Sub1 gene orthologues with strong sequence identity to the Sub1B and Sub1C alleles of cultivated rice. The phylogeny of the Sub1 genes of the domesticated and wild rice suggests that Sub1A arose from duplication of Sub1B. Variation in Sub1B alleles is correlated with the absence or presence of Sub1A. Together, the results indicate that genetic variation at the Sub1 locus is due to gene duplication and divergence that have occurred both prior to and after rice domestication.     

Fine genetic mapping of a gene required for <Emphasis Type="Italic">Rice yellow mottle virus</Emphasis> cell-to-cell movement

The very high resistance to Rice yellow mottle virus observed in the two rice varieties Gigante ( Oryza sativa) and Tog 5681 ( O. glaberrima) is monogenic and recessive. Bulked segregant analysis was carried out to identify AFLP markers linked to the resistance gene. Mapping of PCR-specific markers, CAPS and microsatellite markers on 429 individuals of an IR64 x Gigante F(2) population pinpointed this resistance gene on the long arm of chromosome 4 in a 3.7-cM interval spanned by PCR markers. These markers also flanked the resistance gene of the O. glaberrima accession Tog 5681 and confirmed previous allelism tests. The rarity of this recessive natural resistance was in line with a resistance mechanism model based on point mutations of a host component required for cell-to-cell movement of the virus. Preliminary data on the genetic divergence between the two cultivated rice species in the vicinity of the resistance locus suggested that two different resistance alleles are present in Gigante and Tog 5681. A large set of recombinants is now available to envisage physical mapping and cloning of the gene.     

Variation in the loss of seed dormancy during after-ripening of wild and cultivated rice species

BACKGROUND AND AIMS: The aim of this paper was to verify the variation in the loss of seed dormancy during after-ripening and the interspecific and interpopulation variability in the degree of dormancy of seven wild and two cultivated rice species comprising 21 populations and two cultivars. METHODS: Four wild rice species from South America, Oryza glumaepatula, O. latifolia, O. grandiglumis and O. alta, and two O. sativa cultivars were tested in one experiment. In a second experiment, five wild species, O. punctata, O. eichingeri, O.rufipogon, O. latifolia and O. glumaepatula, and one cultivated species (O. glaberrima) were evaluated. Initial germination tests were performed soon after the seeds were harvested and subsequently at 2-month intervals, for a total of six storage periods in the first experiment and three in the second. All tests were conducted in the dark at a temperature of 27 degrees C. KEY RESULTS: Different patterns of after-ripening among populations within and between species were observed. CONCLUSIONS: The cultivated species (O. sativa and O. glaberrima) and, amongst the wild species, the tetraploids O. latifolia, O. grandiglumis and the diploids O. eichingeri and O. punctata, had weak dormancy, losing it completely 2 months after harvest, while O. rufipogon and O. glumaepatula exhibited pronounced dormancy. The latter showed different patterns of after-ripening between populations indigenous to the Amazon region and those originating in the Paraguay River system. Seeds of Solimoes (Amazon) and Japura origin showed weak dormancy whereas those of Paraguay origin showed deep dormancy. Ecological differences among natural habitats may be involved in such differentiation.     

The plastid chromosome of Atropa belladonna and its comparison with that of Nicotiana tabacum: the role of RNA editing in generating divergence in the process of plant speciation

The nuclear and plastid genomes of the plant cell form a coevolving unit which in interspecific combinations can lead to genetic incompatibility of compartments even between closely related taxa. This phenomenon has been observed for instance in Atropa-Nicotiana cybrids. We have sequenced the plastid chromosome of Atropa belladonna (deadly nightshade), a circular DNA molecule of 156,688 bp, and compared it with the corresponding published sequence of its relative Nicotiana tabacum (tobacco) to understand how divergence at the level of this genome can contribute to nuclear-plastid incompatibilities and to speciation. It appears that (1) regulatory elements, i.e., promoters as well as translational and replicational signal elements, are well conserved between the two species; (2) genes--including introns--are even more highly conserved, with differences residing predominantly in regions of low functional importance; and (3) RNA editotypes differ between the two species, which makes this process an intriguing candidate for causing rapid reproductive isolation of populations.     

Dynamics of the evolution of orthologous and paralogous portions of a complex locus region in two genomes of allopolyploid wheat

Two overlapping bacterial artificial chromosome (BAC) clones from the B genome of the tetraploid wheat Triticum turgidum were identified, each of which contains one of the two high-molecular-weight (HMW) glutenin genes, comprising the complex Glu-B1 locus. The complete sequence (285 506 bp of DNA) of this chromosomal region was determined. The two paralogous x-type ( Glu-1-1 ) and y-type ( Glu-1-2 ) HMW-glutenin genes of the complex Glu-B1 locus were found to be separated by ca. 168 000 bp instead of the 51 000 bp separation previously reported for the orthologous Glu-D1 locus of Aegilops tauschii, the D-genome donor of hexaploid wheat. This difference in intergene spacing is due almost entirely to be the insertion of clusters of nested retrotransposons. Otherwise, the orientation and order of the HMW glutenins and adjacent genes were identical in the two genomes. A comparison of these orthologous regions indicates modes and patterns of sequence divergence, with implications for the overall Triticeae genome structure and evolution. A duplicate globulin gene, found 5' of each HMW-glutenin gene, assists to tentatively define the original duplication event leading to the paralogous x- and y-type HMW-glutenin genes. The intergenic regions of the two loci are composed of different patterns and classes of retrotransposons, indicating that insertion times of these retroelements were after the divergence of the two wheat genomes. In addition, a putative receptor kinase gene near the y-type HMW-glutenin gene at the Glu-B1 locus is likely active as it matches recently reported ESTs from germinating barley endosperm. The presence of four genes represented only in the Triticeae endosperm ESTs suggests an endosperm-specific chromosome domain.     

The genomic bases of morphological divergence and reproductive isolation driven by ecological speciation in Senecio (Asteraceae)

Ecological speciation, driven by adaptation to contrasting environments, provides an attractive opportunity to study the formation of distinct species, and the role of selection and genomic divergence in this process. Here, we focus on a particularly clear‐cut case of ecological speciation to reveal the genomic bases of reproductive isolation and morphological differences between closely related Senecio species, whose recent divergence within the last ~200 000 years was likely driven by the uplift of Mt. Etna (Sicily). These species form a hybrid zone, yet remain morphologically and ecologically distinct, despite active gene exchange. Here, we report a high‐density genetic map of the Senecio genome and map hybrid breakdown to one large and several small quantitative trait loci (QTL). Loci under diversifying selection cluster in three 5 cM regions which are characterized by a significant increase in relative (F_ST), but not absolute (d_XY), interspecific differentiation. They also correspond to some of the regions of greatest marker density, possibly corresponding to ‘cold‐spots’ of recombination, such as centromeres or chromosomal inversions. Morphological QTL for leaf and floral traits overlap these clusters. We also detected three genomic regions with significant transmission ratio distortion (TRD), possibly indicating accumulation of intrinsic genetic incompatibilities between these recently diverged species. One of the TRD regions overlapped with a cluster of high species differentiation, and another overlaps the large QTL for hybrid breakdown, indicating that divergence of these species may have occurred due to a complex interplay of ecological divergence and accumulation of intrinsic genetic incompatibilities.     

Maintenance of genetic and morphological identity in two sibling Syrrhopodon species (Calymperaceae,Bryopsida) despite extensive introgression

Bryophytes are a group of land plants in which the role of hybridization has long been challenged. Using genotyping by sequencing to circumvent the lack of molecular variation at selected loci previously used for phylogeny and morphology, we determine the level of genetic and morphological divergence and reproductive isolation between the sibling Syrrhopodon annotinus and S. simmondsii (Calymperaceae, Bryopsida) that occur in sympatry but in different habitats in lowland Amazonian rainforests. A clear morphological differentiation and a low (0.06), but significant Fst derived from the analysis of 183 single nucleotide polymorphisms were observed between the two species. Conspecific pairs of individuals consistently exhibited higher average kinship coefficients along a gradient of geographic isolation than interspecific pairs. The weak, but significant genetic divergence observed is consistent with growing evidence that ecological specialization can lead to genetic differentiation among bryophyte species. Nevertheless, the spatial genetic structures of the two species were significantly correlated, as evidenced by the significant slope of the Mantel test based on kinship coefficients between pairs of interspecific individuals and the geographic distance separating them. Interspecific pairs of individuals are thus more closely related when they are geographically closer, suggesting that isolation‐by‐distance is stronger than the interspecific reproductive barrier and pointing to interspecific gene flow. We conclude that interspecific introgression, whose role has long been questioned in bryophytes, may take place even in species wherein sporophyte production is scarce due to dioicy, raising the question as to what mechanisms maintain differentiation despite weak reproductive isolation.     

The beta globin gene cluster of the prosimian primate Galago crassicaudatus: nucleotide sequence determination of the 41-kb cluster and comparative sequence analyses.

The nucleotide sequence of the beta globin gene cluster of the prosimian Galago crassicaudatus has been determined. A total sequence spanning 41,101 bp contains and links together previously published sequences of the five galago beta-like globin genes (5'-epsilon-gamma-psi eta-delta-beta-3'). A computer-aided search for middle interspersed repetitive sequences identified 10 LINE (L1) elements, including a 5' truncated repeat that is orthologous to the full-length L1 element found in the human epsilon-gamma intergenic region. SINE elements that were identified included one Alu type I repeat, four Alu type II repeats, and two methionine tRNA-derived Monomer (type III) elements. Alu type II and Monomer sequences are unique to the galago genome. Structural analyses of the cluster sequence reveals that it is relatively A+T rich (about 62%) and regions with high G+C content are associated primarily with globin coding regions. Comparative analyses with the beta globin cluster sequences of human, rabbit, and mouse reveal extensive sequence homologies in their genic regions, but only human, galago, and rabbit sequences share extensive intergenic sequence homologies. Divergence analyses of aligned intergenic and flanking sequences from orthologous human, galago, and rabbit sequences show a gradation in the rate of nucleotide sequence evolution along the cluster where sequences 5' of the epsilon globin gene region show the least sequence divergence and sequences just 5' of the beta globin gene region show the greatest sequence divergence.