Carl J. Schmidt

Multi-platform next-generation sequencing of the domestic Turkey (Meleagris gallopavo): Genome assembly and analysis

The combined application of next-generation sequencing platforms has provided an economical approach to unlocking the potential of the turkey genome.

GoFigure: Automated Gene Ontology ™ annotation

SUMMARY We have developed a web tool to predict Gene Ontology (GO) terms. The tool accepts an input DNA or protein sequence, and uses BLAST to identify homologous sequences in GO annotated databases. A graph is returned to the user via email. AVAILABILITY The tool is freely available at: http://udgenome.ags.udel.edu/frm_go.html/

Three crocodilian genomes reveal ancestral patterns of evolution among archosaurs

INTRODUCTION Crocodilians and birds are the two extant clades of archosaurs, a group that includes the extinct dinosaurs and pterosaurs. Fossils suggest that living crocodilians (alligators, crocodiles, and gharials) have a most recent common ancestor 80 to 100 million years ago. Extant crocodilians are notable for their distinct morphology, limited intraspecific variation, and slow karyotype evolution. Despite their unique biology and phylogenetic position, little is known about genome evolution within crocodilians. Evolutionary rates of tetrapods inferred from DNA sequences anchored by ultraconserved elements. Evolutionary rates among reptiles vary, with especially low rates among extant crocodilians but high rates among squamates. We have reconstructed the genomes of the common ancestor of birds and of all archosaurs (shown in gray silhouette, although the morphology of these species is uncertain). RATIONALE Genome sequences for the American alligator, saltwater crocodile, and Indian gharial—representatives of all three extant crocodilian families—were obtained to facilitate better understanding of the unique biology of this group and provide a context for studying avian genome evolution. Sequence data from these three crocodilians and birds also allow reconstruction of the ancestral archosaurian genome. RESULTS We sequenced shotgun genomic libraries from each species and used a variety of assembly strategies to obtain draft genomes for these three crocodilians. The assembled scaffold N50 was highest for the alligator (508 kilobases). Using a panel of reptile genome sequences, we generated phylogenies that confirm the sister relationship between crocodiles and gharials, the relationship with birds as members of extant Archosauria, and the outgroup status of turtles relative to birds and crocodilians. We also estimated evolutionary rates along branches of the tetrapod phylogeny using two approaches: ultraconserved element–anchored sequences and fourfold degenerate sites within stringently filtered orthologous gene alignments. Both analyses indicate that the rates of base substitution along the crocodilian and turtle lineages are extremely low. Supporting observations were made for transposable element content and for gene family evolution. Analysis of whole-genome alignments across a panel of reptiles and mammals showed that the rate of accumulation of micro-insertions and microdeletions is proportionally lower in crocodilians, consistent with a single underlying cause of a reduced rate of evolutionary change rather than intrinsic differences in base repair machinery. We hypothesize that this single cause may be a consistently longer generation time over the evolutionary history of Crocodylia. Low heterozygosity was observed in each genome, consistent with previous analyses, including the Chinese alligator. Pairwise sequential Markov chain analysis of regional heterozygosity indicates that during glacial cycles of the Pleistocene, each species suffered reductions in effective population size. The reduction was especially strong for the American alligator, whose current range extends farthest into regions of temperate climates. CONCLUSION We used crocodilian, avian, and outgroup genomes to reconstruct 584 megabases of the archosaurian common ancestor genome and the genomes of key ancestral nodes. The estimated accuracy of the archosaurian genome reconstruction is 91% and is higher for conserved regions such as genes. The reconstructed genome can be improved by adding more crocodilian and avian genome assemblies and may provide a unique window to the genomes of extinct organisms such as dinosaurs and pterosaurs. To provide context for the diversification of archosaurs—the group that includes crocodilians, dinosaurs, and birds—we generated draft genomes of three crocodilians: Alligator mississippiensis (the American alligator), Crocodylus porosus (the saltwater crocodile), and Gavialis gangeticus (the Indian gharial). We observed an exceptionally slow rate of genome evolution within crocodilians at all levels, including nucleotide substitutions, indels, transposable element content and movement, gene family evolution, and chromosomal synteny. When placed within the context of related taxa including birds and turtles, this suggests that the common ancestor of all of these taxa also exhibited slow genome evolution and that the comparatively rapid evolution is derived in birds. The data also provided the opportunity to analyze heterozygosity in crocodilians, which indicates a likely reduction in population size for all three taxa through the Pleistocene. Finally, these data combined with newly published bird genomes allowed us to reconstruct the partial genome of the common ancestor of archosaurs, thereby providing a tool to investigate the genetic starting material of crocodilians, birds, and dinosaurs.

Sequencing three crocodilian genomes to illuminate the evolution of archosaurs and amniotes

The International Crocodilian Genomes Working Group (ICGWG) will sequence and assemble the American alligator (Alligator mississippiensis), saltwater crocodile (Crocodylus porosus) and Indian gharial (Gavialis gangeticus) genomes. The status of these projects and our planned analyses are described.

Coordinated international action to accelerate genome-to-phenome with FAANG, the Functional Annotation of Animal Genomes project

We describe the organization of a nascent international effort, the Functional Annotation of Animal Genomes (FAANG) project, whose aim is to produce comprehensive maps of functional elements in the genomes of domesticated animal species.

Comparison of a modern broiler line and a heritage line unselected since the 1950s

Selecting chicken for improved meat production has altered the relative growth of organs in modern broiler lines compared with heritage lines. In this study, we compared the growth and feed efficiency of a heritage line, UIUC, with a modern production line, Ross 708, for 5 wk posthatch. During this period, the BW and feed efficiency of the modern strain was higher than that of the heritage line, indicating that the Ross 708 birds were more efficient than the UIUC birds at converting feed to body mass. The relative growth of the breast, heart, liver, and intestine were also compared during these 5 wk. The breast muscle of the heritage line constituted 9% of the total body mass at 5 wk, whereas in the modern line, the breast muscle was 18% of the total mass of the bird. In contrast, the relative size of the heart decreased after d 14 in the modern line, suggesting that selection for increased breast muscle has translated into relatively less weight of the heart muscle. The liver matured earlier in modern lines, possibly improving nutrient utilization as the birds shift from lipid- to carbohydrate-rich feed. Finally, jejunal and ileal sections of the intestine were 20% longer in the modern line, perhaps allowing for increased nutrient absorption.

THE G PROTEIN GAMMA SUBUNIT : REQUIREMENTS FOR DIMERIZATION WITH BETA SUBUNITS

Guanine nucleotide-binding protein β and γ subunits form a tightly bound complex that can only be separated by denaturation. Assembly of β and γ subunits is a complicated process. The β1 and γ2 subunits can be synthesized in vitro in rabbit reticulocyte lysate and then assembled into dimers, but β1 cannot form βγ dimers when synthesized in a wheat germ extract. In contrast, γ2 translated in either system can dimerize with β1, suggesting that dimerization-competent γ2 can be synthesized without the aid of specific chaperonins or other cofactors. Dimerization-competent γ2 in solution forms an asymmetric particle with a Stokes radius of about 21 ± 0.4 Å (n = 4), s20,w of 0.9 S (range 0.8-1.0 S, n = 2), and frictional ratio of 1.57 (assuming no hydration). To define the part of γ2 that is needed for native βγ dimer formation, a series of N- and C-terminal truncations were generated, synthesized in vitro, and incubated with β1. Dimerization was assessed by stabilization of β1 to tryptic proteolysis. Truncation of up to 13 amino acids at the C terminus did not affect dimerization with β1, whereas removal of 27 amino acids prevented it. Therefore, a region between residues 45 and 59 of γ2 is important for dimerization. Truncation of 15 amino acids from the N terminus greatly diminished the formation of βγ dimers, while removal of 25 amino acids entirely blocked it. Thus, another region important for forming native βγ is near the N terminus. Extension of the N terminus by 12 amino acids that include the influenza virus hemagglutinin epitope did not prevent βγ dimerization. Furthermore, in intact 35S-labeled COS cells, epitope-tagged γ2 coimmunoprecipitates with β and α subunits. The N-terminal epitope tag must lie at the surface of the heterotrimer since it prevents neither heterotrimer formation nor access of the antibody.

Myosin I heavy-chain genes of Acanthamoeba castellanii: cloning of a second gene and evidence for the existence of a third isoform

We have determined the complete sequence and structure of a second myosin I heavy-chain gene from Acanthamoeba castellanii. This gene, which we have named MIL, spans approx. 6kb, is split by 17 introns, encodes a 1147-aa polypeptide, and is transcribed in log-phase cells. The positions of six of the introns are conserved relative to a vertebrate muscle myosin gene. Similar to the previously characterized MIB heavy-chain gene, the deduced MIL heavy-chain aa sequence reveals a 125-kDa protein composed of a myosin globular head domain joined to a novel, approx. 50-kDa C-terminal domain that is rich in glycine, proline and alanine residues. There are differences, however, between MIL and MIB in the sequence organization of their unconventional C-terminal domains. We conclude from this and other data that Acanthamoeba express at least three myosin I heavy-chain isoforms: MIL, plus MIA and MIB, whose purifications have been published previously. Amoeba genomic DNA blots probed with a short, highly conserved sequence whose position is transposed between MIB and MIL indicate that the Acanthamoeba myosin I heavy-chain gene family may actually contain as many as six genes. Finally, we compared the myosin I sequences with those of two related proteins, Drosophila NinaC and the bovine myosin I-like protein, and found that a portion of the unconventional C-terminal domains of the amoeba myosins I and the bovine protein appear to be related.

The chicken gene nomenclature committee report

Comparative genomics is an essential component of the post-genomic era. The chicken genome is the first avian genome to be sequenced and it will serve as a model for other avian species. Moreover, due to its unique evolutionary niche, the chicken genome can be used to understand evolution of functional elements and gene regulation in mammalian species. However comparative biology both within avian species and within amniotes is hampered due to the difficulty of recognising functional orthologs. This problem is compounded as different databases and sequence repositories proliferate and the names they assign to functional elements proliferate along with them. Currently, genes can be published under more than one name and one name sometimes refers to unrelated genes. Standardized gene nomenclature is necessary to facilitate communication between scientists and genomic resources. Moreover, it is important that this nomenclature be based on existing nomenclature efforts where possible to truly facilitate studies between different species. We report here the formation of the Chicken Gene Nomenclature Committee (CGNC), an international and centralized effort to provide standardized nomenclature for chicken genes. The CGNC works in conjunction with public resources such as NCBI and Ensembl and in consultation with existing nomenclature committees for human and mouse. The CGNC will develop standardized nomenclature in consultation with the research community and relies on the support of the research community to ensure that the nomenclature facilitates comparative and genomic studies.

A multi-agent system for the quantitative simulation of biological networks

We apply the multi-agent system (MAS) platform to the task of biological network simulation. In this paper, we describe the simulation of signal transduction (ST) networks using the DECAF [9] MAS architecture. Unlike previous approaches that relied on systems of differential equations (DE), the distributed framework of MAS scales well and allows us to model large, highly interconnected ST pathways. This scalability is achieved by adopting a hybrid strategy that factors macro-level measures, such as reaction rateconstants, to calculate the stochastic kinetics at the level of individual molecules. Thus, by capturing the ST domain at an intermediate level of abstraction, we are able to retain much of the granularity afforded by a purely individual-based approach. The task distribution within a MAS enables us to model certain physical properies, such as diffusion and subcellular compartmentalization, which have proven to be difficult for DE systems. We demonstrate that large-grained agents are well suited to maintaining interal state representations and efficient in computing reactant concentration, both of which are vital considerations in modeling the ST domain. In our system, a molecular species is modeled as an individual agent with hierarchical task network structures to represent self- and externally-initiated reactions. An agents identity is determined by a rule file (one for every participating molecular species) that specifies the reactions it may participate in, as well as its initial concentration. Reactions within the system are actuated by inter-agent communication. We present results from modeling the well-studied epidermal growth factor (EGF) pathway, demonstrating the viability of MAS technologies as a simulation platform for biological networks.