New resources inform study of genome size, content, and organization in nonavian reptiles

Genomic resources for studies of nonavian reptiles have recently improved and will reach a new level of access once the genomes of the painted turtle (Chrysemys picta) and the green anole (Anolis carolinensis) have been published. Eleven speakers gathered for a symposium on reptilian genomics and evolutionary genetics at the 2008 meeting of the Society for Integrative and Comparative Biology in San Antonio, Texas. Presentations described results of reptilian genetic studies concerning molecular evolution, chromosomal evolution, genomic architecture, population dynamics, endocrinology and endocrine disruption, and the evolution of developmental mechanisms. The presented studies took advantage of the recent generation of genetic and genomic tools and resources. Novel findings demonstrated the positive impact made by the improved availability of resources like genome annotations and bacterial artificial chromosomes (BACs). The symposium was timely and important because it provided a vehicle for the dissemination of novel findings that advance the field. Moreover, this meeting fostered the synergistic interaction of the participants as a group, which is anticipated to encourage the funding and creation of further resources such as additional BAC libraries and genomic projects. Novel data have already been collected and studies like those presented in this symposium promise to shape and improve our understanding of overall amniote evolution. Additional reptilian taxa such as the American alligator (Alligator mississippiensis), tuatara (Sphenodon punctatus), and garter snake (Thamnophis sirtalis) should be the foci of future genomic projects. We hope that the following articles in this volume will help promote these efforts by describing the conclusions and the potential that the improvement of genomic resources for nonavian reptiles can continue having in this important area of integrative and comparative biology.     

Molecular systematics of primary reptilian lineages and the tuatara mitochondrial genome

We provide phylogenetic analyses for primary Reptilia lineages including, for the first time, Sphenodon punctatus (tuatara) using data from whole mitochondrial genomes. Our analyses firmly support a sister relationship between Sphenodon and Squamata, which includes lizards and snakes. Using Sphenodon as an outgroup for select squamates, we found evidence indicating a sister relationship, among our study taxa, between Serpentes (represented by Dinodon) and Varanidae. Our analyses support monophyly of Archosauria, and a sister relationship between turtles and archosaurs. This latter relationship is congruent with a growing set of morphological and molecular analyses placing turtles within crown Diapsida and recognizing them as secondarily anapsid (lacking a skull fenestration). Inclusion of Sphenodon, as the only surviving member of Sphenodontia (with fossils from the mid-Triassic), helps to fill a sampling gap within previous analyses of reptilian phylogeny. We also report a unique configuration for the mitochondrial genome of Sphenodon, including two tRNA(Lys) copies and an absence of ND5, tRNA(His), and tRNA(Thr) genes.     

InSatDb: a genomic tool for insect geneticists

Microsatellites show tremendous variation between genomes in terms of their occurrence and composition. Availability of whole genome sequences allows us to study microsatellite characteristics of fully sequenced insect genomes to understand the evolution and biological significance of microsatellites. InSatDb is an insect microsatellite database that provides an interactive interface to query information on microsatellites annotated with size (in base pairs and repeat units), genomic location (exon, intron, up-stream or transposon), nature (perfect or imperfect), and sequence composition (repeat motif and GC%). Here we present a snapshot of the distribution and composition of microsatellites in introns and exons of insect genomes. The data present interesting observations regarding the microsatellite life-cycle and genome flux.     

InSatDb: A Genomic Tool for Insect Geneticists

Microsatellites show tremendous variation between genomes in terms of their occurrence and composition. Availability of whole genome sequences allows us to study microsatellite characteristics of fully sequenced insect genomes to understand the evolution and biological significance of microsatellites. InSatDb is an insect microsatellite database that provides an interactive interface to query information on microsatellites annotated with size (in base pairs and repeat units); genomic location (exon, intron, up-stream or transposon); nature (perfect or imperfect); and sequence composition (repeat motif and GC%). Here, we present a snap shot of the distribution and composition of microsatellites in introns and exons of insect genomes. The data present interesting observations regarding the microsatellite life-cycle and genome flux.     

Spontaneous symmetry breaking in genome evolution

The quest for evolutionary mechanisms providing separation between the coding (exons) and noncoding (introns) parts of genomic DNA remains an important focus of genetics. This work combines an analysis of the most recent achievements of genomics and fundamental concepts of random processes to provide a novel point of view on genome evolution. Exon sizes in sequenced genomes show a lognormal distribution typical of a random Kolmogoroff fractioning process. This implies that the process of intron incretion may be independent of exon size, and therefore could be dependent on intron-exon boundaries. All genomes examined have two distinctive classes of exons, each with different evolutionary histories. In the framework proposed in this article, these two classes of exons can be derived from a hypothetical ancestral genome by (spontaneous) symmetry breaking. We note that one of these exon classes comprises mostly alternatively spliced exons.     

Tuatara (Sphenodon) genomics: BAC library construction, sequence survey, and application to the DMRT gene family

The tuatara (Sphenodon punctatus) is of "extraordinary biological interest" as the most distinctive surviving reptilian lineage (Rhyncocephalia) in the world. To provide a genomic resource for an understanding of genome evolution in reptiles, and as part of a larger project to produce genomic resources for various reptiles (evogen.jgi.doe.gov/second_levels/BACs/our_libraries.html), a large-insert bacterial artificial chromosome (BAC) library from a male tuatara was constructed. The library consists of 215 424 individual clones whose average insert size was empirically determined to be 145 kb, yielding a genomic coverage of approximately 6.3x. A BAC-end sequencing analysis of 121 420 bp of sequence revealed a genomic GC content of 46.8%, among the highest observed thus far for vertebrates, and identified several short interspersed repetitive elements (mammalian interspersed repeat-type repeats) and long interspersed repetitive elements, including chicken repeat 1 element. Finally, as a quality control measure the arrayed library was screened with probes corresponding to 2 conserved noncoding regions of the candidate sex-determining gene DMRT1 and the DM domain of the related DMRT2 gene. A deep coverage contig spanning nearly 300 kb was generated, supporting the deep coverage and utility of the library for exploring tuatara genomics.     

Genome evolution in Reptilia: <Emphasis Type="Italic">in silico</Emphasis> chicken mapping of 12,000 BAC-end sequences from two reptiles and a basal bird

Background

With the publication of the draft chicken genome and the recent production of several BAC clone libraries from non-avian reptiles and birds, it is now possible to undertake more detailed comparative genomic studies in Reptilia. Of interest in particular are the genomic events that transformed the large, repeat-rich genomes of mammals and non-avian reptiles into the minimalist chicken genome. We have used paired BAC end sequences (BESs) from the American alligator (Alligator mississippiensis), painted turtle (Chrysemys picta) and emu (Dromaius novaehollandiae) to investigate patterns of sequence divergence, gene and retroelement content, and microsynteny between these species and chicken.

Results

From a total of 11,967 curated BESs, we successfully mapped 725, 773 and 2597 sequences in alligator, turtle, and emu, respectively, to sites in the draft chicken genome using a stringent BLAST protocol. Most commonly, sequences mapped to a single site in the chicken genome. Of 1675, 1828 and 2936 paired BESs obtained for alligator, turtle, and emu, respectively, a total of 34 (alligator, 2%), 24 (turtle, 1.3%) and 479 (emu, 16.3%) pairs were found to map with high confidence and in the correct orientation and with BAC-sized intermarker distances to single chicken chromosomes, including 25 such paired hits in emu mapping to the chicken Z chromosome. By determining the insert sizes of a subset of BAC clones from these three species, we also found a significant correlation between the intermarker distance in alligator and turtle and in chicken, with slopes as expected on the basis of the ratio of the genome sizes.

Conclusion

Our results suggest that a large number of small-scale chromosomal rearrangements and deletions in the lineage leading to chicken have drastically reduced the number of detected syntenies observed between the chicken and alligator, turtle, and emu genomes and imply that small deletions occurring widely throughout the genomes of reptilian and avian ancestors led to the ~50% reduction in genome size observed in birds compared to reptiles. We have also mapped and identified likely gene regions in hundreds of new BAC clones from these species.

     

Paleogenomic data suggest mammal-like genome size in the ancestral amniote and derived large genome size in amphibians

An unsolved question in evolutionary genomics is whether amniote genomes have been expanding or contracting since the common ancestor of this diverse group. Here, we report on the polarity of amniote genome size evolution using genome size estimates for 14 extinct tetrapod genera from the Paleozoic and early Mesozoic Eras using osteocyte lacunae size as a correlate. We find substantial support for a phylogenetically controlled regression model relating genome size to osteocyte lacunae size (P of slopes <0.01, r2=0.65, phylogenetic signal λ=0.83). Genome size appears to have been homogeneous across Paleozoic crown-tetrapod lineages (average haploid genome size 2.9-3.7 pg) with values similar to those of extant mammals. The differentiation in genome size and underlying architecture among extant tetrapod lineages likely evolved in the Mesozoic and Cenozoic Eras, with expansion in amphibians, contractions along the diapsid lineage, and no directional change within the synapsid lineage leading to mammals.     

A comparative genomics strategy for targeted discovery of single-nucleotide polymorphisms and conserved-noncoding sequences in orphan crops

Completed genome sequences provide templates for the design of genome analysis tools in orphan species lacking sequence information. To demonstrate this principle, we designed 384 PCR primer pairs to conserved exonic regions flanking introns, using Sorghum/Pennisetum expressed sequence tag alignments to the Oryza genome. Conserved-intron scanning primers (CISPs) amplified single-copy loci at 37% to 80% success rates in taxa that sample much of the approximately 50-million years of Poaceae divergence. While the conserved nature of exons fostered cross-taxon amplification, the lesser evolutionary constraints on introns enhanced single-nucleotide polymorphism detection. For example, in eight rice (Oryza sativa) genotypes, polymorphism averaged 12.1 per kb in introns but only 3.6 per kb in exons. Curiously, among 124 CISPs evaluated across Oryza, Sorghum, Pennisetum, Cynodon, Eragrostis, Zea, Triticum, and Hordeum, 23 (18.5%) seemed to be subject to rigid intron size constraints that were independent of per-nucleotide DNA sequence variation. Furthermore, we identified 487 conserved-noncoding sequence motifs in 129 CISP loci. A large CISP set (6,062 primer pairs, amplifying introns from 1,676 genes) designed using an automated pipeline showed generally higher abundance in recombinogenic than in nonrecombinogenic regions of the rice genome, thus providing relatively even distribution along genetic maps. CISPs are an effective means to explore poorly characterized genomes for both DNA polymorphism and noncoding sequence conservation on a genome-wide or candidate gene basis, and also provide anchor points for comparative genomics across a diverse range of species.     

Evolutionary convergence on highly-conserved 3' intron structures in intron-poor eukaryotes and insights into the ancestral eukaryotic genome

The presence of spliceosomal introns in eukaryotes raises a range of questions about genomic evolution. Along with the fundamental mysteries of introns' initial proliferation and persistence, the evolutionary forces acting on intron sequences remain largely mysterious. Intron number varies across species from a few introns per genome to several introns per gene, and the elements of intron sequences directly implicated in splicing vary from degenerate to strict consensus motifs. We report a 50-species comparative genomic study of intron sequences across most eukaryotic groups. We find two broad and striking patterns. First, we find that some highly intron-poor lineages have undergone evolutionary convergence to strong 3' consensus intron structures. This finding holds for both branch point sequence and distance between the branch point and the 3' splice site. Interestingly, this difference appears to exist within the genomes of green alga of the genus Ostreococcus, which exhibit highly constrained intron sequences through most of the intron-poor genome, but not in one much more intron-dense genomic region. Second, we find evidence that ancestral genomes contained highly variable branch point sequences, similar to more complex modern intron-rich eukaryotic lineages. In addition, ancestral structures are likely to have included polyT tails similar to those in metazoans and plants, which we found in a variety of protist lineages. Intriguingly, intron structure evolution appears to be quite different across lineages experiencing different types of genome reduction: whereas lineages with very few introns tend towards highly regular intronic sequences, lineages with very short introns tend towards highly degenerate sequences. Together, these results attest to the complex nature of ancestral eukaryotic splicing, the qualitatively different evolutionary forces acting on intron structures across modern lineages, and the impressive evolutionary malleability of eukaryotic gene structures.     

Genomic Organization of the Human Skeletal Muscle Sodium Channel Gene

Voltage-dependent sodium channels are essential for normal membrane excitability and contractility in adult skeletal muscle. The gene encoding the principal sodium channel α-subunit isoform in human skeletal muscle (SCN4A) has recently been shown to harbor point mutations in certain hereditary forms of periodic paralysis. We have carried out an analysis of the detailed structure of this gene including delineation of intron-exon boundaries by genomic DNA cloning and sequence analysis. The complete coding region of SCN4A is found in 32.5 kb of genomic DNA and consists of 24 exons (54 to > 2.2 kb) and 23 introns (97 bp-4.85 kb). The exon organization of the gene shows no relationship to the predicted functional domains of the channel protein and splice junctions interrupt many of the transmembrane segments. The genomic organization of sodium channels may have been partially conserved during evolution as evidenced by the observation that 10 of the 24 splice junctions in SCN4A are positioned in homologous locations in a putative sodium channel gene in Drosophila (para). The information presented here should be extremely useful both for further identifying sodium channel mutations and for gaining a better understanding of sodium channel evolution.     

A non-adaptationist perspective on evolution of genomic complexity or the continued dethroning of man

A new, non-adaptationist theory of evolution of genomic complexity was recently proposed by Lynch and Conery. This concept holds that increase in complexity seen in eukaryotic genomes is a 'syndrome' caused by increase in genome entropy, which is inevitably triggered by reduction of population size. Here, I discuss the definitions of genomic entropy and complexity and the evidence supporting the entropic theory of genome complexity evolution, including new observations on concordant gain and loss of genes and introns in eukaryotic genomes. I further consider the far-reaching biological and philosophical implications of this theory.     

Isolation and characterization of repetitive DNA sequences from Panax ginseng

Repetitive sequences constitute a significant component of most eukaryotic genomes, and the isolation and characterization of repetitive DNA sequences provide an insight into the organization and evolution of the genome of interest. We report the isolation and characterization of the major classes of repetitive sequences from the genome of Panax ginseng. The isolation of repetitive DNA from P. ginseng was achieved by the reannealing of chemically hydrolyzed (200 bp-1 kb fragments) and heat-denatured genomic DNA to low C(o)t value. The low C(o)t fraction was cloned, and fifty-five P. ginseng clones were identified that contained repetitive sequences. Sequence analysis revealed that the fraction includes repetitive telomeric sequences, species-specific satellite sequences, chloroplast DNA fragments and sequences that are homologous to retrotransposons. Two of the retrotransposon-like sequences are homologous to Ty1/ copia-type retroelements of Zea mays, and six cloned sequences are homologous to various regions of the del retrotransposon of Lilium henryi. The del retrotransposon-like sequences and several novel repetitive DNA sequences from P. ginseng were used to differentiate P. ginseng from P. quinquefolius, and should be useful for evolutionary studies of these disjunct species.     

Analysis of genome organization, composition and microsynteny using 500 kb BAC sequences in chickpea

The small genome size (740 Mb), short life cycle (3 months) and high economic importance as a food crop legume make chickpea (Cicer arietinum L.) an important system for genomics research. Although several genetic linkage maps using various markers and genomic tools have become available, sequencing efforts and their use are limited in chickpea genomic research. In this study, we explored the genome organization of chickpea by sequencing approximately 500 kb from 11 BAC clones (three representing ascochyta blight resistance QTL1 (ABR-QTL1) and eight randomly selected BAC clones). Our analysis revealed that these sequenced chickpea genomic regions have a gene density of one per 9.2 kb, an average gene length of 2,500 bp, an average of 4.7 exons per gene, with an average exon and intron size of 401 and 316 bp, respectively, and approximately 8.6% repetitive elements. Other features analyzed included exon and intron length, number of exons per gene, protein length and %GC content. Although there are reports on high synteny among legume genomes, the microsynteny between the 500 kb chickpea and available Medicago truncatula genomic sequences varied depending on the region analyzed. The GBrowse-based annotation of these BACs is available at http://www.genome.ou.edu/plants_totals.html . We believe that our work provides significant information that supports a chickpea genome sequencing effort in the future.     

Isolation and characterization of genomic clones covering the chicken vitellogenin gene

A series of overlapping recombinant clones, which cover the vitellogenin gene, has been isolated from a phage-lambda linked chicken gene library. The DNA of the overlapping clones spans 28 kb of contiguous DNA sequences in the chicken genome. Electron microscopic analysis of hybrids between vitellogenin mRNA and the genomic clones indicates that the chicken vitellogenin gene has a length of approximately 22 kb, about 3.8 times the size of the mRNA. The mRNA sequence is interrupted by at least 33 intervening sequences (introns). Comparison with the vitellogenin gene A2 from Xenopus laevis (Wahli et al., 1980, Cell 20: 107-117) indicates conservation of the number and length of the exons during evolution. Heteroduplex analysis reveals a short stretch of sequence homology between the genes from chicken and frog.     

Conservation of δ-crystallin gene structure between ducks and chickens

A cloned chicken delta-crystallin cDNA was used to identify two putative delta-crystallin genes in the duck by Southern blot hybridization. A DNA fragment containing most of one of these genes was isolated from a library made in bacteriophage lambda Charon 28A containing genomic DNA from 14-day-old embryonic ducks. Electron microscopy, partial gene sequencing, primer extension analysis using duck mRNA, and comparison with the well-characterized chicken delta-crystallin genes suggest that our cloned duck delta-crystallin gene, like the chicken delta-crystallin genes, is 8-10 kb long and contains 17 exons. Hybridization and sequencing data show great similarity between the homologous 5' untranslated and coding exons of the duck and chicken delta-crystallin genes. Overall, the homologous introns also appear to have approximately 30% sequence similarity, and have been subject to deletion/insertion events. Our partial characterization of duck delta-crystallin gene sequences suggests that this avian and reptilian crystallin family has been conserved during evolution, as have the other crystallin gene families that are expressed in the eye lens.