Large effect quantitative trait loci for salicinoid phenolic glycosides in Populus: Implications for gene discovery

Genomic studies have been used to identify genes underlying many important plant secondary metabolic pathways. However, genes for salicinoid phenolic glycosides (SPGs)—ecologically important compounds with significant commercial, cultural, and medicinal applications—remain largely undescribed. We used a linkage map derived from a full‐sib population of hybrid cottonwoods (Populus spp.) to search for quantitative trait loci (QTL) for the SPGs salicortin and HCH‐salicortin. SSR markers and primer sequences were used to anchor the map to the V3.0 P. trichocarpa genome. We discovered 21 QTL for the two traits, including a major QTL for HCH‐salicortin (R² = .52) that colocated with a QTL for salicortin on chr12. Using the V3.0 Populus genome sequence, we identified 2,983 annotated genes and 1,480 genes of unknown function within our QTL intervals. We note ten candidate genes of interest, including a BAHD‐type acyltransferase that has been potentially linked to PopulusSPGs. Our results complement other recent studies in Populus with implications for gene discovery and the evolution of defensive chemistry in a model genus. To our knowledge, this is the first study to use a full‐sib mapping population to identify QTL intervals and gene lists associated with SPGs.     

蔷薇科果树类受体激酶的研究进展

类受体激酶(receptor like kinase,RLK)参与调控植物几乎所有的生命活动,是植物生长发育和环境适应的“中央处理器”。该文对近年来国内外有关蔷薇科果树RLK基因鉴定、进化特征及其在各器官生长发育、非生物和生物逆境中的作用及调控机制等方面的研究进展进行了综述。蔷薇科果树基因组中存在数目庞大的RLKs,不同树种间的RLK数目和各亚家族成员数目都存在较大差异,而且蔷薇科果树RLK存在极为普遍的部分重复和串联重复现象,是导致家族成员迅速变化的重要原因。有研究发现,一些RLKs调控蔷薇科果树器官发育和对环境的适应性。在器官发育方面,LRR RLK亚家族成员调控根系发育,CrRLK1L、LysM RLK和LRR RLK亚家族部分成员参与调控果实发育,CrRLK1L亚家族成员参与调控花粉管发育,LRR RLK、LysM RLK、L LEC RLK和B Lectin RLK亚家族部分成员调控蔷薇科果树对生物逆境的适应。今后RLK功能研究可侧重于蔷薇科果树特色性状,通过提高目标基因的筛选和验证的效率,加速主效RLKs的筛选进程,并通过筛选主效RLKs诱导方式和加速分子育种进程等途径,将研究成果应用于实际生产。     

Bioinformatics-Based Identification of Candidate Genes from QTLs Associated with Cell Wall Traits in Populus

Quantitative trait locus (QTL) studies are an integral part of plant research and are used to characterize the genetic basis of phenotypic variation observed in structured populations and inform marker-assisted breeding efforts. These QTL intervals can span large physical regions on a chromosome comprising hundreds of genes, thereby hampering candidate gene identification. Genome history, evolution, and expression evidence can be used to narrow the genes in the interval to a smaller list that is manageable for detailed downstream functional genomics characterization. Our primary motivation for the present study was to address the need for a research methodology that identifies candidate genes within a broad QTL interval. Here we present a bioinformatics-based approach for subdividing candidate genes within QTL intervals into alternate groups of high probability candidates. Application of this approach in the context of studying cell wall traits, specifically lignin content and S/G ratios of stem and root in Populus plants, resulted in manageable sets of genes of both known and putative cell wall biosynthetic function. These results provide a roadmap for future experimental work leading to identification of new genes controlling cell wall recalcitrance and, ultimately, in the utility of plant biomass as an energy feedstock.     

QTL‐seq approach identified genomic regions and diagnostic markers for rust and late leaf spot resistance in groundnut (Arachis hypogaea L.)

Rust and late leaf spot (LLS) are the two major foliar fungal diseases in groundnut, and their co‐occurrence leads to significant yield loss in addition to the deterioration of fodder quality. To identify candidate genomic regions controlling resistance to rust and LLS, whole‐genome resequencing (WGRS)‐based approach referred as ‘QTL‐seq’ was deployed. A total of 231.67 Gb raw and 192.10 Gb of clean sequence data were generated through WGRS of resistant parent and the resistant and susceptible bulks for rust and LLS. Sequence analysis of bulks for rust and LLS with reference‐guided resistant parent assembly identified 3136 single‐nucleotide polymorphisms (SNPs) for rust and 66 SNPs for LLS with the read depth of ≥7 in the identified genomic region on pseudomolecule A03. Detailed analysis identified 30 nonsynonymous SNPs affecting 25 candidate genes for rust resistance, while 14 intronic and three synonymous SNPs affecting nine candidate genes for LLS resistance. Subsequently, allele‐specific diagnostic markers were identified for three SNPs for rust resistance and one SNP for LLS resistance. Genotyping of one RIL population (TAG 24 × GPBD 4) with these four diagnostic markers revealed higher phenotypic variation for these two diseases. These results suggest usefulness of QTL‐seq approach in precise and rapid identification of candidate genomic regions and development of diagnostic markers for breeding applications.     

Genome‐wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305‐day milk yield in Holstein dairy cows

Milk production traits, such as 305‐day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome‐wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10^?5). A gene set enrichment analysis with SNP data (GSEA‐SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA‐SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.     

Linking plant genes to insect communities: Identifying the genetic bases of plant traits and community composition

Community genetics aims to understand the effects of intraspecific genetic variation on community composition and diversity, thereby connecting community ecology with evolutionary biology. Thus far, research has shown that plant genetics can underlie variation in the composition of associated communities (e.g., insects, lichen and endophytes), and those communities can therefore be considered as extended phenotypes. This work, however, has been conducted primarily at the plant genotype level and has not identified the key underlying genes. To address this gap, we used genome‐wide association mapping with a population of 445 aspen (Populus tremuloides) genets to identify the genes governing variation in plant traits (defence chemistry, bud phenology, leaf morphology, growth) and insect community composition. We found 49 significant SNP associations in 13 Populus genes that are correlated with chemical defence compounds and insect community traits. Most notably, we identified an early nodulin‐like protein that was associated with insect community diversity and the abundance of interacting foundation species (ants and aphids). These findings support the concept that particular plant traits are the mechanistic link between plant genes and the composition of associated insect communities. In putting the “genes” into “genes to ecosystems ecology”, this work enhances understanding of the molecular genetic mechanisms that underlie plant–insect associations and the consequences thereof for the structure of ecological communities.     

Transmission of important chromosomal regions under selection revealed in rice pedigree breeding programs

Genetic analysis across a whole plant genome based on pedigree information offers considerable potential for enhancing genetic gain from plant breeding programs through quantitative trait loci (QTL) mapping and marker-assisted selection. Here, we report its application for graphically genotyping varieties used in Chinese japonica rice (Oryza sativa L.) pedigree breeding programs. We identified 34 important chromosomal regions from the founder parent that are under selection in the breeding programs, and by comparing donor genomic regions that are under selection with QTL locations of agronomic traits, we found that QTL clustered in important genomic regions, in accordance with association analyses of natural populations and other previous studies. The convergence of genomic regions under selection with QTL locations suggests that donor genomic regions harboring key genes/QTL for important agronomic traits have been selected by plant breeders since the 1950s from the founder rice plants. The results provide better understanding of the effects of selection in breeding programs on the traits of rice cultivars. They also provide potentially valuable information for enhancing rice breeding programs through screening candidate parents for targeted molecular markers, improving crop yield potential and identifying suitable genetic material for use in future breeding programs.     

Quantitative trait locus (QTL) analysis of wood quality traits in <Emphasis Type="Italic">Eucalyptus nitens</Emphasis>

To identify the chromosomal regions affecting wood quality traits, we conducted a genome-wide quantitative trait locus (QTL) analysis of wood quality traits in Eucalyptus nitens. This information is important to exploit the full potential of the impending Eucalyptus genome sequence. A three generational mapping population consisting of 296 progeny trees was used to identify QTL associated with several wood quality traits in E. nitens. Thirty-six QTL positions for cellulose content, pulp yield, lignin content, density, and microfibril angle (MFA) were identified across different linkage groups. On linkage groups (LG)2 and 8, cellulose QTL cluster with pulp yield and extractives QTL while on LG4 and 10 cellulose and pulp yield QTLs cluster together. Similarly, on LG4, 5, and 6 QTL for lignin traits were clustered together. At two positions, QTL for MFA, a physical trait related to wood stiffness, were clustered with QTL for lignin traits. Several cell wall candidate genes were co-located to QTL positions affecting different traits. Comparative QTL analysis with Eucalyptus globulus revealed two common QTL regions for cellulose and pulp yield. The QTL positions identified in this study provide a resource for identifying wood quality genes using the impending Eucalyptus genome sequence. Candidate genes identified in this study through co-location to QTL regions may be useful in association studies.     

Whole‐genome resequencing‐based QTL‐seq identified candidate genes and molecular markers for fresh seed dormancy in groundnut

The subspecies fastigiata of cultivated groundnut lost fresh seed dormancy (FSD) during domestication and human‐made selection. Groundnut varieties lacking FSD experience precocious seed germination during harvest imposing severe losses. Development of easy‐to‐use genetic markers enables early‐generation selection in different molecular breeding approaches. In this context, one recombinant inbred lines (RIL) population (ICGV 00350 × ICGV 97045) segregating for FSD was used for deploying QTL‐seq approach for identification of key genomic regions and candidate genes. Whole‐genome sequencing (WGS) data (87.93 Gbp) were generated and analysed for the dormant parent (ICGV 97045) and two DNA pools (dormant and nondormant). After analysis of resequenced data from the pooled samples with dormant parent (reference genome), we calculated delta‐SNP index and identified a total of 10,759 genomewide high‐confidence SNPs. Two candidate genomic regions spanning 2.4 Mb and 0.74 Mb on the B05 and A09 pseudomolecules, respectively, were identified controlling FSD. Two candidate genes—RING‐H2 finger protein and zeaxanthin epoxidase—were identified in these two regions, which significantly express during seed development and control abscisic acid (ABA) accumulation. QTL‐seq study presented here laid out development of a marker, GMFSD1, which was validated on a diverse panel and could be used in molecular breeding to improve dormancy in groundnut.     

Detection and resolution of genetic loci affecting circadian period in <Emphasis Type="Italic">Brassica oleracea</Emphasis>

Circadian rhythms regulate many aspects of plant growth, fitness and vigour. The components and detailed mechanism of circadian regulation to date have been dissected in the reference species Arabidopsis thaliana. To determine the genetic basis and range of natural allelic variation for intrinsic circadian period in the closest crop relatives, we used an accurate and high throughput data capture system to record rhythmic cotyledon movement in two immortal segregating populations of Brassica oleracea, derived from parent lines representing different crop types. Periods varied between 24.4 and 26.1 h between the parent lines, with transgressive segregation between extreme recombinant lines in both populations of ∼3.5 h. The additive effect of individual QTL identified in each population varied from 0.17 to 0.36 h. QTL detected in one doubled haploid population were verified and the mapping intervals further resolved by determining circadian period in genomic substitution lines derived from the parental lines. Comparative genomic analysis based on collinearity between Brassica and Arabidopsis also allowed identification of candidate orthologous genes known to regulate period in Arabidopsis, that may account for the additive circadian effects of specific QTL. The distinct QTL positions detected in the two populations, and the extent of transgressive segregation suggest that there is likely to be considerable scope for modulating the range of available circadian periods in natural populations and crop species of Brassica. This may provide adaptive advantage for optimising growth and development in different latitudes, seasons or climate conditions.     

Genome‐wide evidence for divergent selection between populations of a major agricultural pathogen

The genetic and environmental homogeneity in agricultural ecosystems is thought to impose strong and uniform selection pressures. However, the impact of this selection on plant pathogen genomes remains largely unknown. We aimed to identify the proportion of the genome and the specific gene functions under positive selection in populations of the fungal wheat pathogen Zymoseptoria tritici. First, we performed genome scans in four field populations that were sampled from different continents and on distinct wheat cultivars to test which genomic regions are under recent selection. Based on extended haplotype homozygosity and composite likelihood ratio tests, we identified 384 and 81 selective sweeps affecting 4% and 0.5% of the 35 Mb core genome, respectively. We found differences both in the number and the position of selective sweeps across the genome between populations. Using a XtX‐based outlier detection approach, we identified 51 extremely divergent genomic regions between the allopatric populations, suggesting that divergent selection led to locally adapted pathogen populations. We performed an outlier detection analysis between two sympatric populations infecting two different wheat cultivars to identify evidence for host‐driven selection. Selective sweep regions harboured genes that are likely to play a role in successfully establishing host infections. We also identified secondary metabolite gene clusters and an enrichment in genes encoding transporter and protein localization functions. The latter gene functions mediate responses to environmental stress, including interactions with the host. The distinct gene functions under selection indicate that both local host genotypes and abiotic factors contributed to local adaptation.     

Genome‐guided investigation of plant natural product biosynthesis

The medicinal plant Madagascar periwinkle, Catharanthus roseus (L.) G. Don, produces hundreds of biologically active monoterpene‐derived indole alkaloid (MIA) metabolites and is the sole source of the potent, expensive anti‐cancer compounds vinblastine and vincristine. Access to a genome sequence would enable insights into the biochemistry, control, and evolution of genes responsible for MIA biosynthesis. However, generation of a near‐complete, scaffolded genome is prohibitive to small research communities due to the expense, time, and expertise required. In this study, we generated a genome assembly for C. roseus that provides a near‐comprehensive representation of the genic space that revealed the genomic context of key points within the MIA biosynthetic pathway including physically clustered genes, tandem gene duplication, expression sub‐functionalization, and putative neo‐functionalization. The genome sequence also facilitated high resolution co‐expression analyses that revealed three distinct clusters of co‐expression within the components of the MIA pathway. Coordinated biosynthesis of precursors and intermediates throughout the pathway appear to be a feature of vinblastine/vincristine biosynthesis. The C. roseus genome also revealed localization of enzyme‐rich genic regions and transporters near known biosynthetic enzymes, highlighting how even a draft genome sequence can empower the study of high‐value specialized metabolites.     

Sap‐feeding insects on forest trees along latitudinal gradients in northern Europe: a climate‐driven patterns

Knowledge of the latitudinal patterns in biotic interactions, and especially in herbivory, is crucial for understanding the mechanisms that govern ecosystem functioning and for predicting their responses to climate change. We used sap‐feeding insects as a model group to test the hypotheses that the strength of plant–herbivore interactions in boreal forests decreases with latitude and that this latitudinal pattern is driven primarily by midsummer temperatures. We used a replicated sampling design and quantitatively collected and identified all sap‐feeding insects from four species of forest trees along five latitudinal gradients (750–1300 km in length, ten sites in each gradient) in northern Europe (59 to 70°N and 10 to 60°E) during 2008–2011. Similar decreases in diversity of sap‐feeding insects with latitude were observed in all gradients during all study years. The sap‐feeder load (i.e. insect biomass per unit of foliar biomass) decreased with latitude in typical summers, but increased in an exceptionally hot summer and was independent of latitude during a warm summer. Analysis of combined data from all sites and years revealed dome‐shaped relationships between the loads of sap‐feeders and midsummer temperatures, peaking at 17 °C in Picea abies, at 19.5 °C in Pinus sylvestris and Betula pubescens and at 22 °C in B. pendula. From these relationships, we predict that the losses of forest trees to sap‐feeders will increase by 0–45% of the current level in southern boreal forests and by 65–210% in subarctic forests with a 1 °C increase in summer temperatures. The observed relationships between temperatures and the loads of sap‐feeders differ between the coniferous and deciduous tree species. We conclude that climate warming will not only increase plant losses to sap‐feeding insects, especially in subarctic forests, but can also alter plant‐plant interactions, thereby affecting both the productivity and the structure of future forest ecosystems.     

Mapping of quantitative trait loci controlling adaptive traits in coastal Douglas fir. III. Quantitative trait loci-by-environment interactions

Quantitative trait loci (QTL) were mapped in the woody perennial Douglas fir (Pseudotsuga menziesii var. menziesii [Mirb.] Franco) for complex traits controlling the timing of growth initiation and growth cessation. QTL were estimated under controlled environmental conditions to identify QTL interactions with photoperiod, moisture stress, winter chilling, and spring temperatures. A three-generation mapping population of 460 cloned progeny was used for genetic mapping and phenotypic evaluations. An all-marker interval mapping method was used for scanning the genome for the presence of QTL and single-factor ANOVA was used for estimating QTL-by-environment interactions. A modest number of QTL were detected per trait, with individual QTL explaining up to 9.5% of the phenotypic variation. Two QTL-by-treatment interactions were found for growth initiation, whereas several QTL-by-treatment interactions were detected among growth cessation traits. This is the first report of QTL interactions with specific environmental signals in forest trees and will assist in the identification of candidate genes controlling these important adaptive traits in perennial plants.