Daniel S. Rokhsar

The Sorghum bicolor genome and the diversification of grasses

Sorghum, an African grass related to sugar cane and maize, is grown for food, feed, fibre and fuel. We present an initial analysis of the ∼730-megabase Sorghum bicolor (L.) Moench genome, placing ∼98% of genes in their chromosomal context using whole-genome shotgun sequence validated by genetic, physical and syntenic information. Genetic recombination is largely confined to about one-third of the sorghum genome with gene order and density similar to those of rice. Retrotransposon accumulation in recombinationally recalcitrant heterochromatin explains the ∼75% larger genome size of sorghum compared with rice. Although gene and repetitive DNA distributions have been preserved since palaeopolyploidization ∼70 million years ago, most duplicated gene sets lost one member before the sorghum–rice divergence. Concerted evolution makes one duplicated chromosomal segment appear to be only a few million years old. About 24% of genes are grass-specific and 7% are sorghum-specific. Recent gene and microRNA duplications may contribute to sorghum’s drought tolerance.

The genome of the choanoflagellate Monosiga brevicollis and the origin of metazoans

Choanoflagellates are the closest known relatives of metazoans. To discover potential molecular mechanisms underlying the evolution of metazoan multicellularity, we sequenced and analysed the genome of the unicellular choanoflagellate Monosiga brevicollis. The genome contains approximately 9,200 intron-rich genes, including a number that encode cell adhesion and signalling protein domains that are otherwise restricted to metazoans. Here we show that the physical linkages among protein domains often differ between M. brevicollis and metazoans, suggesting that abundant domain shuffling followed the separation of the choanoflagellate and metazoan lineages. The completion of the M. brevicollis genome allows us to reconstruct with increasing resolution the genomic changes that accompanied the origin of metazoans.

Genome sequence of the lignocellulose degrading fungus Phanerochaete chrysosporium strain RP78

White rot fungi efficiently degrade lignin, a complex aromatic polymer in wood that is among the most abundant natural materials on earth. These fungi use extracellular oxidative enzymes that are also able to transform related aromatic compounds found in explosive contaminants, pesticides and toxic waste. We have sequenced the 30-million base-pair genome of Phanerochaete chrysosporium strain RP78 using a whole genome shotgun approach. The P. chrysosporium genome reveals an impressive array of genes encoding secreted oxidases, peroxidases and hydrolytic enzymes that cooperate in wood decay. Analysis of the genome data will enhance our understanding of lignocellulose degradation, a pivotal process in the global carbon cycle, and provide a framework for further development of bioprocesses for biomass utilization, organopollutant degradation and fiber bleaching. This genome provides a high quality draft sequence of a basidiomycete, a major fungal phylum that includes important plant and animal pathogens.

A genomewide survey of developmentally relevant genes in Ciona intestinalis. II. Genes for homeobox transcription factors.

Homeobox-containing genes play crucial roles in various developmental processes, including body-plan specification, pattern formation and cell-type specification. The present study searched the draft genome sequence and cDNA/EST database of the basal chordate Ciona intestinalis to identify 83 homeobox-containing genes in this animal. This number of homeobox genes in the Ciona genome is smaller than that in the Caenorhabditis elegans, Drosophila melanogaster, human and mouse genomes. Of the 83 genes, 76 have possible human orthologues and 7 may be unique to Ciona. The ascidian homeobox genes were classified into 11 classes, including Hox class, NK class, Paired class, POU class, LIM class, TALE class, SIX class, Prox class, Cut class, ZFH class and HNF1 class, according to the classification scheme devised for known homeobox genes. As to the Hox cluster, the Ciona genome contains single copies of each of the paralogous groups, suggesting that there is a single Hox cluster, if any, but genes orthologous to Hox7, 8, 9 and 11 were not found in the genome. In addition, loss of genes had occurred independently in the Ciona lineage and was noticed in Gbx of the EHGbox subclass, Sax, NK3, Vax and vent of the NK class, Cart, Og9, Anf and Mix of the Paired class, POU-I, III, V and VI of the POU class, Lhx6/7 of the LIM class, TGIF of the TALE class, Cux and SATB of the Cut class, and ZFH1 of the ZFH class, which might have reduced the number of Ciona homeobox genes. Interestingly, one of the newly identified Ciona intestinalis genes and its vertebrate counterparts constitute a novel subclass of HNF1 class homeobox genes. Furthermore, evidence for the gene structures and expression of 54 of the 83 homeobox genes was provided by analysis of ESTs, suggesting that cDNAs for these 54 genes are available. The present data thus reveal the repertoire of homeodomain-containing transcription factors in the Ciona genome, which will be useful for future research on the development and evolution of chordates.

Genomic analysis of immunity in a Urochordate and the emergence of the vertebrate immune system: “waiting for Godot”

Genome-wide sequence analysis in the invertebrate chordate, Ciona intestinalis, has provided a comprehensive picture of immune-related genes in an organism that occupies a key phylogenetic position in vertebrate evolution. The pivotal genes for adaptive immunity, such as the major histocompatibility complex (MHC) classxa0I and II genes, T-cell receptors, or dimeric immunoglobulin molecules, have not been identified in the Ciona genome. Many genes involved in innate immunity have been identified, including complement components, Toll-like receptors, and the genes involved in intracellular signal transduction of immune responses, and show both expansion and unexpected diversity in comparison with the vertebrates. In addition, a number of genes were identified which predicted integral membrane proteins with extracellular C-type lectin or immunoglobulin domains and intracellular immunoreceptor tyrosine-based inhibitory motifs (ITIMs) and immunoreceptor tyrosine-based activation motifs (ITAMs) (plus their associated signal transduction molecules), suggesting that activating and inhibitory receptors have an MHC-independent function and an early evolutionary origin. A crucial component of vertebrate adaptive immunity is somatic diversification, and the recombination activating genes (RAG) and activation-induced cytidine deaminase (AID) genes responsible for the Generation of diversity are not present in Ciona. However, there are key V regions, the essential feature of an immunoglobulin superfamily VC1-like core, and possible proto-MHC regions scattered throughout the genome waiting for Godot.

The DNA sequence and biology of human chromosome 19

Chromosome 19 has the highest gene density of all human chromosomes, more than double the genome-wide average. The large clustered gene families, corresponding high G + C content, CpG islands and density of repetitive DNA indicate a chromosome rich in biological and evolutionary significance. Here we describe 55.8 million base pairs of highly accurate finished sequence representing 99.9% of the euchromatin portion of the chromosome. Manual curation of gene loci reveals 1,461 protein-coding genes and 321 pseudogenes. Among these are genes directly implicated in mendelian disorders, including familial hypercholesterolaemia and insulin-resistant diabetes. Nearly one-quarter of these genes belong to tandemly arranged families, encompassing more than 25% of the chromosome. Comparative analyses show a fascinating picture of conservation and divergence, revealing large blocks of gene orthology with rodents, scattered regions with more recent gene family expansions and deletions, and segments of coding and non-coding conservation with the distant fish species Takifugu.

Insights into the red algae and eukaryotic evolution from the genome of Porphyra umbilicalis (Bangiophyceae, Rhodophyta)

Significance Fossil evidence shows that red algae (Rhodophyta) are one of the most ancient multicellular lineages. Their ecological, evolutionary, and commercial importance notwithstanding, few red algal nuclear genomes have been sequenced. Our analyses of the Porphyra umbilicalis genome provide insights into how this macrophyte thrives in the stressful intertidal zone and into the basis for its nutritional value as human food. Many of the novel traits (e.g., cytoskeletal organization, calcium signaling pathways) we find encoded in the Porphyra genome are extended to other red algal genomes, and our unexpected findings offer a potential explanation for why the red algae are constrained to small stature relative to other multicellular lineages. Porphyra umbilicalis (laver) belongs to an ancient group of red algae (Bangiophyceae), is harvested for human food, and thrives in the harsh conditions of the upper intertidal zone. Here we present the 87.7-Mbp haploid Porphyra genome (65.8% G + C content, 13,125 gene loci) and elucidate traits that inform our understanding of the biology of red algae as one of the few multicellular eukaryotic lineages. Novel features of the Porphyra genome shared by other red algae relate to the cytoskeleton, calcium signaling, the cell cycle, and stress-tolerance mechanisms including photoprotection. Cytoskeletal motor proteins in Porphyra are restricted to a small set of kinesins that appear to be the only universal cytoskeletal motors within the red algae. Dynein motors are absent, and most red algae, including Porphyra, lack myosin. This surprisingly minimal cytoskeleton offers a potential explanation for why red algal cells and multicellular structures are more limited in size than in most multicellular lineages. Additional discoveries further relating to the stress tolerance of bangiophytes include ancestral enzymes for sulfation of the hydrophilic galactan-rich cell wall, evidence for mannan synthesis that originated before the divergence of green and red algae, and a high capacity for nutrient uptake. Our analyses provide a comprehensive understanding of the red algae, which are both commercially important and have played a major role in the evolution of other algal groups through secondary endosymbioses.

Development of a toolbox to dissect host-endosymbiont interactions and protein trafficking in the trypanosomatid Angomonas deanei

BackgroundBacterial endosymbionts are found across the eukaryotic kingdom and profoundly impacted eukaryote evolution. In many endosymbiotic associations with vertically inherited symbionts, highly complementary metabolic functions encoded by host and endosymbiont genomes indicate integration of metabolic processes between the partner organisms. While endosymbionts were initially expected to exchange only metabolites with their hosts, recent evidence has demonstrated that also host-encoded proteins can be targeted to the bacterial symbionts in various endosymbiotic systems. These proteins seem to participate in regulating symbiont growth and physiology. However, mechanisms required for protein targeting and the specific endosymbiont targets of these trafficked proteins are currently unexplored owing to a lack of molecular tools that enable functional studies of endosymbiotic systems.ResultsHere we show that the trypanosomatid Angomonas deanei, which harbors a β-proteobacterial endosymbiont, is readily amenable to genetic manipulation. Its rapid growth, availability of full genome and transcriptome sequences, ease of transfection, and high frequency of homologous recombination have allowed us to stably integrate transgenes into the A. deanei nuclear genome, efficiently generate null mutants, and elucidate protein localization by heterologous expression of a fluorescent protein fused to various putative targeting signals. Combining these novel tools with proteomic analysis was key for demonstrating the routing of a host-encoded protein to the endosymbiont, suggesting the existence of a specific endosymbiont-sorting machinery in A. deanei.ConclusionsAfter previous reports from plants, insects, and a cercozoan amoeba we found here that also in A. deanei, i.e. a member of a fourth eukaryotic supergroup, host-encoded proteins can be routed to the bacterial endosymbiont. This finding adds further evidence to our view that the targeting of host proteins is a general strategy of eukaryotes to gain control over and interact with a bacterial endosymbiont. The molecular resources reported here establish A. deanei as a time and cost efficient reference system that allows for a rigorous dissection of host-symbiont interactions that have been, and are still being shaped over evolutionary time. We expect this system to greatly enhance our understanding of the biology of endosymbiosis.

Adaptation to heavy-metal contaminated environments proceeds via selection on pre-existing genetic variation

Across a species range, islands of stressful habitats impose similar selection pressures on isolated populations. It is as yet unclear, when populations respond to these selective pressures, the extent to which this results in convergent genetic evolution and whether convergence is due to independent mutations or shared ancestral variation. We address these questions investigating a classic example of adaptation by natural selection - the colonization of plant species to heavy metal contaminated soils. We use field-based reciprocal transplant experiments to demonstrate that mine alleles at a major copper tolerance QTL, Tol1, are strongly selected in the mine environment, but are neutral, or nearly so, in the off-mine environment. To identify scaffolds in genetic linkage with this locus, we assemble the genome of a mine adapted M. guttatus genotype and sequence near isogenic lines (NILs) homozygous for tolerant or non-tolerant alleles at Tol1. We identify genes with differential expression between NILs and differences in allele frequency between independent pairs of mine and off-mine populations to identify Tol1 candidate genes. We identify a single gene, a multicopper oxidase, with large differences in expression between NILs and allele frequency between populations. Furthermore, we find patterns of genetic variation at Tol1, and four additional candidate adaptation loci, are consistent with selection acting upon beneficial haplotypes that predates the existence of the copper mine habitat. We estimate the age of selected Tol1 haplotype to be at least 1700 years old and was at a frequency of 0.4-0.6% in the ancestral population when mining was initiated 150 years ago. These results suggest that adaptation to the mine habitat routinely occurs via selection on ancestral variation, rather than independent de-novo mutations or migration between populations.

Populus resequencing: towards genome-wide association studies

Genome-wide association studies (GWAS) have been used to identify regions of the genome related to various phenotypes in humans, corn, rice and cattle. Successful application of this approach to bioenergy crops such as Populus requires 1) an appropriate mapping population, 2) high-quality phenotypic data and 3) informative genotypic data. With the goal of reducing the recalcitrance of lignocellulosic biomass for economic production of biofuels and understanding basic mechanisms of cell wall formation in Populus we established 4 clonally replicated common gardens experiments each with 1100 native P. trichocarpa genotypes collected from along the northwest coast of the U.S. and Canada. A high-throughput phenotyping pipeline was developed to measure cell wall chemistry, pretreatment response and enzymatic sugar release. Initially 18 genotypes were resequenced to an average 30X depth in order to design a SNP array to test for statistical association using MMAX and PCA methods of testing among ca. 2500 candidate genes. Genetic structure and linkage disequilibrium (LD) was assessed using SSR and SNP markers. Outlying genotypes were excluded from the analyses and estimates of LD were used to design the bead array. Candidate genes were selected based on QTL intervals, expression profiling within developing xylem and expert opinion. MMAX and PCA results revealed similar significant associations for all measured phenotypes and several SNPs within the candidate gene set explain a relatively high degree of the phenotypic variance. As a result, resequencing has continued in order to conduct GWAS in Populus; the complete set of 1100 genotypes will be complete in 2012.