Development of genomic resources for Nothofagus species using next‐generation sequencing data

Using next‐generation sequencing, we developed the first whole‐genome resources for two hybridizing Nothofagus species of the Patagonian forests that crucially lack genomic data, despite their ecological and industrial value. A de novo assembly strategy combining base quality control and optimization of the putative chloroplast gene map yielded ~32 000 contigs from 43% of the reads produced. With 12.5% of assembled reads, we covered ~96% of the chloroplast genome and ~70% of the mitochondrial gene content, providing functional and structural annotations for 112 and 52 genes, respectively. Functional annotation was possible on 15% of the contigs, with ~1750 potentially novel nuclear genes identified for Nothofagus species. We estimated that the new resources (13.41 Mb in total) included ~4000 gene regions representing ~6.5% of the expected genic partition of the genome, the remaining contigs potentially being nongenic DNA. A high‐quality single nucleotide polymorphisms resource was developed by comparing various filtering methods, and preliminary results indicate a strong conservation of cpDNA genomes in contrast to numerous exclusive nuclear polymorphisms in both species. Finally, we characterized 2274 potential simple sequence repeat (SSR) loci, designed primers for 769 of them and validated nine of 29 loci in 42 individuals per species. Nothofagus obliqua had more alleles (4.89) on average than N. nervosa (2.89), 8 SSRs were efficient to discriminate species, and three were successfully transferred in three other Nothofagus species. These resources will greatly help for future inferences of demographic, adaptive and hybridizing events in Nothofagus species, and for conserving and managing natural populations.     

Cunningham's skinks show low genetic connectivity and signatures of divergent selection across its distribution

Establishing corridors of connecting habitat has become a mainstay conservation strategy to maintain gene flow and facilitate climate‐driven range shifts. Yet, little attention has been given to ascertaining the extent to which corridors will benefit philopatric species, which might exhibit localized adaptation. Measures of genetic connectivity and adaptive genetic variation across species’ ranges can help fill this knowledge gap. Here, we characterized the spatial genetic structure of Cunningham's skink (Egernia cunninghami), a philopatric species distributed along Australia's Great Dividing Range, and assessed evidence of localized adaptation. Analysis of 4,274 SNPs from 94 individuals sampled at four localities spanning 500 km and 4° of latitude revealed strong genetic structuring at neutral loci (mean F_ST ± SD = 0.603 ± 0.237) among the localities. Putatively neutral SNPs and those under divergent selection yielded contrasting spatial patterns, with the latter identifying two genetically distinct clusters. Given low genetic connectivity of the four localities, we suggest that the natural movement rate of this species is insufficient to keep pace with spatial shifts to its climate envelope, irrespective of habitat availability. In addition, our finding of localized adaptation highlights the risk of outbreeding depression should the translocation of individuals be adopted as a conservation management strategy.     

Targeted genotyping by sequencing: a new way to genome profile the cat

Targeted GBS is a recent approach for obtaining an effective characterization for hundreds to thousands of markers. The high throughput of next‐generation sequencing technologies, moreover, allows sample multiplexing. The aims of this study were to (i) define a panel of single nucleotide polymorphisms (SNPs) in the cat, (ii) use GBS for profiling 16 cats, and (iii) evaluate the performance with respect to the inference using standard approaches at different coverage thresholds, thereby providing useful information for designing similar experiments. Probes for sequencing 230 variants were designed based on the Felis_catus_8.0. 8.0 genome. The regions comprised anonymous and non‐anonymous SNPs. Sixteen cat samples were analysed, some of which had already been genotyped in a large group of loci and one having been whole‐genome sequenced in the 99_Lives Cat Genome Sequencing Project. The accuracy of the method was assessed by comparing the GBS results with the genotypes already available. Overall, GBS achieved good performance, with 92–96% correct assignments, depending on the coverage threshold used to define the set of trustable genotypes. Analyses confirmed that (i) the reliability of the inference of each genotype depends on the coverage at that locus and (ii) the fraction of target loci whose genotype can be inferred correctly is a function of the total coverage. GBS proves to be a valid alternative to other methods. Data suggested a depth of less than 11× is required for greater than 95% accuracy. However, sequencing depth must be adapted to the total size of the targets to ensure proper genotype inference.     

Genome rearrangements derived from homoeologous recombination following allopolyploidy speciation in coffee

Allopolyploidization is widespread and has played a major role in flowering plant diversification. Genomic changes are common consequences of allopolyploidization, but their mechanisms of occurrence and dynamics over time are still poorly understood. Coffea arabica, a recently formed allotetraploid, was chosen as a model to investigate genetic changes in allopolyploid using an approach that exploits next‐generation sequencing technologies. Genes affected by putative homoeolog loss were inferred by comparing the numbers of single‐nucleotide polymorphisms detected using RNA‐seq in individual accessions of C. arabica, and between accessions of its two diploid progenitor species for common sequence positions. Their physical locations were investigated and clusters of genes exhibiting homoeolog loss were identified. To validate these results, genome sequencing data were generated from one accession of C. arabica and further analyzed. Genomic rearrangements involving homoeologous exchanges appear to occur in C. arabica and to be a major source of genetic diversity. At least 5% of the C. arabica genes were inferred to have undergone homoeolog loss. The detection of a large number of homoeologous exchange events (HEEs) shared by all accessions of C. arabica strongly reinforces the assumption of a single allopolyploidization event. Furthermore, HEEs were specific to one or a few accessions, suggesting that HEE accumulates gradually. Our results provide evidence for the important role of HEE in allopolyploid genome evolution.     

Remarkable life history polymorphism may be evolving under divergent selection in the silverleaf sunflower

Substantial intraspecific variation in life history is rare and potentially a signal of incipient ecological speciation, if variation is driven by geographically heterogenous natural selection. We present the first report of extensive life history polymorphism in Helianthus argophyllus, the silverleaf sunflower, and examine evidence for its evolution by divergent selection. In 18 populations sampled from across the species range and grown in a common garden, most quantitative traits covaried such that individuals could be assigned to two distinct life history syndromes: tall and late flowering with small initial flowerheads, or short and early flowering with larger initial flowerheads. Helianthus argophyllus exhibits regional genetic structure, but this population structure does not closely correspond with patterns of phenotypic variation. The early‐flowering syndrome is primarily observed in populations from coastal barrier islands, while populations from the nearby mainland coast, although geographically and genetically close, are primarily late flowering. Additionally, several traits are more differentiated among regions than expected based on neutral genetic divergence (Q_ST > F_ST), including the first principal component score corresponding with life history syndrome. This discordance between patterns of phenotypic and genetic variation suggests that divergent selection is driving genetic differences in life history across the species range. If so, the silverleaf sunflower may be in early stages of ecological speciation.     

Genome‐wide SNP analysis unveils genetic structure and phylogeographic history of snow sheep (Ovis nivicola) populations inhabiting the Verkhoyansk Mountains and Momsky Ridge (northeastern Siberia)

Insights into the genetic characteristics of a species provide important information for wildlife conservation programs. Here, we used the OvineSNP50 BeadChip developed for domestic sheep to examine population structure and evaluate genetic diversity of snow sheep (Ovis nivicola) inhabiting Verkhoyansk Range and Momsky Ridge. A total of 1,121 polymorphic SNPs were used to test 80 specimens representing five populations, including four populations of the Verkhoyansk Mountain chain: Kharaulakh Ridge–Tiksi Bay (TIK, n = 22), Orulgan Ridge (ORU, n = 22), the central part of Verkhoyansk Range (VER, n = 15), Suntar‐Khayata Ridge (SKH, n = 13), and Momsky Ridge (MOM, n = 8). We showed that the studied populations were genetically structured according to a geographic pattern. Pairwise F_ST values ranged from 0.044 to 0.205. Admixture analysis identified K = 2 as the most likely number of ancestral populations. A Neighbor‐Net tree showed that TIK was an isolated group related to the main network through ORU. TreeMix analysis revealed that TIK and MOM originated from two different ancestral populations and detected gene flow from MOM to ORU. This was supported by the f3 statistic, which showed that ORU is an admixed population with TIK and MOM/SKH heritage. Genetic diversity in the studied groups was increasing southward. Minimum values of observed (Ho) and expected (He) heterozygosity and allelic richness (Ar) were observed in the most northern population—TIK, and maximum values were observed in the most southern population—SKH. Thus, our results revealed clear genetic structure in the studied populations of snow sheep and showed that TIK has a different origin from MOM, SKH, and VER even though they are conventionally considered a single subspecies known as Yakut snow sheep (Ovis nivicola lydekkeri). Most likely, TIK was an isolated group during the Late Pleistocene glaciations of Verkhoyansk Range.     

山地虎耳草和棒腺虎耳草转录组SSR和SNP分析

利用Illumina HiSeqTM2000对山地虎耳草和棒腺虎耳草进行转录组测序,分析和比较其SSR和SNP特征。结果表明:山地虎耳草63 763条Unigene序列中含有4 622个SSR,发生频率为7.25%,有110种重复基元,平均每10.00kb出现一个SSR位点;棒腺虎耳草60 972条Unigene序列中含有4 542个SSR,发生频率为7.45%,有85种重复基元,平均每10.40kb出现一个SSR位点,略低于山地虎耳草。山地虎耳草和棒腺虎耳草转录组序列的SSR优势基元均为三核苷酸重复。2个物种的转录组SSR以5~10次的较低重复次数为主,长度主要集中于12~30bp,具有较高的多态性。山地虎耳草和棒腺虎耳草中分别获得118 424个和112 006个SNP位点,编码区的SNP位点分别占30.40%和28.59%,且在编码SNP中同义突变所占比例(30.27%、28.48%)远高于非同义突变(0.13%、0.11%)。比较发现,2个物种的各项检索结果基本一致,推测与选取的组织部位、组织的发育阶段以及物种的亲缘关系有关。     

Ecological genomics of local adaptation in Cornus florida L. by genotyping by sequencing

Discovering local adaptation, its genetic underpinnings, and environmental drivers is important for conserving forest species. Ecological genomic approaches coupled with next‐generation sequencing are useful means to detect local adaptation and uncover its underlying genetic basis in nonmodel species. We report results from a study on flowering dogwood trees (Cornus florida L.) using genotyping by sequencing (GBS). This species is ecologically important to eastern US forests but is severely threatened by fungal diseases. We analyzed subpopulations in divergent ecological habitats within North Carolina to uncover loci under local selection and associated with environmental–functional traits or disease infection. At this scale, we tested the effect of incorporating additional sequencing before scaling for a broader examination of the entire range. To test for biases of GBS, we sequenced two similarly sampled libraries independently from six populations of three ecological habitats. We obtained environmental–functional traits for each subpopulation to identify associations with genotypes via latent factor mixed modeling (LFMM) and gradient forests analysis. To test whether heterogeneity of abiotic pressures resulted in genetic differentiation indicative of local adaptation, we evaluated F_st per locus while accounting for genetic differentiation between coastal subpopulations and Piedmont‐Mountain subpopulations. Of the 54 candidate loci with sufficient evidence of being under selection among both libraries, 28–39 were Arlequin–BayeScan F_st outliers. For LFMM, 45 candidates were associated with climate (of 54), 30 were associated with soil properties, and four were associated with plant health. Reanalysis of combined libraries showed that 42 candidate loci still showed evidence of being under selection. We conclude environment‐driven selection on specific loci has resulted in local adaptation in response to potassium deficiencies, temperature, precipitation, and (to a marginal extent) disease. High allele turnover along ecological gradients further supports the adaptive significance of loci speculated to be under selection.     

Strategies for determining kinship in wild populations using genetic data

Knowledge of kin relationships between members of wild animal populations has broad application in ecology and evolution research by allowing the investigation of dispersal dynamics, mating systems, inbreeding avoidance, kin recognition, and kin selection as well as aiding the management of endangered populations. However, the assessment of kinship among members of wild animal populations is difficult in the absence of detailed multigenerational pedigrees. Here, we first review the distinction between genetic relatedness and kinship derived from pedigrees and how this makes the identification of kin using genetic data inherently challenging. We then describe useful approaches to kinship classification, such as parentage analysis and sibship reconstruction, and explain how the combined use of marker systems with biparental and uniparental inheritance, demographic information, likelihood analyses, relatedness coefficients, and estimation of misclassification rates can yield reliable classifications of kinship in groups with complex kin structures. We outline alternative approaches for cases in which explicit knowledge of dyadic kinship is not necessary, but indirect inferences about kinship on a group‐ or population‐wide scale suffice, such as whether more highly related dyads are in closer spatial proximity. Although analysis of highly variable microsatellite loci is still the dominant approach for studies on wild populations, we describe how the long‐awaited use of large‐scale single‐nucleotide polymorphism and sequencing data derived from noninvasive low‐quality samples may eventually lead to highly accurate assessments of varying degrees of kinship in wild populations.     

SNP discovery in nonmodel organisms: strand bias and base‐substitution errors reduce conversion rates

Single nucleotide polymorphisms (SNPs) have become the marker of choice for genetic studies in organisms of conservation, commercial or biological interest. Most SNP discovery projects in nonmodel organisms apply a strategy for identifying putative SNPs based on filtering rules that account for random sequencing errors. Here, we analyse data used to develop 4723 novel SNPs for the commercially important deep‐sea fish, orange roughy (Hoplostethus atlanticus), to assess the impact of not accounting for systematic sequencing errors when filtering identified polymorphisms when discovering SNPs. We used SAMtools to identify polymorphisms in a velvet assembly of genomic DNA sequence data from seven individuals. The resulting set of polymorphisms were filtered to minimize ‘bycatch’—polymorphisms caused by sequencing or assembly error. An Illumina Infinium SNP chip was used to genotype a final set of 7714 polymorphisms across 1734 individuals. Five predictors were examined for their effect on the probability of obtaining an assayable SNP: depth of coverage, number of reads that support a variant, polymorphism type (e.g. A/C), strand‐bias and Illumina SNP probe design score. Our results indicate that filtering out systematic sequencing errors could substantially improve the efficiency of SNP discovery. We show that BLASTX can be used as an efficient tool to identify single‐copy genomic regions in the absence of a reference genome. The results have implications for research aiming to identify assayable SNPs and build SNP genotyping assays for nonmodel organisms.     

Identification and characterization of more than 4 million intervarietal SNPs across the group 7 chromosomes of bread wheat

Despite being a major international crop, our understanding of the wheat genome is relatively poor due to its large size and complexity. To gain a greater understanding of wheat genome diversity, we have identified single nucleotide polymorphisms between 16 Australian bread wheat varieties. Whole‐genome shotgun Illumina paired read sequence data were mapped to the draft assemblies of chromosomes 7A, 7B and 7D to identify more than 4 million intervarietal SNPs. SNP density varied between the three genomes, with much greater density observed on the A and B genomes than the D genome. This variation may be a result of substantial gene flow from the tetraploid Triticum turgidum, which possesses A and B genomes, during early co‐cultivation of tetraploid and hexaploid wheat. In addition, we examined SNP density variation along the chromosome syntenic builds and identified genes in low‐density regions which may have been selected during domestication and breeding. This study highlights the impact of evolution and breeding on the bread wheat genome and provides a substantial resource for trait association and crop improvement. All SNP data are publically available on a generic genome browser GBrowse at www.wheatgenome.info .     

A second‐generation genetic linkage map for bighead carp (Aristichthys nobilis) based on microsatellite markers

Bighead carp (Aristichthys nobilis) is an important aquaculture fish worldwide. Genetic linkage maps for the species were previously reported, but map resolution remained to be improved. In this study, a second‐generation genetic linkage map was constructed for bighead carp through a pseudo‐testcross strategy using interspecific hybrids between bighead carp and silver carp. Of the 754 microsatellites genotyped in two interspecific mapping families (with 77 progenies for each family), 659 markers were assigned to 24 linkage groups, which were equal to the chromosome numbers of the haploid genome. The consensus map spanned 1917.3 cM covering 92.8% of the estimated bighead carp genome with an average marker interval of 2.9 cM. The length of linkage groups ranged from 52.2 to 133.5 cM with an average of 79.9 cM. The number of markers per linkage group varied from 11 to 55 with an average of 27.5 per linkage group. Normality tests on interval distances of the map showed a non‐normal marker distribution; however, significant correlation was found between the length of linkage group and the number of markers below the 0.01 significance level (two‐tailed). The length of the female map was 1.12 times that of the male map, and the average recombination ratio of female to male was 1.10:1. Visual inspection showed that distorted markers gathered in some linkage groups and in certain regions of the male and female maps. This well‐defined genetic linkage map will provide a basic framework for further genome mapping of quantitative traits, comparative mapping and marker‐assisted breeding in bighead carp.     

Discovery and characterization of novel genetic markers for use in the management of Lahontan cutthroat trout (Oncorhynchus clarkii henshawi)

The Lahontan cutthroat trout (Oncorhynchus clarkii henshawi) is threatened by habitat destruction, over‐harvest and hybridization with nonnative trout. Currently, three Geographic Management Units (GMUs) are recognized within the taxon. Here, we describe a suite of 68 single‐nucleotide polymorphism (SNP) genetic markers for use in the study and management of Lahontan cutthroat trout and a closely related subspecies, the Paiute cutthroat trout (O. c. seleneris). These include markers variable within the two subspecies (n = 35), diagnostic for the two subspecies (n = 23) and diagnostic for Yellowstone cutthroat trout (O. c. bouvieri) and other closely related subspecies (n = 10). Sixty‐three markers were discovered by Sanger sequencing of 171 EST loci in an ascertainment panel including Lahontan cutthroat trout from four populations representing all GMUs. Five markers were identified in a secondary sequencing effort with a single population of Lahontan cutthroat trout. TaqMan assays were validated on six Lahontan cutthroat trout populations and a diverse panel of other trout. Over 90% of the markers variable in Lahontan cutthroat trout were polymorphic in at least two populations, and 66% were variable within all three GMUs. All Lahontan diagnostic markers were also fixed for the Lahontan allele in Paiute cutthroat trout. Most of the Yellowstone diagnostic markers can also be used for this purpose in other cutthroat trout subspecies. This is the first set of SNP markers to be developed for Lahontan cutthroat trout, and will be an important tool for conservation and management.     

A novel genome‐wide microsatellite resource for species of Eucalyptus with linkage‐to‐physical correspondence on the reference genome sequence

Keystone species in their native ranges, eucalypts, are ecologically and genetically very diverse, growing naturally along extensive latitudinal and altitudinal ranges and variable environments. Besides their ecological importance, eucalypts are also the most widely planted trees for sustainable forestry in the world. We report the development of a novel collection of 535 microsatellites for species of Eucalyptus, 494 designed from ESTs and 41 from genomic libraries. A selected subset of 223 was evaluated for individual identification, parentage testing, and ancestral information content in the two most extensively studied species, Eucalyptus grandis and Eucalyptus globulus. Microsatellites showed high transferability and overlapping allele size range, suggesting they have arisen still in their common ancestor and confirming the extensive genome conservation between these two species. A consensus linkage map with 437 microsatellites, the most comprehensive microsatellite‐only genetic map for Eucalyptus, was built by assembling segregation data from three mapping populations and anchored to the Eucalyptus genome. An overall colinearity between recombination‐based and physical positioning of 84% of the mapped microsatellites was observed, with some ordering discrepancies and sporadic locus duplications, consistent with the recently described whole genome duplication events in Eucalyptus. The linkage map covered 95.2% of the 605.8‐Mbp assembled genome sequence, placing one microsatellite every 1.55 Mbp on average, and an overall estimate of physical to recombination distance of 618 kbp/cM. The genetic parameters estimates together with linkage and physical position data for this large set of microsatellites should assist marker choice for genome‐wide population genetics and comparative mapping in Eucalyptus.     

High indirect fitness benefits for helpers across the nesting cycle in the tropical paper wasp Polistes canadensis

Explaining the evolution of helping behaviour in the eusocial insects where nonreproductive (“worker”) individuals help raise the offspring of other individuals (“queens”) remains one of the most perplexing phenomena in the natural world. Polistes paper wasps are popular study models, as workers retain the ability to reproduce: such totipotency is likely representative of the early stages of social evolution. Polistes is thought to have originated in the tropics, where seasonal constraints on reproductive options are weak and social groups are effectively perennial. Yet, most Polistes research has focused on nontropical species, where seasonality causes family groups to disperse; cofoundresses forming new nests the following spring are often unrelated, leading to the suggestion that direct fitness through nest inheritance is key in the evolution of helping behaviour. Here, we present the first comprehensive genetic study of social structure across the perennial nesting cycle of a tropical Polistes—Polistes canadensis. Using both microsatellites and newly developed single nucleotide polymorphism markers, we show that adult cofoundresses are highly related and that brood production is monopolized by a single female across the nesting cycle. Nonreproductive cofoundresses in tropical Polistes therefore have the potential to gain high indirect fitness benefits as helpers from the outset of group formation, and these benefits persist through the nesting cycle. Direct fitness may have been less important in the origin of Polistes sociality than previously suggested. These findings stress the importance of studying a range of species with diverse life history and ecologies when considering the evolution of reproductive strategies.