Tetraodon genome analysis provides further evidence for whole-genome duplication in the ray-finned fish lineage

A whole-genome duplication in the ray-finned fish lineage has been supported by the analyses of the genome sequence of the Japanese pufferfish, Fugu rubripes. Recently, genome sequence of a second teleost fish, the freshwater pufferfish, Tetraodon nigroviridis, was completed. Comparisons of long-range synteny between the Tetraodon and human genomes provided additional evidence for the whole-genome duplication in the ray-finned fish lineage. In the present study, we conducted phylogenetic analysis of the Tetraodon and human proteins to identify ray-finned fish lineage-specific (‘fish-specific’) duplicate genes in the Tetraodon genome. Our analyses provide evidence for 1087 well defined fish-specific duplicate genes in Tetraodon. We also analyzed the Fugu proteome that was predicted in the recent Fugu genome assembly, and identified 346 duplicate genes in addition to the 425 duplicates previously identified. We estimated the ages of duplicate genes using the molecular clock. The ages of duplicate genes in the two pufferfishes independently support a large-scale gene duplication around 380–400 Myr ago. In addition, a burst of recent gene duplications was evident in the Tetraodon lineage. These findings provide further evidence for a whole-genome duplication early in the evolution of ray-finned fishes, and suggest that independent gene duplications have occurred recently in the Tetraodon lineage.     

Phylogenetic Timing of the Fish-Specific Genome Duplication Correlates with the Diversification of Teleost Fish

For many genes, ray-finned fish (Actinopterygii) have two paralogous copies, where only one ortholog is present in tetrapods. The discovery of an additional, almost-complete set of Hox clusters in teleosts (zebrafish, pufferfish, medaka, and cichlid) but not in basal actinopterygian lineages (Polypterus) led to the formulation of the fish-specific genome duplication hypothesis. The phylogenetic timing of this genome duplication during the evolution of ray-finned fish is unknown, since only a few species of basal fish lineages have been investigated so far. In this study, three nuclear genes (fzd8, sox11, tyrosinase) were sequenced from sturgeons (Acipenseriformes), gars (Semionotiformes), bony tongues (Osteoglossomorpha), and a tenpounder (Elopomorpha). For these three genes, two copies have been described previously teleosts (e.g., zebrafish, pufferfish), but only one orthologous copy is found in tetrapods. Individual gene trees for these three genes and a concatenated dataset support the hypothesis that the fish-specific genome duplication event took place after the split of the Acipenseriformes and the Semionotiformes from the lineage leading to teleost fish but before the divergence of Osteoglossiformes. If these three genes were duplicated during the proposed fish-specific genome duplication event, then this event separates the species-poor early-branching lineages from the species-rich teleost lineage. The additional number of genes resulting from this event might have facilitated the evolutionary radiation and the phenotypic diversification of the teleost fish.[Reviewing Editor: Martin Kreitman]     

Microcolinearity and genome evolution in the AdhA region of diploid and polyploid cotton (Gossypium)

Genome sizes vary by several orders of magnitude, driven by mechanisms such as illegitimate recombination and transposable element proliferation. Prior analysis of the CesA region in two cotton genomes that diverged 5–10 million years ago (Ma), and acquired a twofold difference in genome size, revealed extensive local conservation of genic and intergenic regions, with no evidence of the global genome size difference. The present study extends the comparison to include BAC sequences surrounding the gene encoding alcohol dehydrogenase A ( AdhA ) from four cotton genomes: the two co-resident genomes (A_T and D_T) of the allotetraploid, Gossypium hirsutum , as well as the model diploid progenitors, Gossypium arboreum (A) and Gossypium raimondii (D). In contrast to earlier work, evolution in the AdhA region reflects, in a microcosm, the overall difference in genome size, with a nearly twofold difference in aligned sequence length. Most size differences may be attributed to differential accumulation of retroelements during divergence of the genome diploids from their common ancestor, but in addition there has been a biased accumulation of small deletions, such that those in the smaller D genome are on average twice as large as those in the larger A genome. The data also provide evidence for the global phenomenon of 'genomic downsizing' in polyploids shortly after formation. This in part reflects a higher frequency of small deletions post-polyploidization, and increased illegitimate recombination. In conjunction with previous work, the data here confirm the conclusion that genome size evolution reflects many forces that collectively operate heterogeneously among genomic regions.     

Identification of novel tropomyosin 1 genes of pufferfish (Fugu rubripes) on genomic sequences and tissue distribution of their transcripts

Fugu genome database enabled us to identify two novel tropomyosin 1 (TPM1) genes through in silico data mining and isolation of their corresponding cDNAs in vivo. The duplicate TPM1 genes in Japanese pufferfish Fugu rubripes suggest that additional an ancient segmental duplication or whole genome duplication occurred in fish lineage, which, like many other reported Fugu genes, showed reduction in genomic size in comparison with their human homologue. Computer analysis predicted that the coiled-coil probabilities, that were thought to be the most major function of TPM, were the same between the two TPM1 isoforms. We confirmed that the tissue expression profiles of the two TPM1 genes differed from each other, which implied that changes in expression pattern could fix duplicated TPM1 genes although the two TPM1 isoforms appear to have similar function.     

Phylogenetic determination of the pace of transposable element proliferation in plants: copia and LINE-like elements in Gossypium.

Transposable elements contribute significantly to plant genome evolution in myriad ways, ranging from local insertional mutations to global effects exerted on genome size through accumulation. Differential accumulation and deletion of transposable elements may profoundly affect genome size, even among members of the same genus. One example is that of Gossypium (cotton), where much of the 3-fold genome size variation is due to differential accumulation of one gypsy-like LTR retrotransposon, Gorge3. Copia and non-LTR LINE retrotransposons are also major components of the Gossypium genome, but unlike Gorge3, their extant copy numbers do not correlate with genome size. In the present study, we describe the nature and timing of transposition for copia and LINE retrotransposons in Gossypium. Our findings indicate that copia retrotransposons have been active in each lineage since divergence from a common ancestor, and that they have proliferated in a punctuated manner. However, the evolutionary history of LINEs contrasts markedly with that of the copia retrotransposons. Although LINEs have also been active in each lineage, they have accumulated in a stochastically regular manner, and phylogenetic analysis suggests that extant LINE populations in Gossypium are dominated by ancient insertions. Interestingly, the magnitude of transpositional bursts in each lineage corresponds directly with extant estimated copy number.     

利用SSR标记技术研究棉属A、D染色体组的进化

利用SSR分子标记技术,对棉属A、D染色体二倍体及四倍体代表棉种进行了遗传多样性分析。供试的10个二倍体代表棉种间遗传多态性丰富,分子聚类结果与Fryxell棉属分类结果相同。分子水平上进一步揭示出属于D染色体组的拟似棉与其他D染色体组棉种的相似系数最低,A,D染色体组间相似系数很高,该结果支持拟全民族似棉是D染色体组最原始棉种,棉属不同染色体组是共同起源,单元进化的理论,利用栽培的异源四倍体棉种不太适于研究棉属A、D染色体组的进化。     

Characterization and expression analysis of TERMINAL FLOWER1 homologs from cultivated alloteraploid cotton (Gossypium hirsutum) and its diploid progenitors

The seasonal cycle and persistence of a plant is governed by a combination of the determinate or indeterminate status of shoot and root apical meristems. A perennial plant is one in which the apical meristem of at least one of its shoot axes remains indeterminate beyond the first growth season.TERMINAL FLOWER1 (TFL1) genes play important roles in regulating flowering time, the fate of inflorescence meristem and perenniality. To investigate the role of TFL1-like genes in the determination of the apical meristems in an industrially important crop cultivated for its fibers, we isolated and characterized two TFL1 homologs (TFL1a and TFL1b) from tetraploid cultivated cotton (Gossypium hirsutum) and its diploid progenitors (Gossypium arboreum and Gossypium raimondii). All isolated genes maintain the same exon–intron organization. Their phylogenetic analysis at the amino acid level confirmed that the isolated sequences are TFL1-like genes and collocate in the TFL1 clade of the PEBP protein family. Expression analysis revealed that the genes TFL1a and TFL1b have slightly different expression patterns, suggesting different functional roles in the determination of the meristems. Additionally, promoter analysis by computational methods revealed the presence of common binding motifs in TFL1-like promoters. These are the first reported TFL1-like genes isolated from cotton, the most important crop for the textile industry.     

Restriction fragment length polymorphisms in diploid and allotetraploid Gossypium: assigning the late embryogenesis-abundant (Lea) alloalleles in G. hirsutum

Summary We have determined the copy number of 21 genes in an allotetraploid and several diploid species of cotton by gel and dot blot hybridization with cloned cDNAs. The legumin A, legumin B, and all 18 unique Lea (late embryogenesis-abundant) cDNA sequences isolated from the AD allotetraploid Gossypium hirsutum are present in one copy in A, D, E, and F diploid species and in two copies in G. hirsutum. Gel blot analysis of DNAs digested with EcoRI or BamHI usually detects different sized fragments in A and D diploids. Conservation of these restriction fragment length polymorphisms in G. hirsutum allows most of these fragments to be assigned to their respective subgenomes. Furthermore, both subgenomes in G. hirsutum can be distinguished from those in the interfertile allotetraploid G. barbadense. These results show that physical mapping of both sets of chromosomes in an allotetraploid should be possible by segregation analysis.     

Comparative mapping between<Emphasis Type="Italic"> Medicago sativa</Emphasis> and<Emphasis Type="Italic"> Pisum sativum</Emphasis>

Comparative genome analysis has been performed between alfalfa ( Medicago sativa) and pea ( Pisum sativum), species which represent two closely related tribes of the subfamily Papilionoideae with different basic chromosome numbers. The positions of genes on the most recent linkage map of diploid alfalfa were compared to those of homologous loci on the combined genetic map of pea to analyze the degree of co-linearity between their linkage groups. In addition to using unique genes, analysis of the map positions of multicopy (homologous) genes identified syntenic homologs (characterized by similar positions on the maps) and pinpointed the positions of non-syntenic homologs. The comparison revealed extensive conservation of gene order between alfalfa and pea. However, genetic rearrangements (due to breakage and reunion) were localized which can account for the difference in chromosome number (8 for alfalfa and 7 for pea). Based on these genetic events and our increasing knowledge of the genomic structure of pea, it was concluded that the difference in genome size between the two species (the pea genome is 5- to 10-fold larger than that of alfalfa) is not a consequence of genome duplication in pea. The high degree of synteny observed between pea and Medicago loci makes further map-based cloning of pea genes based on the genome resources now available for M. truncatula a promising strategy.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by W. R. McCombie     

Analysis of repetitive DNA in three species of Gossypium

The rate of reassociation of denatured DNA was determined for two selected diploid species, Gossypium thurberi (D genome) and G. arboreum (A genome), and one allotetraploid species, G. hirsutum (AD genome). The relative genome size and DNA content of the chromosomes of the diploids were A greater than D. Renaturation curves indicated that the differences in genome sizes were due primarily to the repetitive DNA content.     

Rate variation among nuclear genes and the age of polyploidy in Gossypium

Molecular evolutionary rate variation in Gossypium (cotton) was characterized using sequence data for 48 nuclear genes from both genomes of allotetraploid cotton, models of its diploid progenitors, and an outgroup. Substitution rates varied widely among the 48 genes, with silent and replacement substitution levels varying from 0.018 to 0.162 and from 0.000 to 0.073, respectively, in comparisons between orthologous Gossypium and outgroup sequences. However, about 90% of the genes had silent substitution rates spanning a more narrow threefold range. Because there was no evidence of rate heterogeneity among lineages for any gene and because rates were highly correlated in independent tests, evolutionary rate is inferred to be a property of each gene or its genetic milieu rather than the clade to which it belongs. Evidence from approximately 200,000 nucleotides (40,000 per genome) suggests that polyploidy in Gossypium led to a modest enhancement in rates of nucleotide substitution. Phylogenetic analysis for each gene yielded the topology expected from organismal history, indicating an absence of gene conversion or recombination among homoeologs subsequent to allopolyploid formation. Using the mean synonymous substitution rate calculated across the 48 genes, allopolyploid cotton is estimated to have formed circa 1.5 million years ago (MYA), after divergence of the diploid progenitors about 6.7 MYA.     

SNP and haplotype identification of the wheat monomeric α-amylase inhibitor genes

Seventy-three gene sequences encoding monomeric α-amylase inhibitors were characterized from cultivated wheat “Chinese Spring”, group 6 nullisomic-tetrasomic lines of “Chinese Spring” and diploid putative progenitors of common wheat. The monomeric α-amylase inhibitors from the different sources shared very high homology (99.54%). The different α-amylase inhibitors, which were determined by the 24 single nucleotide polymorphisms (SNPs) of their gene sequences, were investigated. A total of 15 haplotypes were defined by sequence alignment, among which 9 haplotypes were found with only one single sequence sample. Haplotype H02 was found to be the main haplotype occurring in 83 WMAI sequence samples, followed by haplotype H11. The median-joining network for the 15 haplotypes of monomeric α-amylase inhibitor gene sequences from hexaploid wheats was star like, and at least two subclusters emerged. Furthermore evidence of homologous recombination was found between the haplotypes. The relationship between nucleotide substitutions and the amino acid changes in WMAI of hexaploid wheats was summarized. It was clear that only five polymorphic sites in the nucleotide sequence of WMAI resulted in amino acid variations, and that should be the reason for different structure and function of inhibitors. However, little evidence could be found that there were WMAI genes in the A genome of hexaploid wheat, whereas it could conclude from our results that the A genome diploid wheat had WMAI genes. The overall information on the monomeric α-amylase inhibitors from wheat and Aegilops strongly support the view that these inhibitors have evolved from a common ancestral gene through duplication and mutation. Ji-Rui Wang and Yu-Ming Wei are contributed equally to this paper.     

Characterization of the Brassinosteroid Insensitive 1 Genes of Cotton

Suppression of brassinosteroid (BR) biosynthesis in cotton ovules by treatment with brassinazole inhibits fiber formation, indicating that BR plays an important role in cotton fiber development. Plant responses to brassinosteroids (BR) are mediated through a plasma membrane-bound leucine-rich repeat (LRR) receptor-like protein kinase known as BRI1. Mutations in the BRI1 genes of several species result in dwarfed plants with reduced sensitivity to BR. A single expressed sequence tag (EST) from cotton with strong sequence similarity to Arabidopsis BRI1 ( AtBRI1 ) was identified in a search of publicly available databases. With this EST as a starting point, full-length cDNAs and genomic coding sequences from upland cotton ( Gossypium hirsutum ) BRI1 ( GhBRI1 ) were obtained and characterized. Ectopic expression of this coding sequence in BR-insensitive Arabidopsis plants resulted in recovery of normal growth indicating that GhBRI1 is a functional homologue of AtBRI1. G. hirsutum is an allotetraploid (AADD) derived from diploid ancestors. Analysis of several GhBRI1 cDNAs showed two distinct sequences indicating the presence of two GhBRI1 genes, denoted GhBRI1-1 and GhBRI1-2. Sequence comparisons between these GhBRI1 coding sequences and those from related A and D genome diploid Gossypium species ( G. arboreum and G. thurberi ) indicated that GhBRI1-1 is likely to the A sub-genome orthologue while GhBRI1-2 is from the D sub-genome.     

Types, levels and patterns of low-copy DNA sequence divergence, and phylogenetic implications, for Gossypium genome types

To explore types, levels and patterns of genetic divergence among diploid Gossypium (cotton) genomes, 780 cDNA, genomic DNA and simple sequence repeat (SSR) loci were re-sequenced in Gossypium herbaceum (A1 genome), G. arboreum (A2), G. raimondii (D5), G. trilobum (D8), G. sturtianum (C1) and an outgroup, Gossypioides kirkii. Divergence among these genomes ranged from 7.32 polymorphic base pairs per 100 between G. kirkii and G. herbaceum (A1) to only 1.44 between G. herbaceum (A1) and G. arboreum (A2). SSR loci are least conserved with 12.71 polymorphic base pairs and 3.77 polymorphic sites per 100 base pairs, whereas expressed sequence tags are most conserved with 3.96 polymorphic base pairs and 2.06 sites. SSR loci also exhibit the highest percentage of 'extended polymorphisms' (spanning multiple consecutive nucleotides). The A genome lineage was particularly rapidly evolving, with the D genome also showing accelerated evolution relative to the C genome. Unexpected asymmetry in mutation rates was found, with much more transition than transversion mutation in the D genome after its divergence from a common ancestor shared with the A genome. This large quantity of orthologous DNA sequence strongly supports a phylogeny in which A-C divergence is more recent than A-D divergence, a subject that is of much importance in view of A-D polyploid formation being key to the evolution of the most productive and finest-quality cottons. Loci that are monomorphic within A or D genome types, but polymorphic between genome types, may be of practical importance for identifying locus-specific DNA markers in tetraploid cottons including leading cultivars.     

Spatial and temporal expression of the two sucrose synthase genes in maize: immunohistological evidence

Summary The DNAs of two diploid species of Gossypium, G. herbaceum var. africanum (A₁ genome) and G. raimondii (D₅ genome), and the allotetraploid species, G. hirsutum (A_h and D_h genomes), were characterized by kinetic analyses of single copy and repetitive sequences. Estimated haploid genome sizes of A₁ and D₅ were 1.04 pg and 0.68 pg, respectively, in approximate agreement with cytological observations that A genome chromosomes are about twice the size of D genome chromosomes. This differences in genome size was accounted for entirely by differences in the major repetitive fraction (0.56 pg versus 0.20 pg), as single copy fractions of the two genomes were essentially identical (0.41 pg for A₁ and 0.43 pg for D₅). Kinetic analyses and thermal denaturation measurements of single copy duplexes from reciprocal intergenomic hybridizations showed considerable sequence similarity between A₁ and D₅ genomes (77% duplex formation with an average thermal depression of 6 °C). Moreover, little sequence divergence was detectable between diploid single copy sequences and their corresponding genomes in the allotetraploid, consistent with previous chromosome pairing observations in interspecific F₁ hybrids.Journal paper No. 4461 of the Arizona Agricultural Experiment Station     

Transgene integration and organization in Cotton (Gossypium hirsutum L.) genome

While genetically modified upland cotton (Gossypium hirsutum L.) varieties are ranked among the most successful genetically modified organisms (GMO), there is little knowledge on transgene integration in the cotton genome, partly because of the difficulty in obtaining large numbers of transgenic plants. In this study, we analyzed 139 independently derived T0 transgenic cotton plants transformed by Agrobacterium tumefaciens strain AGL1 carrying a binary plasmid pPZP-GFP. It was found by PCR that as many as 31% of the plants had integration of vector backbone sequences. Of the 110 plants with good genomic Southern blot results, 37% had integration of a single T-DNA, 24% had two T-DNA copies and 39% had three or more copies. Multiple copies of the T-DNA existed either as repeats in complex loci or unlinked loci. Our further analysis of two T1 populations showed that segregants with a single T-DNA and no vector sequence could be obtained from T0 plants having multiple T-DNA copies and vector sequence. Out of the 57 T-DNA/T-DNA junctions cloned from complex loci, 27 had canonical T-DNA tandem repeats, the rest (30) had deletions to T-DNAs or had inclusion of vector sequences. Overlapping micro-homology was present for most of the T-DNA/T-DNA junctions (38/57). Right border (RB) ends of the T-DNA were precise while most left border (LB) ends (64%) had truncations to internal border sequences. Sequencing of collinear vector integration outside LB in 33 plants gave evidence that collinear vector sequence was determined in agrobacterium culture. Among the 130 plants with characterized flanking sequences, 12% had the transgene integrated into coding sequences, 12% into repetitive sequences, 7% into rDNAs. Interestingly, 7% had the transgene integrated into chloroplast derived sequences. Nucleotide sequence comparison of target sites in cotton genome before and after T-DNA integration revealed overlapping microhomology between target sites and the T-DNA (8/8), deletions to cotton genome in most cases studied (7/8) and some also had filler sequences (3/8). This information on T-DNA integration in cotton will facilitate functional genomic studies and further crop improvement.     

Genetics of allozyme variation in Gossypium arboreum L. and Gossypium herbaceum L. (Malvaceae)

Summary Seed protein extracts from 90 accessions of Gossypium arboreum and 70 accessions of Gossypium herbaceum were electrophoretically analyzed for isozyme variation. Eighteen enzyme systems were resolved, ten of which were polymorphic among accessions. No within accession isozyme variation was observed within these highly inbred lines. A minimum of 24 genes encode the isozymes resolved and data is presented for codominant inheritance at 13 loci. Tests for non-random joint segregation in 63 of the 78 possible two-locus combinations from the 13 characterized loci give evidence for four pairs of linked genes (Lap2/Me1 [r=0.160+/-0.027], Lap2/Pgi1 [r= 0.285+/-0.055], Mdh6/Tpi1 [r= 0.197+/-0.028], and 6Pgd2/6Pgd3[r 0.000]. Numerous presumptive duplicate isozyme loci were observed and these were usually expressed as patterns of nonsegregating heteromultimers within accessions. Single gene expression was also observed at several loci. The observed results are in agreement with those of previous cytological investigations which have proposed the polyploid origin of the diploid Old World Gossypiums.     

Opsins and clusters of sensory G-protein-coupled receptors in the sea urchin genome

Rhodopsin-type G-protein-coupled receptors (GPCRs) contribute the majority of sensory receptors in vertebrates. With 979 members, they form the largest GPCR family in the sequenced sea urchin genome, constituting more than 3% of all predicted genes. The sea urchin genome encodes at least six Opsin proteins. Of these, one rhabdomeric, one ciliary and two G(o)-type Opsins can be assigned to ancient bilaterian Opsin subfamilies. Moreover, we identified four greatly expanded subfamilies of rhodopsin-type GPCRs that we call sea urchin specific rapidly expanded lineages of GPCRs (surreal-GPCRs). Our analysis of two of these groups revealed genomic clustering and single-exon gene structures similar to the most expanded group of vertebrate rhodopsin-type GPCRs, the olfactory receptors. We hypothesize that these genes arose by rapid duplication in the echinoid lineage and act as chemosensory receptors of the animal. In support of this, group B surreal-GPCRs are most prominently expressed in distinct classes of pedicellariae and tube feet of the adult sea urchin, structures that have previously been shown to react to chemical stimuli and to harbor sensory neurons in echinoderms. Notably, these structures also express different opsins, indicating that sea urchins possess an intricate molecular set-up to sense their environment.     

Cryptic repeated genomic recombination during speciation in Gossypium gossypioides

The Mexican cotton Gossypium gossypioides is a perplexing entity, with conflicting morphological, cytogenetic, and molecular evidence of its phylogenetic affinity to other American cottons. We reevaluated the evolutionary history of this enigmatic species using 16.4 kb of DNA sequence. Phylogenetic analyses show that chloroplast DNA (7.3 kb), nuclear ribosomal internal transcribed spacers (ITS; 0.69 kb), and unique nuclear genes (8.4 kb) yield conflicting resolutions for G. gossypioides. Eight low-copy nuclear genes provide a nearly unanimous resolution of G. gossypioides as the basalmost American diploid cotton, whereas cpDNA sequences resolve G. gossypioides deeply nested within the American diploid clade sister to Peruvian G. raimondii, and ITS places G. gossypioides in an African (rather than an American) clade. These data, in conjunction with previous evidence from the repetitive fraction of the genome, implicate a complex history for G. gossypioides possibly involving temporally separated introgression events from genetically divergent cottons that are presently restricted to different hemispheres. Based on repetitive nuclear DNA, it appears that G. gossypioides experienced nuclear introgression from an African species shortly after divergence from the remainder of the American assemblage. More recently, hybridization with a Mexican species may have resulted in cpDNA introgression, and possibly a second round of cryptic nuclear introgression. Gossypium gossypioides provides a striking example of the previously unsuspected chimeric nature of some plant genomes and the resulting phylogenetic complexity produced by multiple historical reticulation events.     

Evidence in favour of ancient octaploidy in the vertebrate genome

Vertebrate genomes are larger than invertebrates and show evidence of extensive gene duplication, including many collinear chromosomal segments. On the basis of this intra-genomic synteny, it has been proposed that two rounds of whole genome duplication (octaploidy) occurred early in the vertebrate lineage. Recently, this early vertebrate octaploidy has been challenged on the basis of gene trees. We report new linkage groups encompassing the matrilin (MATN), syndecan (SDC), Eyes Absent (EYA), HCK kinase and SRC kinase paralogous gene quartets. In contrast to other studies, the sequence trees are weakly supportive of ancient octaploidy. It is concluded that there is no strong evidence against the octaploidy, provided that consecutive genome duplication was rapid.