Demography and selection analysis of the incipient adaptive radiation of a Hawaiian woody species.

Ecological divergence in a species provides a valuable opportunity to study the early stages of speciation. We focused on Metrosideros polymorpha, a unique example of the incipient radiation of woody species, to examine how an ecological divergence continues in the face of gene flow. We analyzed the whole genomes of 70 plants collected throughout the island of Hawaii, which is the youngest island with the highest altitude in the archipelago and encompasses a wide range of environments. The continuous M. polymorpha forest stands on the island of Hawaii were differentiated into three genetic clusters, each of which grows in a distinctive environment and includes substantial genetic and phenotypic diversity. The three genetic clusters showed signatures of selection in genomic regions encompassing genes relevant to environmental adaptations, including genes associated with light utilization, oxidative stress, and leaf senescence, which are likely associated with the ecological differentiation of the species. Our demographic modeling suggested that the glaberrima cluster in wet environments maintained a relatively large population size and two clusters split: polymorpha in the subalpine zone and incana in dry and hot conditions. This ecological divergence possibly began before the species colonized the island of Hawaii. Interestingly, the three clusters recovered genetic connectivity coincidentally with a recent population bottleneck, in line with the weak reproductive isolation observed in the species. This study highlights that the degree of genetic differentiation between ecologically-diverged populations can vary depending on the strength of natural selection in the very early phases of speciation.

A global analysis of matches and mismatches between human genetic and linguistic histories.

Human history is written in both our genes and our languages. The extent to which our biological and linguistic histories are congruent has been the subject of considerable debate, with clear examples of both matches and mismatches. To disentangle the patterns of demographic and cultural transmission, we need a global systematic assessment of matches and mismatches. Here, we assemble a genomic database (GeLaTo, or Genes and Languages Together) specifically curated to investigate genetic and linguistic diversity worldwide. We find that most populations in GeLaTo that speak languages of the same language family (i.e., that descend from the same ancestor language) are also genetically highly similar. However, we also identify nearly 20% mismatches in populations genetically close to linguistically unrelated groups. These mismatches, which occur within the time depth of known linguistic relatedness up to about 10,000 y, are scattered around the world, suggesting that they are a regular outcome in human history. Most mismatches result from populations shifting to the language of a neighboring population that is genetically different because of independent demographic histories. In line with the regularity of such shifts, we find that only half of the language families in GeLaTo are genetically more cohesive than expected under spatial autocorrelations. Moreover, the genetic and linguistic divergence times of population pairs match only rarely, with Indo-European standing out as the family with most matches in our sample. Together, our database and findings pave the way for systematically disentangling demographic and cultural history and for quantifying processes of shifts in language and social identities on a global scale.

UV RESISTANCE LOCUS 8-Mediated UV-B Response Is Required Alongside CRYPTOCHROME 1 for Plant Survival in Sunlight under Field Conditions.

As sessile, photoautotrophic organisms, plants are subjected to fluctuating sunlight that includes potentially detrimental ultraviolet-B (UV-B) radiation. Experiments under controlled conditions have shown that the UV-B photoreceptor UV RESISTANCE LOCUS 8 (UVR8) controls acclimation and tolerance to UV-B in Arabidopsis thaliana; however, its long-term impact on plant fitness under naturally fluctuating environments remain poorly understood. Here, we quantified the survival and reproduction of different Arabidopsis mutant genotypes under diverse field and laboratory conditions. We found that uvr8 mutants produced more fruits than wild type when grown in growth chambers under artificial low-UV-B conditions but not under natural field conditions, indicating a fitness cost in the absence of UV-B stress. Importantly, independent double mutants of UVR8 and the blue light photoreceptor gene CRYPTOCHROME 1 (CRY1) in two genetic backgrounds showed a drastic reduction in fitness in the field. Experiments with UV-B attenuation in the field and with supplemental UV-B in growth chambers demonstrated that UV-B caused the cry1 uvr8 conditional lethal phenotype. Using RNA-seq data of field-grown single and double mutants, we explicitly identified genes showing significant statistical interaction of UVR8 and CRY1 mutations in the presence of UV-B in the field. They were enriched in Gene Ontology categories related to oxidative stress, photoprotection and DNA damage repair in addition to UV-B response. Our study demonstrates the functional importance of the UVR8-mediated response across life stages in natura, which is partially redundant with that of cry1. Moreover, these data provide an integral picture of gene expression associated with plant responses under field conditions.

Hyper-Cryptic radiation of a tropical montane plant lineage.

Species are seen as the fundamental unit of biotic diversity, and thus their delimitation is crucial for defining measures for diversity assessments and studying evolution. Differences between species have traditionally been associated with variation in morphology. And yet, the discovery of cryptic diversity suggests that the evolution of distinct lineages does not necessarily involve morphological differences. Here, we analyze 1,684,987 variant sites and over 4,000 genes for more than 400 samples to show how a tropical montane plant lineage (Geonoma undata species complex) is composed of numerous unrecognized genetic groups that are not morphologically distinct. We find that 11 to 14 clades do not correspond to the three currently recognized species. Most clades are genetically different and geographic distance and topography are the most important factors determining this genetic divergence. The genetic structure of this lineage does not match its morphological variation. Instead, this species complex constitutes the first example of a hyper-cryptic plant radiation in tropical mountains.

Homeolog expression quantification methods for allopolyploids.

Genome duplication with hybridization, or allopolyploidization, occurs in animals, fungi and plants, and is especially common in crop plants. There is an increasing interest in the study of allopolyploids because of advances in polyploid genome assembly; however, the high level of sequence similarity in duplicated gene copies (homeologs) poses many challenges. Here we compared standard RNA-seq expression quantification approaches used currently for diploid species against subgenome-classification approaches which maps reads to each subgenome separately. We examined mapping error using our previous and new RNA-seq data in which a subgenome is experimentally added (synthetic allotetraploid Arabidopsis kamchatica) or reduced (allohexaploid wheat Triticum aestivum versus extracted allotetraploid) as ground truth. The error rates in the two species were very similar. The standard approaches showed higher error rates (>10% using pseudo-alignment with Kallisto) while subgenome-classification approaches showed much lower error rates (<1% using EAGLE-RC, <2% using HomeoRoq). Although downstream analysis may partly mitigate mapping errors, the difference in methods was substantial in hexaploid wheat, where Kallisto appeared to have systematic differences relative to other methods. Only approximately half of the differentially expressed homeologs detected using Kallisto overlapped with those by any other method in wheat. In general, disagreement in low-expression genes was responsible for most of the discordance between methods, which is consistent with known biases in Kallisto. We also observed that there exist uncertainties in genome sequences and annotation which can affect each method differently. Overall, subgenome-classification approaches tend to perform better than standard approaches with EAGLE-RC having the highest precision.

Development of a high-throughput field phenotyping rover optimized for size-limited breeding fields as open-source hardware.

Phenotyping is a critical process in plant breeding, especially when there is an increasing demand for streamlining a selection process in a breeding program. Since manual phenotyping has limited efficiency, high-throughput phenotyping methods are recently popularized owing to progress in sensor and image processing technologies. However, in a size-limited breeding field, which is common in Japan and other Asian countries, it is challenging to introduce large machinery in the field or fly unmanned aerial vehicles over the field. In this study, we developed a ground-based high-throughput field phenotyping rover that could be easily introduced to a field regardless of the scale and location of the field even without special facilities. We also made the field rover open-source hardware, making its system available to public for easy modification, so that anyone can build one for their own use at a low cost. The trial run of the field rover revealed that it allowed the collection of detailed remote-sensing images of plants and quantitative analyses based on the images. The results suggest that the field rover developed in this study could allow efficient phenotyping of plants especially in a small breeding field.

KOMPEITO, an Atypical Arabidopsis Rhomboid-Related Gene, Is Required for Callose Accumulation and Pollen Wall Development.

Fertilization is a key event for sexually reproducing plants. Pollen-stigma adhesion, which is the first step in male-female interaction during fertilization, requires proper pollen wall patterning. Callose, which is a ß-1.3-glucan, is an essential polysaccharide that is required for pollen development and pollen wall formation. Mutations in CALLOSE SYNTHASE 5 (CalS5) disrupt male meiotic callose accumulation; however, how CalS5 activity and callose synthesis are regulated is not fully understood. In this paper, we report the isolation of a kompeito-1 (kom-1) mutant defective in pollen wall patterning and pollen-stigma adhesion in Arabidopsis thaliana. Callose was not accumulated in kom-1 meiocytes or microspores, which was very similar to the cals5 mutant. The KOM gene encoded a member of a subclass of Rhomboid serine protease proteins that lacked active site residues. KOM was localized to the Golgi apparatus, and both KOM and CalS5 genes were highly expressed in meiocytes. A 220 kDa CalS5 protein was detected in wild-type (Col-0) floral buds but was dramatically reduced in kom-1. These results suggested that KOM was required for CalS5 protein accumulation, leading to the regulation of meiocyte-specific callose accumulation and pollen wall formation.

ARPEGGIO: Automated Reproducible Polyploid EpiGenetic GuIdance workflOw.

BACKGROUND: Whole genome duplication (WGD) events are common in the evolutionary history of many living organisms. For decades, researchers have been trying to understand the genetic and epigenetic impact of WGD and its underlying molecular mechanisms. Particular attention was given to allopolyploid study systems, species resulting from an hybridization event accompanied by WGD. Investigating the mechanisms behind the survival of a newly formed allopolyploid highlighted the key role of DNA methylation. With the improvement of high-throughput methods, such as whole genome bisulfite sequencing (WGBS), an opportunity opened to further understand the role of DNA methylation at a larger scale and higher resolution. However, only a few studies have applied WGBS to allopolyploids, which might be due to lack of genomic resources combined with a burdensome data analysis process. To overcome these problems, we developed the Automated Reproducible Polyploid EpiGenetic GuIdance workflOw (ARPEGGIO): the first workflow for the analysis of epigenetic data in polyploids. This workflow analyzes WGBS data from allopolyploid species via the genome assemblies of the allopolyploid's parent species. ARPEGGIO utilizes an updated read classification algorithm (EAGLE-RC), to tackle the challenge of sequence similarity amongst parental genomes. ARPEGGIO offers automation, but more importantly, a complete set of analyses including spot checks starting from raw WGBS data: quality checks, trimming, alignment, methylation extraction, statistical analyses and downstream analyses. A full run of ARPEGGIO outputs a list of genes showing differential methylation. ARPEGGIO was made simple to set up, run and interpret, and its implementation ensures reproducibility by including both package management and containerization. RESULTS: We evaluated ARPEGGIO in two ways. First, we tested EAGLE-RC's performance with publicly available datasets given a ground truth, and we show that EAGLE-RC decreases the error rate by 3 to 4 times compared to standard approaches. Second, using the same initial dataset, we show agreement between ARPEGGIO's output and published results. Compared to other similar workflows, ARPEGGIO is the only one supporting polyploid data. CONCLUSIONS: The goal of ARPEGGIO is to promote, support and improve polyploid research with a reproducible and automated set of analyses in a convenient implementation. ARPEGGIO is available at https://github.com/supermaxiste/ARPEGGIO .

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes.

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Fine-scale empirical data on niche divergence and homeolog expression patterns in an allopolyploid and its diploid progenitor species.

Polyploidization is pervasive in plants, but little is known about the niche divergence of wild allopolyploids (species that harbor polyploid genomes originating from different diploid species) relative to their diploid progenitor species and the gene expression patterns that may underlie such ecological divergence. We conducted a fine-scale empirical study on habitat and gene expression of an allopolyploid and its diploid progenitors. We quantified soil properties and light availability of habitats of an allotetraploid Cardamine flexuosa and its diploid progenitors Cardamine amara and Cardamine hirsuta in two seasons. We analyzed expression patterns of genes and homeologs (homeologous gene copies in allopolyploids) using RNA sequencing. We detected niche divergence between the allopolyploid and its diploid progenitors along water availability gradient at a fine scale: the diploids in opposite extremes and the allopolyploid in a broader range between diploids, with limited overlap with diploids at both ends. Most of the genes whose homeolog expression ratio changed among habitats in C. flexuosa varied spatially and temporally. These findings provide empirical evidence for niche divergence between an allopolyploid and its diploid progenitor species at a fine scale and suggest that divergent expression patterns of homeologs in an allopolyploid may underlie its persistence in diverse habitats.

Neighbor GWAS: incorporating neighbor genotypic identity into genome-wide association studies of field herbivory.

An increasing number of field studies have shown that the phenotype of an individual plant depends not only on its genotype but also on those of neighboring plants; however, this fact is not taken into consideration in genome-wide association studies (GWAS). Based on the Ising model of ferromagnetism, we incorporated neighbor genotypic identity into a regression model, named "Neighbor GWAS". Our simulations showed that the effective range of neighbor effects could be estimated using an observed phenotype when the proportion of phenotypic variation explained (PVE) by neighbor effects peaked. The spatial scale of the first nearest neighbors gave the maximum power to detect the causal variants responsible for neighbor effects, unless their effective range was too broad. However, if the effective range of the neighbor effects was broad and minor allele frequencies were low, there was collinearity between the self and neighbor effects. To suppress the false positive detection of neighbor effects, the fixed effect and variance components involved in the neighbor effects should be tested in comparison with a standard GWAS model. We applied neighbor GWAS to field herbivory data from 199 accessions of Arabidopsis thaliana and found that neighbor effects explained 8% more of the PVE of the observed damage than standard GWAS. The neighbor GWAS method provides a novel tool that could facilitate the analysis of complex traits in spatially structured environments and is available as an R package at CRAN ( https://cran.rproject.org/package=rNeighborGWAS ).

The Current Genomic Landscape of Western South America: Andes, Amazonia, and Pacific Coast.

Studies of Native South American genetic diversity have helped to shed light on the peopling and differentiation of the continent, but available data are sparse for the major ecogeographic domains. These include the Pacific Coast, a potential early migration route; the Andes, home to the most expansive complex societies and to one of the most widely spoken indigenous language families of the continent (Quechua); and Amazonia, with its understudied population structure and rich cultural diversity. Here, we explore the genetic structure of 176 individuals from these three domains, genotyped with the Affymetrix Human Origins array. We infer multiple sources of ancestry within the Native American ancestry component; one with clear predominance on the Coast and in the Andes, and at least two distinct substrates in neighboring Amazonia, including a previously undetected ancestry characteristic of northern Ecuador and Colombia. Amazonian populations are also involved in recent gene-flow with each other and across ecogeographic domains, which does not accord with the traditional view of small, isolated groups. Long-distance genetic connections between speakers of the same language family suggest that indigenous languages here were spread not by cultural contact alone. Finally, Native American populations admixed with post-Columbian European and African sources at different times, with few cases of prolonged isolation. With our results we emphasize the importance of including understudied regions of the continent in high-resolution genetic studies, and we illustrate the potential of SNP chip arrays for informative regional-scale analysis.

Genomic signatures of convergent adaptation to Alpine environments in three Brassicaceae species.

It has long been discussed to what extent related species develop similar genetic mechanisms to adapt to similar environments. Most studies documenting such convergence have either used different lineages within species or surveyed only a limited portion of the genome. Here, we investigated whether similar or different sets of orthologous genes were involved in genetic adaptation of natural populations of three related plant species to similar environmental gradients in the Alps. We used whole-genome pooled population sequencing to study genome-wide SNP variation in 18 natural populations of three Brassicaceae (Arabis alpina, Arabidopsis halleri, and Cardamine resedifolia) from the Swiss Alps. We first de novo assembled draft reference genomes for all three species. We then ran population and landscape genomic analyses with ~3 million SNPs per species to look for shared genomic signatures of selection and adaptation in response to similar environmental gradients acting on these species. Genes with a signature of convergent adaptation were found at significantly higher numbers than expected by chance. The most closely related species pair showed the highest relative over-representation of shared adaptation signatures. Moreover, the identified genes of convergent adaptation were enriched for nonsynonymous mutations, suggesting functional relevance of these genes, even though many of the identified candidate genes have hitherto unknown or poorly described functions based on comparison with Arabidopsis thaliana. We conclude that adaptation to heterogeneous Alpine environments in related species is partly driven by convergent evolution, but that most of the genomic signatures of adaptation remain species-specific.

Fermentation Ability of Gut Microbiota of Wild Japanese Macaques in the Highland and Lowland Yakushima: In Vitro Fermentation Assay and Genetic Analyses.

Wild Japanese macaques (Macaca fuscata Blyth) living in the highland and lowland areas of Yakushima are known to have different diets, with highland individuals consuming more leaves. We aim to clarify whether and how these differences in diet are also reflected by gut microbial composition and fermentation ability. Therefore, we conduct an in vitro fermentation assay using fresh feces from macaques as inoculum and dry leaf powder of Eurya japonica Thunb. as a substrate. Fermentation activity was higher for feces collected in the highland, as evidenced by higher gas and butyric acid production and lower pH. Genetic analysis indicated separation of highland and lowland in terms of both community structure and function of the gut microbiota. Comparison of feces and suspension after fermentation indicated that the community structure changed during fermentation, and the change was larger for lowland samples. Analysis of the 16S rRNA V3-V4 barcoding region of the gut microbiota showed that community structure was clearly clustered between the two areas. Furthermore, metagenomic analysis indicated separation by gene and pathway abundance patterns. Two pathways (glycogen biosynthesis I and D-galacturonate degradation I) were enriched in lowland samples, possibly related to the fruit-eating lifestyle in the lowland. Overall, we demonstrated that the more leaf-eating highland Japanese macaques harbor gut microbiota with higher leaf fermentation ability compared with the more fruit-eating lowland ones. Broad, non-specific taxonomic and functional gut microbiome differences suggest that this pattern may be driven by a complex interplay between many taxa and pathways rather than single functional traits.

Plant trichomes and a single gene GLABRA1 contribute to insect community composition on field-grown Arabidopsis thaliana.

BACKGROUND: Genetic variation in plants alters insect abundance and community structure in the field; however, little is known about the importance of a single gene among diverse plant genotypes. In this context, Arabidopsis trichomes provide an excellent system to discern the roles of natural variation and a key gene, GLABRA1, in shaping insect communities. In this study, we transplanted two independent glabrous mutants (gl1-1 and gl1-2) and 17 natural accessions of Arabidopsis thaliana to two localities in Switzerland and Japan. RESULTS: Fifteen insect species inhabited the plant accessions, with the insect community composition significantly attributed to variations among plant accessions. The total abundance of leaf-chewing herbivores was negatively correlated with trichome density at both field sites, while glucosinolates had variable effects on leaf chewers between the sites. Interestingly, there was a parallel tendency for the abundance of leaf chewers to be higher on gl1-1 and gl1-2 than on their different parental accessions, Ler-1 and Col-0, respectively. Furthermore, the loss of function in the GLABRA1 gene significantly decreased the resistance of plants to the two predominant chewers; flea beetles and turnip sawflies. CONCLUSIONS: Overall, our results indicate that insect community composition significantly varies among A. thaliana accessions across two distant field sites, with GLABRA1 playing a key role in altering the abundance of leaf-chewing herbivores. Given that such a trichome variation is widely observed in Brassicaceae plants, the present study exemplifies the community-wide effect of a single plant gene on crucifer-feeding insects in the field.

Agrobacterium-mediated floral dip transformation of the model polyploid species Arabidopsis kamchatica.

Polyploidization has played an important role in the speciation and diversification of plant species. However, genetic analyses of polyploids are challenging because the vast majority of the model species are diploids. The allotetraploid Arabidopsis kamchatica, which originated through the hybridization of the diploid Arabidopsis halleri and Arabidopsis lyrata, is an emerging model system for studying various aspects of polyploidy. However, a transgenic method that allows the insertion of a gene of interest into A. kamchatica is still lacking. In this study, we investigated the early development of pistils in A. kamchatica and confirmed the formation of open pistils in young flower buds (stages 8-9), which is important for allowing Agrobacterium to access female reproductive tissues. We established a simple Agrobacterium-mediated floral dip transformation method to transform a gene of interest into A. kamchatica by dipping A. kamchatica inflorescences bearing many young flower buds into a 5% sucrose solution containing 0.05% Silwet L-77 and Agrobacterium harboring the gene of interest. We showed that a screenable marker comprising fluorescence-accumulating seed technology with green fluorescent protein was useful for screening the transgenic seeds of two accessions of A. kamchatica subsp. kamchatica and an accession of A. kamchatica subsp. kawasakiana.

Reference-guided de novo assembly approach improves genome reconstruction for related species.

BACKGROUND: The development of next-generation sequencing has made it possible to sequence whole genomes at a relatively low cost. However, de novo genome assemblies remain challenging due to short read length, missing data, repetitive regions, polymorphisms and sequencing errors. As more and more genomes are sequenced, reference-guided assembly approaches can be used to assist the assembly process. However, previous methods mostly focused on the assembly of other genotypes within the same species. We adapted and extended a reference-guided de novo assembly approach, which enables the usage of a related reference sequence to guide the genome assembly. In order to compare and evaluate de novo and our reference-guided de novo assembly approaches, we used a simulated data set of a repetitive and heterozygotic plant genome. RESULTS: The extended reference-guided de novo assembly approach almost always outperforms the corresponding de novo assembly program even when a reference of a different species is used. Similar improvements can be observed in high and low coverage situations. In addition, we show that a single evaluation metric, like the widely used N50 length, is not enough to properly rate assemblies as it not always points to the best assembly evaluated with other criteria. Therefore, we used the summed z-scores of 36 different statistics to evaluate the assemblies. CONCLUSIONS: The combination of reference mapping and de novo assembly provides a powerful tool to improve genome reconstruction by integrating information of a related genome. Our extension of the reference-guided de novo assembly approach enables the application of this strategy not only within but also between related species. Finally, the evaluation of genome assemblies is often not straight forward, as the truth is not known. Thus one should always use a combination of evaluation metrics, which not only try to assess the continuity but also the accuracy of an assembly.

Positional bias in variant calls against draft reference assemblies.

BACKGROUND: Whole genome resequencing projects may implement variant calling using draft reference genomes assembled de novo from short-read libraries. Despite lower quality of such assemblies, they allowed researchers to extend a wide range of population genetic and genome-wide association analyses to non-model species. As the variant calling pipelines are complex and involve many software packages, it is important to understand inherent biases and limitations at each step of the analysis. RESULTS: In this article, we report a positional bias present in variant calling performed against draft reference assemblies constructed from de Bruijn or string overlap graphs. We assessed how frequently variants appeared at each position counted from ends of a contig or scaffold sequence, and discovered unexpectedly high number of variants at the positions related to the length of either k-mers or reads used for the assembly. We detected the bias in both publicly available draft assemblies from Assemblathon 2 competition as well as in the assemblies we generated from our simulated short-read data. Simulations confirmed that the bias causing variants are predominantly false positives induced by reads from spatially distant repeated sequences. The bias is particularly strong in contig assemblies. Scaffolding does not eliminate the bias but tends to mitigate it because of the changes in variants' relative positions and alterations in read alignments. The bias can be effectively reduced by filtering out the variants that reside in repetitive elements. CONCLUSIONS: Draft genome sequences generated by several popular assemblers appear to be susceptible to the positional bias potentially affecting many resequencing projects in non-model species. The bias is inherent to the assembly algorithms and arises from their particular handling of repeated sequences. It is recommended to reduce the bias by filtering especially if higher-quality genome assembly cannot be achieved. Our findings can help other researchers to improve the quality of their variant data sets and reduce artefactual findings in downstream analyses.

Estimating genomic diversity and population differentiation - an empirical comparison of microsatellite and SNP variation in Arabidopsis halleri.

BACKGROUND: Microsatellite markers are widely used for estimating genetic diversity within and differentiation among populations. However, it has rarely been tested whether such estimates are useful proxies for genome-wide patterns of variation and differentiation. Here, we compared microsatellite variation with genome-wide single nucleotide polymorphisms (SNPs) to assess and quantify potential marker-specific biases and derive recommendations for future studies. Overall, we genotyped 180 Arabidopsis halleri individuals from nine populations using 20 microsatellite markers. Twelve of these markers were originally developed for Arabidopsis thaliana (cross-species markers) and eight for A. halleri (species-specific markers). We further characterized 2 million SNPs across the genome with a pooled whole-genome re-sequencing approach (Pool-Seq). RESULTS: Our analyses revealed that estimates of genetic diversity and differentiation derived from cross-species and species-specific microsatellites differed substantially and that expected microsatellite heterozygosity (SSR-H e) was not significantly correlated with genome-wide SNP diversity estimates (SNP-H e and θ Watterson) in A. halleri. Instead, microsatellite allelic richness (A r) was a better proxy for genome-wide SNP diversity. Estimates of genetic differentiation among populations (F ST) based on both marker types were correlated, but microsatellite-based estimates were significantly larger than those from SNPs. Possible causes include the limited number of microsatellite markers used, marker ascertainment bias, as well as the high variance in microsatellite-derived estimates. In contrast, genome-wide SNP data provided unbiased estimates of genetic diversity independent of whether genome- or only exome-wide SNPs were used. Further, we inferred that a few thousand random SNPs are sufficient to reliably estimate genome-wide diversity and to distinguish among populations differing in genetic variation. CONCLUSIONS: We recommend that future analyses of genetic diversity within and differentiation among populations use randomly selected high-throughput sequencing-based SNP data to draw conclusions on genome-wide diversity patterns. In species comparable to A. halleri, a few thousand SNPs are sufficient to achieve this goal.

Conserved but Attenuated Parental Gene Expression in Allopolyploids: Constitutive Zinc Hyperaccumulation in the Allotetraploid Arabidopsis kamchatica.

Allopolyploidization combines parental genomes and often confers broader species distribution. However, little is known about parentally transmitted gene expression underlying quantitative traits following allopolyploidization because of the complexity of polyploid genomes. The allopolyploid species Arabidopsis kamchatica is a natural hybrid of the zinc hyperaccumulator Arabidopsis halleri and of the nonaccumulator Arabidopsis lyrata We found that A. kamchatica retained the ability to hyperaccumulate zinc from A. halleri and grows in soils with both low and high metal content. Hyperaccumulation of zinc by A. kamchatica was reduced to about half of A. halleri, but is 10-fold greater than A. lyrata Homeologs derived from A. halleri had significantly higher levels of expression of genes such as HEAVY METAL ATPASE4 (HMA4), METAL TRANSPORTER PROTEIN1 and other metal ion transporters than those derived from A. lyrata, which suggests cis-regulatory differences. A. kamchatica has on average about half the expression of these genes compared with A. halleri due to fixed heterozygosity inherent in allopolyploids. Zinc treatment significantly changed the ratios of expression of 1% of homeologous pairs, including genes putatively involved in metal homeostasis. Resequencing data showed a significant reduction in genetic diversity over a large genomic region (290 kb) surrounding the HMA4 locus derived from the A. halleri parent compared with the syntenic A. lyrata-derived region, which suggests different evolutionary histories. We also estimated that three A. halleri-derived HMA4 copies are present in A. kamchatica Our findings support a transcriptomic model in which environment-related transcriptional patterns of both parents are conserved but attenuated in the allopolyploids.