High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers.

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Reciprocal expression of MADS-box genes and DNA methylation reconfiguration initiate bisexual cones in spruce.

The naturally occurring bisexual cone of gymnosperms has long been considered a possible intermediate stage in the origin of flowers, but the mechanisms governing bisexual cone formation remain largely elusive. Here, we employed transcriptomic and DNA methylomic analyses, together with hormone measurement, to investigate the molecular mechanisms underlying bisexual cone development in the conifer Picea crassifolia. Our study reveals a "bisexual" expression profile in bisexual cones, especially in expression patterns of B-class, C-class and LEAFY genes, supporting the out of male model. GGM7 could be essential for initiating bisexual cones. DNA methylation reconfiguration in bisexual cones affects the expression of key genes in cone development, including PcDAL12, PcDAL10, PcNEEDLY, and PcHDG5. Auxin likely plays an important role in the development of female structures of bisexual cones. This study unveils the potential mechanisms responsible for bisexual cone formation in conifers and may shed light on the evolution of bisexuality.

Chromosome-level de novo genome assembly of two conifer-parasitic wasps, Megastigmus duclouxiana and Megastigmus sabinae, reveals genomic imprints of adaptation to hosts.

Conifers make up about one third of global forests but are threatened by seed parasitoid wasp species. Many of these wasps belong to the genus Megastigmus, yet little is known about their genomic background. In this study, we provide chromosome-level genome assemblies for two oligophagous conifer parasitoid species of Megastigmus, which represent the first two chromosome-level genomes of the genus. The assembled genomes of Megastigmus duclouxiana and M. sabinae are 878.48 Mb (scaffold N50 of 215.60 Mb) and 812.98 Mb (scaffold N50 of 139.16 Mb), respectively, which are larger than the genome size of most hymenopterans due to the expansion of transposable elements. Expanded gene families highlight the difference in sensory-related genes between the two species, reflecting the difference in their hosts. We further found that these two species have fewer family members but more single-gene duplications than polyphagous congeners in the gene families of ATP-binding cassette transporter (ABC), cytochrome P450 (P450) and olfactory receptors (OR). These findings shed light on the pattern of adaptation to a narrow spectrum of hosts in oligophagous parasitoids. Our findings suggest potential drivers underlying genome evolution and parasitism adaptation, and provide valuable resources for understanding the ecology, genetics and evolution of Megastigmus, as well as for the research and biological control of global conifer forest pests.

Understanding the Potential Gene Regulatory Network of Starch Biosynthesis in Tartary Buckwheat by RNA-Seq.

Starch is a major component of crop grains, and its content affects food quality and taste. Tartary buckwheat is a traditional pseudo-cereal used in food as well as medicine. Starch content, granule morphology, and physicochemical properties have been extensively studied in Tartary buckwheat. However, the complex regulatory network related to its starch biosynthesis needs to be elucidated. Here, we performed RNA-seq analyses using seven Tartary buckwheat varieties differing in starch content and combined the RNA-seq data with starch content by weighted correlation network analysis (WGCNA). As a result, 10,873 differentially expressed genes (DEGs) were identified and were functionally clustered to six hierarchical clusters. Fifteen starch biosynthesis genes had higher expression level in seeds. Four trait-specific modules and 3131 hub genes were identified by WGCNA, with the lightcyan and brown modules positively correlated with starch-related traits. Furthermore, two potential gene regulatory networks were proposed, including the co-expression of FtNAC70, FtPUL, and FtGBSS1-3 in the lightcyan module and FtbHLH5, C3H, FtBE2, FtISA3, FtSS3-5, and FtSS1 in the brown. All the above genes were preferentially expressed in seeds, further suggesting their role in seed starch biosynthesis. These results provide crucial guidance for further research on starch biosynthesis and its regulatory network in Tartary buckwheat.

Inbreeding in Chinese Fir: Insight into the Rare Self-Fertilizing Event from a Genetic View.

Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.) is a fast-growing conifer with great forestation value and prefers outcrossing with high inbreeding depression effect. Previously, we captured a special Chinese fir parent clone named as 'cx569' that lacks early inbreeding depression. In view of the fact that very little has been published about the rare self-fertilizing event in Chinese fir from a genetic view, herein, we conduct an SSR-based study on the variation of open- and self-pollinated offspring of this parent to gain a view of the rare self-fertilizing event. The results indicated that genetic diversity of self-pollinated offspring was significantly reduced by half (Ho: 0.302, vs. 0.595, p = 0.001; He: 0.274 vs. 0.512, p = 0.002) when compared to an open-pollinated set. Self-pollinated offspring also had significantly positive FIS values (FIS = 0.057, p = 0.034) with a much higher proportion of common allele (20.59% vs. 0), reflecting their heterozygote deficiency. Clustering analysis further indicated a separation of the self- and opened- pollinated groups, implying a natural preference of outcrossing for cx569. However, the cx569 still had 6% acceptance for selfing. When accepted 100% for its own pollen, the cx569 led to a genetically unique selfing group. Additionally, this selfing group seemed to be consistently homozygous at seven particular loci. These findings gave us more genetic clues to gain insight into the rare self-fertilizing event in conifer (Chinese fir).

Stress-responsive retrotransposable elements in conifers.

Conifers are important in many forest ecosystems. They have a long generation time and are immobile; therefore, they require considerable plasticity to adapt to environmental stresses. Moreover, conifers have a large genome, a high proportion of which is occupied by repetitive elements. Retrotransposons are the most highly represented repetitive elements in conifers whose whole-genome sequences have been examined. These retrotransposons are usually silenced, to maintain genome integrity; however, some are activated by environmental stress. The insertion of retrotransposons into genic regions is associated with phenotypic and genetic diversity. The large number and high diversity of retrotransposons in conifer genomes suggest that they play a role in adaptation to the environment. In this review, progress in research on the roles of retrotransposons in the stress responses of conifers is reviewed, and potential future work is discussed.

The western redcedar genome reveals low genetic diversity in a self-compatible conifer.

We assembled the 9.8-Gbp genome of western redcedar (WRC; Thuja plicata), an ecologically and economically important conifer species of the Cupressaceae. The genome assembly, derived from a uniquely inbred tree produced through five generations of self-fertilization (selfing), was determined to be 86% complete by BUSCO analysis, one of the most complete genome assemblies for a conifer. Population genomic analysis revealed WRC to be one of the most genetically depauperate wild plant species, with an effective population size of approximately 300 and no significant genetic differentiation across its geographic range. Nucleotide diversity, π, is low for a continuous tree species, with many loci showing zero diversity, and the ratio of π at zero- to fourfold degenerate sites is relatively high (approximately 0.33), suggestive of weak purifying selection. Using an array of genetic lines derived from up to five generations of selfing, we explored the relationship between genetic diversity and mating system. Although overall heterozygosity was found to decline faster than expected during selfing, heterozygosity persisted at many loci, and nearly 100 loci were found to deviate from expectations of genetic drift, suggestive of associative overdominance. Nonreference alleles at such loci often harbor deleterious mutations and are rare in natural populations, implying that balanced polymorphisms are maintained by linkage to dominant beneficial alleles. This may account for how WRC remains responsive to natural and artificial selection, despite low genetic diversity.

Comparative transcriptomics of Pinus massoniana organs provides insights on terpene biosynthesis regulation.

Terpenoids are the most important natural products collected from conifer species. However, the molecular mechanisms and core factors underlying terpenoid biosynthesis in Pinus massoniana remain unclear. To clarify these mechanisms, this study aimed to identify potential genes that might participate in the terpenoid biosynthesis of P. massoniana. In this study, single molecule real-time (SMRT) sequencing and expression analysis were used to confirm the expression patterns of genes involved in the cones, immature needles, mature needles, immature branches, and mature branches of P. massoniana. A total of 31,331 lncRNAs and 71,240 mRNAs were identified from these organs, and the greatest number of differentially expressed genes (DEGs) was discovered between needles and branches. Weighted gene coexpression network analysis (WGCNA) classified all expressed genes into nine typical modules with 11 kinds of transcription factors (TFs), namely, AP2-ERF, ARF, AUX-IAA, C2H2, Dof, F-box, SBP, WRKY, bHLH, bZIP, and GRAS, and seven kinds of functional genes, namely, ABC transporter, cellulose synthase (CesA), leucine-rich repeats (LRR), cytochrome P450 (CYT P450), pathogenesis-related protein (PR), terpene synthase (TPS), and chlorophyllase enzyme. A molecular network was constructed for hub genes, TFs, and functional genes in three modules. The potential function of eight candidate genes, including PmbHLH2, PmERF1, PmRGA, PmGAI, PmbZIP1, PmLOB1, PmMADS1, and PmMYB1, was validated through correlation analysis between terpenoid contents and expression levels, subcellular localization, and transcriptional activation activity, which provides us with probable regulators of terpenoid biosynthesis in conifers.

Demographic modelling helps track the rapid and recent divergence of a conifer species pair from Central Mexico.

Secondary contact of recently diverged species may have several outcomes, ranging from rampant hybridization to reinforced reproductive isolation. In plants, selfing tolerance and disjunct reproductive phenology may lead to reproductive isolation at contact zones. However, they may also evolve under both allopatric or parapatric frameworks and originate from adaptive and/or neutral forces. Inferring the historical demography of diverging taxa is thus a crucial step to identify factors that may have led to putative reproductive isolation. We explored various competing demographypotheses to account for the rapid divergence of a fir species complex (Abies flinckii-A. religiosa) distributed in "sky-islands" across central Mexico (i.e., along the Trans-Mexican Volcanic Belt; TMVB). Despite co-occurring in two independent sympatric regions (west and centre), these taxa rarely interbreed because of disjunct reproductive phenologies. We genotyped 1147 single nucleotide polymorphisms, generated by GBS (genotyping by sequencing), across 23 populations, and compared multiple scenarios based on the geological history of the TMVB. The best-fitting model revealed one of the most rapid and complete speciation cases for a conifer species-pair, dating back to ~1.2 million years ago. Coupled with the lack of support for stepwise colonization, our coalescent inferences point to an early cessation of interspecific gene flow under parapatric speciation; ancestral gene flow during divergence was asymmetrical (mostly from western firs into A. religiosa) and exclusive to the most ancient (i.e., central) contact zone. Factors promoting rapid reproductive isolation should be explored in other slowly evolving species complexes as they may account for the large tropical and subtropical diversity.

Intercontinental Diversity of Caballeronia Gut Symbionts in the Conifer Pest Bug Leptoglossus occidentalis.

Many stinkbugs in the superfamily Coreoidea (Hemiptera: Heteroptera) develop crypts in the posterior midgut, harboring Caballeronia (Burkholderia) symbionts. These symbionts form a monophyletic group in Burkholderia sensu lato, called the "stinkbug-associated beneficial and environmental (SBE)" group, recently reclassified as the new genus Caballeronia. SBE symbionts are separated into the subclades SBE-α and SBE-ß. Previous studies suggested a regional effect on the symbiont infection pattern; Japanese and American bug species are more likely to be associated with SBE-α, while European bug species are almost exclusively associated with SBE-ß. However, since only a few insect species have been investigated, it remains unclear whether region-specific infection is general. We herein investigated Caballeronia gut symbionts in diverse Japanese, European, and North American populations of a cosmopolitan species, the Western conifer seed bug Leptoglossus occidentalis (Coreoidea: Coreidae). A mole-cular phylogenetic ana-lysis of the 16S rRNA gene demonstrated that SBE-ß was the most dominant in all populations. Notably, SBE-α was rarely detected in any region, while a third clade, the "Coreoidea clade" occupied one fourth of the tested populations. Although aposymbiotic bugs showed high mortality, SBE-α- and SBE-ß-inoculated insects both showed high survival rates; however, a competition assay demonstrated that SBE-ß outcompeted SBE-α in the midgut crypts of L. occidentalis. These results strongly suggest that symbiont specificity in the Leptoglossus-Caballeronia symbiotic association is influenced by the host rather than geography, while the geographic distribution of symbionts may be more important in other bugs.

Plastid Phylogenomics and Plastomic Diversity of the Extant Lycophytes.

Although extant lycophytes represent the most ancient surviving lineage of early vascular plants, their plastomic diversity has long been neglected. The ancient evolutionary history and distinct genetic diversity patterns of the three lycophyte families, each with its own characteristics, provide an ideal opportunity to investigate the interfamilial relationships of lycophytes and their associated patterns of evolution. To compensate for the lack of data on Lycopodiaceae, we sequenced and assembled 14 new plastid genomes (plastomes). Combined with other lycophyte plastomes available online, we reconstructed the phylogenetic relationships of the extant lycophytes based on 93 plastomes. We analyzed, traced, and compared the plastomic diversity and divergence of the three lycophyte families (Isoëtaceae, Lycopodiaceae, and Selaginellaceae) in terms of plastomic diversity by comparing their plastome sizes, GC contents, substitution rates, structural rearrangements, divergence times, ancestral states, RNA editings, and gene losses. Comparative analysis of plastid phylogenomics and plastomic diversity of three lycophyte families will set a foundation for further studies in biology and evolution in lycophytes and therefore in vascular plants.

Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes.

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

Club-mosses (Diphasiastrum, Lycopodiaceae) from the Far East - Introgression and possible cryptic speciation.

Hybridization occurs often in the genus Diphasiastrum (Lycopodiaceae), which corroborates reports for the two other recognized lycophyte families, Isoëtaceae and Selaginellaceae. Here we investigate the case of D. alpinum and D. sitchense from the Russian Far East (Kamchatka). Their hybrid, D. × takedae, was morphologically recognizable in 16 out of 22 accessions showing molecular signatures of hybridization; the remaining accessions displayed the morphology of either D. alpinum (3) or D. sitchense (3). We sequenced markers for chloroplast microsatellites (cp, 175 accessions from Kamchatka) and for the two nuclear markers RPB and LFY (175 and 152 accessions). A selection of 42 accessions, including all hybrid accessions, was analysed via genotyping by sequencing (GBS). We found multiple, but apparently uniparental hybridization, clearly characterized by a deviating group of haplotypes for D. sitchense and all hybrids. All accessions showing molecular signatures of hybridization in nuclear markers revealed the parental haplotype of D. sitchense, however only the LFY marker differentiated between the parent species. GBS, including 69,819 quality-filtered single nucleotid polymorphisms, unambiguously identified the hybrids and revealed introgression to occur. Most of the hybrids were F1, but three turned out to be backcrosses with D. alpinum (one) and with D. sitchense (two). These observations are in contrast to prior findings on three European species and their intermediates where all three hybrids turned out to be independent F1 crosses without evidence of recent backcrossing. In this study, backcrossing was detected, which indicates a limited fertility of the hybrid taxon D. × takedae. A comparison of accessions of Kamchatkian D. alpinum with plants from Europe indicated possible cryptic speciation. Accessions from the Far East had (i) a lower DNA content (7.0 vs. 7.5 pg/2C), (ii) different prevailing cp haplotypes, and (iii) RPB genotypes, and (iv) a clearly different SNP pattern in GBS. Diphasiastrum sitchense and the similar D. nikoënse, for the latter additional accessions from Japan were investigated, appeared as forms of one diverse species, sharing genotypes in both nuclear markers, although chloroplast haplotypes and DNA content show slight variations.

Evolution of mitochondrial RNA editing in extant gymnosperms.

To unveil the evolution of mitochondrial RNA editing in gymnosperms, we characterized mitochondrial genomes (mitogenomes), plastid genomes, RNA editing sites, and pentatricopeptide repeat (PPR) proteins from 10 key taxa representing four of the five extant gymnosperm clades. The assembled mitogenomes vary in gene content due to massive gene losses in Gnetum and Conifer II clades. Mitochondrial gene expression levels also vary according to protein function, with the most highly expressed genes involved in the respiratory complex. We identified 9132 mitochondrial C-to-U editing sites, as well as 2846 P-class and 8530 PLS-class PPR proteins. Regains of editing sites were demonstrated in Conifer II rps3 transcripts whose corresponding mitogenomic sequences lack introns due to retroprocessing. Our analyses reveal that non-synonymous editing is efficient and results in more codons encoding hydrophobic amino acids. In contrast, synonymous editing, although performed with variable efficiency, can increase the number of U-ending codons that are preferentially utilized in gymnosperm mitochondria. The inferred loss-to-gain ratio of mitochondrial editing sites in gymnosperms is 2.1:1, of which losses of non-synonymous editing are mainly due to genomic C-to-T substitutions. However, such substitutions only explain a small fraction of synonymous editing site losses, indicating distinct evolutionary mechanisms. We show that gymnosperms have experienced multiple lineage-specific duplications in PLS-class PPR proteins. These duplications likely contribute to accumulated RNA editing sites, as a mechanistic correlation between RNA editing and PLS-class PPR proteins is statistically supported.

Transitions of foliar mycobiota community and transcriptome in response to pathogenic conifer needle interactions.

Profiling the host-mycobiota interactions in healthy vs. diseased forest ecosystems helps understand the dynamics of understudied yet increasingly important threats to forest health that are emerging due to climate change. We analyzed the structural and functional changes of the mycobiota and the responses of Pinus contorta in the Lophodermella needle cast pathosystem through metabarcoding and metatranscriptomics. When needles transitioned from asymptomatic to symptomatic, dysbiosis of the mycobiota occurred, but with an enrichment of Lophodermella pathogens. Many pathogenicity-related genes were highly expressed by the mycobiota at the necrotrophic phase, showing an active pathogen response that are absent in asymptomatic needles. This study also revealed that Lophodermella spp. are members of a healthy needle mycobiota that have latent lifestyles suggesting that other pine needle pathogens may have similar biology. Interestingly, Pinus contorta upregulated defense genes in healthy needles, indicating response to fungal recognition, while a variety of biotic and abiotic stresses genes were activated in diseased needles. Further investigation to elucidate the possible antagonistic interplay of other biotic members leading to disease progression and/or suppression is warranted. This study provides insights into microbial interactions in non-model pathosystems and contributes to the development of new forest management strategies against emerging latent pathogens.

Phylogeny, character evolution, and systematics of the fern family Ophioglossaceae based on Sanger sequence data, plastomes, and morphology.

Adder's tongue ferns or Ophioglossaceae are best known among evolutionary biologists and botanists for their highest chromosome count of any known organisms, the presence of sporophores, and simple morphology. Previous studies recovered and strongly supported the monophyly of the family and the two multi-generic subfamilies, Botrychioideae and Ophioglossoideae, but the relationships among these and two other subfamilies (Helminthostachyoideae and Mankyuoideae) are not well resolved preventing us from understanding the character evolution. The monophyly of and the relationships in the species-rich genus, Ophioglossum, have not well been understood. In this study, new phylogenetic trees are reconstructed based on four datasets: Sanger sequences of eight plastid markers of 184 accessions, 22 plastomes (12 are new), 29 morphological characters, and combined Sanger and morphological data. Our major results include: (1) the relationships among the four subfamilies are well resolved and strongly supported in Bayesian and parsimony analyses based on plastomes: Mankyua is sister to the rest, followed by Ophioglossoideae which are sister to Helminthostachys + Botrychioideae; (2) Sanger data, plastomes, and combined Sanger and morphological data recovered and strongly supported the monophyly of Ophioglossum in its current circumscription (sensu lato; s.l.) in Bayesian and/or parsimony analyses; (3) within Ophioglossum s.l., four deeply diverged clades are identified and the relationships among the four clades are well resolved; (4) evolution of 34 morphological characters is analyzed in the context of the new phylogeny, among which shape of rhizomes, germination time of spores, shape of early gametophytes, and a number of other characters are found to contain interesting phylogenetic signal; and (5) based on the new phylogeny and character evolution, we propose a new classification of Ophioglossaceae in which the currently circumscribed Ophioglossum is divided into four genera including three new ones: Goswamia, Haukia, and Whittieria considering their molecular, morphological, ecological, and biogeographical distinctiveness.

Incomplete lineage sorting and local extinction shaped the complex evolutionary history of the Paleogene relict conifer genus, Chamaecyparis (Cupressaceae).

Inferring accurate biogeographic history of plant taxa with an East Asia (EA)-North America (NA) is usually hindered by conflicting phylogenies and a poor fossil record. The current distribution of Chamaecyparis (false cypress; Cupressaceae) with four species in EA, and one each in western and eastern NA, and its relatively rich fossil record, make it an excellent model for studying the EA-NA disjunction. Here we reconstruct phylogenomic relationships within Chamaecyparis using > 1400 homologous nuclear and 61 plastid genes. Our phylogenomic analyses using concatenated and coalescent approaches revealed strong cytonuclear discordance and conflicting topologies between nuclear gene trees. Incomplete lineage sorting (ILS) and hybridization are possible explanations of conflict; however, our coalescent analyses and simulations suggest that ILS is the major contributor to the observed phylogenetic discrepancies. Based on a well-resolved species tree and four fossil calibrations, the crown lineage of Chamaecyparis is estimated to have originated in the upper Cretaceous, followed by diversification events in the early and middle Paleogene. Ancestral area reconstructions suggest that Chamaecyparis had an ancestral range spanning both EA and NA. Fossil records further indicate that this genus is a relict of the "boreotropical" flora, and that local extinctions of European species were caused by global cooling. Overall, our results unravel a complex evolutionary history of a Paleogene relict conifer genus, which may have involved ILS, hybridization and the extinction of local species.

A novel cupulate seed plant, Xadzigacalix quatsinoensis gen. et sp. nov., provides new insight into the Mesozoic radiation of gymnosperms.

PREMISE: Anatomically preserved evidence for a novel clade of gymnosperms emphasizes diversity of seed plants immediately prior to the appearance of angiosperm fossils in the paleontological record. METHODS: Cupulate seeds from the Early Cretaceous Apple Bay locality (Vancouver Island) are described from serial cellulose acetate peels and three-dimensional reconstruction. Phylogenetic context is assessed through the comparative analysis of gymnosperm seed producing fructifications and maximum parsimony analysis of a revised morphological data set for seed plant phylogeny. RESULTS: Xadzigacalix quatsinoensis gen. et sp. nov. is characterized by an orthotropous ovule with an elongated micropyle and complex integument, enclosed within a radial cupule. The micropylar canal is elongated; and the nucellus extends into the micropyle to seal the post pollination ovule. Except at the apex of the micropyle, the seed is completely enclosed by a parenchymatous cupule with ca. 20 axially elongated secretory ducts. The cupulate seed is produced upon a triangular woody stele, consisting of a parenchymatous pith surrounded by radially aligned tracheids. The stele produces three short terete traces that terminate within the base of the cupule as transfusion tissue at the seed chalaza. CONCLUSIONS: Organography, vascularization, nature of the integument and nucellus, and configuration of the micropylar canal distinguish Xadzigacalix quatsinoensis from all other gymnosperm clades. Cladistic analyses suggest the new plant may have affinities with gnetophytes or angiosperms. These results are complemented with a critical re-evaluation of ovulate structures for Mesozoic gymnosperms, providing new insight into plant diversity immediately antecedent to the explosive diversification of flowering plants.

Genome-Wide Prediction of Transcription Start Sites in Conifers.

The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.