Development and characterization of microsatellite markers for the poplar rust fungi Melampsora medusae and Melampsora larici‐populina

Melampsora species and their hybrids are obligate parasites of Populus and Salix species worldwide. The increasing interest in Populus and Salix for biomass and fibre production necessitates genetic markers for population studies of Melampsora spp. Libraries enriched for simple sequence repeats were used to develop five microsatellite markers for Melampsora medusae and Melampsora larici‐populina. The variation detected by these markers will be valuable for phylogenetic and population genetic studies, substantiate putative hybrids, and deployment of resistant poplar clones.     

A set of polymorphic EST‐derived markers for Picea species

A pooled DNA method was used to produce fully informative EST (expressed sequence tag)‐derived markers for the Picea genus. Nine markers were produced from 10 cDNA identified as candidates for cold tolerance or embryogenesis. Indels and SNPs (single nucleotide polymorphisms) were characterized from sequences obtained from pools of 10 individuals for each of the three species: Picea glauca (white spruce), Picea mariana (black spruce) and Picea abies (Norway spruce). Indels were present in 28% of the sequences and SNPs with a frequency greater than 10% were present on average in 1.2% of the positions.     

Genetic architecture of qualitative and quantitative Melampsora larici-populina leaf rust resistance in hybrid poplar: genetic mapping and QTL detection

In order to elucidate the genetic control of resistance to Melampsora larici-populina leaf rust in hybrid poplars, a Populus deltoides x P. trichocarpa F(1) progeny was analysed for qualitative and quantitative rust resistances. This progeny was evaluated for three components of quantitative resistance (latent period, uredinia number and uredinia size) to seven M. larici-populina strains in controlled conditions, and for one component of field susceptibility (rust colonization on the most infected leaf). One qualitative resistance locus inherited from P. deltoides, R(1), was localized on the genetic map. It segregates 1 : 1 in the F(1) progeny and is effective against four of the studied strains. QTL analysis was performed separately on R(1) and r(1) genotype subsets. An additional detection was conducted on the entire F(1) progeny for the three strains able to overcome R(1) and for MAX2. A total of nine QTLs were detected. Two had large, broad-spectrum effects. One (R(US)) is inherited from the P. trichocarpa parent; the other is inherited from P. deltoides and colocalized with R(1). Seven QTLs had only limited and specific effects. Significant interaction effects were detected mainly between the two major QTLs. Implications of these results for durable resistance breeding strategies, and possible benefits from the Populus genome sequence, are discussed.     

Sequencing primers and SNPs for rapidly evolving reproductive loci in endangered ibex and their kin (Bovidae,Capra spp.)

Rapidly evolving genes (e.g. candidate selected loci) are of increasing interest to molecular ecologists and conservation geneticists. Here, we report primers for five regions from three independent nuclear reproductive genes that reliably generate polymorphic sequences across the widespread wild goats of the Capra ibex species group and likely many other species of bovids. From three to nine single‐nucleotide polymorphisms (SNPs) were identified in each gene region among C. ibex subspecies. Average numbers of SNPs per 1000 bp across all five gene regions was 15.0, with a high of 25.3 in the ZP3 exons 3 and 4 sequence and a low of 6.1 in the TNP1 sequence.     

Isolation and characterization of 15 microsatellite loci in the poplar rust fungus,Melampsora larici‐populina,and cross‐amplification in related species

We developed 15 microsatellite loci in the poplar rust fungus, Melampsora larici‐populina, using two enrichment protocols. Polymorphism of each locus was assessed on a panel of 30 isolates, comprising three subpanels (world, regional and local scales). Thirteen loci were polymorphic with three to eight alleles detected. The 15 loci were also tested on five related Melampsora species, M. allii‐populina, M. medusae f. sp. deltoidae, M. larici‐tremulae, M. rostrupii and M. pinitorqua, and partial or global cross‐amplification events were detected.     

Screening for signatures of selection of Tianzhu white yak using genome‐wide re‐sequencing

The Tianzhu white yak, a domestic yak indigenous to the Qilian Mountains, migrated inland from the Qinghai‐Tibet Plateau. Specific ecological and long‐term artificial selection influenced the evolution of its pure white coat and physiological characteristics. Therefore, it is not only a natural population that represents a genomic selective region of environmental adaptability but is also an animal model for studying the pigmentation of the yak coat. A total of 24 261 829 variants, including 22 445 252 SNPs, were obtained from 29 yaks by genome‐wide re‐sequencing. According to the results of a selective sweep analysis of Tianzhu white yak in comparison to Tibetan yaks, nine candidate genes under selection in Tianzhu white yak were identified by combining π, Tajima's D, πA/πB and F_ST statistics, with threshold standards of 5%. These genes include PDCD1, NUP210, ABCG8, NEU4, LOC102287650, D2HGDH, COL4A1, RTP5 and HDAC11. Five of the nine genes were classified into 12 molecular signaling pathways, and most of these signaling pathways are involved in environmental information processing, organismal systems and metabolism. A majority of these genes has not been implicated in previous studies of yak coat color and high‐altitude animals. Our findings are helpful not only for explaining the molecular mechanism of yak coat pigmentation but also for exploring the genetic changes in Tianzhu white yak due to environmental adaptation.     

Development of EST‐derived microsatellite markers for Arabidopsis lyrata subspecies petraea (L.)

A set of expressed sequence tag–simple sequence repeat (EST‐SSR) loci has been developed for Arabidopsis lyrata ssp. petraea. From 768 root cDNA clones, 126 microsatellites, including di‐, tri‐, tetra‐ and pentanucleotide repeat motifs were identified and primers were designed to 24 EST‐SSRs. Eleven loci were subsequently screened on 150 individuals sampled from five natural populations, which revealed three to nine alleles per locus (mean 5.36) and expected heterozygosity (H_E) estimates ranging from 0.046 to 0.698. Significant deviations from random mating were observed at 10 EST‐SSR loci, likely due to inbreeding (global F_IS = 0.151) and population structure (global F_ST = 0.246).     

Signatures of diversifying selection at EST‐SSR loci and association with climate in natural Eucalyptus populations

Understanding the environmental parameters that drive adaptation among populations is important in predicting how species may respond to global climatic changes and how gene pools might be managed to conserve adaptive genetic diversity. Here, we used Bayesian F_ST outlier tests and allele–climate association analyses to reveal two Eucalyptus EST‐SSR loci as strong candidates for diversifying selection in natural populations of a southwestern Australian forest tree, Eucalyptus gomphocephala (Myrtaceae). The Eucalyptus homolog of a CONSTANS‐like gene was an F_ST outlier, and allelic variation showed significant latitudinal clinal associations with annual and winter solar radiation, potential evaporation, summer precipitation and aridity. A second F_ST outlier locus, homologous to quinone oxidoreductase, was significantly associated with measures of temperature range, high summer temperature and summer solar radiation, with important implications for predicting the effect of temperature on natural populations in the context of climate change. We complemented these data with investigations into neutral population genetic structure and diversity throughout the species range. This study provides an investigation into selection signatures at gene‐homologous EST‐SSRs in natural Eucalyptus populations, and contributes to our understanding of the relationship between climate and adaptive genetic variation, informing the conservation of both putatively neutral and adaptive components of genetic diversity.     

Association of dopamine‐ and serotonin‐related genes with canine aggression

Human‐directed canine aggression was studied using 50 aggressive and 81 non‐aggressive dogs. We examined 62 single nucleotide polymorphisms (SNPs) occurring in or in the close vicinity of 16 neurotransmitter‐related genes. Allelic associations with aggression were identified for DRD1, HTR1D, HTR2C and SLC6A1. Risk or protective haplotypes for aggressive behaviour based on 2–5 SNPs were identified. The frequency of aggressive dogs varied significantly between the haplotypes within loci and the odds ratios of aggression in dogs with risk haplotypes compared with protective haplotypes varied from 4.4 (HTR2C) to 9.0 (SLC6A1). A risk haplotype across the neurotransmitter receptor gene HTR1D harboured a non‐synonymous SNP with a potential effect on protein function. We identified no haplotypes in complete association with the recorded phenotypes, supporting a complex inheritance of aggression.     

Evolutionary and ecological genomics of non‐model plants

Abstract Dissecting evolutionary dynamics of ecologically important traits is a long‐term challenge for biologists. Attempts to understand natural variation and molecular mechanisms have motivated a move from laboratory model systems to non‐model systems in diverse natural environments. Next generation sequencing methods, along with an expansion of genomic resources and tools, have fostered new links between diverse disciplines, including molecular biology, evolution, ecology, and genomics. Great progress has been made in a few non‐model wild plants, such as Arabidopsis relatives, monkey flowers, and wild sunflowers. Until recently, the lack of comprehensive genomic information has limited evolutionary and ecological studies to larger QTL (quantitative trait locus) regions rather than single gene resolution, and has hindered recognition of general patterns of natural variation and local adaptation. Further efforts in accumulating genomic data and developing bioinformatic and biostatistical tools are now poised to move this field forward. Integrative national and international collaborations and research communities are needed to facilitate development in the field of evolutionary and ecological genomics.     

Pathogenic divergence of Central European and Australian populations of Blumeria graminis f. sp. hordei

Hosts and pathogens have adapted their response to each other through genetic changes that have arisen during the course of their co‐evolution. In developed countries the longevity of varieties is often short; new varieties frequently possess novel genes with specific resistance to pathogens. The latter must adapt to the resistance genes to maintain pathogenicity. To study this adaptation, 50 Central European and 50 Australian isolates of Blumeria graminis f. sp. hordei (Bgh) were tested on 50 barley differential varieties with different specific resistance genes. All the Central European isolates differed from each other in their virulence combinations and belonged to 50 various pathotypes, whereas Australian isolates comprised 37 pathotypes. None of the pathotypes detected in Central Europe was identical or similar to any of those in Australia. This can be attributed to the much higher number of virulences in Central European isolates that developed over a long period of contact with a range of host varieties containing specific resistance genes. This has led to a gradual divergence of the Australian and the European Bgh populations. In Europe, unlike Australia, new specific resistance genes are still widely used in breeding barley varieties and the divergence of both populations will continue.     

Interacting elevated CO2 and tropospheric O3 predisposes aspen (Populus tremuloides Michx.) to infection by rust (Melampsora medusae f. sp. tremuloidae)

We investigated the interaction of elevated CO₂ and/or (Ozone) O₃ on the occurrence and severity of aspen leaf rust (Melampsora medusae Thuem. f. sp. tremuloidae) on trembling aspen (Populus tremuloides Michx.). Furthermore, we examined the role of changes in leaf surface properties induced by elevated CO₂ and/or O₃ in this host–pathogen interaction. Three‐ to five‐fold increases in levels of rust infection index were found in 2 consecutive years following growing‐season‐long exposures with either O₃ alone or CO₂ + O₃ depending on aspen clone. Examination of leaf surface properties (wax appearance, wax amount, wax chemical composition, leaf surface and wettability) suggested significant effects by O₃ and CO₂ + O₃. We conclude that elevated O₃ is altering aspen leaf surfaces in such a way that it is likely predisposing the plants to increased infection by aspen leaf rust.     

Identification of pathogenicity‐related genes in Fusarium oxysporum f. sp. cepae

Pathogenic isolates of Fusarium oxysporum, distinguished as formae speciales (f. spp.) on the basis of their host specificity, cause crown rots, root rots and vascular wilts on many important crops worldwide. Fusarium oxysporum f. sp. cepae (FOC) is particularly problematic to onion growers worldwide and is increasing in prevalence in the UK. We characterized 31 F. oxysporum isolates collected from UK onions using pathogenicity tests, sequencing of housekeeping genes and identification of effectors. In onion seedling and bulb tests, 21 isolates were pathogenic and 10 were non‐pathogenic. The molecular characterization of these isolates, and 21 additional isolates comprising other f. spp. and different Fusarium species, was carried out by sequencing three housekeeping genes. A concatenated tree separated the F. oxysporum isolates into six clades, but did not distinguish between pathogenic and non‐pathogenic isolates. Ten putative effectors were identified within FOC, including seven Secreted In Xylem (SIX) genes first reported in F. oxysporum f. sp. lycopersici. Two highly homologous proteins with signal peptides and RxLR motifs (CRX1/CRX2) and a gene with no previously characterized domains (C5) were also identified. The presence/absence of nine of these genes was strongly related to pathogenicity against onion and all were shown to be expressed in planta. Different SIX gene complements were identified in other f. spp., but none were identified in three other Fusarium species from onion. Although the FOC SIX genes had a high level of homology with other f. spp., there were clear differences in sequences which were unique to FOC, whereas CRX1 and C5 genes appear to be largely FOC specific.     

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.     

Signatures of selection in mammalian clock genes with coding trinucleotide repeats: Implications for studying the genomics of high‐pace adaptation

Climate change is predicted to affect the reproductive ecology of wildlife; however, we have yet to understand if and how species can adapt to the rapid pace of change. Clock genes are functional genes likely critical for adaptation to shifting seasonal conditions through shifts in timing cues. Many of these genes contain coding trinucleotide repeats, which offer the potential for higher rates of change than single nucleotide polymorphisms (SNPs) at coding sites, and, thus, may translate to faster rates of adaptation in changing environments. We characterized repeats in 22 clock genes across all annotated mammal species and evaluated the potential for selection on repeat motifs in three clock genes (NR1D1, CLOCK, and PER1) in three congeneric species pairs with different latitudinal range limits: Canada lynx and bobcat (Lynx canadensis and L. rufus), northern and southern flying squirrels (Glaucomys sabrinus and G. volans), and white‐footed and deer mouse (Peromyscus leucopus and P. maniculatus). Signatures of positive selection were found in both the interspecific comparison of Canada lynx and bobcat, and intraspecific analyses in Canada lynx. Northern and southern flying squirrels showed differing frequencies at common CLOCK alleles and a signature of balancing selection. Regional excess homozygosity was found in the deer mouse at PER1 suggesting disruptive selection, and further analyses suggested balancing selection in the white‐footed mouse. These preliminary signatures of selection and the presence of trinucleotide repeats within many clock genes warrant further consideration of the importance of candidate gene motifs for adaptation to climate change.     

Two different CC‐NBS‐LRR genes are required for Lr10‐mediated leaf rust resistance in tetraploid and hexaploid wheat

Comparative study of disease resistance genes in crop plants and their relatives provides insight on resistance gene function, evolution and diversity. Here, we studied the allelic diversity of the Lr10 leaf rust resistance gene, a CC‐NBS‐LRR coding gene originally isolated from hexaploid wheat, in 20 diploid and tetraploid wheat lines. Besides a gene in the tetraploid wheat variety ‘Altar’ that is identical to the hexaploid wheat Lr10, two additional, functional resistance alleles showing sequence diversity were identified by virus‐induced gene silencing in tetraploid wheat lines. In contrast to most described NBS‐LRR proteins, the N‐terminal CC domain of LR10 was found to be under strong diversifying selection. A second NBS‐LRR gene at the Lr10 locus, RGA2, was shown through silencing to be essential for Lr10 function. Interestingly, RGA2 showed much less sequence diversity than Lr10. These data demonstrate allelic diversity of functional genes at the Lr10 locus in tetraploid wheat, and these new genes can now be analyzed for agronomic relevance. Lr10‐based resistance is highly unusual both in its dependence on two, only distantly, related CC‐NBS‐LRR proteins, as well as in the pattern of diversifying selection in the N‐terminal domain. This indicates a new and complex molecular mechanism of pathogen detection and signal transduction.     

A new genomic tool for walnut (Juglans regia L.): development and validation of the high‐density Axiom™ J. regia 700K SNP genotyping array

Over the last 20 years, global production of Persian walnut (Juglans regia L.) has grown enormously, likely reflecting increased consumption due to its numerous benefits to human health. However, advances in genome‐wide association (GWA) studies and genomic selection (GS) for agronomically important traits in walnut remain limited due to the lack of powerful genomic tools. Here, we present the development and validation of a high‐density 700K single nucleotide polymorphism (SNP) array in Persian walnut. Over 609K high‐quality SNPs have been thoroughly selected from a set of 9.6 m genome‐wide variants, previously identified from the high‐depth re‐sequencing of 27 founders of the Walnut Improvement Program (WIP) of University of California, Davis. To validate the effectiveness of the array, we genotyped a collection of 1284 walnut trees, including 1167 progeny of 48 WIP families and 26 walnut cultivars. More than half of the SNPs (55.7%) fell in the highest quality class of ‘Poly High Resolution’ (PHR) polymorphisms, which were used to assess the WIP pedigree integrity. We identified 151 new parent‐offspring relationships, all confirmed with the Mendelian inheritance test. In addition, we explored the genetic variability among cultivars of different origin, revealing how the varieties from Europe and California were differentiated from Asian accessions. Both the reconstruction of the WIP pedigree and population structure analysis confirmed the effectiveness of the Applied Biosystems? Axiom? J. regia 700K SNP array, which initiates a novel genomic and advanced phase in walnut genetics and breeding.