A BIL Population Derived from G. hirsutum and G. barbadense Provides a Resource for Cotton Genetics and Breeding

To provide a resource for cotton genetics and breeding, an interspecific hybridization between Gossypium hirsutum cv. Emian22 and G. barbadense acc. 3–79 was made. A population of 54 BILs (backcross inbred lines, BC₁F₈) was developed with the aim of transferring G. barbadense genes into G. hirsutum in order to genetically analyze these genes’ function in a G. hirsutum background and create new germplasms for breeding. Preliminary investigation of the morphological traits showed that the BILs had diverse variations in plant architecture, seed size, and fuzz color; the related traits of yield and fiber quality evaluated in 4 environments also showed abundant phenotypic variation. In order to explore the molecular diversity of the BIL population, 446 SSR markers selected at an average genetic distance of 10 cM from our interspecific linkage map were used to genotype the BIL population. A total of 393 polymorphic loci accounting for 84.4% MAF (major allele frequency) > 0.05 and 922 allele loci were detected, and the Shannon diversity index (I) was 0.417 per locus. The average introgression segment length was 16.24 cM, and an average of 29.53 segments were introgressed in each BIL line with an average background recovery of 79.8%. QTL mapping revealed 58 QTL associated with fiber quality and yield traits, and 47 favored alleles derived from the donor parent were discovered. This study demonstrated that the interspecific BIL population was enriched with much phenotypic and molecular variation which could be a resource for cotton genetics and breeding.     

Admixture mapping of quantitative traits in Populus hybrid zones: power and limitations

Uncovering the genetic architecture of species differences is of central importance for understanding the origin and maintenance of biological diversity. Admixture mapping can be used to identify the number and effect sizes of genes that contribute to the divergence of ecologically important traits, even in taxa that are not amenable to laboratory crosses because of their long generation time or other limitations. Here, we apply admixture mapping to naturally occurring hybrids between two ecologically divergent Populus species. We map quantitative trait loci for eight leaf morphological traits using 77 mapped microsatellite markers from all 19 chromosomes of Populus. We apply multivariate linear regression analysis allowing the modeling of additive and non-additive gene action and identify several candidate genomic regions associated with leaf morphology using an information-theoretic approach. We perform simulation studies to assess the power and limitations of admixture mapping of quantitative traits in natural hybrid populations for a variety of genetic architectures and modes of gene action. Our results indicate that (1) admixture mapping has considerable power to identify the genetic architecture of species differences if sample sizes and marker densities are sufficiently high, (2) modeling of non-additive gene action can help to elucidate the discrepancy between genotype and phenotype sometimes seen in interspecific hybrids, and (3) the genetic architecture of leaf morphological traits in the studied Populus species involves complementary and overdominant gene action, providing the basis for rapid adaptation of these ecologically important forest trees.     

Geographic genetic structure of Iberian columbines (gen. <Emphasis Type="Italic">Aquilegia</Emphasis>)

Southern European columbines (genus Aquilegia) are involved in active processes of diversification, and the Iberian Peninsula offers a privileged observatory to witness the process. Studies on Iberian columbines have provided significant advances on species diversification, but we still lack a complete perspective of the genetic diversification in the Iberian scenario. This work explores how genetic diversity of the genus Aquilegia is geographically structured across the Iberian Peninsula. We used Bayesian clustering methods, principal coordinates analyses, and NJ phenograms to assess the genetic relationships among 285 individuals from 62 locations and detect the main lineages. Genetic diversity of Iberian columbines consists of five geographically structured lineages, corresponding to different Iberian taxa. Differentiation among lineages shows particularly complex admixture patterns at Northeast and highly homogeneous toward Northwest and Southeast. This geographic genetic structure suggests the existence of incomplete lineage sorting and interspecific hybridization as could be expected in recent processes of diversification under the influence of quaternary postglacial migrations. This scenario is consistent with what is proposed by the most recent studies on European and Iberian columbines, which point to geographic isolation and divergent selection by habitat specialization as the main diversification drivers of the Iberian Aquilegia complex.     

Pathway position constrains the evolution of an ecologically important pathway in aspens (Populus tremula L.)

Many ecological interactions of aspens and their relatives (Populus spp.) are affected by products of the phenylpropanoid pathway synthesizing condensed tannins (CTs), whose production involves trade‐offs with other ecologically important compounds and with growth. Genes of this pathway are candidates for investigating the role of selection on ecologically important, polygenic traits. We analysed sequences from 25 genes representing 10 steps of the CT synthesis pathway, which produces CTs used in defence and lignins used for growth, in 12 individuals of European aspen (Populus tremula). We compared these to homologs from P. trichocarpa, to a control set of 77 P. tremula genes, to genome‐wide resequencing data and to RNA‐seq expression levels, in order to identify signatures of selection distinct from those of demography. In Populus, pathway position exerts a strong influence on the evolution of these genes. Nonsynonymous diversity, divergence and allele frequency shifts (Tajima's D) were much lower than for synonymous measures. Expression levels were higher, and the direction of selection more negative, for upstream genes than for those downstream. Selective constraints act with increasing intensity on upstream genes, despite the presence of multiple paralogs in most gene families. Pleiotropy, expression level, flux control and codon bias appear to interact in determining levels and patterns of variation in genes of this pathway, whose products mediate a wide array of ecological interactions for this widely distributed species.     

Composition and functional analysis of low-molecular-weight glutenin alleles with Aroona near-isogenic lines of bread wheat

Background

The polyphenolic products of the phenylpropanoid pathway, including proanthocyanidins, anthocyanins and flavonols, possess antioxidant properties that may provide health benefits. To investigate the genetic architecture of control of their biosynthesis in apple fruit, various polyphenolic compounds were quantified in progeny from a 'Royal Gala' × 'Braeburn' apple population segregating for antioxidant content, using ultra high performance liquid chromatography of extracts derived from fruit cortex and skin.

Results

Construction of genetic maps for 'Royal Gala' and 'Braeburn' enabled detection of 79 quantitative trait loci (QTL) for content of 17 fruit polyphenolic compounds. Seven QTL clusters were stable across two years of harvest and included QTLs for content of flavanols, flavonols, anthocyanins and hydroxycinnamic acids. Alignment of the parental genetic maps with the apple whole genome sequence in silico enabled screening for co-segregation with the QTLs of a range of candidate genes coding for enzymes in the polyphenolic biosynthetic pathway. This co-location was confirmed by genetic mapping of markers derived from the gene sequences. Leucoanthocyanidin reductase (LAR1) co-located with a QTL cluster for the fruit flavanols catechin, epicatechin, procyanidin dimer and five unknown procyanidin oligomers identified near the top of linkage group (LG) 16, while hydroxy cinnamate/quinate transferase (HCT/HQT) co-located with a QTL for chlorogenic acid concentration mapping near the bottom of LG 17.

Conclusion

We conclude that LAR1 and HCT/HQT are likely to influence the concentration of these compounds in apple fruit and provide useful allele-specific markers for marker assisted selection of trees bearing fruit with healthy attributes.     

Natural phenological variation in aspen (Populus tremula): the SwAsp collection

The genus Populus is currently the main model system for genetic, genomic, and physiological research in trees. Phenotypic variation in aspen (Populus tremula) populations growing in different environments across Sweden is expected to reflect genetic variation that is important for local adaptation. To analyze such natural phenotypic and genetic variation, the Swedish Aspen (SwAsp) Collection was established. Trees were taken from 12 different populations across Sweden, from 56° to 66° latitude north and planted in two common gardens in Ekebo (55.9°N) and Sävar (63.4°N). Data related to phenological and growth traits were collected during the second year of growth. Some traits like the date of bud set and leaf area duration showed strong clinal variation patterns with latitude in both field trials, but the date of bud flush did not change along a latitudinal cline. The phenological traits showed moderate within-populations heritabilities, although growth traits showed weaker clinal patterns and lower heritabilities than the phenological traits. This research forms the starting point for the development of the SwAsp collection, a resource facilitating analysis of the natural genetic variation in aspen, the elucidation of the structure and dynamics of aspen populations, and the future identification of the genes controlling adaptive traits using association mapping of selected candidate genes.     

MicroRNA expression changes following synthesis of three full-sib <Emphasis Type="Italic">Populus</Emphasis> triploid populations with different heterozygosities

Key message

Through high-throughput sequencing, we compared the relative expression levels of miRNA in three full-sib Populus triploid populations with that in their parents and one diploid hybrid population. We found similar numbers of miRNAs differentially expressed between the parents and the four progeny hybrid populations. In addition, unbalanced parental expression level dominance of miRNAs were found in the three allotriploid and interspecific hybrid populations, which may reprogram gene expression networks and contribute to the growth of Populus hybrids. These results indicated that hybridization has a great impact on the miRNA expression variation in the newly synthesized Populus triploid and diploid hybrid populations. However, we also found no significant differences in miRNA expression among one diploid and three triploid hybrid populations, hinting that miRNA abundances do not increase with the genome content. No dosage effect of miRNA expression could lead to dosage-dependent negative effects on target genes and their downstream pathway in polyploids. We speculate that polyploids may gain advantages from the slight decrease in miRNA regulation, suggesting an important molecular mechanism of polyploid advantage.

Abstract

Hybridization with three types of induced 2n gametes transmitted different parental heterozygosities has been proven as an efficient method for Populus triploid production. Several researches have shown that miRNA could be non-additively expressed in allopolyploids. However, it is still unclear whether the non-additively expressed miRNAs result from the effect of hybridization or polyploidization, and whether a dose response to the additional genomic content exists for the expression of miRNA. Toward this end, through high-throughput sequencing, we compared the expression levels of miRNA in three full-sib Populus triploid populations with that in their parents and one interspecific hybrid population. We found similar numbers of miRNAs differentially expressed between the parents and the four progeny hybrid populations. Unbalanced parental expression level dominance of miRNAs were found in the three triploid and diploid hybrid populations, which may reprogram gene expression networks and affect the growth of Populus hybrids. These results indicated that hybridization has a great impact on the miRNA expression variation in the newly synthesized Populus triploid and diploid hybrid populations. However, we also found no significant differences in miRNA expression among the three triploid populations and the diploid hybrid population. No dosage effect of miRNA expression could lead to dosage-dependent negative effects on target genes and their downstream pathway in polyploids. We speculate that polyploids may gain advantages from the decrease in miRNA negative regulation, suggesting an important molecular mechanism of polyploid advantage.

     

Diverging pathways: Differential benzenoid and phenylpropanoid volatile production in Phlox subulata L. cultivars

Previous research into the genetic mechanisms of benzenoid and phenylpropanoid volatile biosynthesis has suggested the potential for metabolic flux, in which phenylalanine substrate dedicated to one pathway branch (i.e., benzenoid production) could alter volatile production in other pathways (i.e., phenylpropanoid production). However, little research has been conducted in planta to verify the validity of this hypothesis. We examined the emission rates of representative benzenoid and phenylpropanoid volatiles from seven cultivars of Phlox subulata L. to determine if cultivars had metabolic flux differences in terms of these two compound categories. Cultivars that produced large quantities of methyl benzoate and benzaldehyde were found to emit little or no phenylacetaldehyde and 2-phenylethanol, and vice versa. These results suggest that P. subulata cultivars experience phenylalanine substrate flux directed toward one pathway and away from the alternate branch. Such a pattern may be the result of differential selection pressures, in which gene expression has been altered to direct flux away from either benzenoid or phenylpropanoid production. Moreover, if these patterns hold true in wild populations, metabolic flux may lead to differential pollinator behavior and further phenotypic evolution. Future research using molecular tools could verify the role of metabolic flux in determining scent phenotypes and pinpoint the exact nature of the genetic mutations leading to phenotypic differences in odor.     

QTLs for Woolly Poplar Aphid (Phloeomyzus passerinii L.) Resistance Detected in an Inter-Specific Populus deltoides x P. nigra Mapping Population

The genus Populus represents one of the most economically important groups of forest trees. It is composed by approximately 30 species used for wood and non-wood products, phytoremediation and biomass. Poplar is subjected to several biological and environmental threats although, compared to annual crops, we know far less about the genetic bases of biotic stress resistance. Woolly poplar aphid (Phloeomyzus passerinii) is considered a main pest of cultivated poplars in European and American countries. In this work we present two high density linkage maps in poplar obtained by a genotyping by sequencing (GBS) approach and the identification of QTLs involved in Ph. passerinii resistance. A total of 5,667 polymorphic markers (5,606 SNPs and 61 SSRs) identified on expressed sequences have been used to genotype 131 plants of an F1 population P ×canadensis obtained by an interspecific mate between Populus deltoides (resistant to woolly poplar aphid) and Populus nigra (susceptible to woolly poplar aphid). The two linkage maps, obtained following the two-way pseudo-testcross mapping strategy, have been used to investigate the genetic bases of woolly poplar aphid resistance. One major QTL and two QTLs with minor effects (mapped on LGV, LGXVI and LG XIX) explaining the 65.8% of the genetic variance observed in the progeny in response to Ph. passerinii attack were found. The high density coverage of functional markers allowed the identification of three genes belonging to disease resistance pathway as putative candidates for P. deltoides resistance to woolly poplar aphid. This work is the first report on genetic of woolly poplar aphid genetic resistance and the resistant loci associated markers identified represent a valuable tool in resistance poplar breeding programs.     

Genetic Mapping of MAPK-Mediated Complex Traits Across S. cerevisiae

Signaling pathways enable cells to sense and respond to their environment. Many cellular signaling strategies are conserved from fungi to humans, yet their activity and phenotypic consequences can vary extensively among individuals within a species. A systematic assessment of the impact of naturally occurring genetic variation on signaling pathways remains to be conducted. In S. cerevisiae, both response and resistance to stressors that activate signaling pathways differ between diverse isolates. Here, we present a quantitative trait locus (QTL) mapping approach that enables us to identify genetic variants underlying such phenotypic differences across the genetic and phenotypic diversity of S. cerevisiae. Using a Round-robin cross between twelve diverse strains, we identified QTL that influence phenotypes critically dependent on MAPK signaling cascades. Genetic variants under these QTL fall within MAPK signaling networks themselves as well as other interconnected signaling pathways. Finally, we demonstrate how the mapping results from multiple strain background can be leveraged to narrow the search space of causal genetic variants.     

Intra- and interspecific hybridization in invasive Siberian elm

Hybridization creates unique allele combinations which can facilitate the evolution of invasiveness. Frequent interspecific hybridization between the Siberian elm, Ulmus pumila, and native elm species has been detected in the Midwestern United States, Italy and Spain. However, Ulmus pumila also occurs in the western United States and Argentina, regions where no native elm species capable of hybridizing with it occurs. We examined whether inter- or intraspecific hybridization could be detected in these regions. Nuclear markers and the program STRUCTURE helped detect interspecific hybridization and determine the population genetic structure in both the native and the two non-native ranges. Chloroplast markers identified sources of introduction into these two non-native ranges. No significant interspecific hybridization was detected between U. pumila and U. rubra in the western United States or between U. pumila and U. minor in Argentina and vice versa. However, the genetic findings supported the presence of intraspecific hybridization and high levels of genetic diversity in both non-native ranges. The evidence presented for intraspecific hybridization in the current study, combined with reports of interspecific hybridization from previous studies, identifies elm as a genus where both inter- and intraspecific hybridization may occur and help maintain high levels of genetic diversity potentially associated with invasiveness.     

Genetic investigations on the caffeine and chlorogenic acid relationship in an interspecific cross between Coffea liberica dewevrei and C. pseudozanguebariae

The objective of the paper was to identify the number of major loci explaining caffeine content in coffee seeds. Investigations were based on previously published results: (1) Caffeine binds to chlorogenic acids in a 1:1 molecular ratio; (2) Between species, the caffeine content is correlated to the chlorogenic acid content; (3) Only a part of chlorogenic acids is bound to caffeine. Especially, the content ratio between caffeine and chlorogenic acids varied between species. For identifying the number of major loci, a quantitative trait locus (QTL) approach was carried out using an interspecific cross between two highly differentiated species—Coffea liberica dewevrei and Coffea pseudozanguebariae, the latter being a caffeine-free species. As main finding, two QTLs, i.e., RCQ1 and CQA1, were identified allowing us to explain up to 97 % of the caffeine content variance. RCQ1 explained variation of the caffeine/chlorogenic acid ratio and was genetically independent of the second QTL. The latter explained the part of the caffeine content which was dependent on the chlorogenic acid content. The findings also confirmed that only a part of chlorogenic acids were trapped by caffeine, as in wild species.     

Novel nuclear-cytoplasmic interaction in wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis>) induces vigorous plants

Interspecific hybridization can be considered an accelerator of evolution, otherwise a slow process, solely dependent on mutation and recombination. Upon interspecific hybridization, several novel interactions between nuclear and cytoplasmic genomes emerge which provide additional sources of diversity. The magnitude and essence of intergenomic interactions between nuclear and cytoplasmic genomes remain unknown due to the direction of many crosses. This study was conducted to address the role of nuclear-cytoplasmic interactions as a source of variation upon hybridization. Wheat (Triticum aestivum) alloplasmic lines carrying the cytoplasm of Aegilops mutica along with an integrated approach utilizing comparative quantitative trait locus (QTL) and epigenome analysis were used to dissect this interaction. The results indicate that cytoplasmic genomes can modify the magnitude of QTL controlling certain physiological traits such as dry matter weight. Furthermore, methylation profiling analysis detected eight polymorphic regions affected by the cytoplasm type. In general, these results indicate that novel nuclear-cytoplasmic interactions can potentially trigger an epigenetic modification cascade in nuclear genes which eventually change the genetic network controlling physiological traits. These modified genetic networks can serve as new sources of variation to accelerate the evolutionary process. Furthermore, this variation can synthetically be produced by breeders in their programs to develop epigenomic-segregating lines.     

Hybridization in contact zone between temperate European pine species

Hybridization studies are important to advance our understanding of the interspecific gene flow and its evolutionary consequences in closely related species. Hybridization and admixture patterns were assessed in a contact zone and reference populations of European pine species using sequence data from 26 nuclear genes and a species-diagnostic cpDNA marker. Reference populations formed three distinct genetic clusters comprising Pinus sylvestris, Pinus mugo/Pinus uliginosa, and Pinus uncinata. Evidence of population structure was found only in P. uliginosa. Based on phenotypic characteristics and molecular data, we identified five groups of individuals in the contact zone in Poland, comprising forms of the parental species and intermediates that were most probably the result of interspecific crosses. A combination of nuclear gene sequence data and a diagnostic organelle marker were used to show that hybridization is frequent in the contact zone and results in hybrid trees with distinct phenotypic identity. The influence of selection in maintaining hybrid phenotypes in environments unsuited to parental species was inferred from nucleotide polymorphism data. A lack of admixture in reference populations suggests that hybridization has not occurred during post-glacial migration and so the contact zone represents a distinct, active example of ongoing evolution. Pine populations in this zone will be a valuable system for studying the genetic basis of hybrid advantage in environmental conditions untypical of pure parental species.     

Variable levels of introgression between the endangered Podarcis carbonelli and highly divergent congeneric species

Recent empirical studies have demonstrated that speciation with gene flow is more common than previously thought. From a conservation perspective, the potential negative effects of hybridization raise concerns on the genetic integrity of endangered species. However, introgressive hybridization has also been growingly recognized as a source of diversity and new advantageous alleles. Carbonell’s wall lizard (Podarcis carbonelli) is an endangered species whose distribution overlaps with four other congeneric species. Our goal here was to determine whether P. carbonelli is completely reproductively isolated from its congeners and to evaluate the relevance of hybridization and interspecific gene flow for developing a conservation plan. We used restriction site associated DNA (RAD) sequencing to discover SNPs in samples from four contact zones between P. carbonelli and four other species. Principal component analysis, multilocus genotype assignment and interspecific heterozygosity suggest incomplete reproductive isolation and ongoing gene flow between species. However, hybridization dynamics vary across all pairs, suggesting complex interactions between multiple intrinsic and extrinsic barriers. Despite seemingly ubiquitous interspecific gene flow, we found evidence of strong reproductive isolation across most contact zones. Instead, indirect effects of hybridization like waste of reproductive effort in small isolated populations may be more problematic. Our results highlight the need to further evaluate the consequences of introgression for P. carbonelli, both on a geographic and genomic level and included in a comprehensive and urgently needed conservation plan. Besides, those findings will add important insights on the potential effects of hybridization and introgression for endangered species.Subject terms: Population genetics, Speciation     

A peroxisomal β-oxidative pathway contributes to the formation of C6–C1 aromatic volatiles in poplar

Benzenoids (C₆–C₁ aromatic compounds) play important roles in plant defense and are often produced upon herbivory. Black cottonwood (Populus trichocarpa) produces a variety of volatile and nonvolatile benzenoids involved in various defense responses. However, their biosynthesis in poplar is mainly unresolved. We showed feeding of the poplar leaf beetle (Chrysomela populi) on P. trichocarpa leaves led to increased emission of the benzenoid volatiles benzaldehyde, benzylalcohol, and benzyl benzoate. The accumulation of salicinoids, a group of nonvolatile phenolic defense glycosides composed in part of benzenoid units, was hardly affected by beetle herbivory. In planta labeling experiments revealed that volatile and nonvolatile poplar benzenoids are produced from cinnamic acid (C₆–C₃). The biosynthesis of C₆–C₁ aromatic compounds from cinnamic acid has been described in petunia (Petunia hybrida) flowers where the pathway includes a peroxisomal-localized chain shortening sequence, involving cinnamate-CoA ligase (CNL), cinnamoyl-CoA hydratase/dehydrogenase (CHD), and 3-ketoacyl-CoA thiolase (KAT). Sequence and phylogenetic analysis enabled the identification of small CNL, CHD, and KAT gene families in P. trichocarpa. Heterologous expression of the candidate genes in Escherichia coli and characterization of purified proteins in vitro revealed enzymatic activities similar to those described in petunia flowers. RNA interference-mediated knockdown of the CNL subfamily in gray poplar (Populus x canescens) resulted in decreased emission of C₆–C₁ aromatic volatiles upon herbivory, while constitutively accumulating salicinoids were not affected. This indicates the peroxisomal β-oxidative pathway participates in the formation of volatile benzenoids. The chain shortening steps for salicinoids, however, likely employ an alternative pathway.

A three-step peroxisomal β-oxidative pathway mediates the shortening of the propyl side chain of cinnamic acid and contributes to herbivore-induced aromatic volatile formation in poplar.     

Joint linkage QTL analyses for partial resistance to Phytophthora sojae in soybean using six nested inbred populations with heterogeneous conditions

Partial resistance to Phytophthora sojae in soybean is controlled by multiple quantitative trait loci (QTL). With traditional QTL mapping approaches, power to detect such QTL, frequently of small effect, can be limited by population size. Joint linkage QTL analysis of nested recombinant inbred line (RIL) populations provides improved power to detect QTL through increased population size, recombination, and allelic diversity. However, uniform development and phenotyping of multiple RIL populations can prove difficult. In this study, the effectiveness of joint linkage QTL analysis was evaluated on combinations of two to six nested RIL populations differing in inbreeding generation, phenotypic assay method, and/or marker set used in genotyping. In comparison to linkage analysis in a single population, identification of QTL by joint linkage analysis was only minimally affected by different phenotypic methods used among populations once phenotypic data were standardized. In contrast, genotyping of populations with only partially overlapping sets of markers had a marked negative effect on QTL detection by joint linkage analysis. In total, 16 genetic regions with QTL for partial resistance against P. sojae were identified, including four novel QTL on chromosomes 4, 9, 12, and 16, as well as significant genotype-by-isolate interactions. Resistance alleles from PI 427106 or PI 427105B contributed to a major QTL on chromosome 18, explaining 10–45 % of the phenotypic variance. This case study provides guidance on the application of joint linkage QTL analysis of data collected from populations with heterogeneous assay conditions and a genetic framework for partial resistance to P. sojae.