Population genetic structure,genetic diversity,and natural history of the South American species of Nothofagus subgenus Lophozonia (Nothofagaceae) inferred from nuclear microsatellite data

The effect of glaciation on the levels and patterns of genetic variation has been well studied in the Northern Hemisphere. However, although glaciation has undoubtedly shaped the genetic structure of plants in the Southern Hemisphere, fewer studies have characterized the effect, and almost none of them using microsatellites. Particularly, complex patterns of genetic structure might be expected in areas such as the Andes, where both latitudinal and altitudinal glacial advance and retreat have molded modern plant communities. We therefore studied the population genetics of three closely related, hybridizing species of Nothofagus (N. obliqua, N. alpina, and N. glauca, all of subgenus Lophozonia; Nothofagaceae) from Chile. To estimate population genetic parameters and infer the influence of the last ice age on the spatial and genetic distribution of these species, we examined and analyzed genetic variability at seven polymorphic microsatellite DNA loci in 640 individuals from 40 populations covering most of the ranges of these species in Chile. Populations showed no significant inbreeding and exhibited relatively high levels of genetic diversity (H_E = 0.502–0.662) and slight, but significant, genetic structure (R_ST = 8.7–16.0%). However, in N. obliqua, the small amount of genetic structure was spatially organized into three well‐defined latitudinal groups. Our data may also suggest some introgression of N. alpina genes into N. obliqua in the northern populations. These results allowed us to reconstruct the influence of the last ice age on the genetic structure of these species, suggesting several centers of genetic diversity for N. obliqua and N. alpina, in agreement with the multiple refugia hypothesis.     

The Ficus erecta genome aids Ceratocystis canker resistance breeding in common fig (F. carica)

Ficus erecta, a wild relative of the common fig (F. carica), is a donor of Ceratocystis canker resistance in fig breeding programmes. Interspecific hybridization followed by recurrent backcrossing is an effective method to transfer the resistance trait from wild to cultivated fig. However, this process is time consuming and labour intensive for trees, especially for gynodioecious plants such as fig. In this study, genome resources were developed for F. erecta to facilitate fig breeding programmes. The genome sequence of F. erecta was determined using single‐molecule real‐time sequencing technology. The resultant assembly spanned 331.6 Mb with 538 contigs and an N50 length of 1.9 Mb, from which 51 806 high‐confidence genes were predicted. Pseudomolecule sequences corresponding to the chromosomes of F. erecta were established with a genetic map based on single nucleotide polymorphisms from double‐digest restriction‐site‐associated DNA sequencing. Subsequent linkage analysis and whole‐genome resequencing identified a candidate gene for the Ceratocystis canker resistance trait. Genome‐wide genotyping analysis enabled the selection of female lines that possessed resistance and effective elimination of the donor genome from the progeny. The genome resources provided in this study will accelerate and enhance disease‐resistance breeding programmes in fig.     

Identification of chloroplast genome loci suitable for high‐resolution phylogeographic studies of Colocasia esculenta (L.) Schott (Araceae) and closely related taxa

Recently, we reported the chloroplast genome‐wide association of oligonucleotide repeats, indels and nucleotide substitutions in aroid chloroplast genomes. We hypothesized that the distribution of oligonucleotide repeat sequences in a single representative genome can be used to identify mutational hotspots and loci suitable for population genetic, phylogenetic and phylogeographic studies. Using information on the location of oligonucleotide repeats in the chloroplast genome of taro (Colocasia esculenta), we designed 30 primer pairs to amplify and sequence polymorphic loci. The primers have been tested in a range of intra‐specific to intergeneric comparisons, including ten taro samples (Colocasia esculenta) from diverse geographical locations, four other Colocasia species (C. affinis, C. fallax, C. formosana, C. gigantea) and three other aroid genera (represented by Remusatia vivipara, Alocasia brisbanensis and Amorphophallus konjac). Multiple sequence alignments for the intra‐specific comparison revealed nucleotide substitutions (point mutations) at all 30 loci and microsatellite polymorphisms at 14 loci. The primer pairs reported here reveal levels of genetic variation suitable for high‐resolution phylogeographic and evolutionary studies of taro and other closely related aroids. Our results confirm that information on repeat distribution can be used to identify loci suitable for such studies, and we expect that this approach can be used in other plant groups.     

Chromosome‐level genome assembly of an important pine defoliator,Dendrolimus punctatus (Lepidoptera; Lasiocampidae)

Dendrolimus spp. are important destructive pests of conifer forests, and Dendrolimus punctatus Walker (Lepidoptera; Lasiocampidae) is the most widely distributed Dendrolimus species. During periodic outbreaks, this species is said to make “fire without smoke” because large areas of pine forest can be quickly and heavily damaged. Yet, little is known about the molecular mechanisms that underlie the unique ecological characteristics of this forest insect. Here, we combined Pacific Biosciences (PacBio) RSII single‐molecule long reads and high‐throughput chromosome conformation capture (Hi‐C) genomics‐linked reads to produce a high‐quality, chromosome‐level reference genome for D. punctatus. The final assembly was 614 Mb with contig and scaffold N50 values of 1.39 and 22.15 Mb, respectively, and 96.96% of the contigs anchored onto 30 chromosomes. Based on the prediction, this genome contained 17,593 protein‐coding genes and 56.16% repetitive sequences. Phylogenetic analyses indicated that D. punctatus diverged from the common ancestor of Hyphantria cunea, Spodoptera litura and Thaumetopoea pityocampa ~ 108.91 million years ago. Many gene families that were expanded in the D. punctatus genome were significantly enriched for the xenobiotic biodegradation system, especially the cytochrome P450 gene family. This high‐quality, chromosome‐level reference genome will be a valuable resource for understanding mechanisms of D. punctatus outbreak and host resistance adaption. Because this is the first Lasiocampidae insect genome to be sequenced, it also will serve as a reference for further comparative genomics.     

A unigene set for European beech (Fagus sylvatica L.) and its use to decipher the molecular mechanisms involved in dormancy regulation

Systematic sequencing is the method of choice for generating genomic resources for molecular marker development and candidate gene identification in nonmodel species. We generated 47 357 Sanger ESTs and 2.2M Roche‐454 reads from five cDNA libraries for European beech (Fagus sylvatica L.). This tree species of high ecological and economic value in Europe is among the most representative trees of deciduous broadleaf forests. The sequences generated were assembled into 21 057 contigs with MIRA software. Functional annotations were obtained for 85% of these contigs, from the proteomes of four plant species, Swissprot accessions and the Gene Ontology database. We were able to identify 28 079 in silico SNPs for future marker development. Moreover, RNAseq and qPCR approaches identified genes and gene networks regulated differentially between two critical phenological stages preceding vegetative bud burst (the quiescent and swelling buds stages). According to climatic model‐based projection, some European beech populations may be endangered, particularly at the southern and eastern edges of the European distribution range, which are strongly affected by current climate change. This first genomic resource for the genus Fagus should facilitate the identification of key genes for beech adaptation and management strategies for preserving beech adaptability.     

Development of chloroplast genomic resources for Cynara

In this study, new chloroplast (cp) resources were developed for the genus Cynara, using whole cp genomes from 20 genotypes, by means of high‐throughput sequencing technologies. Our target species included seven globe artichokes, two cultivated cardoons, eight wild artichokes, and three other wild Cynara species (C. baetica, C. cornigera and C. syriaca). One complete cp genome was isolated using short reads from a whole‐genome sequencing project, while the others were obtained by means of long‐range PCR, for which primer pairs are provided here. A de novo assembly strategy combined with a reference‐based assembly allowed us to reconstruct each cp genome. Comparative analyses among the newly sequenced genotypes and two additional Cynara cp genomes (‘Brindisino’ artichoke and C. humilis) retrieved from public databases revealed 126 parsimony informative characters and 258 singletons in Cynara, for a total of 384 variable characters. Thirty‐nine SSR loci and 34 other INDEL events were detected. After data analysis, 37 primer pairs for SSR amplification were designed, and these molecular markers were subsequently validated in our Cynara genotypes. Phylogenetic analysis based on all cp variable characters provided the best resolution when compared to what was observed using only parsimony informative characters, or only short ‘variable’ cp regions. The evaluation of the molecular resources obtained from this study led us to support the ‘super‐barcode’ theory and consider the total cp sequence of Cynara as a reliable and valuable molecular marker for exploring species diversity and examining variation below the species level.     

Phylogenetic position of the coral symbiont Ostreobium (Ulvophyceae) inferred from chloroplast genome data

The green algal genus Ostreobium is an important symbiont of corals, playing roles in reef decalcification and providing photosynthates to the coral during bleaching events. A chloroplast genome of a cultured strain of Ostreobium was available, but low taxon sampling and Ostreobium's early‐branching nature left doubt about its phylogenetic position. Here, we generate and describe chloroplast genomes from four Ostreobium strains as well as Avrainvillea mazei and Neomeris sp., strategically sampled early‐branching lineages in the Bryopsidales and Dasycladales respectively. At 80,584 bp, the chloroplast genome of Ostreobium sp. HV05042 is the most compact yet found in the Ulvophyceae. The Avrainvillea chloroplast genome is ~94 kbp and contains introns in infA and cysT that have nearly complete sequence identity except for an open reading frame (ORF) in infA that is not present in cysT. In line with other bryopsidalean species, it also contains regions with possibly bacteria‐derived ORFs. The Neomeris data did not assemble into a canonical circular chloroplast genome but a large number of contigs containing fragments of chloroplast genes and showing evidence of long introns and intergenic regions, and the Neomeris chloroplast genome size was estimated to exceed 1.87 Mb. Chloroplast phylogenomics and 18S nrDNA data showed strong support for the Ostreobium lineage being sister to the remaining Bryopsidales. There were differences in branch support when outgroups were varied, but the overall support for the placement of Ostreobium was strong. These results permitted us to validate two suborders and introduce a third, the Ostreobineae.     

Decoding the oak genome: public release of sequence data,assembly, annotation and publication strategies

The 1.5 Gbp/2C genome of pedunculate oak (Quercus robur) has been sequenced. A strategy was established for dealing with the challenges imposed by the sequencing of such a large, complex and highly heterozygous genome by a whole‐genome shotgun (WGS) approach, without the use of costly and time‐consuming methods, such as fosmid or BAC clone‐based hierarchical sequencing methods. The sequencing strategy combined short and long reads. Over 49 million reads provided by Roche 454 GS‐FLX technology were assembled into contigs and combined with shorter Illumina sequence reads from paired‐end and mate‐pair libraries of different insert sizes, to build scaffolds. Errors were corrected and gaps filled with Illumina paired‐end reads and contaminants detected, resulting in a total of 17 910 scaffolds (>2 kb) corresponding to 1.34 Gb. Fifty per cent of the assembly was accounted for by 1468 scaffolds (N50 of 260 kb). Initial comparison with the phylogenetically related Prunus persica gene model indicated that genes for 84.6% of the proteins present in peach (mean protein coverage of 90.5%) were present in our assembly. The second and third steps in this project are genome annotation and the assignment of scaffolds to the oak genetic linkage map. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement, the oak genome data have been released into public sequence repositories in advance of publication. In this presubmission paper, the oak genome consortium describes its principal lines of work and future directions for analyses of the nature, function and evolution of the oak genome.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

The complete chloroplast genome of Gracilariopsis lemaneiformis (Rhodophyta) gives new insight into the evolution of family Gracilariaceae

The complete chloroplast genome of Gracilariopsis lemaneiformis was recovered from a Next Generation Sequencing data set. Without quadripartite structure, this chloroplast genome (183,013 bp, 27.40% GC content) contains 202 protein‐coding genes, 34 tRNA genes, 3 rRNA genes, and 1 tmRNA gene. Synteny analysis showed plasmid incorporation regions in chloroplast genomes of three species of family Gracilariaceae and in Grateloupia taiwanensis of family Halymeniaceae. Combined with reported red algal plasmid sequences in nuclear and mitochondrial genomes, we postulated that red algal plasmids may have played an important role in ancient horizontal gene transfer among nuclear, chloroplast, and mitochondrial genomes. Substitution rate analysis showed that purifying selective forces maintaining stability of protein‐coding genes of nine red algal chloroplast genomes over long periods must be strong and that the forces acting on gene groups and single genes of nine red algal chloroplast genomes were similar and consistent. The divergence of Gp. lemaneiformis occurred ~447.98 million years ago (Mya), close to the divergence time of genus Pyropia and Porphyra (443.62 Mya).     

Resistance gene enrichment sequencing (RenSeq) enables reannotation of the NB‐LRR gene family from sequenced plant genomes and rapid mapping of resistance loci in segregating populations

RenSeq is a NB‐LRR (nucleotide binding‐site leucine‐rich repeat) gene‐targeted, Resistance gene enrichment and sequencing method that enables discovery and annotation of pathogen resistance gene family members in plant genome sequences. We successfully applied RenSeq to the sequenced potato Solanum tuberosum clone DM, and increased the number of identified NB‐LRRs from 438 to 755. The majority of these identified R gene loci reside in poorly or previously unannotated regions of the genome. Sequence and positional details on the 12 chromosomes have been established for 704 NB‐LRRs and can be accessed through a genome browser that we provide. We compared these NB‐LRR genes and the corresponding oligonucleotide baits with the highest sequence similarity and demonstrated that ~80% sequence identity is sufficient for enrichment. Analysis of the sequenced tomato S. lycopersicum ‘Heinz 1706’ extended the NB‐LRR complement to 394 loci. We further describe a methodology that applies RenSeq to rapidly identify molecular markers that co‐segregate with a pathogen resistance trait of interest. In two independent segregating populations involving the wild Solanum species S. berthaultii (Rpi‐ber2) and S. ruiz‐ceballosii (Rpi‐rzc1), we were able to apply RenSeq successfully to identify markers that co‐segregate with resistance towards the late blight pathogen Phytophthora infestans. These SNP identification workflows were designed as easy‐to‐adapt Galaxy pipelines.     

Chromosome‐level genome assembly of Triplophysa tibetana,a fish adapted to the harsh high‐altitude environment of the Tibetan Plateau

Triplophysa is an endemic fish genus of the Tibetan Plateau in China. Triplophysa tibetana, which lives at a recorded altitude of ~4,000 m and plays an important role in the highland aquatic ecosystem, serves as an excellent model for investigating high‐altitude environmental adaptation. However, evolutionary and conservation studies of T. tibetana have been limited by scarce genomic resources for the genus Triplophysa. In the present study, we applied PacBio sequencing and the Hi‐C technique to assemble the T. tibetana genome. A 652‐Mb genome with 1,325 contigs with an N50 length of 3.1 Mb was obtained. The 1,137 contigs were further assembled into 25 chromosomes, representing 98.7% and 80.47% of all contigs at the base and sequence number level, respectively. Approximately 260 Mb of sequence, accounting for ~39.8% of the genome, was identified as repetitive elements. DNA transposons (16.3%), long interspersed nuclear elements (12.4%) and long terminal repeats (11.0%) were the most repetitive types. In total, 24,372 protein‐coding genes were predicted in the genome, and ~95% of the genes were functionally annotated via a search in public databases. Using whole genome sequence information, we found that T. tibetana diverged from its common ancestor with Danio rerio ~121.4 million years ago. The high‐quality genome assembled in this work not only provides a valuable genomic resource for future population and conservation studies of T. tibetana, but it also lays a solid foundation for further investigation into the mechanisms of environmental adaptation of endemic fishes in the Tibetan Plateau.     

Flip‐flop organization in the chloroplast genome of Capsosiphon fulvescens (Ulvophyceae,Chlorophyta)

To better understand organelle genome evolution of the ulvophycean green alga Capsosiphon fulvescens, we sequenced and characterized its complete chloroplast genome. The circular chloroplast genome was 111,561 bp in length with 31.3% GC content that contained 108 genes including 77 protein‐coding genes, two copies of rRNA operons, and 27 tRNAs. In this analysis, we found the two types of isoform, called heteroplasmy, were likely caused by a flip‐flop organization. The flip‐flop mechanism may have caused structural variation and gene conversion in the chloroplast genome of C. fulvescens. In a phylogenetic analysis based on all available ulvophycean chloroplast genome data, including a new C. fulvescens genome, we found three major conflicting signals for C. fulvescens and its sister taxon Pseudoneochloris marina within 70 individual genes: (i) monophyly with Ulotrichales, (ii) monophyly with Ulvales, and (iii) monophyly with the clade of Ulotrichales and Ulvales. Although the 70‐gene concatenated phylogeny supported monophyly with Ulvales for both species, these complex phylogenetic signals of individual genes need further investigations using a data‐rich approach (i.e., organelle genome data) from broader taxon sampling.     

Construction and comparison of three reference‐quality genome assemblies for soybean

We report reference‐quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single‐nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan‐gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40–42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.