Drivers of genomic landscapes of differentiation across a Populus divergence gradient.

Speciation, the continuous process by which new species form, is often investigated by looking at the variation of nucleotide diversity and differentiation across the genome (hereafter genomic landscapes). A key challenge lies in how to determine the main evolutionary forces at play shaping these patterns. One promising strategy, albeit little used to date, is to comparatively investigate these genomic landscapes as progression through time by using a series of species pairs along a divergence gradient. Here, we resequenced 201 whole-genomes from eight closely related Populus species, with pairs of species at different stages along the divergence gradient to learn more about speciation processes. Using population structure and ancestry analyses, we document extensive introgression between some species pairs, especially those with parapatric distributions. We further investigate genomic landscapes, focusing on within-species (i.e. nucleotide diversity and recombination rate) and among-species (i.e. relative and absolute divergence) summary statistics of diversity and divergence. We observe relatively conserved patterns of genomic divergence across species pairs. Independent of the stage across the divergence gradient, we find support for signatures of linked selection (i.e. the interaction between natural selection and genetic linkage) in shaping these genomic landscapes, along with gene flow and standing genetic variation. We highlight the importance of investigating genomic patterns on multiple species across a divergence gradient and discuss prospects to better understand the evolutionary forces shaping the genomic landscapes of diversity and differentiation.

Adaptive Introgression Facilitates Adaptation to High Latitudes in European Aspen (Populus tremula L.).

Understanding local adaptation has become a key research area given the ongoing climate challenge and the concomitant requirement to conserve genetic resources. Perennial plants, such as forest trees, are good models to study local adaptation given their wide geographic distribution, largely outcrossing mating systems, and demographic histories. We evaluated signatures of local adaptation in European aspen (Populus tremula) across Europe by means of whole-genome resequencing of a collection of 411 individual trees. We dissected admixture patterns between aspen lineages and observed a strong genomic mosaicism in Scandinavian trees, evidencing different colonization trajectories into the peninsula from Russia, Central and Western Europe. As a consequence of the secondary contacts between populations after the last glacial maximum, we detected an adaptive introgression event in a genome region of ∼500 kb in chromosome 10, harboring a large-effect locus that has previously been shown to contribute to adaptation to the short growing seasons characteristic of Northern Scandinavia. Demographic simulations and ancestry inference suggest an Eastern origin-probably Russian-of the adaptive Nordic allele which nowadays is present in a homozygous state at the north of Scandinavia. The strength of introgression and positive selection signatures in this region is a unique feature in the genome. Furthermore, we detected signals of balancing selection, shared across regional populations, that highlight the importance of standing variation as a primary source of alleles that facilitate local adaptation. Our results, therefore, emphasize the importance of migration-selection balance underlying the genetic architecture of key adaptive quantitative traits.

The avocado genome informs deep angiosperm phylogeny, highlights introgressive hybridization, and reveals pathogen-influenced gene space adaptation.

The avocado, Persea americana, is a fruit crop of immense importance to Mexican agriculture with an increasing demand worldwide. Avocado lies in the anciently diverged magnoliid clade of angiosperms, which has a controversial phylogenetic position relative to eudicots and monocots. We sequenced the nuclear genomes of the Mexican avocado race, P. americana var. drymifolia, and the most commercially popular hybrid cultivar, Hass, and anchored the latter to chromosomes using a genetic map. Resequencing of Guatemalan and West Indian varieties revealed that ∼39% of the Hass genome represents Guatemalan source regions introgressed into a Mexican race background. Some introgressed blocks are extremely large, consistent with the recent origin of the cultivar. The avocado lineage experienced 2 lineage-specific polyploidy events during its evolutionary history. Although gene-tree/species-tree phylogenomic results are inconclusive, syntenic ortholog distances to other species place avocado as sister to the enormous monocot and eudicot lineages combined. Duplicate genes descending from polyploidy augmented the transcription factor diversity of avocado, while tandem duplicates enhanced the secondary metabolism of the species. Phenylpropanoid biosynthesis, known to be elicited by Colletotrichum (anthracnose) pathogen infection in avocado, is one enriched function among tandems. Furthermore, transcriptome data show that tandem duplicates are significantly up- and down-regulated in response to anthracnose infection, whereas polyploid duplicates are not, supporting the general view that collections of tandem duplicates contribute evolutionarily recent "tuning knobs" in the genome adaptive landscapes of given species.

The genetic basis of adaptation in phenology in an introduced population of Black Cottonwood (Populus trichocarpa, Torr. & Gray).

BACKGROUND: Entering and exiting winter dormancy present important trade-offs between growth and survival at northern latitudes. Many forest trees display local adaptation across latitude in traits associated with these phenology transitions. Transfers of a species outside its native range introduce the species to novel combinations of environmental conditions potentially requiring different combinations of alleles to optimize growth and survival. In this study, we performed genome wide association analyses and a selection scan in a P. trichocarpa mapping population derived from crossings between clones collected across the native range and introduced into Sweden. GWAS analyses were performed using phenotypic data collected across two field seasons and in a controlled phytotron experiment. RESULTS: We uncovered 584 putative candidate genes associated with spring and autumn phenology traits as well as with growth. Many regions harboring variation significantly associated with the initiation of leaf shed and leaf autumn coloring appeared to have been evolving under positive selection in the native environments of P. trichocarpa. A comparison between the candidate genes identified with results from earlier GWAS analyses performed in the native environment found a smaller overlap for spring phenology traits than for autumn phenology traits, aligning well with earlier observations that spring phenology transitions have a more complex genetic basis than autumn phenology transitions. CONCLUSIONS: In a small and structured introduced population of P. trichocarpa, we find complex genetic architectures underlying all phenology and growth traits, and identify multiple putative candidate genes despite the limitations of the study population.

Antarease-like Zn-metalloproteases are ubiquitous in the venom of different scorpion genera.

BACKGROUND: The venoms of several scorpion species have long been associated with pancreatitis in animal models and humans. Antarease, a Zn-metalloprotease from Tityus serrulatus, is able to penetrate intact pancreatic tissue and disrupts the normal vesicular traffic necessary for secretion, so it could play a relevant role in the onset of acute pancreatitis. METHODS: The cDNA libraries from five different scorpion species were screened for antarease homologs with specific primers. The amplified PCR products were cloned and sequenced. A structural model was constructed to assess the functionality of the putative metalloproteases. A phylogenetic analysis was performed to identify clustering patterns of these venom components. RESULTS: Antarease-like sequences were amplified from all the screened cDNA libraries. The complete sequence of the antarease from T. serrulatus was obtained. The structural model of the putative antarease from Tityus trivittatus shows that it may adopt a catalytically active conformation, sharing relevant structural elements with previously reported metalloproteases of the ADAM family. The phylogenetic analysis reveals that the reported sequences cluster in groups that correlate with the geographical localization of the respective species. CONCLUSIONS: Antareases are ubiquitous to a broad range of scorpion species, where they could be catalytically active enzymes. These molecules can be used to describe the evolution of scorpion venoms under different ecogeographic constrains. GENERAL SIGNIFICANCE: For the first time the complete sequence of the antareases is reported. It is demonstrated that antareases are common in the venom of different scorpion species. They are now proposed as targets for antivenom therapies.

Synthesis of 6-substituted 1-oxoindanoyl isoleucine conjugates and modeling studies with the COI1-JAZ co-receptor complex of lima bean.

The conjugates of 6-substituted 1-oxoindanoyl carboxylic acids with L-isoleucine are mimics of the plant hormone (+)-7-iso-JA-L-Ile (3) that controls and regulates secondary metabolism and stress responses. In order to generate ligands that can be used as hormone-like compounds possessing different biological activities, an efficient and short synthesis of 6-bromo-1-oxoindane-4-carboxylic acid opens a general route to 6-Br-1-oxoindanoyl L-isoleucine conjugate (Br-In-L-Ile) (9a) as a key intermediate for several bioactive 6-halogen-In-L-Ile analogs (7a, 8a, 10a). The 6-ethynyl-In-L-Ile analog (11a) might be a valuable tool to localize macromolecular receptor molecules by click-chemistry. The activities of In-Ile derivatives were evaluated by assays inducing the release of volatile organic compounds (VOCs) in lima bean (Phaseolus lunatus). Each compound showed slightly different VOC induction patterns. To correlate such differences with structural features, modeling studies of In-Ile derivatives with COI-JAZa/b/c co-receptors of P. lunatus were performed. The modeling profits from the rigid backbone of the 1-oxoindanonoyl conjugates, which allows only well defined interactions with the receptor complex.

Dissecting the genetic basis of drought responses in common bean using natural variation.

The common bean (Phaseolus vulgaris L) is the most important legume for human consumption, contributing 30% of the total daily protein intake in developing countries. A major limitation for its cultivation is drought, which causes more than 60% of the annual losses. Among physiological adaptations to drought, delaying senescence and extending the photosynthetic capacity can improve crop productivity. This strategy is known as functional "stay-green" (SG) and has been discussed as a goal in plant breeding to alleviate the loss of yield under water scarcity conditions. The genetic components behind SG traits have been explored specially in cereals, but they are to date poorly studied in the common bean. For this, we screened 71 common bean cultivars belonging to the three most important gene-pools, Mesoamerica, Andes and Europe, selected to cover the natural variation of the species. Phenotyping experiments under terminal drought during long-days in greenhouse conditions, identified six photoperiod insensitive cultivars of European origin with a clear SG phenotype. Using SNP data produced from whole genome re-sequencing data, we obtained 10 variants significantly associated to the SG phenotype on chromosomes 1, 3, 7, 8, 9 and 10 that are in close proximity to gene models with functional annotations related to hormone signaling and anti-oxidant production. Calculating pairwise FST between subgroups of cultivars divided according to their drought response (susceptibility, escape, recovery or SG), we identified up to 29 genomic windows accounting for 1,45Mb that differentiate SG cultivars; these signals were especially strong on chromosomes 1, 5 and 10. Within these windows, we found genes directly involved in photosynthetic processes and trehalose synthesis. Altogether, these signals represent good targets for further characterization and highlight the multigenic nature of the SG response in legumes.

Comparative Study of Pine Reference Genomes Reveals Transposable Element Interconnected Gene Networks.

Sequencing the giga-genomes of several pine species has enabled comparative genomic analyses of these outcrossing tree species. Previous studies have revealed the wide distribution and extraordinary diversity of transposable elements (TEs) that occupy the large intergenic spaces in conifer genomes. In this study, we analyzed the distribution of TEs in gene regions of the assembled genomes of Pinus taeda and Pinus lambertiana using high-performance computing resources. The quality of draft genomes and the genome annotation have significant consequences for the investigation of TEs and these aspects are discussed. Several TE families frequently inserted into genes or their flanks were identified in both species' genomes. Potentially important sequence motifs were identified in TEs that could bind additional regulatory factors, promoting gene network formation with faster or enhanced transcription initiation. Node genes that contain many TEs were observed in multiple potential transposable element-associated networks. This study demonstrated the increased accumulation of TEs in the introns of stress-responsive genes of pines and suggests the possibility of rewiring them into responsive networks and sub-networks interconnected with node genes containing multiple TEs. Many such regulatory influences could lead to the adaptive environmental response clines that are characteristic of naturally spread pine populations.

The advantage of parallel selection of domestication genes to accelerate crop improvement.

A recent study identifies a locus controlling seed dormancy - a key trait of the 'domestication syndrome' - that has been selected for in parallel across multiple crop families.

Genomic history of the origin and domestication of common bean unveils its closest sister species.

BACKGROUND: Modern civilization depends on only a few plant species for its nourishment. These crops were derived via several thousands of years of human selection that transformed wild ancestors into high-yielding domesticated descendants. Among cultivated plants, common bean (Phaseolus vulgaris L.) is the most important grain legume. Yet, our understanding of the origins and concurrent shaping of the genome of this crop plant is limited. RESULTS: We sequenced the genomes of 29 accessions representing 12 Phaseolus species. Single nucleotide polymorphism-based phylogenomic analyses, using both the nuclear and chloroplast genomes, allowed us to detect a speciation event, a finding further supported by metabolite profiling. In addition, we identified ~1200 protein coding genes (PCGs) and ~100 long non-coding RNAs with domestication-associated haplotypes. Finally, we describe asymmetric introgression events occurring among common bean subpopulations in Mesoamerica and across hemispheres. CONCLUSIONS: We uncover an unpredicted speciation event in the tropical Andes that gave rise to a sibling species, formerly considered the "wild ancestor" of P. vulgaris, which diverged before the split of the Mesoamerican and Andean P. vulgaris gene pools. Further, we identify haplotypes strongly associated with genes underlying the emergence of domestication traits. Our findings also reveal the capacity of a predominantly autogamous plant to outcross and fix loci from different populations, even from distant species, which led to the acquisition by domesticated beans of adaptive traits from wild relatives. The occurrence of such adaptive introgressions should be exploited to accelerate breeding programs in the near future.

Genome and transcriptome analysis of the Mesoamerican common bean and the role of gene duplications in establishing tissue and temporal specialization of genes.

BACKGROUND: Legumes are the third largest family of angiosperms and the second most important crop class. Legume genomes have been shaped by extensive large-scale gene duplications, including an approximately 58 million year old whole genome duplication shared by most crop legumes. RESULTS: We report the genome and the transcription atlas of coding and non-coding genes of a Mesoamerican genotype of common bean (Phaseolus vulgaris L., BAT93). Using a comprehensive phylogenomics analysis, we assessed the past and recent evolution of common bean, and traced the diversification of patterns of gene expression following duplication. We find that successive rounds of gene duplications in legumes have shaped tissue and developmental expression, leading to increased levels of specialization in larger gene families. We also find that many long non-coding RNAs are preferentially expressed in germ-line-related tissues (pods and seeds), suggesting that they play a significant role in fruit development. Our results also suggest that most bean-specific gene family expansions, including resistance gene clusters, predate the split of the Mesoamerican and Andean gene pools. CONCLUSIONS: The genome and transcriptome data herein generated for a Mesoamerican genotype represent a counterpart to the genomic resources already available for the Andean gene pool. Altogether, this information will allow the genetic dissection of the characters involved in the domestication and adaptation of the crop, and their further implementation in breeding strategies for this important crop.

Global transcriptome analysis of the scorpion Centruroides noxius: new toxin families and evolutionary insights from an ancestral scorpion species.

Scorpion venoms have been studied for decades, leading to the identification of hundreds of different toxins with medical and pharmacological implications. However, little emphasis has been given to the description of these arthropods from cellular and evolutionary perspectives. In this report, we describe a transcriptomic analysis of the Mexican scorpion Centruroides noxius Hoffmann, performed with a pyrosequencing platform. Three independent sequencing experiments were carried out, each including three different cDNA libraries constructed from RNA extracted from the whole body of the scorpion after telson removal, and from the venom gland before and after venom extraction. Over three million reads were obtained and assembled in almost 19000 isogroups. Within the telson-specific sequences, 72 isogroups (0.4% of total unique transcripts) were found to be similar to toxins previously reported in other scorpion species, spiders and sea anemones. The annotation pipeline also revealed the presence of important elements of the small non-coding RNA processing machinery, as well as microRNA candidates. A phylogenomic analysis of concatenated essential genes evidenced differential evolution rates in this species, particularly in ribosomal proteins and proteasome components. Additionally, statistical comparison of transcript abundance before and after venom extraction showed that 3% and 2% of the assembled isogroups had higher expression levels in the active and replenishing gland, respectively. Thus, our sequencing and annotation strategies provide a general view of the cellular and molecular processes that take place in these arthropods, allowed the discovery of new pharmacological and biotechnological targets and uncovered several regulatory and metabolic responses behind the assembly of the scorpion venom. The results obtained in this report represent the first high-throughput study that thoroughly describes the universe of genes that are expressed in the scorpion Centruroides noxius Hoffmann, a highly relevant organism from medical and evolutionary perspectives.