The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible.

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

Brassica rapa Domestication: Untangling Wild and Feral Forms and Convergence of Crop Morphotypes.

The study of domestication contributes to our knowledge of evolution and crop genetic resources. Human selection has shaped wild Brassica rapa into diverse turnip, leafy, and oilseed crops. Despite its worldwide economic importance and potential as a model for understanding diversification under domestication, insights into the number of domestication events and initial crop(s) domesticated in B. rapa have been limited due to a lack of clarity about the wild or feral status of conspecific noncrop relatives. To address this gap and reconstruct the domestication history of B. rapa, we analyzed 68,468 genotyping-by-sequencing-derived single nucleotide polymorphisms for 416 samples in the largest diversity panel of domesticated and weedy B. rapa to date. To further understand the center of origin, we modeled the potential range of wild B. rapa during the mid-Holocene. Our analyses of genetic diversity across B. rapa morphotypes suggest that noncrop samples from the Caucasus, Siberia, and Italy may be truly wild, whereas those occurring in the Americas and much of Europe are feral. Clustering, tree-based analyses, and parameterized demographic inference further indicate that turnips were likely the first crop type domesticated, from which leafy types in East Asia and Europe were selected from distinct lineages. These findings clarify the domestication history and nature of wild crop genetic resources for B. rapa, which provides the first step toward investigating cases of possible parallel selection, the domestication and feralization syndrome, and novel germplasm for Brassica crop improvement.

Genes derived from ancient polyploidy have higher genetic diversity and are associated with domestication in Brassica rapa.

Many crops are polyploid or have a polyploid ancestry. Recent phylogenetic analyses have found that polyploidy often preceded the domestication of crop plants. One explanation for this observation is that increased genetic diversity following polyploidy may have been important during the strong artificial selection that occurs during domestication. In order to test the connection between domestication and polyploidy, we identified and examined candidate genes associated with the domestication of the diverse crop varieties of Brassica rapa. Like all 'diploid' flowering plants, B. rapa has a diploidized paleopolyploid genome and experienced many rounds of whole genome duplication (WGD). We analyzed transcriptome data of more than 100 cultivated B. rapa accessions. Using a combination of approaches, we identified > 3000 candidate genes associated with the domestication of four major B. rapa crop varieties. Consistent with our expectation, we found that the candidate genes were significantly enriched with genes derived from the Brassiceae mesohexaploidy. We also observed that paleologs were significantly more diverse than non-paleologs. Our analyses find evidence for that genetic diversity derived from ancient polyploidy played a key role in the domestication of B. rapa and provide support for its importance in the success of modern agriculture.

Genomic inferences of domestication events are corroborated by written records in Brassica rapa.

Demographic modelling is often used with population genomic data to infer the relationships and ages among populations. However, relatively few analyses are able to validate these inferences with independent data. Here, we leverage written records that describe distinct Brassica rapa crops to corroborate demographic models of domestication. Brassica rapa crops are renowned for their outstanding morphological diversity, but the relationships and order of domestication remain unclear. We generated genomewide SNPs from 126 accessions collected globally using high-throughput transcriptome data. Analyses of more than 31,000 SNPs across the B. rapa genome revealed evidence for five distinct genetic groups and supported a European-Central Asian origin of B. rapa crops. Our results supported the traditionally recognized South Asian and East Asian B. rapa groups with evidence that pak choi, Chinese cabbage and yellow sarson are likely monophyletic groups. In contrast, the oil-type B. rapa subsp. oleifera and brown sarson were polyphyletic. We also found no evidence to support the contention that rapini is the wild type or the earliest domesticated subspecies of B. rapa. Demographic analyses suggested that B. rapa was introduced to Asia 2,400-4,100 years ago, and that Chinese cabbage originated 1,200-2,100 years ago via admixture of pak choi and European-Central Asian B. rapa. We also inferred significantly different levels of founder effect among the B. rapa subspecies. Written records from antiquity that document these crops are consistent with these inferences. The concordance between our age estimates of domestication events with historical records provides unique support for our demographic inferences.

QTL Analysis for Resistance to Blast Disease in U.S. Weedy Rice.

Understanding the genetic architecture of adaptation is of great importance in evolutionary biology. U.S. weedy rice is well adapted to the local conditions in U.S. rice fields. Rice blast disease is one of the most destructive diseases of cultivated rice worldwide. However, information about resistance to blast in weedy rice is limited. Here, we evaluated the disease reactions of 60 U.S. weedy rice accessions with 14 blast races, and investigated the quantitative trait loci (QTL) associated with blast resistance in two major ecotypes of U.S. weedy rice. Our results revealed that U.S. weedy rice exhibited a broad resistance spectrum. Using genotyping by sequencing, we identified 28 resistance QTL in two U.S. weedy rice ecotypes. The resistance QTL with relatively large and small effects suggest that U.S. weedy rice groups have adapted to blast disease using two methods, both major resistance (R) genes and QTL. Three genomic loci shared by some of the resistance QTL indicated that these loci may contribute to no-race-specific resistance in weedy rice. Comparing with known blast disease R genes, we found that the R genes at these resistance QTL are novel, suggesting that U.S. weedy rice is a potential source of novel blast R genes for resistant breeding.

More than one way to evolve a weed: parallel evolution of US weedy rice through independent genetic mechanisms.

Many different crop species were selected for a common suite of 'domestication traits', which facilitates their use for studies of parallel evolution. Within domesticated rice (Oryza sativa), there has also been independent evolution of weedy strains from different cultivated varieties. This makes it possible to examine the genetic basis of parallel weed evolution and the extent to which this process occurs through shared genetic mechanisms. We performed comparative QTL mapping of weediness traits using two recombinant inbred line populations derived from crosses between an indica crop variety and representatives of each of the two independently evolved weed strains found in US rice fields, strawhull (S) and blackhull awned (B). Genotyping-by-sequencing provided dense marker coverage for linkage map construction (average marker interval <0.25 cM), with 6016 and 13 730 SNPs mapped in F5 lines of the S and B populations, respectively. For some weediness traits (awn length, hull pigmentation and pericarp pigmentation), QTL mapping and sequencing of underlying candidate genes confirmed that trait variation was largely attributable to individual loci. However, for more complex quantitative traits (including heading date, panicle length and seed shattering), we found multiple QTL, with little evidence of shared genetic bases between the S and B populations or across previous studies of weedy rice. Candidate gene sequencing revealed causal genetic bases for 8 of 27 total mapped QTL. Together these findings suggest that despite the genetic bottleneck that occurred during rice domestication, there is ample genetic variation in this crop to allow agricultural weed evolution through multiple genetic mechanisms.

A strong 'filter' effect of the East China Sea land bridge for East Asia's temperate plant species: inferences from molecular phylogeography and ecological niche modelling of Platycrater arguta (Hydrangeaceae).

BACKGROUND: In East Asia, an increasing number of studies on temperate forest tree species find evidence for migration and gene exchange across the East China Sea (ECS) land bridge up until the last glacial maximum (LGM). However, it is less clear when and how lineages diverged in this region, whether in full isolation or in the face of post-divergence gene flow. Here, we investigate the effects of Quaternary changes in climate and sea level on the evolutionary and demographic history of Platycrater arguta, a rare temperate understorey shrub with disjunct distributions in East China (var. sinensis) and South Japan (var. arguta). Molecular data were obtained from 14 P. arguta populations to infer current patterns of molecular structure and diversity in relation to past (Last Interglacial and Last Glacial Maximum) and present distributions based on ecological niche modelling (ENM). A coalescent-based isolation-with-migration (IM) model was used to estimate lineage divergence times and population demographic parameters. RESULTS: Combining information from nuclear/chloroplast sequence data with nuclear microsatellites, our IM analyses identify the two varieties as genetically distinct units that evolved in strict allopatry since the mid-Pleistocene, c. 0.89 (0.51-1.2) Ma. Together with Bayesian Skyeline Plots, our data further suggest that both lineages experienced post-divergence demographic growth, followed by refugial isolation, divergence, and in the case of var. arguta post-glacial admixture. However, past species distribution modelling indicates that the species' overall distribution has not greatly changed over the last glacial cycles. CONCLUSIONS: Our findings highlight the important influence of ancient sea-level changes on the diversification of East Asia's temperate flora. Implicitly, they challenge the notion of general temperate forest expansion across the ECS land bridge, demonstrating instead its 'filter' effect owing to an unsuitable environment for certain species and their biological (e.g., recruitment) properties.

Phylogeny and evolutionary history of glycogen synthase kinase 3/SHAGGY-like kinase genes in land plants.

BACKGROUND: GSK3 (glycogen synthase kinase 3) genes encode signal transduction proteins with roles in a variety of biological processes in eukaryotes. In contrast to the low copy numbers observed in animals, GSK3 genes have expanded into a multi-gene family in land plants (embryophytes), and have also evolved functions in diverse plant specific processes, including floral development in angiosperms. However, despite previous efforts, the phylogeny of land plant GSK3 genes is currently unclear. Here, we analyze genes from a representative sample of phylogenetically pivotal taxa, including basal angiosperms, gymnosperms, and monilophytes, to reconstruct the evolutionary history and functional diversification of the GSK3 gene family in land plants. RESULTS: Maximum Likelihood phylogenetic analyses resolve a gene tree with four major gene duplication events that coincide with the emergence of novel land plant clades. The single GSK3 gene inherited from the ancestor of land plants was first duplicated along the ancestral branch to extant vascular plants, and three subsequent duplications produced three GSK3 loci in the ancestor of euphyllophytes, four in the ancestor of seed plants, and at least five in the ancestor of angiosperms. A single gene in the Amborella trichopoda genome may be the sole survivor of a sixth GSK3 locus that originated in the ancestor of extant angiosperms. Homologs of two Arabidopsis GSK3 genes with genetically confirmed roles in floral development, AtSK11 and AtSK12, exhibit floral preferential expression in several basal angiosperms, suggesting evolutionary conservation of their floral functions. Members of other gene lineages appear to have independently evolved roles in plant reproductive tissues in individual taxa. CONCLUSIONS: Our phylogenetic analyses provide the most detailed reconstruction of GSK3 gene evolution in land plants to date and offer new insights into the origins, relationships, and functions of family members. Notably, the diversity of this "green" branch of the gene family has increased in concert with the increasing morphological and physiological complexity of land plant life forms. Expression data for seed plants indicate that the functions of GSK3 genes have also diversified during evolutionary time.

Molecular data and ecological niche modelling reveal a highly dynamic evolutionary history of the East Asian Tertiary relict Cercidiphyllum (Cercidiphyllaceae).

East Asia's temperate deciduous forests served as sanctuary for Tertiary relict trees, but their ages and response to past climate change remain largely unknown. To address this issue, we elucidated the evolutionary and population demographic history of Cercdiphyllum, comprising species in China/Japan (Cercdiphyllum japonicum) and central Japan (Cercdiphyllum magnificum). Fifty-three populations were genotyped using chloroplast and ribosomal DNA sequences and microsatellite loci to assess molecular structure and diversity in relation to past (Last Glacial Maximum) and present distributions based on ecological niche modelling. Late Tertiary climate cooling was reflected in a relatively recent speciation event, dated at the Mio-/Pliocene boundary. During glacials, the warm-temperate C. japonicum experienced massive habitat losses in some areas (north-central China/north Japan) but increases in others (southwest/-east China, East China Sea landbridge, south Japan). In China, the Sichuan Basin and/or the middle-Yangtze were source areas of postglacial northward recolonization; in Japan, this may have been facilitated through introgressive hybridization with the cool-temperate C. magnificum. Our findings challenge the notion of relative evolutionary and demographic stability of Tertiary relict trees, and may serve as a guideline for assessing the impact of Neogene climate change on the evolution and distribution of East Asian temperate plants.

Climate oscillation during the Quaternary associated with landscape heterogeneity promoted allopatric lineage divergence of a temperate tree Kalopanax septemlobus (Araliaceae) in East Asia.

We investigated the biogeographic history of Kalopanax septemlobus, one of the most widespread temperate tree species in East Asia, using a combined phylogeographic and palaeodistribution modelling approach. Range-wide genetic differentiation at nuclear microsatellites (G'(ST) = 0.709; 2205 samples genotyped at five loci) and chloroplast DNA (G(ST) = 0.697; 576 samples sequenced for 2055 bp at three fragments) was high. A major phylogeographic break in Central China corresponded with those of other temperate species and the spatial delineation of the two temperate forest subkingdoms of East Asia, consistent with the forests having been isolated within both East and West China for multiple glacial-interglacial cycles. Evidence for multiple glacial refugia was found in most of its current range in China, South Japan and the southernmost part of the Korean Peninsula. In contrast, lineage admixture and absence of private alleles and haplotypes in Hokkaido and the northern Korean Peninsula support a postglacial origin of northernmost populations. Although palaeodistribution modelling predicted suitable climate across a land-bridge extending from South Japan to East China during the Last Glacial Maximum, the genetic differentiation of regional populations indicated a limited role of the exposed sea floor as a dispersal corridor at that time. Overall, this study provides evidence that differential impacts of Quaternary climate oscillation associated with landscape heterogeneity have shaped the genetic structure of a wide-ranging temperate tree in East Asia.

Development of 12 microsatellite markers for Platycrater arguta (Hydrangeaceae) endemic to East Asia.

PREMISE OF THE STUDY: Microsatellite markers were characterized in Platycrater arguta, a rare endemic shrub from eastern China and southern Japan, to investigate its population genetic structure, phylogeographic patterns, and breeding system. METHODS AND RESULTS: Using both the Fast Isolation by AFLP of Sequences COntaining repeats (FIASCO) protocol and the compound microsatellite marker technique, we developed 12 microsatellite markers. All primers showed polymorphism when assessed in 43 individuals from two populations in eastern China and southern Japan. Overall, the number of alleles ranged from five to 21, with an average of 14.91 alleles per locus. CONCLUSIONS: These markers can be used in further studies on population genetics, phylogeography, and breeding system of P. arguta.

Genetic effects of recent habitat fragmentation in the Thousand-Island Lake region of southeast China on the distylous herb Hedyotis chrysotricha (Rubiaceae).

PREMISE OF THE STUDY: Known-age artificial-lake islands provide ideal model systems to elucidate the genetic and evolutionary consequences of anthropogenic habitat fragmentation on very recent time scales. Here, we studied a distylous herb, Hedyotis chrysotricha (Rubiaceae), in the artificially created Thousand-Island Lake (TIL) region of southeast China to explore the genetic consequences of islanding for this species. • METHODS: Seven microsatellite loci were used to genotype 384 individuals of H. chrysotricha from 18 populations to estimate genetic diversity, population structure, and demographic parameters. • KEY RESULTS: Island populations had significantly lower mean genetic diversity than those from the western/eastern mainland (e.g., H(E) = 0.381 vs. 0.461) and also displayed higher mean subdivision (F(ST) = 0.12 vs. 0.042/0.051). BayesAss analyses indicated moderate levels of migration rates among most populations, whereas Bottleneck did not provide strong evidence for such effects. In consequence, 2MOD strongly favored a gene flow-drift model over a pure drift model in the study area, but concomitantly revealed a relatively greater influence of drift in the island populations as evidenced by their significantly higher probabilities of allelic coancestry (F = 0.184 vs. 0.085). • CONCLUSIONS: The observed genetic patterns in H. chrysotricha indicate that recent anthropogenic habitat fragmentation in the TIL region can lead to significant loss of genetic diversity in isolated fragments (islands) due to ongoing drift. By contrast, patterns of random mating, gene flow, and population connectivity have not greatly been modified yet, possibly owing to the species' fruit (seed) dispersal capabilities providing resilience in the face of habitat fragmentation.

Transcriptomic profiles in pulmonary arterial hypertension associate with disease severity and identify novel candidate genes.

Using RNAseq, we identified a 61 gene-based circulating transcriptomic profile most correlated with four indices of pulmonary arterial hypertension severity. In an independent dataset, 13/61 (21%) genes were differentially expressed in lung tissues of pulmonary arterial hypertension cases versus controls, highlighting potentially novel candidate genes involved in pulmonary arterial hypertension development.

Call of the wild rice: Oryza rufipogon shapes weedy rice evolution in Southeast Asia.

Agricultural weeds serve as productive models for studying the genetic basis of rapid adaptation, with weed-adaptive traits potentially evolving multiple times independently in geographically distinct but environmentally similar agroecosystems. Weedy relatives of domesticated crops can be especially interesting systems because of the potential for weed-adaptive alleles to originate through multiple mechanisms, including introgression from cultivated and/or wild relatives, standing genetic variation, and de novo mutations. Weedy rice populations have evolved multiple times through dedomestication from cultivated rice. Much of the genomic work to date in weedy rice has focused on populations that exist outside the range of the wild crop progenitor. In this study, we use genome-wide SNPs generated through genotyping-by-sequencing to compare the evolution of weedy rice in regions outside the range of wild rice (North America, South Korea) and populations in Southeast Asia, where wild rice populations are present. We find evidence for adaptive introgression of wild rice alleles into weedy rice populations in Southeast Asia, with the relative contributions of wild and cultivated rice alleles varying across the genome. In addition, gene regions underlying several weed-adaptive traits are dominated by genomic contributions from wild rice. Genome-wide nucleotide diversity is also much higher in Southeast Asian weeds than in North American and South Korean weeds. Besides reflecting introgression from wild rice, this difference in diversity likely reflects genetic contributions from diverse cultivated landraces that may have served as the progenitors of these weedy populations. These important differences in weedy rice evolution in regions with and without wild rice could inform region-specific management strategies for weed control.

Transcriptome and organellar sequencing highlights the complex origin and diversification of allotetraploid Brassica napus.

Brassica napus, an allotetraploid crop, is hypothesized to be a hybrid from unknown varieties of Brassica rapa and Brassica oleracea. Despite the economic importance of B. napus, much is unresolved regarding its phylogenomic relationships, genetic structure, and diversification. Here we conduct a comprehensive study among diverse accessions from 183 B. napus (including rapeseed, rutabaga, and Siberian kale), 112 B. rapa, and 62 B. oleracea and its wild relatives. Using RNA-seq of B. napus accessions, we define the genetic diversity and sub-genome variance of six genetic clusters. Nuclear and organellar phylogenies for B. napus and its progenitors reveal varying patterns of inheritance and post-formation introgression. We discern regions with signatures of selective sweeps and detect 8,187 differentially expressed genes with implications for B. napus diversification. This study highlights the complex origin and evolution of B. napus providing insights that can further facilitate B. napus breeding and germplasm preservation.

Empirical Comparisons of Different Statistical Models To Identify and Validate Kernel Row Number-Associated Variants from Structured Multi-parent Mapping Populations of Maize.

Advances in next generation sequencing technologies and statistical approaches enable genome-wide dissection of phenotypic traits via genome-wide association studies (GWAS). Although multiple statistical approaches for conducting GWAS are available, the power and cross-validation rates of many approaches have been mostly tested using simulated data. Empirical comparisons of single variant (SV) and multi-variant (MV) GWAS approaches have not been conducted to test if a single approach or a combination of SV and MV is effective, through identification and cross-validation of trait-associated loci. In this study, kernel row number (KRN) data were collected from a set of 6,230 entries derived from the Nested Association Mapping (NAM) population and related populations. Three different types of GWAS analyses were performed: 1) single-variant (SV), 2) stepwise regression (STR) and 3) a Bayesian-based multi-variant (BMV) model. Using SV, STR, and BMV models, 257, 300, and 442 KRN-associated variants (KAVs) were identified in the initial GWAS analyses. Of these, 231 KAVs were subjected to genetic validation using three unrelated populations that were not included in the initial GWAS. Genetic validation results suggest that the three GWAS approaches are complementary. Interestingly, KAVs in low recombination regions were more likely to exhibit associations in independent populations than KAVs in recombinationally active regions, probably as a consequence of linkage disequilibrium. The KAVs identified in this study have the potential to enhance our understanding of the genetic basis of ear development.

A targeted enrichment strategy for massively parallel sequencing of angiosperm plastid genomes.

PREMISE OF THE STUDY: We explored a targeted enrichment strategy to facilitate rapid and low-cost next-generation sequencing (NGS) of numerous complete plastid genomes from across the phylogenetic breadth of angiosperms. • METHODS AND RESULTS: A custom RNA probe set including the complete sequences of 22 previously sequenced eudicot plastomes was designed to facilitate hybridization-based targeted enrichment of eudicot plastid genomes. Using this probe set and an Agilent SureSelect targeted enrichment kit, we conducted an enrichment experiment including 24 angiosperms (22 eudicots, two monocots), which were subsequently sequenced on a single lane of the Illumina GAIIx with single-end, 100-bp reads. This approach yielded nearly complete to complete plastid genomes with exceptionally high coverage (mean coverage: 717×), even for the two monocots. • CONCLUSIONS: Our enrichment experiment was highly successful even though many aspects of the capture process employed were suboptimal. Hence, significant improvements to this methodology are feasible. With this general approach and probe set, it should be possible to sequence more than 300 essentially complete plastid genomes in a single Illumina GAIIx lane (achieving ∼50× mean coverage). However, given the complications of pooling numerous samples for multiplex sequencing and the limited number of barcodes (e.g., 96) available in commercial kits, we recommend 96 samples as a current practical maximum for multiplex plastome sequencing. This high-throughput approach should facilitate large-scale plastid genome sequencing at any level of phylogenetic diversity in angiosperms.