Highlights d Chromosome-level assembly and methylome of the largest gymnosperm genome so far d Continuous expansion and slow removal of transposons cause conifer huge genome d Large genes with ultra-long introns tend to be expressed at higher levels d Distinctive reproductive evolutionary trajectory compared to angiosperms

As the consequence of complex interactions between different parts of an organ, shape can be used as a predictor of structural-functional relationships implicated in changing environments. Despite such importance, however, it is no surprise that little is known about the genetic detail involved in shape variation, because no approach is currently available for mapping quantitative trait loci (QTLs) that control shape. Here, we address this problem by developing a statistical model that integrates the principle of shape analysis into a mixture-model-based likelihood formulated for QTL mapping. One state-of-theart approach for shape analysis is to identify and analyze the polar coordinates of anatomical landmarks on a shape measured in terms of radii from the centroid to the contour at regular intervals. A procrustes analysis is used to align shapes to filter out position, scale and rotation effects on shape variation. To the end, the accurate and quantitative representation of a shape is produced with aligned radius-centroid-contour (RCC) curves, that is, a function of radial angle at the centroid. The high dimensionality of the RCC data, crucial for a comprehensive description of the geometric feature of a shape, is reduced by principal component (PC) analysis, and the resulting PC axes are treated as phenotypic traits, allowing specific QTLs for global and local shape variability to be mapped, respectively. The usefulness and utilization of the new model for shape mapping in practice are validated by analyzing a mapping data collected from a natural population of poplar, Populus szechuanica var tibetica, and identifying several QTLs for leaf shape in this species. The model provides a powerful tool to compute which genes determine biological shape in plants, animals and humans.

BackgroundPopulus euphratica is a representative model woody plant species for studying resistance to abiotic stresses such as drought and salt. Salt stress is one of the most common environmental factors that affect plant growth and development. MicroRNAs (miRNAs) are small, noncoding RNAs that have important regulatory functions in plant growth, development, and response to abiotic stress.ResultsTo investigate the miRNAs involved in the salt-stress response, we constructed four small cDNA libraries from P. euphratica plantlets treated with or without salt (300 mM NaCl, 3 days) in either the root or leaf. Using high-throughput sequencing to identify miRNAs, we found 164 conserved miRNAs belonging to 44 families. Of these, 136 novel miRNAs were from the leaf, and 128 novel miRNAs were from the root. In response to salt stress, 95 miRNAs belonging to 46 conserved miRNAs families changed significantly, with 56 miRNAs upregulated and 39 miRNAs downregulated in the leaf. A comparison of the leaf and root tissues revealed 155 miRNAs belonging to 63 families with significantly altered expression, including 84 upregulated and 71 downregulated miRNAs. Furthermore, 479 target genes in the root and 541 targets of novel miRNAs in the leaf were predicted, and functional information was annotated using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases.ConclusionsThis study provides a novel visual field for understanding the regulatory roles of miRNAs in response to salt stress in Populus.

BackgroundThe Tibetan poplar (Populus szechuanica var. tibetica Schneid), which is distributed at altitudes of 2,000-4,500 m above sea level, is an ecologically important species of the Qinghai-Tibet Plateau and adjacent areas. However, the genetic adaptations responsible for its ability to cope with the harsh environment remain unknown.ResultsIn this study, a total of 24 expressed sequence tag microsatellite (EST-SSR) markers were used to evaluate the genetic diversity and population structure of Tibetan poplars along an altitude gradient. The 172 individuals were of genotypes from low-, medium- and high-altitude populations, and 126 alleles were identified. The expected heterozygosity (HE) value ranged from 0.475 to 0.488 with the highest value found in low-altitude populations and the lowest in high-altitude populations. Genetic variation was low among populations, indicating a limited influence of altitude on microsatellite variation. Low genetic differentiation and high levels of gene flow were detected both between and within the populations along the altitude gradient. An analysis of molecular variance (AMOVA) showed that 6.38% of the total molecular variance was attributed to diversity between populations, while 93.62% variance was associated with differences within populations. There was no clear correlation between genetic variation and altitude, and a Mantel test between genetic distance and altitude resulted in a coefficient of association of r = 0.001, indicating virtually no correlation.ConclusionMicrosatellite genotyping results showing genetic diversity and low differentiation suggest that extensive gene flow may have counteracted local adaptations imposed by differences in altitude. The genetic analyses carried out in this study provide new insight for conservation and optimization of future arboriculture.

To understand the genetic and molecular mechanisms underlying floral development in Populus tomentosa, we isolated PtLFY, a LEAFY homolog, from a P. tomentosa floral bud cDNA library. DNA gel blot analysis showed that PtLFY is present as a single copy in the genomes of both male and female individuals of P. tomentosa. The genomic copy is composed of three exons and two introns. Relative expression levels of PtLFY in tissues of P. tomentosa were estimated by RT-PCR; our results revealed that PtLFY mRNA is highly abundant in roots and both male and female floral buds. A low level of gene expression was detected in stems and vegetative buds, and no PtLFY-specific transcripts were detected in leaves. PtLFY expression patterns were analyzed during the development of both male and female floral buds in P. tomentosa via real-time quantitative RT-PCR. Continuous, stable and high-level expression of PtLFY-specific mRNA was detected in both male and female floral buds from September 13th to February 25th, but the level of PtLFY transcripts detected in male floral buds was considerably higher than in female floral buds. Our results also showed an inverted repeat PtLFY fragment (PtLFY-IR) effectively blocked flowering of transgenic tobacco plants, and that this effect appeared to be due to post-transcriptional silencing of the endogenous tobacco LFY homologs NFL1 and NFL2.

No abstract

How genes interact with the environment to shape phenotypic variation and evolution is a fundamental question intriguing to biologists from various fields. Existing linear models built on single genes are inadequate to reveal the complexity of genotype-environment (G-E) interactions. Here, we develop a conceptual model for mechanistically dissecting G-E interplay by integrating previously disconnected theories and methods. Under this integration, evolutionary game theory, developmental modularity theory, and a variable selection method allow us to reconstruct environment-induced, maximally informative, sparse, and casual multilayer genetic networks. We design and conduct two mapping experiments by using a desert-adapted tree species to validate the biological application of the model proposed. The model identifies previously uncharacterized molecular mechanisms that mediate trees' response to saline stress. Our model provides a tool to comprehend the genetic architecture of trait variation and evolution and trace the information flow of each gene toward phenotypes within omnigenic networks.

SUMMARYThe coordination of shoots and roots is critical for plants to adapt to changing environments by fine-tuning energy production in leaves and the availability of water and nutrients from roots. To understand the genetic architecture of how these two organs covary during developmental ontogeny, we conducted a mapping experiment using Euphrates poplar (Populus euphratica), a so-called hero tree able to grow in the desert. We geminated intraspecific F 1 seeds of Euphrates Poplar individually in a tube to obtain a total of 370 seedlings, whose shoot and taproot lengths were measured repeatedly during the early stage of growth. By fitting a growth equation, we estimated asymptotic growth, relative growth rate, the timing of inflection point and duration of linear growth for both shoot and taproot growth. Treating these heterochronic parameters as phenotypes, a univariate mapping model detected 19 heterochronic quantitative trait loci (hQTLs), of which 15 mediate the forms of shoot growth and four mediate taproot growth. A bivariate mapping model identified 11 pleiotropic hQTLs that determine the covariation of shoot and taproot growth. Most QTLs detected reside within the region of candidate genes with various functions, thus confirming their roles in the biochemical processes underlying plant growth.