Local Adaptation in the Flowering-Time Gene Network of Balsam Poplar, Populus balsamifera L

Identifying the signature and targets of local adaptation is an increasingly important goal in empirical population genetics. Using data from 443 balsam poplar Populus balsamifera trees sampled from 31 populations, we tested for evidence of geographically variable selection shaping diversity at 27 homologues of the Arabidopsis flowering-time network. These genes are implicated in the control of seasonal phenology, an important determinant of fitness. Using 335 candidate and 412 reference single nucleotide polymorphisms (SNPs), we tested for evidence of local adaptation by searching for elevated population differentiation using F(ST)-based outlier analyses implemented in BayeScan or a Hierarchical Model in Arelquin and by testing for significant associations between allele frequency and environmental variables using BAYENV. A total of 46 SNPs from 14 candidate genes had signatures of local adaptation-either significantly greater population differentiation or significant covariance with one or more environmental variable relative to reference SNP distributions. Only 11 SNPs from two genes exhibited both elevated population differentiation and covariance with one or more environmental variables. Several genes including the abscisic acid gene ABI1B and the circadian clock genes ELF3 and GI5 harbored a large number of SNPs with signatures of local adaptation-with SNPs in GI5 strongly covarying with both latitude and precipitation and SNPs in ABI1B strongly covarying with temperature. In contrast to several other systems, we find little evidence that photoreceptors, including phytochromes, play an important role in local adaptation. Our results additionally show that detecting local adaptation is sensitive to the analytical approaches used and that model-based significance thresholds should be viewed with caution.     

Insertional mutagenesis in <Emphasis Type="Italic">Populus</Emphasis>: relevance and feasibility

The recent sequencing of the first tree genome, that of the black cottonwood (Populus trichocarpa), opens a new chapter in tree functional genomics. While the completion of the genome is a milestone, mobilizing this significant resource for better understanding the growth and development of woody perennials will be an even greater undertaking in the years to come. In other model organisms, a critical tool for high-throughput analysis of gene function has been the generation of large mutagenized populations. Some mutagenesis technologies and approaches cannot be applied to trees because of their typically outcrossing breeding systems, high heterozygosity, large body size, and delayed flowering. In contrast, gene-tagging approaches that use insertional mutagenesis to create dominant phenotypes are ideally suited for trees and, especially, Populus. Both activation tagging and enhancer trap programs have been successful in identifying new genes important to tree development. The generation of genome-wide insertional mutant populations, which provide direct functional links between genes and phenotypes, should help to integrate in silico analyses of gene and protein expression, association studies of natural genetic polymorphism, and phenotypic analyses of adaptation and development.

Coalescent methods for estimating species trees from phylogenomic data

Genome-scale sequence data have become increasingly available in the phylogenetic studies for understanding the evolutionary histories of species. However, it is challenging to develop probabilistic models to account for heterogeneity of phylogenomic data. The multispecies coalescent model describes gene trees as independent random variables generated from a coalescence process occurring along the lineages of the species tree. Since the multispecies coalescent model allows gene trees to vary across genes, coalescent-based methods have been popularly used to account for heterogeneous gene trees in phylogenomic data analysis. In this paper, we summarize and evaluate the performance of coalescent-based methods for estimating species trees from genome-scale sequence data. We investigate the effects of deep coalescence and mutation on the performance of species tree estimation methods. We found that the coalescent-based methods perform well in estimating species trees for a large number of genes, regardless of the degree of deep coalescence and mutation. The performance of the coalescent methods is negatively correlated with the lengths of internal branches of the species tree.     

Gene disruption by biolistic transformation in serotype D strains of Cryptococcus neoformans

Gene disruption by biolistic transformation in serotype D strains of Cryptococcus neoformans. Fungal Genetics and Biology 29, 38-48. Cryptococcus neoformans is an opportunistic fungal pathogen with a defined sexual cycle and well-developed genetic and molecular approaches. Two different transformation systems have been developed, and a number of genes have been disrupted by homologous recombination. However, the frequency of homologous recombination achieved by these approaches has differed dramatically between strains of the A and D serotypes. Transformation by electroporation in serotype D strains results in homologous recombination at frequencies of 1/1000 to 1/100,000, whereas transformation by the biolistic method has resulted in gene disruption at frequencies between 2 and 50% in serotype A strains. We find that gene disruption by homologous recombination can be achieved in the congenic serotype D strain series by biolistic transformation with frequencies of approximately 1 to 4%. By this approach, we have readily disrupted the genes encoding a MAPK homolog (CPK1), the calcineurin A catalytic subunit (CNA1), and a G protein alpha subunit (GPA1). By physical and genetic methods, we show that these mutations result from targeted recombination events without ectopic integrations. Because genetic approaches can be applied in the congenic serotype D strains, our observations represent a significant advance in molecular approaches to understand the physiology and virulence of this important human pathogen.     

Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits

A central challenge in evolutionary biology is to identify genes underlying ecologically important traits and describe the fitness consequences of naturally occurring variation at these loci. To address this goal, several novel approaches have been developed, including 'population genomics,' where a large number of molecular markers are scored in individuals from different environments with the goal of identifying markers showing unusual patterns of variation, potentially due to selection at linked sites. Such approaches are appealing because of (1) the increasing ease of generating large numbers of genetic markers, (2) the ability to scan the genome without measuring phenotypes and (3) the simplicity of sampling individuals without knowledge of their breeding history. Although such approaches are inherently applicable to non-model systems, to date these studies have been limited in their ability to uncover functionally relevant genes. By contrast, quantitative genetics has a rich history, and more recently, quantitative trait locus (QTL) mapping has had some success in identifying genes underlying ecologically relevant variation even in novel systems. QTL mapping, however, requires (1) genetic markers that specifically differentiate parental forms, (2) a focus on a particular measurable phenotype and (3) controlled breeding and maintenance of large numbers of progeny. Here we present current advances and suggest future directions that take advantage of population genomics and quantitative genetic approaches - in both model and non-model systems. Specifically, we discuss advantages and limitations of each method and argue that a combination of the two provides a powerful approach to uncovering the molecular mechanisms responsible for adaptation.     

Restoration of threatened species: a noble cause for transgenic trees

Some of the first applications of transgenic trees in North America may be for the conservation or restoration of threatened forest trees that have been devastated by fungal pathogens or insect pests. In some cases, where resistance has yet to be found in the natural population of a tree species, incorporating genes from other organisms may offer the only hope for restoration. In others, transgenics may play a role as part of an integrated approach, along with conventional breeding or biocontrol agents. American chestnut (Castanea dentata) was wiped out as a canopy species by a fungal disease accidentally introduced into the United States around 1900. Similarly, American elm (Ulmus americana) virtually disappeared as a favored street tree from Northeastern U.S. cities after the introduction of the Dutch elm disease fungus in the 1940s. In both cases, progress has been made toward restoration via conventional techniques such as selection and propagation of tolerant cultivars (American elm) or breeding with a related resistant species (American chestnut). Recently, progress has also been made with development of systems for engineering antifungal candidate genes into these “heritage trees.” An Agrobacterium-leaf disk system has been used to produce transgenic American elm trees engineered with an antimicrobial peptide gene that may enhance resistance to Dutch elm disease. Two gene transfer systems have been developed for American chestnut using Agrobacterium-mediated transformation of embryogenic cultures, setting the stage for the first tests of potential antifungal genes for their ability to confer resistance to the chestnut blight fungus. Despite the promise of transgenic approaches for restoration of these heritage trees, a number of technical, environmental, economic, and ethical questions remain to be addressed before such trees can be deployed, and the debate around these questions may be quite different from that associated with transgenic trees developed for other purposes.     

Alteration of soil rhizosphere communities following genetic transformation of white spruce

The application of plant genetic manipulations to agriculture and forestry with the aim of alleviating insect damage through Bacillus thuringiensis transformation could lead to a significant reduction in the release of pesticides into the environment. However, many groups have come forward with very valid and important questions related to potentially adverse effects, and it is crucial to assess and better understand the impact that this technology might have on ecosystems. In this study, we analyzed rhizosphere soil samples collected from the first B. thuringiensis-transformed trees [with insertion of the CryIA(b) toxin-encoding gene] grown in Canada (Val-Cartier, QC, Canada) as part of an ecological impact assessment project. Using a robust amplified rRNA gene restriction analysis approach coupled with 16S rRNA gene sequencing, the rhizosphere-inhabiting microbial communities of white spruce (Picea glauca) genetically modified by biolistic insertion of the cryIA(b), uidA (beta-glucuronidase), and nptII genes were compared with the microbial communities associated with non-genetically modified counterparts and with trees in which only the genetic marker genes uidA and nptII have been inserted. Analysis of 1,728 rhizosphere bacterial clones (576 clones per treatment) using a Cramér-von Mises statistic analysis combined with a Monte Carlo comparison clearly indicated that there was a statistically significant difference (P < 0.05) between the microbial communities inhabiting the rhizospheres of trees carrying the cryIA(b), uidA, and nptII transgenes, trees carrying only the uidA and nptII transgenes, and control trees. Clear rhizosphere microbial community alterations due to B. thuringiensis tree genetic modification have to our knowledge never been described previously and open the door to interesting questions related to B. thuringiensis genetic transformation and also to the impact of commonly used uidA and nptII genetic marker genes.     

Genetic engineering of trees: progress and new horizons

Genetic engineering of trees to improve productivity, wood quality, and resistance to biotic and abiotic stresses has been the primary goal of the forest biotechnology community for decades. We review the extensive progress in these areas and their current status with respect to commercial applications. Examples include novel methods for lignin modification, solutions for long-standing problems related to pathogen resistance, modifications to flowering onset and fertility, and drought and freeze tolerance. There have been numerous successful greenhouse and field demonstrations of genetically engineered trees, but commercial application has been severely limited by social and technical considerations. Key social factors are costly and uncertain regulatory hurdles and sweeping market barriers in the form of forest certification systems that disallow genetically modified trees. These factors limit and, in many cases, preclude field research and commercial adoption. Another challenge is the high cost and uncertainty in transformation efficiency that is needed to apply genetic engineering and gene editing methods to most species and genotypes of commercial importance. Recent advances in developmental gene-based transformation systems and gene editing, if combined with regulatory and certification system reform, could provide the foundation for genetic engineering to become a significant tool for coping with the increasing environmental and biological stresses on planted and wild forests.     

High efficiency poplar transformation

With the completion of the poplar tree genome database, Populus species have become one of the most useful model systems for the study of woody plant biology. Populus tremuloides (quaking aspen) is the most wide-spread tree species in North America, and its rapid growth generates the most abundant wood-based biomass out of any other plant species. To study such beneficial traits, there is a need for easier and more efficient transformation procedures that will allow the study of large numbers of tree genes. We have developed transformation procedures that are suitable for high-throughput format transformations using either Agrobacterium tumefaciens to produce transformed trees or Agrobacterium rhizogenes to generate hairy roots. Our method uses Agrobacterium inoculated aspen seedling hypocotyls followed by direct thidiazuron (TDZ)-mediated shoot regeneration on selective media. Transformation was verified through β-glucuronidase (GUS) reporter gene expression in all tree tissues, PCR amplification of appropriate vector products from isolated genomic DNA, and northern hybridization of incorporated and expressed transgenes. The hairy root protocol follows the same inoculation procedures and was tested using GUS reporter gene integration and antibiotic selection. The benefit of these procedures is that they are simple and efficient, requiring no maintenance of starting materials and allowing fully formed transgenic trees (or hairy roots) to be generated in only 3–4 months, rather than the 6–12 months required by more traditional methods. Likewise, the fact that the protocols are amenable to high-throughput formats makes them better suited for large-scale functional genomics studies in poplars.     

Adaptation to climate through flowering phenology: a case study in Medicago truncatula

Local climatic conditions likely constitute an important selective pressure on genes underlying important fitness‐related traits such as flowering time, and in many species, flowering phenology and climatic gradients strongly covary. To test whether climate shapes the genetic variation on flowering time genes and to identify candidate flowering genes involved in the adaptation to environmental heterogeneity, we used a large Medicago truncatula core collection to examine the association between nucleotide polymorphisms at 224 candidate genes and both climate variables and flowering phenotypes. Unlike genome‐wide studies, candidate gene approaches are expected to enrich for the number of meaningful trait associations because they specifically target genes that are known to affect the trait of interest. We found that flowering time mediates adaptation to climatic conditions mainly by variation at genes located upstream in the flowering pathways, close to the environmental stimuli. Variables related to the annual precipitation regime reflected selective constraints on flowering time genes better than the other variables tested (temperature, altitude, latitude or longitude). By comparing phenotype and climate associations, we identified 12 flowering genes as the most promising candidates responsible for phenological adaptation to climate. Four of these genes were located in the known flowering time QTL region on chromosome 7. However, climate and flowering associations also highlighted largely distinct gene sets, suggesting different genetic architectures for adaptation to climate and flowering onset.     

Gene and enhancer trap tagging of vascular-expressed genes in poplar trees

We report a gene discovery system for poplar trees based on gene and enhancer traps. Gene and enhancer trap vectors carrying the beta-glucuronidase (GUS) reporter gene were inserted into the poplar genome via Agrobacterium tumefaciens transformation, where they reveal the expression pattern of genes at or near the insertion sites. Because GUS expression phenotypes are dominant and are scored in primary transformants, this system does not require rounds of sexual recombination, a typical barrier to developmental genetic studies in trees. Gene and enhancer trap lines defining genes expressed during primary and secondary vascular development were identified and characterized. Collectively, the vascular gene expression patterns revealed that approximately 40% of genes expressed in leaves were expressed exclusively in the veins, indicating that a large set of genes is required for vascular development and function. Also, significant overlap was found between the sets of genes responsible for development and function of secondary vascular tissues of stems and primary vascular tissues in other organs of the plant, likely reflecting the common evolutionary origin of these tissues. Chromosomal DNA flanking insertion sites was amplified by thermal asymmetric interlaced PCR and sequenced and used to identify insertion sites by reference to the nascent Populus trichocarpa genome sequence. Extension of the system was demonstrated through isolation of full-length cDNAs for five genes of interest, including a new class of vascular-expressed gene tagged by enhancer trap line cET-1-pop1-145. Poplar gene and enhancer traps provide a new resource that allows plant biologists to directly reference the poplar genome sequence and identify novel genes of interest in forest biology.     

Discovering regulatory and signalling circuits in molecular interaction networks

MOTIVATION: In model organisms such as yeast, large databases of protein-protein and protein-DNA interactions have become an extremely important resource for the study of protein function, evolution, and gene regulatory dynamics. In this paper we demonstrate that by integrating these interactions with widely-available mRNA expression data, it is possible to generate concrete hypotheses for the underlying mechanisms governing the observed changes in gene expression. To perform this integration systematically and at large scale, we introduce an approach for screening a molecular interaction network to identify active subnetworks, i.e., connected regions of the network that show significant changes in expression over particular subsets of conditions. The method we present here combines a rigorous statistical measure for scoring subnetworks with a search algorithm for identifying subnetworks with high score. RESULTS: We evaluated our procedure on a small network of 332 genes and 362 interactions and a large network of 4160 genes containing all 7462 protein-protein and protein-DNA interactions in the yeast public databases. In the case of the small network, we identified five significant subnetworks that covered 41 out of 77 (53%) of all significant changes in expression. Both network analyses returned several top-scoring subnetworks with good correspondence to known regulatory mechanisms in the literature. These results demonstrate how large-scale genomic approaches may be used to uncover signalling and regulatory pathways in a systematic, integrative fashion.     

The dprA gene is required for natural transformation of Helicobacter pylori

Genetic recombination in Helicobacter pylori is believed to be involved in host adaptation of this gastric pathogen and uptake of DNA by natural transformation can result in changes in virulence factors as well as antigenic variation. To elucidate the mechanisms involved in natural transformation we tested two genes with homology to known competence genes (dprA and traG) for their role in this process. Insertion mutants in these genes were constructed in two different H. pylori strains and their competence by natural transformation was compared to the wild-type. Mutation of the traG homolog did not reduce competence. Mutation of the dprA gene, however, severely impaired natural transformation both with plasmid and chromosomal DNA. Our data indicate that dprA and comB3 are essential parts of a common pathway for chromosomal and plasmid transformation.     

Integrative biology and the developing limb bud

The identification or selective construction of mutations within genes has allowed researchers to explore the downstream effects of gene disruption. Although these approaches have been successful, a limitation in our assessment of the consequences of conditional changes, and thereby our understanding of roles or function of genes, limits the degree to which we examine the effects of our manipulations. It is also clear that linear associations are incorrect models for describing development, and newer methods now give us an opportunity to practice an integrative biology. In our attempts to explore the consequences of Hoxa13 disruption in mice and humans, it has become clear that a better understanding of the consequences of gene alteration may be achievable by taking a broader approach with a long-term view. Fundamental questions regarding Hox gene function in vertebrates, including those related to the number of target genes; the degree of overlap of target gene regulation among paralogs; the magnitude of modulation exerted; and the identity of genes that are activated versus repressed need to be explored if a more thorough mechanistic understanding is to be achieved. To begin to address these questions, we undertook a comprehensive analysis of the expression of genes within developing limb buds of mice, and here we present some of our preliminary results. Our efforts will further (1) the exploration of the broader genetic relationships of expressed genes, (2) the determination of parallels or variations in target usage for a given gene in different tissues and between different organisms, (3) the evaluation of limb patterning mechanisms in other animal model systems, and (4) the exploration of gene expression hierarchies regulated by HOX proteins in developmental systems.     

Lentivirus-delivered stable gene silencing by RNAi in primary cells

Genome-wide genetic approaches have proven useful for examining pathways of biological significance in model organisms such as Saccharomyces cerevisiae, Drosophila melanogastor, and Caenorhabditis elegans, but similar techniques have proven difficult to apply to mammalian systems. Although manipulation of the murine genome has led to identification of genes and their function, this approach is laborious, expensive, and often leads to lethal phenotypes. RNA interference (RNAi) is an evolutionarily conserved process of gene silencing that has become a powerful tool for investigating gene function by reverse genetics. Here we describe the delivery of cassettes expressing hairpin RNA targeting green fluorescent protein (GFP) using Moloney leukemia virus-based and lentivirus-based retroviral vectors. Both transformed cell lines and primary dendritic cells, normally refractory to transfection-based gene transfer, demonstrated stable silencing of targeted genes, including the tumor suppressor gene TP53 in normal human fibroblasts. This report demonstrates that both Moloney leukemia virus and lentivirus vector-mediated expression of RNAi can achieve effective, stable gene silencing in diverse biological systems and will assist in elucidating gene functions in numerous cell types including primary cells.     

Comparative genomics of two closely related unicellular thermo-acidophilic red algae, Galdieria sulphuraria and Cyanidioschyzon merolae, reveals the molecular basis of the metabolic flexibility of Galdieria sulphuraria and significant differences in carbohydrate metabolism of both algae

Unicellular algae serve as models for the study and discovery of metabolic pathways, for the functional dissection of cell biological processes such as organellar division and cell motility, and for the identification of novel genes and gene functions. The recent completion of several algal genome sequences and expressed sequence tag collections and the establishment of nuclear and organellar transformation methods has opened the way for functional genomics approaches using algal model systems. The thermo-acidophilic unicellular red alga Galdieria sulphuraria represents a particularly interesting species for a genomics approach owing to its extraordinary metabolic versatility such as heterotrophic and mixotrophic growth on more than 50 different carbon sources and its adaptation to hot acidic environments. However, the ab initio prediction of genes required for unknown metabolic pathways from genome sequences is not trivial. A compelling strategy for gene identification is the comparison of similarly sized genomes of related organisms with different physiologies. Using this approach, candidate genes were identified that are critical to the metabolic versatility of Galdieria. Expressed sequence tags and high-throughput genomic sequence reads covering >70% of the G. sulphuraria genome were compared to the genome of the unicellular, obligate photoautotrophic red alga Cyanidioschyzon merolae. More than 30% of the Galdieria sequences did not relate to any of the Cyanidioschyzon genes. A closer inspection of these sequences revealed a large number of membrane transporters and enzymes of carbohydrate metabolism that are unique to Galdieria. Based on these data, it is proposed that genes involved in the uptake of reduced carbon compounds and enzymes involved in their metabolism are crucial to the metabolic flexibility of G. sulphuraria.     

Quikgene: a gene synthesis method integrated with ligation-free cloning

Gene synthesis is a convenient tool that is widely used to make genes for a variety of purposes. All current protocols essentially take inside-out approaches to assemble complete genes using DNA oligonucleotides or intermediate fragments. Here we present an efficient method that integrates gene synthesis and cloning into one step. Our method, which is evolved from QuikChange mutagenesis, can modify, extend, or even de novo synthesize relatively large genes. The genes are inserted directly into vectors without ligations or subcloning. We de novo synthesized a 600-bp gene through multiple steps of polymerase chain reaction (PCR) directly into a bacterial expression vector. This outside-in gene synthesis method is called Quikgene. Furthermore, we have defined an overlap region of a minimum of nine nucleotides in insertion primers that is sufficient enough to circularize PCR products for efficient transformation, allowing one to significantly reduce the lengths of primers. Taken together, our protocol greatly extends the current length limit for QuikChange insertion. More importantly, it combines gene synthesis and cloning into one step. It has potential applications for high-throughput structural genomics.