Genome sequencing and population genomic analyses provide insights into the adaptive landscape of silver birch

Population structure of Betula albosinensis and Betula platyphylla: evidence for hybridization and a cryptic lineage

BACKGROUND AND AIMS

Differences in local abundance and ploidy level are predicted to impact the direction of introgression between species. Here, we tested these hypotheses on populations of Betula albosinensis (red birch) and Betula platyphylla (white birch) which were thought to differ in ploidy level, the former being tetraploid and the latter diploid.

METHODS

We sampled 391 birch individuals from nine localities in China, and classified them into species based on leaf morphology. Twelve nuclear microsatellite markers were genotyped in each sample, and analysed using principal coordinates analysis and STRUCTURE software. We compared the effects of two different methods of scoring polyploid genotypes on population genetic analyses. We analysed the effect of ploidy, local species abundance and latitude on levels of introgression between the species.

KEY RESULTS

Leaf morphology divided our samples into red and white birch, but genetic analyses unexpectedly revealed two groups within red birch, one of which was tetraploid, as expected, but the other of which appeared to have diploid microsatellite genotypes. Five individuals were identified as early-generation hybrids or backcrosses between white birch and red birch and five were identified between red birch and 'diploid' red birch. Cline analysis showed that levels of admixture were not significantly correlated with latitude. Estimated genetic differentiation among species was not significantly different between determined tetraploid and undetermined tetraploid genotypes.

CONCLUSIONS

Limited hybridization and gene flow have occurred between red birch and white birch. Relative species abundance and ploidy level do not impact the direction of introgression between them, as genetic admixture is roughly symmetrical. We unexpectedly found populations of apparently diploid red birch and this taxon may be a progenitor of allotetraploid red birch populations. Incomplete lineage sorting may explain patterns of genetic admixture between apparently diploid and allotetraploid red birch.

Genome assembly of Nannochloropsis oceanica provides evidence of host nucleus overthrow by the symbiont nucleus during speciation

The species of the genus Nannochloropsis are unique in their maintenance of a nucleus-plastid continuum throughout their cell cycle, non-motility and asexual reproduction. These characteristics should have been endorsed in their gene assemblages (genomes). Here we show that N. oceanica has a genome of 29.3 Mb consisting of 32 pseudochromosomes and containing 7,330 protein-coding genes; and the host nucleus may have been overthrown by an ancient red alga symbiont nucleus during speciation through secondary endosymbiosis. In addition, N. oceanica has lost its flagella and abilities to undergo meiosis and sexual reproduction, and adopted a genome reduction strategy during speciation. We propose that N. oceanica emerged through the active fusion of a host protist and a photosynthesizing ancient red alga and the symbiont nucleus became dominant over the host nucleus while the chloroplast was wrapped by two layers of endoplasmic reticulum. Our findings evidenced an alternative speciation pathway of eukaryotes.

Loss of GLK1 transcription factor function reveals new insights in chlorophyll biosynthesis and chloroplast development

Birch (Betula platyphylla × B. pendula) is an important tree for landscaping due to its attractive white bark and straight trunk. In this study, we characterized a T-DNA yellow-green leaf mutant, yl. We identified six insertion sites (ISs) in the mutant by genome resequencing and found a 40-kb deletion containing BpGLK1 around IS2 on chromosome 2. Complementation experiments with the yl mutant and repression of BpGLK1 in wild-type plants confirmed that BpGLK1 was responsible for the mutated phenotype. Physiological and ultrastructural analyses showed that the leaves of the yl mutant and BpGLK1-repression lines had decreased chlorophyll content and defective chloroplast development compared to the wild-type. Furthermore, the loss function of BpGLK1 also affected photosynthesis in leaves. Transcriptomics, proteomics, and ChIP-PCR analysis revealed that BpGLK1 directly interacted with the promoter of genes related to antenna proteins, chlorophyll biosynthesis, and photosystem subunit synthesis, and regulated their expression. Overall, our research not only provides new insights into the mechanism of chloroplast development and chlorophyll biosynthesis regulated by BpGLK1, but also provides new transgenic birch varieties with various levels of yellowing leaves by repressing BpGLK1 expression.

Tissue-specific study across the stem reveals the chemistry and transcriptome dynamics of birch bark

Tree bark is a highly specialized array of tissues that plays important roles in plant protection and development. Bark tissues develop from two lateral meristems; the phellogen (cork cambium) produces the outermost stem-environment barrier called the periderm, while the vascular cambium contributes with phloem tissues. Although bark is diverse in terms of tissues, functions and species, it remains understudied at higher resolution. We dissected the stem of silver birch (Betula pendula) into eight major tissue types, and characterized these by a combined transcriptomics and metabolomics approach. We further analyzed the varying bark types within the Betulaceae family. The two meristems had a distinct contribution to the stem transcriptomic landscape. Furthermore, inter- and intraspecies analyses illustrated the unique molecular profile of the phellem. We identified multiple tissue-specific metabolic pathways, such as the mevalonate/betulin biosynthesis pathway, that displayed differential evolution within the Betulaceae. A detailed analysis of suberin and betulin biosynthesis pathways identified a set of underlying regulators and highlighted the important role of local, small-scale gene duplication events in the evolution of metabolic pathways. This work reveals the transcriptome and metabolic diversity among bark tissues and provides insights to its development and evolution, as well as its biotechnological applications.

The genomes of pecan and Chinese hickory provide insights into Carya evolution and nut nutrition

BACKGROUND

Pecan (Carya illinoinensis) and Chinese hickory (C. cathayensis) are important commercially cultivated nut trees in the genus Carya (Juglandaceae), with high nutritional value and substantial health benefits.

RESULTS

We obtained >187.22 and 178.87 gigabases of sequence, and ∼288× and 248× genome coverage, to a pecan cultivar ("Pawnee") and a domesticated Chinese hickory landrace (ZAFU-1), respectively. The total assembly size is 651.31 megabases (Mb) for pecan and 706.43 Mb for Chinese hickory. Two genome duplication events before the divergence from walnut were found in these species. Gene family analysis highlighted key genes in biotic and abiotic tolerance, oil, polyphenols, essential amino acids, and B vitamins. Further analyses of reduced-coverage genome sequences of 16 Carya and 2 Juglans species provide additional phylogenetic perspective on crop wild relatives.

CONCLUSIONS

Cooperative characterization of these valuable resources provides a window to their evolutionary development and a valuable foundation for future crop improvement.

Exome sequencing highlights the role of wild-relative introgression in shaping the adaptive landscape of the wheat genome

Introgression is a potential source of beneficial genetic diversity. The contribution of introgression to adaptive evolution and improvement of wheat as it was disseminated worldwide remains unknown. We used targeted re-sequencing of 890 diverse accessions of hexaploid and tetraploid wheat to identify wild-relative introgression. Introgression, and selection for improvement and environmental adaptation, each reduced deleterious allele burden. Introgression increased diversity genome wide and in regions harboring major agronomic genes, and contributed alleles explaining a substantial proportion of phenotypic variation. These results suggest that historic gene flow from wild relatives made a substantial contribution to the adaptive diversity of modern bread wheat.

Comprehensive characterization of T-DNA integration induced chromosomal rearrangement in a birch T-DNA mutant

Background

Integration of T-DNA into plant genomes via Agrobacterium may interrupt gene structure and generate numerous mutants. The T-DNA caused mutants are valuable materials for understanding T-DNA integration model in plant research. T-DNA integration in plants is complex and still largely unknown. In this work, we reported that multiple T-DNA fragments caused chromosomal translocation and deletion in a birch (Betula platyphylla × B. pendula) T-DNA mutant yl.

Results

We performed PacBio genome resequencing for yl and the result revealed that two ends of a T-DNA can be integrated into plant genome independently because the two ends can be linked to different chromosomes and cause chromosomal translocation. We also found that these T-DNA were connected into tandem fragment regardless of direction before integrating into plant genome. In addition, the integration of T-DNA in yl genome also caused several chromosomal fragments deletion. We then summarized three cases for T-DNA integration model in the yl genome. (1) A T-DNA fragment is linked to the two ends of a double-stranded break (DSB); (2) Only one end of a T-DNA fragment is linked to a DSB; (3) A T-DNA fragment is linked to the ends of different DSBs. All the observations in the yl genome supported the DSB repair model.

Conclusions

In this study, we showed a comprehensive genome analysis of a T-DNA mutant and provide a new insight into T-DNA integration in plants. These findings would be helpful for the analysis of T-DNA mutants with special phenotypes.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5636-y) contains supplementary material, which is available to authorized users.

Kermit: linkage map guided long read assembly

BACKGROUND

With long reads getting even longer and cheaper, large scale sequencing projects can be accomplished without short reads at an affordable cost. Due to the high error rates and less mature tools, de novo assembly of long reads is still challenging and often results in a large collection of contigs. Dense linkage maps are collections of markers whose location on the genome is approximately known. Therefore they provide long range information that has the potential to greatly aid in de novo assembly. Previously linkage maps have been used to detect misassemblies and to manually order contigs. However, no fully automated tools exist to incorporate linkage maps in assembly but instead large amounts of manual labour is needed to order the contigs into chromosomes.

RESULTS

We formulate the genome assembly problem in the presence of linkage maps and present the first method for guided genome assembly using linkage maps. Our method is based on an additional cleaning step added to the assembly. We show that it can simplify the underlying assembly graph, resulting in more contiguous assemblies and reducing the amount of misassemblies when compared to de novo assembly.

CONCLUSIONS

We present the first method to integrate linkage maps directly into genome assembly. With a modest increase in runtime, our method improves contiguity and correctness of genome assembly.

Pollens destroy respiratory epithelial cell anchors and drive alphaherpesvirus infection

Pollens are well-known triggers of respiratory allergies and asthma. The pollen burden in today's ambient air is constantly increasing due to rising climate change and air pollution. How pollens interact with the respiratory mucosa remains largely unknown due to a lack of representative model systems. We here demonstrate how pollen proteases of Kentucky bluegrass, white birch and hazel selectively destroy integrity and anchorage of columnar respiratory epithelial cells, but not of basal cells, in both ex vivo respiratory mucosal explants and in vitro primary equine respiratory epithelial cells (EREC). In turn, this pollen protease-induced damage to respiratory epithelial cell anchorage resulted in increased infection by the host-specific and ancestral alphaherpesvirus equine herpesvirus type 1 (EHV1). Pollen proteases of all three plant species were characterized by zymography and those of white birch were fully identified for the first time as serine proteases of the subtilase family and meiotic prophase aminopeptidase 1 using mass spectrometry-based proteomics. Together, our findings demonstrate that pollen proteases selectively and irreversibly damage integrity and anchorage of columnar respiratory epithelial cells. In turn, alphaherpesviruses benefit from this partial loss-of-barrier function, resulting in increased infection of the respiratory epithelium.

High-throughput sequencing data and the impact of plant gene annotation quality

The use of draft genomes of different species and re-sequencing of accessions and populations are now common tools for plant biology research. The de novo assembled draft genomes make it possible to identify pivotal divergence points in the plant lineage and provide an opportunity to investigate the genomic basis and timing of biological innovations by inferring orthologs between species. Furthermore, re-sequencing facilitates the mapping and subsequent molecular characterization of causative loci for traits, such as those for plant stress tolerance and development. In both cases high-quality gene annotation-the identification of protein-coding regions, gene promoters, and 5'- and 3'-untranslated regions-is critical for investigation of gene function. Annotations are constantly improving but automated gene annotations still require manual curation and experimental validation. This is particularly important for genes with large introns, genes located in regions rich with transposable elements or repeats, large gene families, and segmentally duplicated genes. In this opinion paper, we highlight the impact of annotation quality on evolutionary analyses, genome-wide association studies, and the identification of orthologous genes in plants. Furthermore, we predict that incorporating accurate information from manual curation into databases will dramatically improve the performance of automated gene predictors.

Mechanistic insights into the evolution of DUF26-containing proteins in land plants

Large protein families are a prominent feature of plant genomes and their size variation is a key element for adaptation. However, gene and genome duplications pose difficulties for functional characterization and translational research. Here we infer the evolutionary history of the DOMAIN OF UNKNOWN FUNCTION (DUF) 26-containing proteins. The DUF26 emerged in secreted proteins. Domain duplications and rearrangements led to the appearance of CYSTEINE-RICH RECEPTOR-LIKE PROTEIN KINASES (CRKs) and PLASMODESMATA-LOCALIZED PROTEINS (PDLPs). The DUF26 is land plant-specific but structural analyses of PDLP ectodomains revealed strong similarity to fungal lectins and thus may constitute a group of plant carbohydrate-binding proteins. CRKs expanded through tandem duplications and preferential retention of duplicates following whole genome duplications, whereas PDLPs evolved according to the dosage balance hypothesis. We propose that new gene families mainly expand through small-scale duplications, while fractionation and genetic drift after whole genome multiplications drive families towards dosage balance.

Physiological and Transcriptome Analysis of a Yellow-Green Leaf Mutant in Birch (Betula platyphylla × B. Pendula)

Chlorophyll (Chl)-deficient mutants are ideal materials for the study of Chl biosynthesis, chloroplast development, and photosynthesis. Although the genes encoding key enzymes related to Chl biosynthesis have been well-characterized in herbaceous plants, rice (Oryza sativa L.), Arabidopsis (Arabidopsis thaliana), and maize (Zea mays L.), yellow-green leaf mutants have not yet been fully studied in tree species. In this work, we explored the molecular mechanism of the leaf color formation in a yellow-green leaf mutant (yl). We investigated the differentially expressed genes (DEGs) between yl and control plants (wild type birch (WT) and BpCCR1 overexpression line 11, (C11)) by transcriptome sequencing. Approximately 1163 genes (874 down-regulated and 289 up-regulated) and 930 genes (755 down-regulated and 175 up-regulated) were found to be differentially expressed in yl compared with WT and C11, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis for DEGs revealed that photosynthesis antenna proteins represent the most significant enriched pathway. The expressions of photosynthesis antenna proteins are crucial to the leaf color formation in yl. We also found that Chl accumulate, leaf anatomical structure, photosynthesis, and growth were affected in yl. Taken together, our results not only provide the difference of phenomenal, physiological, and gene expression characteristics in leaves between yl mutant and control plants, but also provide a new insight into the mutation underlying the chlorotic leaf phenotype in birch.

The red bayberry genome and genetic basis of sex determination

Morella rubra, red bayberry, is an economically important fruit tree in south China. Here, we assembled the first high-quality genome for both a female and a male individual of red bayberry. The genome size was 313-Mb, and 90% sequences were assembled into eight pseudo chromosome molecules, with 32 493 predicted genes. By whole-genome comparison between the female and male and association analysis with sequences of bulked and individual DNA samples from female and male, a 59-Kb region determining female was identified and located on distal end of pseudochromosome 8, which contains abundant transposable element and seven putative genes, four of them are related to sex floral development. This 59-Kb female-specific region was likely to be derived from duplication and rearrangement of paralogous genes and retained non-recombinant in the female-specific region. Sex-specific molecular markers developed from candidate genes co-segregated with sex in a genetically diverse female and male germplasm. We propose sex determination follow the ZW model of female heterogamety. The genome sequence of red bayberry provides a valuable resource for plant sex chromosome evolution and also provides important insights for molecular biology, genetics and modern breeding in Myricaceae family.

Branching out: new insights into the genetic regulation of shoot architecture in trees

Directional growth in all plants involves both phototropic and gravitropic responses. Accordingly, mechanisms controlling shoot architecture throughout the plant kingdom are likely similar. However, as forms vary between species due in part to gene copy number and functional divergence, some aspects of how plants predetermine and regulate architecture can differ. This is especially true when comparing annual herbaceous species (e.g. model plants) to woody perennials such as trees. In the past decade, inexpensive genomic sequencing and technological advances enabled gene discovery and functional analyses in trees. This led to the identification of genes associated with tree shoot architecture control. Here, we present recent discoveries on the regulation of shoot architectures for which causative genes have been identified, including dwarf, weeping, columnar, and pillar growth habits. We also discuss potential applications of these findings.

Genomic effects of population collapse in a critically endangered ironwood tree Ostrya rehderiana

Increased human activity and climate change are driving numerous tree species to endangered status, and in the worst cases extinction. Here we examine the genomic signatures of the critically endangered ironwood tree Ostrya rehderiana and its widespread congener O. chinensis. Both species have similar demographic histories prior to the Last Glacial Maximum (LGM); however, the effective population size of O. rehderiana continued to decrease through the last 10,000 years, whereas O. chinensis recovered to Pre-LGM numbers. O. rehderiana accumulated more deleterious mutations, but purged more severely deleterious recessive variations than in O. chinensis. This purging and the gradually reduced inbreeding depression together may have mitigated extinction and contributed to the possible future survival of the outcrossing O. rehderiana. Our findings provide critical insights into the evolutionary history of population collapse and the potential for future recovery of the endangered trees.

Complete chloroplast genome sequence of Betula platyphylla: gene organization, RNA editing, and comparative and phylogenetic analyses

BACKGROUND

Betula platyphylla is a common tree species in northern China that has high economic and medicinal value. Our laboratory has been devoted to genome research on B. platyphylla for approximately 10 years. As primary organelle genomes, the complete genome sequences of chloroplasts are important to study the divergence of species, RNA editing and phylogeny. In this study, we sequenced and analyzed the complete chloroplast (cp) genome sequence of B. platyphylla.

RESULTS

The complete cp genome of B. platyphylla was 160,518 bp in length, which included a pair of inverted repeats (IRs) of 26,056 bp that separated a large single copy (LSC) region of 89,397 bp and a small single copy (SSC) region of 19,009 bp. The annotation contained a total of 129 genes, including 84 protein-coding genes, 37 tRNA genes and 8 rRNA genes. There were 3 genes using alternative initiation codons. Comparative genomics showed that the sequence of the Fagales species cp genome was relatively conserved, but there were still some high variation regions that could be used as molecular markers. The IR expansion event of B. platyphylla resulted in larger cp genomes and rps19 pseudogene formation. The simple sequence repeat (SSR) analysis showed that there were 105 SSRs in the cp genome of B. platyphylla. RNA editing sites recognition indicated that at least 80 RNA editing events occurred in the cp genome. Most of the substitutions were C to U, while a small proportion of them were not. In particular, three editing loci on the rRNA were converted to more than two other bases that had never been reported. For synonymous conversion, most of them increased the relative synonymous codon usage (RSCU) value of the codons. The phylogenetic analysis suggested that B. platyphylla had a closer evolutionary relationship with B. pendula than B. nana.

CONCLUSIONS

In this study, we not only obtained and annotated the complete cp genome sequence of B. platyphylla, but we also identified new RNA editing sites and predicted the phylogenetic relationships among Fagales species. These findings will facilitate genomic, genetic engineering and phylogenetic studies of this important species.

Extreme copy number variation at a tRNA ligase gene affecting phenology and fitness in yellow monkeyflowers

Copy number variation (CNV) is a major part of the genetic diversity segregating within populations, but remains poorly understood relative to single nucleotide variation. Here, we report on a tRNA ligase gene (Migut.N02091; RLG1a) exhibiting unprecedented, and fitness-relevant, CNV within an annual population of the yellow monkeyflower Mimulus guttatus. RLG1a variation was associated with multiple traits in pooled population sequencing (PoolSeq) scans of phenotypic and phenological cohorts. Resequencing of inbred lines revealed intermediate-frequency three-copy variants of RLG1a (trip+; 5/35 = 14%), and trip+ lines exhibited elevated RLG1a expression under multiple conditions. trip+ carriers, in addition to being over-represented in late-flowering and large-flowered PoolSeq populations, flowered later under stressful conditions in a greenhouse experiment (p < 0.05). In wild population samples, we discovered an additional rare RLG1a variant (high+) that carries 250-300 copies of RLG1a totalling ~5.7 Mb (20-40% of a chromosome). In the progeny of a high+ carrier, Mendelian segregation of diagnostic alleles and qPCR-based copy counts indicate that high+ is a single tandem array unlinked to the single-copy RLG1a locus. In the wild, high+ carriers had highest fitness in two particularly dry and/or hot years (2015 and 2017; both p < 0.01), while single-copy individuals were twice as fecund as either CNV type in a lush year (2016: p < 0.005). Our results demonstrate fluctuating selection on CNVs affecting phenological traits in a wild population, suggest that plant tRNA ligases mediate stress-responsive life-history traits, and introduce a novel system for investigating the molecular mechanisms of gene amplification.

Complex Haplotypes of GSTM1 Gene Deletions Harbor Signatures of a Selective Sweep in East Asian Populations

The deletion of the metabolizing Glutathione S-transferase Mu 1 (GSTM1) gene has been associated with multiple cancers, metabolic and autoimmune disorders, as well as drug response. It is unusually common, with allele frequency reaching up to 75% in some human populations. Such high allele frequency of a derived allele with apparent impact on an otherwise conserved gene is a rare phenomenon. To investigate the evolutionary history of this locus, we analyzed 310 genomes using population genetics tools. Our analysis revealed a surprising lack of linkage disequilibrium between the deletion and the flanking single nucleotide variants in this locus. Tests that measure extended homozygosity and rapid change in allele frequency revealed signatures of an incomplete sweep in the locus. Using empirical approaches, we identified the Tanuki haplogroup, which carries the GSTM1 deletion and is found in approximately 70% of East Asian chromosomes. This haplogroup has rapidly increased in frequency in East Asian populations, contributing to a high population differentiation among continental human groups. We showed that extended homozygosity and population differentiation for this haplogroup is incompatible with simulated neutral expectations in East Asian populations. In parallel, we found that the Tanuki haplogroup is significantly associated with the expression levels of other GSTM genes. Collectively, our results suggest that standing variation in this locus has likely undergone an incomplete sweep in East Asia with regulatory impact on multiple GSTM genes. Our study provides the necessary framework for further studies to elucidate the evolutionary reasons that maintain disease-susceptibility variants in the GSTM1 locus.

BARCOSEL: a tool for selecting an optimal barcode set for high-throughput sequencing

BACKGROUND

Current high-throughput sequencing platforms provide capacity to sequence multiple samples in parallel. Different samples are labeled by attaching a short sample specific nucleotide sequence, barcode, to each DNA molecule prior pooling them into a mix containing a number of libraries to be sequenced simultaneously. After sequencing, the samples are binned by identifying the barcode sequence within each sequence read. In order to tolerate sequencing errors, barcodes should be sufficiently apart from each other in sequence space. An additional constraint due to both nucleotide usage and basecalling accuracy is that the proportion of different nucleotides should be in balance in each barcode position. The number of samples to be mixed in each sequencing run may vary and this introduces a problem how to select the best subset of available barcodes at sequencing core facility for each sequencing run. There are plenty of tools available for de novo barcode design, but they are not suitable for subset selection.

RESULTS

We have developed a tool which can be used for three different tasks: 1) selecting an optimal barcode set from a larger set of candidates, 2) checking the compatibility of user-defined set of barcodes, e.g. whether two or more libraries with existing barcodes can be combined in a single sequencing pool, and 3) augmenting an existing set of barcodes. In our approach the selection process is formulated as a minimization problem. We define the cost function and a set of constraints and use integer programming to solve the resulting combinatorial problem. Based on the desired number of barcodes to be selected and the set of candidate sequences given by user, the necessary constraints are automatically generated and the optimal solution can be found. The method is implemented in C programming language and web interface is available at http://ekhidna2.biocenter.helsinki.fi/barcosel .

CONCLUSIONS

Increasing capacity of sequencing platforms raises the challenge of mixing barcodes. Our method allows the user to select a given number of barcodes among the larger existing barcode set so that both sequencing errors are tolerated and the nucleotide balance is optimized. The tool is easy to access via web browser.

Linking genome signatures of selection and adaptation in non-model plants: exploring potential and limitations in the angiosperm Amborella

Selective sweeps may be caused by environmental conditions that select for a gene function or trait at one locus, causing reduced variation at neighboring sites due to linkage, with specific non-selected variants being swept along with the selected variant. For many species, genomic and environmental data are available to test hypotheses that environmental conditions are correlated with selected regions. Most genomic studies relating selection to environment use model organisms or crop species; typically, these studies have genomic data from large numbers of individuals and extensive environmental data. Here, we review studies associating selective sweeps with environment and consider the impediments to successful application of these methods to non-model species. We present an initial investigation into linking genomic regions of selection to environmental conditions in the narrowly distributed, non-model plant Amborella trichopoda (Amborellaceae), the sister species to all other living flowering plants and one of over 2500 plant species endemic to New Caledonia.

Ethylene-Related Gene Expression Networks in Wood Formation

Thickening of tree stems is the result of secondary growth, accomplished by the meristematic activity of the vascular cambium. Secondary growth of the stem entails developmental cascades resulting in the formation of secondary phloem outwards and secondary xylem (i.e., wood) inwards of the stem. Signaling and transcriptional reprogramming by the phytohormone ethylene modifies cambial growth and cell differentiation, but the molecular link between ethylene and secondary growth remains unknown. We addressed this shortcoming by analyzing expression profiles and co-expression networks of ethylene pathway genes using the AspWood transcriptome database which covers all stages of secondary growth in aspen (Populus tremula) stems. ACC synthase expression suggests that the ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) is synthesized during xylem expansion and xylem cell maturation. Ethylene-mediated transcriptional reprogramming occurs during all stages of secondary growth, as deduced from AspWood expression profiles of ethylene-responsive genes. A network centrality analysis of the AspWood dataset identified EIN3D and 11 ERFs as hubs. No overlap was found between the co-expressed genes of the EIN3 and ERF hubs, suggesting target diversification and hence independent roles for these transcription factor families during normal wood formation. The EIN3D hub was part of a large co-expression gene module, which contained 16 transcription factors, among them several new candidates that have not been earlier connected to wood formation and a VND-INTERACTING 2 (VNI2) homolog. We experimentally demonstrated Populus EIN3D function in ethylene signaling in Arabidopsis thaliana. The ERF hubs ERF118 and ERF119 were connected on the basis of their expression pattern and gene co-expression module composition to xylem cell expansion and secondary cell wall formation, respectively. We hereby establish data resources for ethylene-responsive genes and potential targets for EIN3D and ERF transcription factors in Populus stem tissues, which can help to understand the range of ethylene targeted biological processes during secondary growth.

FAR1-RELATED SEQUENCE (FRS) and FRS-RELATED FACTOR (FRF) Family Proteins in Arabidopsis Growth and Development

Transposable elements make important contributions to adaptation and evolution of their host genomes. The well-characterized transposase-derived transcription factor FAR-RED ELONGATED HYPOCOTYLS3 (FHY3) and its homologue FAR-RED IMPAIRED RESPONSE1 (FAR1) have crucial functions in plant growth and development. In addition, FHY3 and FAR1 are the founding members of the FRS (FAR1-RELATED SEQUENCE) and FRF (FRS-RELATED FACTOR) families, which are conserved among land plants. Although the coding sequences of many putative FRS and FRF orthologs have been found in various clades of angiosperms, their physiological functions remain elusive. Here, we summarize recent progress toward characterizing the molecular mechanisms of FHY3 and FAR1, as well as other FRS-FRF family proteins, examining their roles in regulating plant growth and development. This review also suggests future directions for further functional characterization of other FRS-FRF family proteins in plants.

Evolution as an ecosystem process: insights from genomics

Evolution is a fundamental ecosystem process. The study of genomic variation of organisms can not only improve our understanding of evolutionary processes, but also of contemporary and future ecosystem dynamics. We argue that integrative research between the fields of genomics and ecosystem ecology could generate new insights. Specifically, studies of biodiversity and ecosystem functioning, evolutionary rescue, and eco-evolutionary dynamics could all benefit from information about variation in genome structure and the genetic architecture of traits, whereas genomic studies could benefit from information about the ecological context of evolutionary dynamics. We propose new ways to help link research on functional genomic diversity with (reciprocal) interactions between phenotypic evolution and ecosystem change. Despite numerous challenges, we anticipate that the wealth of genomic data being collected on natural populations will improve our understanding of ecosystems.