Robert M. Stupar

Selection-free zinc-finger-nuclease engineering by context-dependent assembly (CoDA)

Engineered zinc-finger nucleases (ZFNs) enable targeted genome modification. Here we describe context-dependent assembly (CoDA), a platform for engineering ZFNs using only standard cloning techniques or custom DNA synthesis. Using CoDA-generated ZFNs, we rapidly altered 20 genes in Danio rerio, Arabidopsis thaliana and Glycine max. The simplicity and efficacy of CoDA will enable broad adoption of ZFN technology and make possible large-scale projects focused on multigene pathways or genome-wide alterations.

Cis-transcriptional variation in maize inbred lines B73 and Mo17 leads to additive expression patterns in the F1 hybrid.

Microarray analysis of gene expression patterns in immature ear, seedling, and embryo tissues from the maize inbred lines B73 and Mo17 identified numerous genes with variable expression. Some genes had detectable expression in only one of the two inbreds; most of these genes were detected in the genomic DNA of both inbreds, indicating that the expression differences are likely caused by differential regulation rather than by differences in gene content. Gene expression was also monitored in the reciprocal F1 hybrids B73 × Mo17 and Mo17 × B73. The reciprocal F1 hybrid lines did not display parental effects on gene expression levels. Approximately 80% of the differentially expressed genes displayed additive expression patterns in the hybrids relative to the inbred parents. The ∼20% of genes that display nonadditive expression patterns tend to be expressed at levels within the parental range, with minimal evidence for novel expression levels greater than the high parent or less than the low parent. Analysis of allele-specific expression patterns in the hybrid suggested that intraspecific variation in gene expression levels is largely attributable to cis-regulatory variation in maize. Collectively, our data suggest that allelic cis-regulatory variation between B73 and Mo17 dictates maintenance of inbred allelic expression levels in the F1 hybrid, resulting in additive expression patterns.

Complex mtDNA constitutes an approximate 620-kb insertion on Arabidopsis thaliana chromosome 2: Implication of potential sequencing errors caused by large-unit repeats

Previously conducted sequence analysis of Arabidopsis thaliana (ecotype Columbia-0) reported an insertion of 270-kb mtDNA into the pericentric region on the short arm of chromosome 2. DNA fiber-based fluorescence in situ hybridization analyses reveal that the mtDNA insert is 618 ± 42 kb, ≈2.3 times greater than that determined by contig assembly and sequencing analysis. Portions of the mitochondrial genome previously believed to be absent were identified within the insert. Sections of the mtDNA are repeated throughout the insert. The cytological data illustrate that DNA contig assembly by using bacterial artificial chromosomes tends to produce a minimal clone path by skipping over duplicated regions, thereby resulting in sequencing errors. We demonstrate that fiber-fluorescence in situ hybridization is a powerful technique to analyze large repetitive regions in the higher eukaryotic genomes and is a valuable complement to ongoing large genome sequencing projects.

Reciprocal Silencing, Transcriptional Bias and Functional Divergence of Homeologs in Polyploid Cotton (Gossypium)

Polyploidy is an important force in the evolution of flowering plants. Genomic merger and doubling induce an extensive array of genomic effects, including immediate and long-term alterations in the expression of duplicate genes (“homeologs”). Here we employed a novel high-resolution, genome-specific, mass-spectrometry technology and a well-established phylogenetic framework to investigate relative expression levels of each homeolog for 63 gene pairs in 24 tissues in naturally occurring allopolyploid cotton (Gossypium L.), a synthetic allopolyploid of the same genomic composition, and models of the diploid progenitor species. Results from a total of 2177 successful expression assays permitted us to determine the extent of expression evolution accompanying genomic merger of divergent diploid parents, genome doubling, and genomic coevolution in a common nucleus subsequent to polyploid formation. We demonstrate that 40% of homeologs are transcriptionally biased in at least one stage of cotton development, that genome merger per se has a large effect on relative expression of homeologs, and that the majority of these alterations are caused by cis-regulatory divergence between the diploid progenitors. We describe the scope of transcriptional subfunctionalization and 15 cases of probable neofunctionalization among 8 tissues. To our knowledge, this study represents the first characterization of transcriptional neofunctionalization in an allopolyploid. These results provide a novel temporal perspective on expression evolution of duplicate genomes and add to our understanding of the importance of polyploidy in plants.

Targeted Mutagenesis of Duplicated Genes in Soybean with Zinc-Finger Nucleases

We performed targeted mutagenesis of a transgene and nine endogenous soybean (Glycine max) genes using zinc-finger nucleases (ZFNs). A suite of ZFNs were engineered by the recently described context-dependent assembly platform—a rapid, open-source method for generating zinc-finger arrays. Specific ZFNs targeting DICER-LIKE (DCL) genes and other genes involved in RNA silencing were cloned into a vector under an estrogen-inducible promoter. A hairy-root transformation system was employed to investigate the efficiency of ZFN mutagenesis at each target locus. Transgenic roots exhibited somatic mutations localized at the ZFN target sites for seven out of nine targeted genes. We next introduced a ZFN into soybean via whole-plant transformation and generated independent mutations in the paralogous genes DCL4a and DCL4b. The dcl4b mutation showed efficient heritable transmission of the ZFN-induced mutation in the subsequent generation. These findings indicate that ZFN-based mutagenesis provides an efficient method for making mutations in duplicate genes that are otherwise difficult to study due to redundancy. We also developed a publicly accessible Web-based tool to identify sites suitable for engineering context-dependent assembly ZFNs in the soybean genome.

Phenotypic and Transcriptomic Changes Associated With Potato Autopolyploidization

Polyploidy is remarkably common in the plant kingdom and polyploidization is a major driving force for plant genome evolution. Polyploids may contain genomes from different parental species (allopolyploidy) or include multiple sets of the same genome (autopolyploidy). Genetic and epigenetic changes associated with allopolyploidization have been a major research subject in recent years. However, we know little about the genetic impact imposed by autopolyploidization. We developed a synthetic autopolyploid series in potato (Solanum phureja) that includes one monoploid (1x) clone, two diploid (2x) clones, and one tetraploid (4x) clone. Cell size and organ thickness were positively correlated with the ploidy level. However, the 2x plants were generally the most vigorous and the 1x plants exhibited less vigor compared to the 2x and 4x individuals. We analyzed the transcriptomic variation associated with this autopolyploid series using a potato cDNA microarray containing ∼9000 genes. Statistically significant expression changes were observed among the ploidies for ∼10% of the genes in both leaflet and root tip tissues. However, most changes were associated with the monoploid and were within the twofold level. Thus, alteration of ploidy caused subtle expression changes of a substantial percentage of genes in the potato genome. We demonstrated that there are few genes, if any, whose expression is linearly correlated with the ploidy and can be dramatically changed because of ploidy alteration.

Structural variants in the soybean genome localize to clusters of biotic stress response genes

Genome-wide structural and gene content variations are hypothesized to drive important phenotypic variation within a species. Structural and gene content variations were assessed among four soybean (Glycine max) genotypes using array hybridization and targeted resequencing. Many chromosomes exhibited relatively low rates of structural variation (SV) among genotypes. However, several regions exhibited both copy number and presence-absence variation, the most prominent found on chromosomes 3, 6, 7, 16, and 18. Interestingly, the regions most enriched for SV were specifically localized to gene-rich regions that harbor clustered multigene families. The most abundant classes of gene families associated with these regions were the nucleotide-binding and receptor-like protein classes, both of which are important for plant biotic defense. The colocalization of SV with plant defense response signal transduction pathways provides insight into the mechanisms of soybean resistance gene evolution and may inform the development of new approaches to resistance gene cloning.

Gene expression analyses in maize inbreds and hybrids with varying levels of heterosis

BackgroundHeterosis is the superior performance of F1 hybrid progeny relative to the parental phenotypes. Maize exhibits heterosis for a wide range of traits, however the magnitude of heterosis is highly variable depending on the choice of parents and the trait(s) measured. We have used expression profiling to determine whether the level, or types, of non-additive gene expression vary in maize hybrids with different levels of genetic diversity or heterosis.ResultsWe observed that the distributions of better parent heterosis among a series of 25 maize hybrids generally do not exhibit significant correlations between different traits. Expression profiling analyses for six of these hybrids, chosen to represent diversity in genotypes and heterosis responses, revealed a correlation between genetic diversity and transcriptional variation. The majority of differentially expressed genes in each of the six different hybrids exhibited additive expression patterns, and ~25% exhibited statistically significant non-additive expression profiles. Among the non-additive profiles, ~80% exhibited hybrid expression levels between the parental levels, ~20% exhibited hybrid expression levels at the parental levels and ~1% exhibited hybrid levels outside the parental range.ConclusionWe have found that maize inbred genetic diversity is correlated with transcriptional variation. However, sampling of seedling tissues indicated that the frequencies of additive and non-additive expression patterns are very similar across a range of hybrid lines. These findings suggest that heterosis is probably not a consequence of higher levels of additive or non-additive expression, but may be related to transcriptional variation between parents. The lack of correlation between better parent heterosis levels for different traits suggests that transcriptional diversity at specific sets of genes may influence heterosis for different traits.

Phenotypic and genomic analyses of a fast neutron mutant population resource in soybean

Mutagenized populations have become indispensable resources for introducing variation and studying gene function in plant genomics research. In this study, fast neutron (FN) radiation was used to induce deletion mutations in the soybean (Glycine max) genome. Approximately 120,000 soybean seeds were exposed to FN radiation doses of up to 32 Gray units to develop over 23,000 independent M2 lines. Here, we demonstrate the utility of this population for phenotypic screening and associated genomic characterization of striking and agronomically important traits. Plant variation was cataloged for seed composition, maturity, morphology, pigmentation, and nodulation traits. Mutants that showed significant increases or decreases in seed protein and oil content across multiple generations and environments were identified. The application of comparative genomic hybridization (CGH) to lesion-induced mutants for deletion mapping was validated on a midoleate x-ray mutant, M23, with a known FAD2-1A (for fatty acid desaturase) gene deletion. Using CGH, a subset of mutants was characterized, revealing deletion regions and candidate genes associated with phenotypes of interest. Exome resequencing and sequencing of PCR products confirmed FN-induced deletions detected by CGH. Beyond characterization of soybean FN mutants, this study demonstrates the utility of CGH, exome sequence capture, and next-generation sequencing approaches for analyses of mutant plant genomes. We present this FN mutant soybean population as a valuable public resource for future genetic screens and functional genomics research.

The composition and origins of genomic variation among individuals of the soybean reference cultivar Williams 82

Soybean (Glycine max) is a self-pollinating species that has relatively low nucleotide polymorphism rates compared with other crop species. Despite the low rate of nucleotide polymorphisms, a wide range of heritable phenotypic variation exists. There is even evidence for heritable phenotypic variation among individuals within some cultivars. Williams 82, the soybean cultivar used to produce the reference genome sequence, was derived from backcrossing a Phytophthora root rot resistance locus from the donor parent Kingwa into the recurrent parent Williams. To explore the genetic basis of intracultivar variation, we investigated the nucleotide, structural, and gene content variation of different Williams 82 individuals. Williams 82 individuals exhibited variation in the number and size of introgressed Kingwa loci. In these regions of genomic heterogeneity, the reference Williams 82 genome sequence consists of a mosaic of Williams and Kingwa haplotypes. Genomic structural variation between Williams and Kingwa was maintained between the Williams 82 individuals within the regions of heterogeneity. Additionally, the regions of heterogeneity exhibited gene content differences between Williams 82 individuals. These findings show that genetic heterogeneity in Williams 82 primarily originated from the differential segregation of polymorphic chromosomal regions following the backcross and single-seed descent generations of the breeding process. We conclude that soybean haplotypes can possess a high rate of structural and gene content variation, and the impact of intracultivar genetic heterogeneity may be significant. This detailed characterization will be useful for interpreting soybean genomic data sets and highlights important considerations for research communities that are developing or utilizing a reference genome sequence.