Intraspecific variation in reproductive barriers between two closely related Arabidopsis sister species

Reproductive isolation (RI) is a critical component of speciation and varies strongly in timing and strength among different sister taxa, depending on, for example the geography of speciation and divergence time. However, these factors may also produce variation in timing and strength among populations within species. Here we tested for variation in the expression of RI among replicate population pairs between the sister taxa Arabidopsis lyrata subsp. lyrata and A. arenicola. While the former is predominantly outcrossing, the latter is predominantly selfing. We focused on intrinsic prezygotic and postzygotic RI as both species occur largely in allopatry. We assessed RI by performing within-population crosses and interspecific between-population crosses, and by raising offspring. RI was generally high between all interspecific population pairs, but it varied in timing and strength depending on population history. Prezygotic isolation was strongest between the closest-related population pair, while early postzygotic isolation was high for all other population pairs. Furthermore, the timing and strength of RI depended strongly on cross direction. Our study provides empirical support that reproductive barriers between species are highly variable among population pairs and asymmetric within population pairs, and this variation seems to follow patterns typically described across species pairs.     

Effects of gene expression on molecular evolution in Arabidopsis thaliana and Arabidopsis lyrata

We analyzed the complete genome sequence of Arabidopsis thaliana and sequence data from 83 genes in the outcrossing A. lyrata, to better understand the role of gene expression on the strength of natural selection on synonymous and replacement sites in Arabidopsis. From data on tRNA gene abundance, we find a good concordance between codon preferences and the relative abundance of isoaccepting tRNAs in the complete A. thaliana genome, consistent with models of translational selection. Both EST-based and new quantitative measures of gene expression (MPSS) suggest that codon preferences derived from information on tRNA abundance are more strongly associated with gene expression than those obtained from multivariate analysis, which provides further support for the hypothesis that codon bias in Arabidopsis is under selection mediated by tRNA abundance. Consistent with previous results, analysis of protein evolution reveals a significant correlation between gene expression level and amino acid substitution rate. Analysis by MPSS estimates of gene expression suggests that this effect is primarily the result of a correlation between the number of tissues in which a gene is expressed and the rate of amino acid substitution, which indicates that the degree of tissue specialization may be an important determinant of the rate of protein evolution in Arabidopsis.     

Gene conversion and the evolution of three leucine-rich repeat gene families in Arabidopsis thaliana

The high number of duplicated genes in plant genomes provides a potential template for gene conversion and unequal crossing-over. Within a gene family these two processes can render all members homogeneous or generate diversity by reassorting variants among paralogs. The latter is especially feasible in families where gene diversity confers a selective advantage and thus conversion events are likely to be retained. Consequently, the most complete record of gene conversion is expected to be most evident in gene families commonly subjected to positive selection. Here, we describe the extent and characteristics of gene conversion and unequal crossing-over in the coding and noncoding regions of nucleotide-binding site leucine-rich repeat (NBS-LRR), receptor-like kinases (RLK), and receptor-like proteins (RLP) in the plant Arabidopsis thaliana. Members of these three gene families are associated with disease resistance and their pathogen-recognition domain is a documented target of positive selection. Our bioinformatic approach to study the major family features that may influence gene conversion revealed that in these families there is a significant association between the occurrence of gene conversion and high levels of sequence similarity, close physical clustering, gene orientation, and recombination rate. We discuss these results in the context of the overlap between gene conversion and positive selection during the evolutionary expansion of the NBS-LRR, RLK, and RLP gene families.     

Gene family evolution in green plants with emphasis on the origination and evolution of Arabidopsis thaliana genes

Gene family size variation is an important mechanism that shapes the natural variation for adaptation in various species. Despite its importance, the pattern of gene family size variation in green plants is still not well understood. In particular, the evolutionary pattern of genes and gene families remains unknown in the model plant Arabidopsis thaliana in the context of green plants. In this study, eight representative genomes of green plants are sampled to study gene family evolution and characterize the origination of A. thaliana genes, respectively. Four important insights gained are that: (i) the rate of gene gains and losses is about 0.001359 per gene every million years, similar to the rate in yeast, Drosophila, and mammals; (ii) some gene families evolved rapidly with extreme expansions or contractions, and 2745 gene families present in all the eight species represent the ‘core’ proteome of green plants; (iii) 70% of A. thaliana genes could be traced back to 450 million years ago; and (iv) intriguingly, A. thaliana genes with early origination are under stronger purifying selection and more conserved. In summary, the present study provides genome‐wide insights into evolutionary history and mechanisms of genes and gene families in green plants and especially in A. thaliana.     

Independent origins of self-compatibility in Arabidopsis thaliana

The evolution from outcrossing based on self-incompatibility (SI) to a selfing system is one of the most prevalent transitions in flowering plants. It has been suggested that the loss of SI in Arabidopsis thaliana is associated with pseudogene formation at the SCR male component of the S locus. Recent work, however, suggests that alternative alleles with large deletions at the S locus are also present and may be responsible for the evolution of self-compatibility in this species. We demonstrate that most of these deletion alleles are evolutionarily derived from an S haplotype (haplogroups A) that already possessed the SCR pseudogene. This haplotype and its deletion variants are nearly fixed in Europe. Together with previous transgenic data, these results suggest that the pseudogenization of PsiSCR1 gene changed the SI phenotype in the majority of A. thaliana accessions, and was a critical step in the evolution of selfing in this species. Two other haplogroups (B and C) were also identified, the former of which contains a novel and possibly functional SCR allele. In contrast to haplogroups A, these two haplogroups are found primarily in Africa and Asia. These results suggest that self-compatibility, which appears to be fixed in this species, arose multiple times with different genetic bases, and indicates that a species-specific trait is associated with parallel evolution at the molecular level.     

Nad+-dependent isocitrate dehydrogenase from Arabidopsis thaliana. Characterization of two closely related subunits

Two cDNA clones which appear to encode different subunits of NAD⁺-dependent isocitrate dehydrogenase (IDH; EC 1.1.1.41) were identified by homology searches from the Arabidopsis EST database. These cDNA clones were obtained and sequenced; both encoded full-length messages and displayed 82.7% nucleotide sequence identity over the coding region. The deduced amino acid sequences revealed preprotein lengths of 367 residues, with an amino acid identity of 86.1%. Genomic Southern blot analysis showed distinct single-copy genes for both IDH subunits. Both IDH subunits were expressed as recombinant proteins in Escherichia coli, and polyclonal antibodies were raised to each subunit. The Arabidopsis cDNA clones were expressed in Saccharomyces cerevisiae mutants which were deficient in either one or both of the yeast NAD⁺-dependent IDH subunits. The Arabidopsis cDNA clones failed to complement the yeast mutations; although both IDH-I and IDH-II were expressed at detectable levels, neither protein was imported into the mitochondria.     

Establishing gene function by mutagenesis in Arabidopsis thaliana

The nuclear genome of Arabidopsis thaliana was sequenced to near completion a few years ago, and ahead lies the challenge of understanding its meaning and discerning its potential. How many genes are there? What are they? What do they do? Computer algorithms combined with genome array technologies have proven efficient in addressing the first two questions as shown in a recent report ( Yamada et al., 2003 ). However, assessing the function of every gene in every cell will require years of careful analyses of the phenotypes caused by mutations in each gene. Current progress in generating large numbers of molecular markers and near‐saturation insertion mutant collections has immensely facilitated functional genomics studies in Arabidopsis. In this review, we focus on how gene function can be revealed through the analysis of mutants by either forward or reverse genetics. These mutants generally fall into two distinct classes. The first class typically includes point mutations or small deletions derived from chemical or fast neutron mutagenesis whereas the second class includes insertions of transferred‐DNA or transposon elements. We describe the current methods that are used to identify the gene corresponding to these mutations, which can then be used as a probe to further dissect its function.     

Expression of the Arabidopsis thaliana invertase gene family

Two new loci have been found to be clustered with five other genes for the nitrate assimilation pathway in the Chlamydomonas reinhardtii genome. One gene, located close to the 3′-end of the high-affinity nitrate transporter (HANT) gene Nrt2;2, corresponds to the nitrite reductase (NiR) structural gene Nii1. This is supported by a number of experimental findings: (i) NiR-deficient mutants have lost Nii1 gene expression; (ii) Nii1 mRNA accumulation is co-regulated with the expression of other structural genes of the nitrate assimilation pathway; (iii) nitrite (nitrate) utilization ability is recovered in the NiR mutants by functional complementation with a wild-type Nii1 gene; (iv) the elucidated NII1 amino acid sequence is highly similar to that of the cyanobacterial and higher-plant enzyme, and contains the predicted domains for plastidic ferredoxin-NiRs. Thus, the mutant phenotype and the mRNA sequence and expression of the Nii1 gene have been unequivocally related. Accumulation of mRNA for the second locus identified, Lde1 (light-dependent expression), was not regulated by nitrogen, but like nitrate-assimilation clustered genes, its expression was down-regulated in the dark. Received: 27 November 1997 / Accepted: 19 January 1998     

Molecular analysis of the myosin gene family in Arabidopsis thaliana

Myosin is believed to act as the molecular motor for many actin-based motility processes in eukaryotes. It is becoming apparent that a single species may possess multiple myosin isoforms, and at least seven distinct classes of myosin have been identified from studies of animals, fungi, and protozoans. The complexity of the myosin heavy-chain gene family in higher plants was investigated by isolating and characterizing myosin genomic and cDNA clones from Arabidopsis thaliana. Six myosin-like genes were identified from three polymerase chain reaction (PCR) products (PCR1, PCR11, PCR43) and three cDNA clones (ATM2, MYA2, MYA3). Sequence comparisons of the deduced head domains suggest that these myosins are members of two major classes. Analysis of the overall structure of the ATM2 and MYA2 myosins shows that they are similar to the previously-identified ATM1 and MYA1 myosins, respectively. The MYA3 appears to possess a novel tail domain, with five IQ repeats, a six-member imperfect repeat, and a segment of unique sequence. Northern blot analyses indicate that some of the Arabidopsis myosin genes are preferentially expressed in different plant organs. Combined with previous studies, these results show that the Arabidopsis genome contains at least eight myosin-like genes representing two distinct classes.     

Solution structure of thioredoxin h1 from Arabidopsis thaliana

Present in virtually every species, thioredoxins catalyze disulfide/dithiol exchange with various substrate proteins. While the human genome contains a single thioredoxin gene, plant thioredoxins are a complex protein family. A total of 19 different thioredoxin genes in six subfamilies has emerged from analysis of the Arabidopsis thaliana genome. Some function specifically in mitochondrial and chloroplast redox signaling processes, but target substrates for a group of eight thioredoxin proteins comprising the h subfamily are largely uncharacterized. In the course of a structural genomics effort directed at the recently completed A. thaliana genome, we determined the structure of thioredoxin h1 (At3g51030.1) in the oxidized state. The structure, defined by 1637 NMR-derived distance and torsion angle constraints, displays the conserved thioredoxin fold, consisting of a five-stranded beta-sheet flanked by four helices. Redox-dependent chemical shift perturbations mapped primarily to the conserved WCGPC active-site sequence and other nearby residues, but distant regions of the C-terminal helix were also affected by reduction of the active-site disulfide. Comparisons of the oxidized A. thaliana thioredoxin h1 structure with an h-type thioredoxin from poplar in the reduced state revealed structural differences in the C-terminal helix but no major changes in the active site conformation.     

Expression of three members of the calcium-dependent protein kinase gene family in Arabidopsis thaliana

Calcium-dependent protein kinases (CDPKs) belong to a unique family of enzymes containing a single polypeptide chain with a kinase domain at the amino terminus and a putative calcium-binding EF hands structure at the carboxyl terminus. From Arabidopsis thaliana, we have cloned three distinct cDNA sequences encoding CDPKs, which were designated as atcdpk6, atcdpk9 and atcdpk19. The full-length cDNA sequences for atcdpk6, atcdpk9 and atcdpk19 encode proteins with a molecular weight of 59343, 55376 and 59947, respectively. Recombinant atCDPK6 and atCDPK9 proteins were fully active as kinases whose activities were induced by Ca²⁺. Biochemical studies suggested the presence of an autoinhibitory domain in the junction between the kinase domain and the EF hands structure. Serial deletion of the four EF hands of atCDPK6 demonstrated that the integrity of the four EF hands was crucial to the Ca²⁺ response. All the three atcdpk genes were ubiquitously expressed in the plant as demonstrated by RNA gel blot experiments. Comparison of the genomic sequences suggested that the three cdpk genes have evolved differently. Using antibodies against atCDPK6 and atCDPK9 for immunohistochemical experiments, CDPKs were found to be expressed in specific cell types in a temporally and developmentally regulated manner.     

Genetic differentiation and adaptive evolution at reproductive loci in incipient Drosophila species

Accessory gland proteins (Acps) are part of the seminal fluid of male Drosophila flies. Some Acps have exceptionally high evolutionary rates and evolve under positive selection. Proper interactions between Acps and female reproductive molecules are essential for fertilization. These observations lead to suggestions that fast evolving Acps could be involved in speciation by promoting reproductive incompatibilities between emerging species. To test this hypothesis, we used population genetics data for three sibling species: D. mayaguana, D. parisiena and D. straubae. The latter two species are morphologically very similar and show only incipient reproductive isolation. This system allowed us to examine Acp evolution at different time frames with respect to speciation and reproductive isolation. Comparing data of 14 Acp loci with data obtained for other genomic regions, we found that some Acps show extraordinarily high levels of divergence between D. mayaguana and its two sister species D. parisiena and D. straubae. This divergence was likely driven by adaptive evolution at several loci. No fixed nucleotide differences were found between D. parisiena and D. straubae, however. Nevertheless, some Acp loci did show significant differentiation between these species associated with signs of positive selection; these loci may be involved in this early phase of the speciation process.     

GAGE基因家族的分子进化特征的研究

GAGE基因通常表达于睾丸组织和部分恶性肿瘤组织中, 被认为可能是理想的癌症诊断的标记和治疗的靶位。我们对GAGE基因家族作了生物信息学分析, 发现它们在X染色体上串联成簇排列, 为灵长类所独有, 各拷贝序列趋异度很低。在人类有15个以上的拷贝, 在黑猩猩和猕猴分别有3个和4个。对GAGE基因家族构建进化树, 并估算了复制事件发生的时间, 结果显示在近400万年内陆续发生。用两种方法计算了GAGE各拷贝间的Ka/Ks值, 结果为显著大于1, 表明该基因家族受到正选择作用。这些结果提示该基因可能与灵长类的特征有关, 其在进化上的地位和在配子发育和肿瘤发生中承担的功能值得深入研究。     

Identification of a tRNA isopentenyltransferase gene from Arabidopsis thaliana

The tRNA of most organisms contain modified adenines called cytokinins. Situated next to the anticodon, they have been shown to influence translational fidelity and efficiency. The enzyme that synthesizes cytokinins on pre-tRNA, tRNA isopentenyltransferase (EC 2.5.1.8), has been studied in micro-organisms like Escherichia coli and Saccharomyces cerevisiae, and the corresponding genes have been cloned. We here report the first cloning and functional characterization of a homologous gene from a plant, Arabidopsis thaliana. Expression in S. cerevisiae showed that the gene can complement the anti-suppressor phenotype of a mutant that lacks MOD5, the intrinsic tRNA isopentenyltransferase gene. This was accompanied by the reintroduction of isopentenyladenosine in the tRNA. The Arabidopsis gene is constitutively expressed in seedling tissues.