Identification and chromosomal location of four subfamilies of the rubisco small subunit genes in common wheat

Summary Three different 3 noncoding sequences of wheat rubisco small subunit (SSU) genes (RbcS) were used as probes to identify the gene members of different RbcS subfamilies in the common wheat cultivar Chinese Spring (CS). All genes of the wheat RbcS multigene family were previously assigned to the long arm of homoeologous group 5 and to the short arm of homoeologous group 2 chromosomes of cv CS. Extracted DNA from various aneuploids of these homoeologous groups was digested with four restriction enzymes and hybridized with three different 3 noncoding sequences of wheat SSU clones. All RbcS genes located on the long arm of homoeologous group 5 chromosomes were found to comprise a single subfamily, while those located on the short arm of group 2 comprised three subfamilies. Each of the ancestral diploid genomes A, B, and D has at least one representative gene in each subfamily, suggesting that the divergence into subfamilies preceded the differentiation into species. This divergence of the RbcS genes, which is presumably accompanied by a similar divergence in the 5 region, may lead to differential expression of various subfamilies in different tissues and in different developmental stages, in response to different environmental conditions. Moreover, members of one subfamily that belong to different genomes may have diverged also in the coding sequence and, consequently, code for distinguishable SSU. It is assumed that such utilization of the RbcS multigene family increases the adaptability and phenotypic plasticity of common wheat over its diploid progenitors.     

Selective Pressures on <Emphasis Type="Italic">Drosophila</Emphasis> Chemosensory Receptor Genes

The evolution and patterns of selection of genes encoding 10 Drosophila odorant receptors (Or) and the sex pheromone receptor Gr68a were investigated by comparing orthologous sequences across five to eight ecologically diverse species of Drosophila. Using maximum likelihood estimates of dN/dS ratios we show that all 11 genes sampled are under purifying selection, indicating functional constraint. Four of these genes (Or33c, Or42a, Or85e, and Gr68a) may be under positive selection, and if so, there is good evidence that 12 specific amino acid sites may be under positive selection. All of these sites are predicted to be located either in loop regions or just inside membrane spanning regions, and interestingly one of the two sites in Gr68a is in a similar position to a previously described polymorphism in Gr5a that causes a shift in sensitivity to its ligand trehalose. For three Ors, possible evidence for positive selection was detected along a lineage. These include Or22a in the lineage leading to D. mauritiana and Or22b in the lineage leading to D. simulans. This is of interest in light of previous data showing a change in ligand response profile for these species in the sensory neuron (ab3A) which expresses both Or22a and Or22b in D. melanogaster. In summary, while the main chemosensory function and/or structural integrity of these 10 Or genes and Gr68a are evolutionarily preserved, positive selection appears to be acting on some of these genes, at specific sites and along certain lineages, and provides testable hypotheses for further functional experimentation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. [Reviewing Editor: Dr. David Pollock]     

Evolutionary dynamics of the immunoglobulin heavy chain variable region genes in vertebrates

Immunoglobulin heavy chains are polypeptides encoded by four genes: variable (IGHV), joining (IGHJ), diversity (IGHD), and constant (IGHC) region genes. The number of IGHV genes varies from species to species. To understand the evolution of the IGHV multigene family, we identified and analyzed the IGHV sequences from 16 vertebrate species. The results show that the numbers of functional and nonfunctional IGHV genes among different species are positively correlated. The number of IGHV genes is relatively stable in teleosts, but the intragenomic sequence variation is generally higher in teleosts than in tetrapods. The IGHV genes in tetrapods can be classified into three phylogenetic clans (I, II, and III). The clan III and/or II genes are relatively abundant, whereas clan I genes exist in small numbers or are absent in most species. The genomic organization of clan I, II, and III IGHV genes varies considerably among species, but the entire IGHV locus seems to be conserved in the subtelomeric or near-centromeric region of chromosome. The presence or absence of specific IGHV clan members and the lineage-specific expansion and contraction of IGHV genes indicate that the IGHV locus continues to evolve in a species-specific manner. Our results suggest that the evolution of IGHV multigene family is more complex than previously thought and that several factors may act synergistically for the development of antibody repertoire. Electronic supplementary materials The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Chromosome assignment of four photosynthesis-related genes and their variability in wheat species

Copy numbers of four photosynthesis-related genes, PhyA, Ppc, RbcS and Lhcb1 ^*1, in wheat genomes were estimated by slot-blot analysis, and these genes were assigned to the chromosome arms of common wheat by Southern hybridization of DNA from an aneuploid series of the cultivar Chinese Spring. The copy number of PhyA was estimated to be one locus per haploid genome, and this gene was assigned to chromosomes 4AL, 4BS and 4DS. The Ppc gene showed a low copy number of small multigenes, and was located on the short arm of homoeologous group 3 chromosomes and the long arm of chromosomes of homoeologous group 7. RbcS consisted of a multigene family, with approximately 100 copies in the common wheat genome, and was located on the short arm of group 2 chromosomes and the long arm of group 5 chromosomes. Lhcb1 ^*1 also consisted of a multigene family with about 50 copies in common wheat. Only a limited number of restriction fragments (approximately 15%) were used to determine the locations of members of this family on the long arm of group 1 chromosomes owing to the multiplicity of DNA bands. The variability of hybridized bands with the four genes was less in polyploids, but was more in the case of multigene families. RFLP analysis of polyploid wheats and their presumed ancestors was carried out with probes of the oat PhyA gene, the maize Ppc gene, the wheat RbcS gene and the wheat Lhcb1 ^*1 gene. The RFLP patterns of common wheat most closely resembled those of T. Dicoccum (Emmer wheat), T. urartu (A genome), Ae. speltoides (S genome) and Ae. squarrosa (D genome). Diversification of genes in the wheat complex appear to have occurred mainly at the diploid level. Based on RFLP patterns, B and S genomes were clustered into two major groups. The fragment numbers per genome were reduced in proportion to the increase of ploidy level for all four genes, suggesting that some mechanism(s) might operate to restrict, and so keep to a minimum, the gene numbers in the polyploid genomes. However, the RbcS genes, located on 2BS, were more conserved (double dosage), indicating that the above mechanism(s) does not operate equally on individual genes.     

Functional and molecular evolution of olfactory neurons and receptors for aliphatic esters across the Drosophila genus

Insect olfactory receptor (Or) genes are large, rapidly evolving gene families of considerable interest for evolutionary studies. They determine the responses of sensory neurons which mediate critical behaviours and ecological adaptations. We investigated the evolution across the genus Drosophila of a subfamily of Or genes largely responsible for the perception of ecologically relevant aliphatic esters; products of yeast fermentation and fruits. Odour responses were recorded from eight classes of olfactory receptor neurons known to express this Or subfamily in D. melanogaster and from homologous sensilla in seven other species. Despite the fact that these species have diverged over an estimated 40 million years, we find that odour specificity is largely maintained in seven of the eight species. In contrast, we observe extensive changes in most neurons of the outgroup species D. virilis, and in two neurons across the entire genus. Some neurons show small shifts in specificity, whilst some dramatic changes correlate with gene duplication or loss. An olfactory receptor neuron response similarity tree did not match an Or sequence similarity tree, but by aligning Or proteins of likely functional equivalence we identify residues that may be relevant for odour specificity. This will inform future structure–function studies of Drosophila Ors.     

Phylogenetic and evolutionary relationships and developmental expression patterns of the zebrafish <Emphasis Type="Italic">twist</Emphasis> gene family

Four members of the twist gene family (twist1a, 1b, 2, and 3) are found in the zebrafish, and they are thought to have arisen through three rounds of gene duplication, two of which occurred prior to the tetrapod-fish split. Phylogenetic analysis groups most of the vertebrate Twist1 peptides into clade I, except for the Twist1b proteins of the acanthopterygian fish (medaka, pufferfish, stickleback), which clustered within clade III. Paralogies and orthologies among the zebrafish, medaka, and human twist genes were determined using comparative synteny analysis of the chromosomal regions flanking these genes. Comparative nucleotide substitution analyses also revealed a faster rate of nucleotide mutation/substitution in the acanthopterygian twist1b compared to the zebrafish twist1b, thus accounting for their anomalous phylogenetic clustering. We also observed minimal expression overlap among the four twist genes, suggesting that despite their significant peptide similarity, their regulatory controls have diverged considerably, with minimal functional redundancy between them. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evolution of the <Emphasis Type="Italic">Groucho</Emphasis>/<Emphasis Type="Italic">Tle</Emphasis> gene family: gene organization and duplication events

The Groucho/Tle family of corepressor proteins has important roles in development and in adult tissue in both Protostomes and Deuterostomes. In Drosophila, a single member of this family has been identified. Unlike in Protostomes, most Deutrostomes contain more than two full-length Tle genes. In this study, I analyse the genomic organization and phylogenetic relationship between the long and short forms of the Groucho/Tle family members in Chordata. The genomic location and sequence similarities suggest that Aes/Grg5 and Tle6/Grg6 arose from duplication of the Tle2 gene; each evolved independently and acquired new functions as negative regulators of the other Tle proteins. Based on these data, a model for Groucho/Tle gene evolution is proposed. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Localization of Shaw-related K+ channel genes on mouse and human chromosomes

Four related genes, Shaker, Shab, Shaw, and Shal, encode voltage-gated K⁺ channels in Drosophila. Multigene subfamilies corresponding to each of these Drosophila genes have been identified in rodents and primates; this suggests that the four genes are older than the common ancestor of present-day insects and mammals and that the expansion of each into a family occurred before the divergence of rodents and primates.In order to define these evolutionary relationships more precisely and to facilitate the search for mammalian candidate K⁺ channel gene mutations, we have mapped members of the Shaw-homologous gene family in humans and mice. Fluorescence in situ hybridization analysis of human metaphase chromosomes mapped KCNC2 (KShIIIA, KV3.2) and KCNC3 (KShIIID, KV3.3) to Chromosome (Chr) 19q13.3-q13.4. Inheritance patterns of DNA restriction fragment length variants in recombinant inbred strains of mice placed the homologous mouse genes on distal Chr 10 near Ms15-8 and Mdm-1. The mouse Kcnc1 (KShIIIB, NGK2-KV4, KV3.1) gene mapped to Chr 7 near Tam-1.These results are consistent with the hypothesis that the generation of the mammalian KCNC gene family included both duplication events to generate family members in tandem arrays (KCNC2, KCNC3) and dispersion of family members to unlinked chromosomal sites (KCNC1). The KNCN2 and KCNC3 genes define a new synteny group between humans and mice.     

Chromosomal localizations of mouse Fgf2 and Fgf5 genes

The mouse Fgf2 and Fgf5 genes, two members of the FGF family whose human homologs are on the long arm of chromosome (chr) 4, were localized by using in situ chromosomal hybridization, to mouse Chrs 3 and 5, respectively.     

The genomic organization of HeT-A retroposons inDrosophila melanogaster

Members of theDrosophila HeT-A family of transposable elements are LINE-like retroposons that are found at telomeres and in centric heterochromatin. We recently characterized an active HeT-A element that had transposed to a broken chromosome end fewer than mine generations before it was isolated. The sequence arerangement of this element, called 9D4, most likely represents the organization of an actively transposing member of the HeT-A family. Here we assess the degree of divergence among members of the HeT-A family and test a model of telomere length maintenance based on HeT-A transposition. The region containing the single open reading frame of this element appears to be more highly conserved than the non-coding regions. The HeT-A element has been implicated in theDrosophila telomere elongation process, because frequent transpositions to chromosome ends are sufficient to counter-balance nucleotide loss due to incomplete DNA replication. The proposed elongation model and the hypothetical mechanism of HeT-A transposition predict a predominant orientation of HeT-A elements with their oligo (A) tails facing proximally at chromosome ends, as well as the existence of irregular tandem arrays of HeT-A elements at chromosome ends resulting from transposition of new HeT-A elements onto chromosome ends with existing elements. Twenty-nine different HeT-A fragments were isolated from directional libraries that were enriched in terminal DNA fragments. Sequence analyses of these fragments and comparisons with the organization of the HeT-A element, 9D4, fit these two predictions and support the model ofDrosophila telomere elongation by transposition of HeT-A elements.     

A novel locus (RP33) for autosomal dominant retinitis pigmentosa mapping to chromosomal region 2cen-q12.1

Retinitis pigmentosa (RP) is a heterogeneous group of progressive degenerative disorders of the retina with a strong genetic component. Here, we report the clinical and genetic findings in a Chinese family in which autosomal dominant RP (adRP) was inherited by 13 affected members in four generations. Using a genome-wide linkage screening approach, a novel disease locus (RP33) was assigned to the long arm of chromosome 2. A maximum multi-point LOD score of 4.69 was reached at marker D2S2222 in 2q11.2. Meiotic recombination events in affected members placed RP33 in a 15.5 cM region between D2S329 and D2S2229. From meiotic recombinations in two unaffected members RP33 was further refined to a 4.8 cM (9.5 Mb) interval flanked by D2S2159 and D2S1343 in chromosomal region 2cen-q12.1. No disease-associated mutations were detected in the candidate genes SEMA4C, CNGA3 or HNK1ST from within the region. MERTK, a known disease gene for autosomal recessive RP located close to RP33 was similarly excluded. Clinically, the family presented relatively late onset of night blindness, gradually decreased visual acuity, progressive loss of peripheral visual field and typical RP fundus changes in the mid-periphery of the retina. In conclusion, a novel locus for adRP has been assigned to chromosomal region 2cen-q12.1, which in the present kindred was associated with a relatively late onset form of the disease.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.Chen Zhao and Shasha Lu have contributed equally to this study     

Birth and Death of Genes and Functions in the β-Esterase Cluster of <Emphasis Type="Italic">Drosophila</Emphasis>

Here we analyze the molecular evolution of the β-esterase gene cluster in the Drosophila genus using the recently released genome sequences of 12 Drosophila species. Molecular evolution in this small cluster is noteworthy because it contains contrasting examples of the types and stages of loss of gene function. Specifically, missing orthologs, pseudogenes, and null alleles are all inferred. Phylogenetic analyses also suggest a minimum of 9 gene gain–loss events; however, the exact number and age of these events is confounded by interparalog recombination. A previous enigma, in which allozyme loci were mapped to β-esterase genes that lacked catalytically essential amino acids, was resolved through the identification of neighbouring genes that contain the canonical catalytic residues and thus presumably encode the mapped allozymes. The originally identified genes are evolving with selective constraint, suggesting that they have a “noncatalytic” function. Curiously, 3 of the 4 paralogous β-esterase genes in the D. ananassae genome sequence have single inactivating (frame-shift or nonsense) mutations. To determine whether these putatively inactivating mutations were fixed, we sequenced other D. ananassae alleles of these four loci. We did not find any of the 3 inactivating mutations of the sequenced strain in 12 other strains; however, other inactivating mutations were observed in the same 3 genes. This is reminiscent of the high frequency of null alleles observed in one of the β-esterase genes (Est7/EstP) of D. melanogaster. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

An integrated approach for the comparative analysis of a multigene family: The nicotianamine synthase genes of barley

Recent genomic projects reveal that about half of the gene repertoire in plant genomes is made up by multigene families. In this paper, a set of structural and phylogenetic analyses have been applied to compare the differently sized nicotianamine synthase (NAS) gene families in barley and rice. Nicotianamine acts as a chelator of iron and other heavy metals and plays a key role in uptake, phloem transport and cytoplasmic distribution of iron, challenging efforts for the breeding of iron-efficient crop plants. Nine barley NAS genes have been mapped, and co-linearity of flanking genes in barley and rice was determined. The combined analyses reveal that the NAS multigene family members in barley originated through at least one duplication event that occurred before the divergence of rice and barley. Additional duplications appear to have occurred within each of the species. Although we detected no evidence for positive selection of recently duplicated genes within species, codon-based tests revealed evidence for positive selection having contributed to the divergence of some amino acids. The integrated comparative and phylogenetic analysis improved our current view of NAS gene family evolution, might facilitate the functional characterization of individual members and is applicable to other multigene families. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Genome-wide Comparative Analysis of Annexin Superfamily in Plants

Most annexins are calcium-dependent, phospholipid-binding proteins with suggested functions in response to environmental stresses and signaling during plant growth and development. They have previously been identified and characterized in Arabidopsis and rice, and constitute a multigene family in plants. In this study, we performed a comparative analysis of annexin gene families in the sequenced genomes of Viridiplantae ranging from unicellular green algae to multicellular plants, and identified 149 genes. Phylogenetic studies of these deduced annexins classified them into nine different arbitrary groups. The occurrence and distribution of bona fide type II calcium binding sites within the four annexin domains were found to be different in each of these groups. Analysis of chromosomal distribution of annexin genes in rice, Arabidopsis and poplar revealed their localization on various chromosomes with some members also found on duplicated chromosomal segments leading to gene family expansion. Analysis of gene structure suggests sequential or differential loss of introns during the evolution of land plant annexin genes. Intron positions and phases are well conserved in annexin genes from representative genomes ranging from Physcomitrella to higher plants. The occurrence of alternative motifs such as K/R/HGD was found to be overlapping or at the mutated regions of the type II calcium binding sites indicating potential functional divergence in certain plant annexins. This study provides a basis for further functional analysis and characterization of annexin multigene families in the plant lineage.     

Preferential transposition of aDrosophila P element to the corresponding region of the homologous chromosome

Genetic and physical analyses have demonstrated an intimate interaction or pairing of homologous chromosomes in the nuclei of manyDrosophila cell types. Experiments were performed to determine whether P elements transposing from a given chromosome to its homolog would preferentially insert in the region corresponding to the donor site, perhaps due to such a proximity. AP[lacZ;ry ⁺] element at thecactus locus (35F) on the second chromosome was mobilized and 96 insertions on the homolog were recovered. The distribution of these new insertions was determined by recombination mapping and molecular analysis, and compared with a control set of 93 second-chromosome insertions originating from theX chromosome. A nearly threefold preference was observed for re-insertion in a region of two to three number divisions aroundcactus on the homolog. However, none of these local insertions was actually within 50 kb of the site atcactus corresponding to the starting site. This is in marked contrast to the previously described phenomenon of intrachromosomal local transposition, where the majority of local transpositions are within 10 kb. The data suggest that the mechanisms for intrachromosomal and interchromosomal local transposition are distinct, and are consistent with a model for interchromosomal local transposition involving proximity of homologous chromosomal regions in the nuclei of the germline cells.     

Evolutionary genetics ofDrosophila esterases

Over 30 carboxylester hydrolases have been identified inD. melanogaster. Most are classified as acetyl, carboxyl or cholinesterases. Sequence similarities among most of the carboxyl and all the cholinesterases so far characterised fromD. melanogaster and other eukaryotes justify recognition of a carboxyl/cholinesterase multigene family. This family shows minimal sequence similarities with other esterases but crystallographic data for a few non-drosophilid enzymes show that the family shares a distinctive overall structure with some other carboxyl and aryl esterases, so they are all put in one superfamily of / hydrolases. Fifteen esterase genes have been mapped inD. melanogaster and twelve are clustered at two chromosomal sites. The constitution of each cluster varies acrossDrosophila species but two carboxyl esterases in one cluster are sufficiently conserved that their homologues can be identified among enzymes conferring insecticide resistance in other Diptera. Sequence differences between two other esterases, the EST6 carboxyl esterase and acetylcholinesterase, have been interpreted against the consensus super-secondary structure for the carboxyl/cholinesterase multigene family; their sequence differences are widely dispersed across the structure and include substantial divergence in substrate binding sites and the active site gorge. This also applies when EST6 is compared across species where differences in its expression indicate a difference in function. However, comparisons within and among species where EST6 expression is conserved show that many aspects of the predicted super-secondary structure are tightly conserved. Two notable exceptions are a pair of polymorphisms in the substrate binding site of the enzyme inD. melanogaster. These polymorphisms are associated with differences in substrate interactions_vitro} and demographic data indicate that the alternative forms are not selectively equivalentin vivo.     

Clustering, haplotype diversity and locations of MIC-3: a unique root-specific defense-related gene family in Upland cotton (Gossypium hirsutum L.)

MIC-3 is a recently identified gene family shown to exhibit increased root-specific expression following nematode infection of cotton plants that are resistant to root-knot nematode. Here, we cloned and sequenced MIC-3 genes from selected diploid and tetraploid cotton species to reveal sequence differences at the molecular level and identify chromosomal locations of MIC-3 genes in Gossypium species. Detailed sequence analysis and phylogenetic clustering of MIC-3 genes indicated the presence of multiple MIC-3 gene members in Gossypium species. Haplotypes of a MIC-3 gene family member were discovered by comparative analysis among consensus sequences across genotypes within an individual clade in the phylogram to overcome the problem of duplicated loci in the tetraploid cotton. Deficiency tests of the SNPs delimited six A_t-genome members of the MIC-3 family clustered to chromosome arm 4sh, and one D_t-genome member to chromosome 19. Clustering was confirmed by long-PCR amplification of the intergenic regions using A_t-genome-specific MIC-3 primer pairs. The clustered distribution may have been favored by selection for responsiveness to evolving disease and/or pest pressures, because large variants of the MIC-3 gene family may have been recovered from small physical areas by recombination. This could give a buffer against selection pressure from a broad range of pest and pathogens in the future. To our knowledge, these are the first results on the evolution of clustering and genome-specific haplotype members of a unique cotton gene family associated with resistant response against a major pathogen.