Genome-wide patterns of nucleotide polymorphism in domesticated rice

Biases in ARG-Based Inference of Historical Population Size in Populations Experiencing Selection

Inferring the demographic history of populations provides fundamental insights into species dynamics and is essential for developing a null model to accurately study selective processes. However, background selection and selective sweeps can produce genomic signatures at linked sites that mimic or mask signals associated with historical population size change. While the theoretical biases introduced by the linked effects of selection have been well established, it is unclear whether ancestral recombination graph (ARG)-based approaches to demographic inference in typical empirical analyses are susceptible to misinference due to these effects. To address this, we developed highly realistic forward simulations of human and Drosophila melanogaster populations, including empirically estimated variability of gene density, mutation rates, recombination rates, purifying, and positive selection, across different historical demographic scenarios, to broadly assess the impact of selection on demographic inference using a genealogy-based approach. Our results indicate that the linked effects of selection minimally impact demographic inference for human populations, although it could cause misinference in populations with similar genome architecture and population parameters experiencing more frequent recurrent sweeps. We found that accurate demographic inference of D. melanogaster populations by ARG-based methods is compromised by the presence of pervasive background selection alone, leading to spurious inferences of recent population expansion, which may be further worsened by recurrent sweeps, depending on the proportion and strength of beneficial mutations. Caution and additional testing with species-specific simulations are needed when inferring population history with non-human populations using ARG-based approaches to avoid misinference due to the linked effects of selection.

Analyzing the response of rice to tefuryltrione herbicide: Haplotype variation and evolutionary dynamics of the HIS1 gene

Weeds pose multifaceted challenges in rice cultivation, leading to substantial economic losses through reduced yield and poor grain quality. Harnessing the natural genetic diversity in germplasm collections becomes crucial for identifying novel herbicide resistance loci in crops. A comprehensive analysis was conducted on 475 rice accessions from the KRICE depository, assessing their response to TFT (tefuryltrione) and probing the underlying HIS1 (HPPD INHIBITOR SENSITIVE 1) genotypic variations. The HIS1 gene, responsible for detoxifying benzobicyclon (BBC) and imparting broad-spectrum herbicide resistance, holds significant promise in rice breeding. This study explores the genetic landscape of HIS1 within Korean rice collection (KRICE), aiming to unveil genetic variations, haplotype diversity, and evolutionary relationships across diverse rice ecotypes. The indica ecotype showed the highest nucleotide diversity, while the wild and temperate japonica groups exhibited low diversity, hinting at selective sweeps and possible population expansion. Negative Tajima's D values in temperate japonica and wild groups indicate an excess of low-frequency mutations, potentially resulting from selective sweeps. In contrast, with positive Tajima's D values, admixture, indica, and aus groups suggest balancing selection. Furthermore, haplotype analysis uncovered 42 distinct haplotypes within KRICE, with four shared haplotypes between cultivated and wild accessions, four specific to cultivated accessions, and 34 specific to wild types. Phenotypic assessments of these haplotypes revealed that three haplotypes, viz., Hap_1 (predominant in japonica), Hap_2 (predominant in indica), and Hap_3 (specific to indica), displayed significant differences from aus-specific Hap_4 and indica-specific Hap_5. This study offers insights into genetic diversity, selective pressures, and ecotype-specific responses, ultimately paving the way for developing HPPD-inhibiting herbicide-resistant rice cultivars.

Signatures of Adaptation and Purifying Selection in Highland Populations of Dasiphora fruticosa

High mountains harbor a considerable proportion of biodiversity, but we know little about how diverse plants adapt to the harsh environment. Here we finished a high-quality genome assembly for Dasiphora fruticosa, an ecologically important plant distributed in the Qinghai-Tibetan Plateau and lowland of the Northern Hemisphere, and resequenced 592 natural individuals to address how this horticulture plant adapts to highland. Demographic analysis revealed D. fruticosa underwent a bottleneck after Naynayxungla Glaciation. Selective sweep analysis of two pairs of lowland and highland populations identified 63 shared genes related to cell wall organization or biogenesis, cellular component organization, and dwarfism, suggesting parallel adaptation to highland habitats. Most importantly, we found that stronger purging of estimated genetic load due to inbreeding in highland populations apparently contributed to their adaptation to the highest mountain. Our results revealed how plants could tolerate the extreme plateau, which could provide potential insights for species conservation and crop breeding.

Generation of aroma in three-line hybrid rice through CRISPR/Cas9 editing of BETAINE ALDEHYDE DEHYDROGENASE2 (OsBADH2)

Aroma or fragrance in rice is a genetically controlled trait; Its high appreciation by consumers increases the rice market price. Previous studies have revealed that the rice aroma is controlled by a specific gene called BETAINE ALDEHYDE DEHYDROGENASE (OsBADH2), and mutation of this gene leads to the accumulation of an aromatic substance 2-acetyl-1-pyrroline (2-AP). The use of genetic engineering to produce aroma in commercial and cultivated hybrids is a contemporary need for molecular breeding. The current study reports the generation of aroma in the three-line hybrid restorer line Shu-Hui-313 (SH313). We created knock-out (KO) lines of OsBADH2 through the CRISPR/Cas9. The analysis of KO lines revealed a significantly increased content of 2AP in the grains compared with the control. However, other phenotypic traits (plant height, seed setting rate, and 1000-grain weight) were significantly decreased. These KO lines were crossed with a non-aromatic three-line hybrid rice male sterile line (Rong-7-A) to produce Rong-7-You-626 (R7Y626), R7Y627 and R7Y628. The measurement of 2-AP revealed significantly increased contents in these cross combinations. We compared the content of 2-AP in tissues at the booting stage. Data revealed that young spike stalk base contained the highest content of 2-AP and can be used for identification (by simple chewing) of aromatic lines under field conditions. In conclusion, our dataset offers a genetic source and illustrates the generation of aroma in non-aromatic hybrids, and outlines a straightforward identification under field conditions.

Detection of locally adapted genomic regions in wild rice (Oryza rufipogon) using environmental association analysis

Oryza rufipogon is the wild progenitor of cultivated rice Oryza sativa and exhibits high levels of genetic diversity across its distribution, making it a useful resource for the identification of abiotic stress-tolerant varieties and genes that could limit future climate-changed-induced yield losses. To investigate local adaptation in O. rufipogon, we analyzed single nucleotide polymorphism (SNP) data from a panel of 286 samples located across a diverse range of climates. Environmental association analysis (EAA), a genome-wide association study (GWAS)-based method, was used and revealed 15 regions of the genome significantly associated with various climate factors. Genes within these environmentally associated regions have putative functions in abiotic stress response, phytohormone signaling, and the control of flowering time. This provides an insight into potential local adaptation in O. rufipogon and reveals possible locally adaptive genes that may provide opportunities for breeding novel rice varieties with climate change-resilient phenotypes.

Gene expression differentials driven by mass rearing and artificial selection in black soldier fly colonies

The micro-evolutionary forces that shape genetic diversity during domestication have been assessed in many plant and animal systems. However, the impact of these processes on gene expression, and consequent functional adaptation to artificial environments, remains under-investigated. In this study, whole-transcriptome dynamics associated with the early stages of domestication of the black soldier fly (BSF), Hermetia illucens, were assessed. Differential gene expression (DGE) was evaluated in relation to (i) generational time within the cultured environment (F2 vs. F3), and (ii) two selection strategies [no artificial selective pressure (NS); and selection for greater larval mass (SEL)]. RNA-seq was conducted on 5th instar BSF larvae (n = 36), representing equal proportions of the NS (F2 = 9; F3 = 9) and SEL (F2 = 9; F3 = 9) groups. A multidimensional scaling plot revealed greater gene expression variability within the NS and F2 subgroups, while the SEL group clustered separately with lower levels of variation. Comparisons between generations revealed 898 differentially expressed genes (DEGs; FDR-corrected p < 0.05), while between selection strategies, 213 DEGs were observed (FDR-corrected p < 0.05). Enrichment analyses revealed that metabolic, developmental, and defence response processes were over-expressed in the comparison between F2 and F3 larvae, while metabolic processes were the main differentiating factor between NS and SEL lines. This illustrates the functional adaptations that occur in BSF colonies across generations due to mass rearing; as well as highlighting genic dynamics associated with artificial selection for production traits that might inform future selective breeding strategies.

Genome-Wide Analysis of microRNAs and Their Target Genes in Dongxiang Wild Rice (Oryza rufipogon Griff.) Responding to Salt Stress

Rice (Oryza sativa) is a staple food for more than half of the world's population, and its production is critical for global food security. Moreover, rice yield decreases when exposed to abiotic stresses, such as salinity, which is one of the most detrimental factors for rice production. According to recent trends, as global temperatures continue to rise due to climate change, more rice fields may become saltier. Dongxiang wild rice (Oryza rufipogon Griff., DXWR) is a progenitor of cultivated rice and has a high tolerance to salt stress, making it useful for studying the regulatory mechanisms of salt stress tolerance. However, the regulatory mechanism of miRNA-mediated salt stress response in DXWR remains unclear. In this study, miRNA sequencing was performed to identify miRNAs and their putative target genes in response to salt stress in order to better understand the roles of miRNAs in DXWR salt stress tolerance. A total of 874 known and 476 novel miRNAs were identified, and the expression levels of 164 miRNAs were found to be significantly altered under salt stress. The stem-loop quantitative real-time PCR (qRT-PCR) expression levels of randomly selected miRNAs were largely consistent with the miRNA sequencing results, suggesting that the sequencing results were reliable. The gene ontology (GO) analysis indicated that the predicted target genes of salt-responsive miRNAs were involved in diverse biological pathways of stress tolerance. This study contributes to our understanding of DXWR salt tolerance mechanisms regulated by miRNAs and may ultimately improve salt tolerance in cultivated rice breeding using genetic methods in the future.

Discovery and Validation of Grain Shape Loci in U.S. Rice Germplasm Through Haplotype Characterization

Rice grain shape is a major determinant of rice market value and the end-use. We mapped quantitative trait loci (QTL) for grain shape traits in a bi-parental recombinant inbred line population (Trenasse/Jupiter) and discovered two major grain length QTLs—qGL3.1 and qGL7.1. Previously, a major grain shape gene GS3 was reported in the qGL3.1 region and grain length gene GL7 was reported to be encompassing qGL7.1 locus. The re-sequencing SNP data on the International Rice Research Institute (IRRI) 3K Rice Genome Project (RGP) panel were obtained from the IRRI SNP-Seek database for both genes and haplotype diversity was characterized for each gene in this diverse panel. United States rice germplasm was not well represented in the IRRI 3K RGP database. Therefore, a minimum SNP set was identified for each gene that could differentiate all the characterized haplotypes. These haplotypes in the 3K RGP panel were screened across 323 elite U.S. genotypes using the minimum SNP set. The screening of haplotypes and phenotype association confirmed the role of GS3 under qGL3.1. However, screening of the GL7 haplotypes in the U.S. germplasm panel showed that GL7 did not play a role in qGL7.1, and in addition, GL7.1 did not segregate in the Trenasse/Jupiter RIL population. This concluded that qGL7.1 is a novel QTL discovered on chr7 for grain shape in the Trenasse/Jupiter RIL population. A high-throughput KASP-based SNP marker for each locus (GS3 and qGL7.1) was identified and validated in elite U.S. rice germplasm to be used in an applied rice breeding program.

Identification of quantitative trait loci for tillering, root, and shoot biomass at the maximum tillering stage in rice

Tillering and plant biomass are key determinants of rice crop productivity. Tillering at the vegetative stage is associated with weed competition, nutrient uptake, and methane emissions. However, little information is available on quantitative trait loci (QTLs) associated with tiller number (qTN), root biomass (qRB), and shoot biomass (qSB) at the active tillering stage which occurs approximately 6 weeks after planting. Here, we mapped tiller and biomass QTLs with ~ 250 recombinant inbred lines derived from a ‘Francis’ by ‘Rondo’ cross using data collected at the maximum tillering stage from two years of greenhouse study, and further compared these QTLs with those mapped at the harvest stage from a field study. Across these three studies, we discovered six qTNs, two qRBs, and three qSBs. Multiple linear regression further indicated that qTN1-2, qTN3-3, qTN4-1, qRB3-1, and qRB5-1 were significant at the maximum tillering stage while qTN3-2 was detected only at the harvest stage. Moreover, qTN3-1 was consistently significant across different developmental stages and growing environments. The genes identified from the peak target qTN regions included a carotenoid metabolism enzyme, a MYB transcription factor, a CBS domain-containing protein, a SAC3/GANP family protein, a TIFY motif containing protein, and an ABC transporter protein. Two genes in the qRB peak target regions included an expressed protein and a WRKY gene. This knowledge of the QTLs, associated markers, candidate genes, and germplasm resources with high TN, RB and SB is of value to rice cultivar improvement programs.

Global genomic analyses of wheat powdery mildew reveal association of pathogen spread with historical human migration and trade

The fungus Blumeria graminis f. sp. tritici causes wheat powdery mildew disease. Here, we study its spread and evolution by analyzing a global sample of 172 mildew genomes. Our analyses show that B.g. tritici emerged in the Fertile Crescent during wheat domestication. After it spread throughout Eurasia, colonization brought it to America, where it hybridized with unknown grass mildew species. Recent trade brought USA strains to Japan, and European strains to China. In both places, they hybridized with local ancestral strains. Thus, although mildew spreads by wind regionally, our results indicate that humans drove its global spread throughout history and that mildew rapidly evolved through hybridization.

The fungus Blumeria graminis f. sp. tritici causes wheat powdery mildew disease. Here, Sotiropoulos et al. analyze a global sample of 172 mildew genomes, providing evidence that humans drove global spread of the pathogen throughout history and that mildew rapidly evolved through hybridization with local fungal strains.

Draft genome sequence of Indian mulberry (Morus indica) provides a resource for functional and translational genomics

Mulberry is an important crop plant for the sericulture industry. Here, we report high-quality genome sequence of a cultivated Indian mulberry (Morus indica cv K2) obtained by combining data from four different technologies, including Illumina, single-molecule real-time sequencing, chromosome conformation capture and optical mapping, with a gene completeness of 96.5%. Based on the genome sequence, we identified 49.2% of repetitive DNA and 27,435 high-confidence protein-coding genes with >90% of them supported by transcript evidence. A comparative analysis with other plant genomes identified 4.8% of species-specific genes in the M. indica genome. Transcriptome profiling revealed tissue-specific and differential expression across multiple accessions of ~4.7% and 2-5% of protein-coding genes, respectively, implicated in diverse biological processes. Whole genome resequencing of 21 accessions/species revealed ~2.5 million single nucleotide polymorphisms and ~ 0.2 million insertions/deletions. These data and results provide a comprehensive resource to accelerate the genomics research in mulberry for its improvement.

Optimisation of root traits to provide enhanced ecosystem services in agricultural systems: A focus on cover crops

Roots are the interface between the plant and the soil and play a central role in multiple ecosystem processes. With intensification of agricultural practices, rhizosphere processes are being disrupted and are causing degradation of the physical, chemical and biotic properties of soil. However, cover crops, a group of plants that provide ecosystem services, can be utilised during fallow periods or used as an intercrop to restore soil health. The effectiveness of ecosystem services provided by cover crops varies widely as very little breeding has occurred in these species. Improvement of ecosystem service performance is rarely considered as a breeding trait due to the complexities and challenges of belowground evaluation. Advancements in root phenotyping and genetic tools are critical in accelerating ecosystem service improvement in cover crops. In this study, we provide an overview of the range of belowground ecosystem services provided by cover crop roots: (1) soil structural remediation, (2) capture of soil resources and (3) maintenance of the rhizosphere and building of organic matter content. Based on the ecosystem services described, we outline current and promising phenotyping technologies and breeding strategies in cover crops that can enhance agricultural sustainability through improvement of root traits.

Highly diverged lineages of Saccharomyces paradoxus in temperate to subtropical climate zones in China

The wild yeast Saccharomyces paradoxus has become a new model in ecology and evolutionary biology. Different lineages of S. paradoxus have been recognized across the world, but the distribution and genetic diversity of the species remain unknown in China, where the origin of its sibling species S. cerevisiae lies. In this study, we investigated the ecological and geographic distribution of S. paradoxus through an extensive field survey in China and performed population genomic analysis on a set of S. paradoxus strains, including 27 strains, representing different geographic and ecological origins within China, and 59 strains representing all the known lineages of the species recognized in the other regions of the world so far. We found two distinct lineages of S. paradoxus in China. The majority of the Chinese strains studied belong to the Far East lineage, and six strains belong to a novel highly diverged lineage. The distribution of these two lineages overlaps ecologically and geographically in temperate to subtropical climate zones in China. With the addition of the new China lineage, the Eurasian population of S. paradoxus exhibits higher genetic diversity than the American population. We observed more possible lineage-specific introgression events from the Eurasian lineages than from the American lineages. Our results expand the knowledge on ecology, genetic diversity, biogeography, and evolution of S. paradoxus.

Evaluation of gene flow and genetic diversity in rice accessions across Kurdistan region-iraq using SSR markers

Background

In recent years, farmers have complained that the only way to obtain seeds is to select plants that show good performance under local climate conditions in the region. This study aimed to investigate the diversity of rice accessions grown in the region to build a breeding program.

Methods and results

A total of 62 accessions of rice from farmers and research stations were collected from the Kurdistan region, including short-grain and long-grain types, for molecular genetics and diversity analysis. In this study, 37 polymorphic simple sequence repeat (SSR) markers were selected with several molecular genetics software programs. The results show that these SSR markers are very effective for this investigation, generating a total of 152 observed alleles (Na), 75.166 Effective number of alleles (Ne) and an average of 4.1 and 2.03 alleles per locus, respectively. The average polymorphic information content (PIC) per locus was recorded as 0.404. The research presented here confirms two subpopulations, japonica (C1 and C2) and indica (C3), based on molecular genetics data analysis. Analysis of molecular variance revealed that the 72% variance was due to the variation among populations and 28% within the population.

Conclusions

Altogether, these results indicate that there is very low gene flow. These results show the importance of the study of genetic diversity and relationships for starting breeding and improvement programs for rice in the Kurdistan region.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11033-021-06920-x.

Genome-wide association studies of yield-related traits in high-latitude japonica rice

BACKGROUND

Heilongjiang Province is a high-quality japonica rice cultivation area in China. One in ten bowls of Chinese rice is produced here. Increasing yield is one of the main aims of rice production in this area. However, yield is a complex quantitative trait composed of many factors. The purpose of this study was to determine how many genetic loci are associated with yield-related traits. Genome-wide association studies (GWAS) were performed on 450 accessions collected from northeast Asia, including Russia, Korea, Japan and Heilongjiang Province of China. These accessions consist of elite varieties and landraces introduced into Heilongjiang Province decade ago.

RESULTS

After resequencing of the 450 accessions, 189,019 single nucleotide polymorphisms (SNPs) were used for association studies by two different models, a general linear model (GLM) and a mixed linear model (MLM), examining four traits: days to heading (DH), plant height (PH), panicle weight (PW) and tiller number (TI). Over 25 SNPs were found to be associated with each trait. Among them, 22 SNPs were selected to identify candidate genes, and 2, 8, 1 and 11 SNPs were found to be located in 3' UTR region, intron region, coding region and intergenic region, respectively.

CONCLUSIONS

All SNPs detected in this research may become candidates for further fine mapping and may be used in the molecular breeding of high-latitude rice.

Balancing selection and wild gene pool contribute to resistance in global rice germplasm against planthopper

Interactions and co-evolution between plants and herbivorous insects are critically important in agriculture. Brown planthopper (BPH) is the most severe insect of rice, and the biotypes adapt to feed on different rice genotypes. Here, we present genomics analyses on 1,520 global rice germplasms for resistance to three BPH biotypes. Genome-wide association studies identified 3,502 single nucleotide polymorphisms (SNPs) and 59 loci associated with BPH resistance in rice. We cloned a previously unidentified gene Bph37 that confers resistance to BPH. The associated loci showed high nucleotide diversity. Genome-wide scans for trans-species polymorphisms revealed ancient balancing selection at the loci. The secondarily evolved insect biotypes II and III exhibited significantly higher virulence and overcame more rice varieties than the primary biotype I. In response, more SNPs and loci evolved in rice for resistance to biotypes II and III. Notably, three exceptional large regions with high SNP density and resistance-associated loci on chromosomes 4 and 6 appear distinct between the resistant and susceptible rice varieties. Surprisingly, these regions in resistant rice might have been retained from wild species Oryza nivara. Our findings expand the understanding of long-term interactions between rice and BPH and provide resistance genes and germplasm resources for breeding durable BPH-resistant rice varieties.

Population genomics of apricots unravels domestication history and adaptive events

Among crop fruit trees, the apricot (Prunus armeniaca) provides an excellent model to study divergence and adaptation processes. Here, we obtain nearly 600 Armeniaca apricot genomes and four high-quality assemblies anchored on genetic maps. Chinese and European apricots form two differentiated gene pools with high genetic diversity, resulting from independent domestication events from distinct wild Central Asian populations, and with subsequent gene flow. A relatively low proportion of the genome is affected by selection. Different genomic regions show footprints of selection in European and Chinese cultivated apricots, despite convergent phenotypic traits, with predicted functions in both groups involved in the perennial life cycle, fruit quality and disease resistance. Selection footprints appear more abundant in European apricots, with a hotspot on chromosome 4, while admixture is more pervasive in Chinese cultivated apricots. Our study provides clues to the biology of selected traits and targets for fruit tree research and breeding.

Genetic differentiation and restricted gene flow in rice landraces from Yunnan, China: effects of isolation-by-distance and isolation-by-environment

BACKGROUND

Understanding and identifying the factors responsible for genetic differentiation is of fundamental importance for efficient utilization and conservation of traditional rice landraces. In this study, we examined the spatial genetic differentiation of 594 individuals sampled from 28 locations in Yunnan Province, China, covering a wide geographic distribution and diverse growing conditions. All 594 accessions were studied using ten unlinked target genes and 48 microsatellite loci, and the representative 108 accessions from the whole collection were sampled for resequencing.

RESULTS

The genetic diversity of rice landraces was quite different geographically and exhibited a geographical decline from south to north in Yunnan, China. Population structure revealed that the rice landraces could be clearly differentiated into japonica and indica groups, respectively. In each group, the rice accessions could be further differentiated corresponded to their geographic locations, including three subgroups from northern, southern and middle locations. We found more obvious internal geographic structure in the japonica group than in the indica group. In the japonica group, we found that genetic and phenotypic differentiation were strongly related to geographical distance, suggesting a pattern of isolation by distance (IBD); this relationship remained highly significant when we controlled for environmental effects, where the likelihood of gene flow is inversely proportional to the distance between locations. Moreover, the gene flow also followed patterns of isolation by environment (IBE) whereby gene flow rates are higher in similar environments. We detected 314 and 216 regions had been differentially selected between Jap-N and Jap-S, Ind-N and Ind-S, respectively, and thus referred to as selection signatures for different geographic subgroups. We also observed a number of significant and interesting associations between loci and environmental factors, which implies adaptation to local environment.

CONCLUSIONS

Our findings highlight the influence of geographical isolation and environmental heterogeneity on the pattern of the gene flow, and demonstrate that both geographical isolation and environment drives adaptive divergence play dominant roles in the genetic differentiation of the rice landraces in Yunnan, China as a result of limited dispersal.

Genetic loci mediating circadian clock output plasticity and crop productivity under barley domestication

Circadian clock rhythms are shown to be intertwined with crop adaptation. To realize the adaptive value of changes in these rhythms under crop domestication and improvement, there is a need to compare the genetics of clock and yield traits. We compared circadian clock rhythmicity based on Chl leaf fluorescence and transcriptomics among wild ancestors, landraces, and breeding lines of barley under optimal and high temperatures. We conducted a genome scan to identify pleiotropic loci regulating the clock and field phenotypes. We also compared the allelic diversity in wild and cultivated barley to test for selective sweeps. We found significant loss of thermal plasticity in circadian rhythms under domestication. However, transcriptome analysis indicated that this loss was only for output genes and that temperature compensation in the core clock machinery was maintained. Drivers of the circadian clock (DOC) loci were identified via genome-wide association study. Notably, these loci also modified growth and reproductive outputs in the field. Diversity analysis indicated selective sweep in these pleiotropic DOC loci. These results indicate a selection against thermal clock plasticity under barley domestication and improvement and highlight the importance of identifying genes underlying for understanding the biochemical basis of crop adaptation to changing environments.

Inferring the Demographic History of Inbred Species from Genome-Wide SNP Frequency Data

Demographic inference using the site frequency spectrum (SFS) is a common way to understand historical events affecting genetic variation. However, most methods for estimating demography from the SFS assume random mating within populations, precluding these types of analyses in inbred populations. To address this issue, we developed a model for the expected SFS that includes inbreeding by parameterizing individual genotypes using beta-binomial distributions. We then take the convolution of these genotype probabilities to calculate the expected frequency of biallelic variants in the population. Using simulations, we evaluated the model's ability to coestimate demography and inbreeding using one- and two-population models across a range of inbreeding levels. We also applied our method to two empirical examples, American pumas (Puma concolor) and domesticated cabbage (Brassica oleracea var. capitata), inferring models both with and without inbreeding to compare parameter estimates and model fit. Our simulations showed that we are able to accurately coestimate demographic parameters and inbreeding even for highly inbred populations (F = 0.9). In contrast, failing to include inbreeding generally resulted in inaccurate parameter estimates in simulated data and led to poor model fit in our empirical analyses. These results show that inbreeding can have a strong effect on demographic inference, a pattern that was especially noticeable for parameters involving changes in population size. Given the importance of these estimates for informing practices in conservation, agriculture, and elsewhere, our method provides an important advancement for accurately estimating the demographic histories of these species.

Soybean (Glycine max) Haplotype Map (GmHapMap): a universal resource for soybean translational and functional genomics

Here, we describe a worldwide haplotype map for soybean (GmHapMap) constructed using whole-genome sequence data for 1007 Glycine max accessions and yielding 14.9 million variants as well as 4.3 M tag single-nucleotide polymorphisms (SNPs). When sampling random subsets of these accessions, the number of variants and tag SNPs plateaued beyond approximately 800 and 600 accessions, respectively. This suggests extensive coverage of diversity within the cultivated soybean. GmHapMap variants were imputed onto 21 618 previously genotyped accessions with up to 96% success for common alleles. A local association analysis was performed with the imputed data using markers located in a 1-Mb region known to contribute to seed oil content and enabled us to identify a candidate causal SNP residing in the NPC1 gene. We determined gene-centric haplotypes (407 867 GCHs) for the 55 589 genes and showed that such haplotypes can help to identify alleles that differ in the resulting phenotype. Finally, we predicted 18 031 putative loss-of-function (LOF) mutations in 10 662 genes and illustrated how such a resource can be used to explore gene function. The GmHapMap provides a unique worldwide resource for applied soybean genomics and breeding.

Genetics and Genomics of African Rice (Oryza glaberrima Steud) Domestication

African rice (Oryza glaberrima Steud) is one of the two independently domesticated rice species, the other one being Asian rice (Oryza sativa L.). Despite major progress being made in understanding the evolutionary and domestication history of African rice, key outstanding issues remain controversial. There appears to be an underlying difficulty in identifying the domestication centre and number of times the crop has been domesticated. Advances in genomics have provided unprecedented opportunities for understanding the genetic architecture of domestication related traits. For most of the domestication traits, the underlying genes and mutations have been identified. Comparative analysis of domestication genes between Asian and African rice has revealed that the two species went through an independent but convergent evolution process. The genetic and developmental basis of some of the domestic traits are conserved not only between Asian and African rice but also with other domesticated crop species. Analysis of genome data and its interpretation is emerging as a major challenge facing studies of domestication in African rice as key studies continue giving contradictory findings and conclusions. Insights obtained on the domestication of this species are vital for guiding crop improvement efforts.

The Road to Sorghum Domestication: Evidence From Nucleotide Diversity and Gene Expression Patterns

Native African cereals (sorghum, millets) ensure food security to millions of low-income people from low fertility and drought-prone regions of Africa and Asia. In spite of their agronomic importance, the genetic bases of their phenotype and adaptations are still not well-understood. Here we focus on Sorghum bicolor, which is the fifth cereal worldwide for grain production and constitutes the staple food for around 500 million people. We leverage transcriptomic resources to address the adaptive consequences of the domestication process. Gene expression and nucleotide variability were analyzed in 11 domesticated and nine wild accessions. We documented a downregulation of expression and a reduction of diversity both in nucleotide polymorphism (30%) and gene expression levels (18%) in domesticated sorghum. These findings at the genome-wide level support the occurrence of a global reduction of diversity during the domestication process, although several genes also showed patterns consistent with the action of selection. Nine hundred and forty-nine genes were significantly differentially expressed between wild and domesticated gene pools. Their functional annotation points to metabolic pathways most likely contributing to the sorghum domestication syndrome, such as photosynthesis and auxin metabolism. Coexpression network analyzes revealed 21 clusters of genes sharing similar expression patterns. Four clusters (totaling 2,449 genes) were significantly enriched in differentially expressed genes between the wild and domesticated pools and two were also enriched in domestication and improvement genes previously identified in sorghum. These findings reinforce the evidence that the combined and intricated effects of the domestication and improvement processes do not only affect the behaviors of a few genes but led to a large rewiring of the transcriptome. Overall, these analyzes pave the way toward the identification of key domestication genes valuable for genetic resources characterization and breeding purposes.

The History and Diversity of Rice Domestication as Resolved From 1464 Complete Plastid Genomes

The plastid is an essential organelle in autotrophic plant cells, descending from free-living cyanobacteria and acquired by early eukaryotic cells through endosymbiosis roughly one billion years ago. It contained a streamlined genome (plastome) that is uniparentally inherited and non-recombinant, which makes it an ideal tool for resolving the origin and diversity of plant species and populations. In the present study, a large dataset was amassed by de novo assembling plastomes from 295 common wild rice (Oryza rufipogon Griff.) and 1135 Asian cultivated rice (Oryza sativa L.) accessions, supplemented with 34 plastomes from other Oryza species. From this dataset, the phylogenetic relationships and biogeographic history of O. rufipogon and O. sativa were reconstructed. Our results revealed two major maternal lineages across the two species, which further diverged into nine well supported genetic clusters. Among them, the Or-wj-I/II/III and Or-wi-I/II genetic clusters were shared with cultivated (percentage for each cluster ranging 54.9%∼99.3%) and wild rice accessions. Molecular dating, phylogeographic analyses and reconstruction of population historical dynamics indicated an earlier origin of the Or-wj-I/II genetic clusters from East Asian with at least two population expansions, and later origins of other genetic clusters from multiple regions with one or more population expansions. These results supported a single origin of japonica rice (mainly in Or-wj-I/II) and multiple origins of indica rice (in all five clusters) for the history of rice domestication. The massive plastomic data set presented here provides an important resource for understanding the history and evolution of rice domestication as well as a genomic resources for use in future breeding and conservation efforts.

Evolutionary Genomics of Structural Variation in Asian Rice (Oryza sativa) Domestication

Structural variants (SVs) are a largely unstudied feature of plant genome evolution, despite the fact that SVs contribute substantially to phenotypes. In this study, we discovered SVs across a population sample of 347 high-coverage, resequenced genomes of Asian rice (Oryza sativa) and its wild ancestor (O. rufipogon). In addition to this short-read data set, we also inferred SVs from whole-genome assemblies and long-read data. Comparisons among data sets revealed different features of genome variability. For example, genome alignment identified a large (∼4.3 Mb) inversion in indica rice varieties relative to japonica varieties, and long-read analyses suggest that ∼9% of genes from the outgroup (O. longistaminata) are hemizygous. We focused, however, on the resequencing sample to investigate the population genomics of SVs. Clustering analyses with SVs recapitulated the rice cultivar groups that were also inferred from SNPs. However, the site-frequency spectrum of each SV type-which included inversions, duplications, deletions, translocations, and mobile element insertions-was skewed toward lower frequency variants than synonymous SNPs, suggesting that SVs may be predominantly deleterious. Among transposable elements, SINE and mariner insertions were found at especially low frequency. We also used SVs to study domestication by contrasting between rice and O. rufipogon. Cultivated genomes contained ∼25% more derived SVs and mobile element insertions than O. rufipogon, indicating that SVs contribute to the cost of domestication in rice. Peaks of SV divergence were enriched for known domestication genes, but we also detected hundreds of genes gained and lost during domestication, some of which were enriched for traits of agronomic interest.

Genetic Diversity of Landraces and Improved Varieties of Rice (Oryza sativa L.) in Taiwan

BACKGROUND

Rice, the most important crop in Asia, has been cultivated in Taiwan for more than 5000 years. The landraces preserved by indigenous peoples and brought by immigrants from China hundreds of years ago exhibit large variation in morphology, implying that they comprise rich genetic resources. Breeding goals according to the preferences of farmers, consumers and government policies also alter gene pools and genetic diversity of improved varieties. To unveil how genetic diversity is affected by natural, farmers', and breeders' selections is crucial for germplasm conservation and crop improvement.

RESULTS

A diversity panel of 148 rice accessions, including 47 cultivars and 59 landraces from Taiwan and 42 accessions from other countries, were genotyped by using 75 molecular markers that revealed an average of 12.7 alleles per locus with mean polymorphism information content of 0.72. These accessions could be grouped into five subpopulations corresponding to wild rice, japonica landraces, indica landraces, indica cultivars, and japonica cultivars. The genetic diversity within subpopulations was: wild rices > landraces > cultivars; and indica rice > japonica rice. Despite having less variation among cultivars, japonica landraces had greater genetic variation than indica landraces because the majority of Taiwanese japonica landraces preserved by indigenous peoples were classified as tropical japonica. Two major clusters of indica landraces were formed by phylogenetic analysis, in accordance with immigration from two origins. Genetic erosion had occurred in later japonica varieties due to a narrow selection of germplasm being incorporated into breeding programs for premium grain quality. Genetic differentiation between early and late cultivars was significant in japonica (F_ST = 0.3751) but not in indica (F_ST = 0.0045), indicating effects of different breeding goals on modern germplasm. Indigenous landraces with unique intermediate and admixed genetic backgrounds were untapped, representing valuable resources for rice breeding.

CONCLUSIONS

The genetic diversity of improved rice varieties has been substantially shaped by breeding goals, leading to differentiation between indica and japonica cultivars. Taiwanese landraces with different origins possess various and unique genetic backgrounds. Taiwanese rice germplasm provides diverse genetic variation for association mapping to unveil useful genes and is a precious genetic reservoir for rice improvement.

Genomic inference of complex domestication histories in three Solanaceae species

Domestication is a human-induced selection process that imprints the genomes of domesticated populations over a short evolutionary time scale and that occurs in a given demographic context. Reconstructing historical gene flow, effective population size changes and their timing is therefore of fundamental interest to understand how plant demography and human selection jointly shape genomic divergence during domestication. Yet, the comparison under a single statistical framework of independent domestication histories across different crop species has been little evaluated so far. Thus, it is unclear whether domestication leads to convergent demographic changes that similarly affect crop genomes. To address this question, we used existing and new transcriptome data on three crop species of Solanaceae (eggplant, pepper and tomato), together with their close wild relatives. We fitted twelve demographic models of increasing complexity on the unfolded joint allele frequency spectrum for each wild/crop pair, and we found evidence for both shared and species-specific demographic processes between species. A convergent history of domestication with gene flow was inferred for all three species, along with evidence of strong reduction in the effective population size during the cultivation stage of tomato and pepper. The absence of any reduction in size of the crop in eggplant stands out from the classical view of the domestication process; as does the existence of a "protracted period" of management before cultivation. Our results also suggest divergent management strategies of modern cultivars among species as their current demography substantially differs. Finally, the timing of domestication is species-specific and supported by the few historical records available.

The shared genetic architecture of schizophrenia, bipolar disorder and lifespan

Psychiatric disorders such as Schizophrenia (SCZ) and Bipolar Disorder (BD) represent an evolutionary paradox, as they exhibit strong negative effects on fitness, such as decreased fecundity and early mortality, yet they persist at a worldwide prevalence of approximately 1%. Molecular mechanisms affecting lifespan, which may be widely common among complex diseases with fitness effects, can be studied by the integrated analysis of data from genome-wide association studies (GWAS) of human longevity together with any disease of interest. Here, we report the first of such studies, focusing on the genetic overlap-pleiotropy-between two psychiatric disorders with shortened lifespan, SCZ and BD, and human parental lifespan (PLS) as a surrogate of life expectancy. Our results are twofold: first, we demonstrate extensive polygenic overlap between SCZ and PLS and to a lesser extent between BD and PLS. Second, we identified novel loci shared between PLS and SCZ (n = 39), and BD (n = 8). Whereas most of the identified SCZ (66%) and BD (62%) pleiotropic risk alleles were associated with reduced lifespan, we also detected some antagonistic protective alleles associated to shorter lifespans. In fact, top-associated SNPs with SCZ seems to explain longevity variance explained (LVE) better than many other life-threatening diseases, including Type 2 diabetes and most cancers, probably due to a high overlap with smoking-related pathways. Overall, our study provides evidence of a genetic burden driven through premature mortality among people with SCZ, which can have profound implications for understanding, and potentially treating, the mortality gap associated with this psychiatric disorder.

Comparisons of Natural and Cultivated Populations of Corydalis yanhusuo Indicate Divergent Patterns of Genetic and Epigenetic Variation

Epigenetic variation may contribute to traits that are important in domestication, but how patterns of genetic and epigenetic variation differ between cultivated and wild plants remains poorly understood. In particular, we know little about how selection may shape epigenetic variation in natural and cultivated populations. In this study, we investigated 11 natural populations and 6 major cultivated populations using amplified fragment length polymorphism (AFLP) and methylation-sensitive AFLP (MS-AFLP or MSAP) markers to identify patterns of genetic and epigenetic diversity among Corydalis yanhusuo populations. We further explored correlations among genetic, epigenetic, alkaloidal, and climatic factors in natural and cultivated C. yanhusuo. We found support for a single origin for all cultivated populations, from a natural population which was differentiated from the other natural populations. The magnitude of F _ST based on AFLP was significantly correlated with that for MSAP in pairwise comparisons in both natural and cultivated populations, suggesting a relationship between genetic and epigenetic variation in C. yanhusuo. This relationship was further supported by dbRDA (distance-based redundancy analyses) where some of the epigenetic variation could be explained by genetic variation in natural and cultivated populations. Genetic variation was slightly higher in natural than cultivated populations, and exceeded epigenetic variation in both types of populations. However, epigenetic differentiation exceeded that of genetic differentiation among cultivated populations, while the reverse was observed among natural populations. The differences between wild and cultivated plants may be partly due to processes inherent to cultivation and in particular the differences in mode of reproduction. The importance of epigenetic compared to genetic modifications is thought to vary depending on reproductive strategies, and C. yanhusuo usually reproduces sexually in natural environments, while the cultivated C. yanhusuo are propagated clonally. In addition, alkaloid content of C. yanhusuo varied across cultivated populations, and alkaloid content was significantly correlated to climatic variation, but also to genetic (6.89%) and even more so to epigenetic (14.09%) variation in cultivated populations. Our study demonstrates that epigenetic variation could be important in cultivation of C. yanhusuo and serve as a source of variation for response to environmental conditions.

Characterization of the transcriptional divergence between the subspecies of cultivated rice (Oryza sativa)

Background

Cultivated rice consists of two subspecies, Indica and Japonica, that exhibit well-characterized differences at the morphological and genetic levels. However, the differences between these subspecies at the transcriptome level remains largely unexamined. Here, we provide a comprehensive characterization of transcriptome divergence and cis-regulatory variation within rice using transcriptome data from 91 accessions from a rice diversity panel (RDP1).

Results

The transcriptomes of the two subspecies of rice are highly divergent. Japonica have significantly lower expression and genetic diversity relative to Indica, which is likely a consequence of a population bottleneck during Japonica domestication. We leveraged high-density genotypic data and transcript levels to identify cis-regulatory variants that may explain the genetic divergence between the subspecies. We identified significantly more eQTL that were specific to the Indica subspecies compared to Japonica, suggesting that the observed differences in expression and genetic variability also extends to cis-regulatory variation.

Conclusions

Using RNA sequencing data for 91diverse rice accessions and high-density genotypic data, we show that the two species are highly divergent with respect to gene expression levels, as well as the genetic regulation of expression. The data generated by this study provide, to date, the largest collection of genome-wide transcriptional levels for rice, and provides a community resource to accelerate functional genomic studies in rice.

Advances in genome-wide association studies of complex traits in rice

Genome-wide association studies (GWAS), genetic surveys of the whole genome to detect variants associated with a trait in natural populations, are a powerful approach for dissecting complex traits. This genetic mapping approach has been applied in rice over the last 10 years. During the last decade, GWAS was used to identify the loci underlying tens of rice traits, and several important genes were detected in GWAS and further confirmed in follow-up functional experiments. In this review, we present an overview of the whole process in a typical GWAS, including population design, genotyping, phenotyping and analysis methods. Recent advances in rice GWAS are also provided, including several examples of the functional characterization of candidate genes. The possible breakthroughs of rice GWAS in the next decade are discussed with regard to their application in breeding, the consideration of epistatic interactions and in-depth functional annotations of DNA elements and genetic variants throughout the rice genome.

Natural Sequence Variations and Combinations of GNP1 and NAL1 Determine the Grain Number per Panicle in Rice

BACKGROUND

The grain number per panicle (GNP), which is one of three grain yield components, is an important trait for the genetic improvement of rice. Although the NAL1 and GNP1 genes regulating the rice GNP and grain yield have been cloned, their allelic diversity, functional differences in rice germplasms, and effects of their combination on GNP and grain yield remain unclear.

RESULTS

Based on DNA sequences of these two genes in 198 cultivated rice (Oryza sativa) and 8-10 wild rice (Oryza rufipogon) germplasms, 16 and 14 haplotypes were identified for NAL1 and GNP1, respectively. The NAL1 gene had the strongest effects on GNP in indica (xian) and japonica (geng) subpopulations. In contrast, GNP1 had no significant effects in the geng subpopulation and was rare in the xian background, in which the superior GNP1 allele (GNP1-6) was detected in only 4.0% of the 198 germplasms. Compared with the transgenic lines with GNP1 or NAL1, the transgenic lines with both genes had a higher GNP (15.5%-25.4% and 11.6%-15.9% higher, respectively) and grain yield (5.7%-9.0% and 8.3%-12.3% higher, respectively) across 3 years. The two genes combined in the introgression lines in Lemont background resulted in especially favorable effects on the GNP.

CONCLUSIONS

Our results indicated that the GNP1 and NAL1 exhibited obvious differentiation and their combinations can significantly increase the grain yield in geng rice cultivars. These observations provide insights into the molecular basis of the GNP and may be useful for rice breeding of high yield potential by pyramiding GNP1 and NAL1.

The evolution of metabolism: How to test evolutionary hypotheses at the genomic level

The origin of primordial metabolism and its expansion to form the metabolic networks extant today represent excellent systems to study the impact of natural selection and the potential adaptive role of novel compounds. Here we present the current hypotheses made on the origin of life and ancestral metabolism and present the theories and mechanisms by which the large chemical diversity of plants might have emerged along evolution. In particular, we provide a survey of statistical methods that can be used to detect signatures of selection at the gene and population level, and discuss potential and limits of these methods for investigating patterns of molecular adaptation in plant metabolism.

The strength and pattern of natural selection on gene expression in rice

Levels of gene expression underpin organismal phenotypes^1,2, but the nature of selection that acts on gene expression and its role in adaptive evolution remain unknown^1,2. Here we assayed gene expression in rice (Oryza sativa)³, and used phenotypic selection analysis to estimate the type and strength of selection on the levels of more than 15,000 transcripts^4,5. Variation in most transcripts appears (nearly) neutral or under very weak stabilizing selection in wet paddy conditions (with median standardized selection differentials near zero), but selection is stronger under drought conditions. Overall, more transcripts are conditionally neutral (2.83%) than are antagonistically pleiotropic⁶ (0.04%), and transcripts that display lower levels of expression and stochastic noise^7-9 and higher levels of plasticity⁹ are under stronger selection. Selection strength was further weakly negatively associated with levels of cis-regulation and network connectivity⁹. Our multivariate analysis suggests that selection acts on the expression of photosynthesis genes^4,5, but that the efficacy of selection is genetically constrained under drought conditions¹⁰. Drought selected for earlier flowering^11,12 and a higher expression of OsMADS18 (Os07g0605200), which encodes a MADS-box transcription factor and is a known regulator of early flowering¹³-marking this gene as a drought-escape gene^11,12. The ability to estimate selection strengths provides insights into how selection can shape molecular traits at the core of gene action.

An inferred fitness consequence map of the rice genome

The extent to which sequence variation impacts plant fitness is poorly understood. High-resolution maps detailing the constraint acting on the genome, especially in regulatory sites, would be beneficial as functional annotation of noncoding sequences remains sparse. Here, we present a fitness consequence (fitCons) map for rice (Oryza sativa). We inferred fitCons scores (ρ) for 246 inferred genome classes derived from nine functional genomic and epigenomic datasets, including chromatin accessibility, messenger RNA/small RNA transcription, DNA methylation, histone modifications and engaged RNA polymerase activity. These were integrated with genome-wide polymorphism and divergence data from 1,477 rice accessions and 11 reference genome sequences in the Oryzeae. We found ρ to be multimodal, with ~9% of the rice genome falling into classes where more than half of the bases would probably have a fitness consequence if mutated. Around 2% of the rice genome showed evidence of weak negative selection, frequently at candidate regulatory sites, including a novel set of 1,000 potentially active enhancer elements. This fitCons map provides perspective on the evolutionary forces associated with genome diversity, aids in genome annotation and can guide crop breeding programs.

Eggplant Domestication: Pervasive Gene Flow, Feralization, and Transcriptomic Divergence

In the context of food security, examining the genomics of domestication will help identify genes underlying adaptive and economically important phenotypes, for example, larger fruit, improved taste, and loss of agronomically inferior phenotypes.  Examination of genome-scale single nucleotide polymorphisms demonstrates the relationships between wild ancestors of eggplant (Solanum melongena L.), confirming that Solanum insanum L. is the wild progenitor. This species is split roughly into an Eastern (Malaysian, Thai, and Vietnamese) and Western (Indian, Madagascan, and Sri Lankan) group, with domesticates derived from the former. Additional "wild" accessions from India appear to be feral escapes, derived multiple times from domesticated varieties through admixture. Accessions with small egg-shaped fruit are generally found intermixed with East Asian Solanum insanum confirming they are primitive relative to the large-fruited domesticates.  Comparative transcriptomics was used to track the loci under selection. Sequence analysis revealed a genetic bottleneck reducing variation by almost 50% in the primitive accessions relative to the wild species and a further 10% in the landraces. We also show evidence for selection on genes with a role in response to wounding and apoptosis.  Genes showing a significant difference in expression between wild and primitive or between primitive and landrace genepools were mostly (>75%) downregulated in the derived populations and enriched for gene ontologies related to defense, flowering, signaling, and response to biotic and abiotic stimuli.  This work reveals genomic changes involved in crop domestication and improvement, and the population genetics work explains why defining the eggplant domestication trajectory has been so challenging.

Intraspecific diversification of the crop wild relative Brassica cretica Lam. using demographic model selection

BACKGROUND

Crop wild relatives (CWRs) contain genetic diversity, representing an invaluable resource for crop improvement. Many of their traits have the potential to help crops to adapt to changing conditions that they experience due to climate change. An impressive global effort for the conservation of various CWR will facilitate their use in crop breeding for food security. The genus Brassica is listed in Annex I of the International Treaty on Plant Genetic Resources for Food and Agriculture. Brassica oleracea (or wild cabbage), a species native to southern and western Europe, has become established as an important human food crop plant because of its large reserves stored over the winter in its leaves. Brassica cretica Lam. (Bc) is a CWR in the brassica group and B. cretica subsp. nivea (Bcn) has been suggested as a separate subspecies. The species Bc has been proposed as a potential gene donor to brassica crops, including broccoli, cabbage, cauliflower, oilseed rape, etc. RESULTS: We sequenced genomes of four Bc individuals, including two Bcn and two Bc. Demographic analysis based on our whole-genome sequence data suggests that populations of Bc are not isolated. Classification of the Bc into distinct subspecies is not supported by the data. Using only the non-coding part of the data (thus, the parts of the genome that has evolved nearly neutrally), we find the gene flow between different Bc population is recent and its genomic diversity is high.

CONCLUSIONS

Despite predictions on the disruptive effect of gene flow in adaptation, when selection is not strong enough to prevent the loss of locally adapted alleles, studies show that gene flow can promote adaptation, that local adaptations can be maintained despite high gene flow, and that genetic architecture plays a fundamental role in the origin and maintenance of local adaptation with gene flow. Thus, in the genomic era it is important to link the selected demographic models with the underlying processes of genomic variation because, if this variation is largely selectively neutral, we cannot assume that a diverse population of crop wild relatives will necessarily exhibit the wide-ranging adaptive diversity required for further crop improvement.

Detecting Positive Selection in Populations Using Genetic Data

High-throughput genomic sequencing allows to disentangle the evolutionary forces acting in populations. Among evolutionary forces, positive selection has received a lot of attention because it is related to the adaptation of populations in their environments, both biotic and abiotic. Positive selection, also known as Darwinian selection, occurs when an allele is favored by natural selection. The frequency of the favored allele increases in the population and, due to genetic hitchhiking, neighboring linked variation diminishes, creating so-called selective sweeps. Such a process leaves traces in genomes that can be detected in a future time point. Detecting traces of positive selection in genomes is achieved by searching for signatures introduced by selective sweeps, such as regions of reduced variation, a specific shift of the site frequency spectrum, and particular linkage disequilibrium (LD) patterns in the region. A variety of approaches can be used for detecting selective sweeps, ranging from simple implementations that compute summary statistics to more advanced statistical approaches, e.g., Bayesian approaches, maximum-likelihood-based methods, and machine learning methods. In this chapter, we discuss selective sweep detection methodologies on the basis of their capacity to analyze whole genomes or just subgenomic regions, and on the specific polymorphism patterns they exploit as selective sweep signatures. We also summarize the results of comparisons among five open-source software releases (SweeD, SweepFinder, SweepFinder2, OmegaPlus, and RAiSD) regarding sensitivity, specificity, and execution times. Furthermore, we test and discuss machine learning methods and present a thorough performance analysis. In equilibrium neutral models or mild bottlenecks, most methods are able to detect selective sweeps accurately. Methods and tools that rely on linkage disequilibrium (LD) rather than single SNPs exhibit higher true positive rates than the site frequency spectrum (SFS)-based methods under the model of a single sweep or recurrent hitchhiking. However, their false positive rate is elevated when a misspecified demographic model is used to build the distribution of the statistic under the null hypothesis. Both LD and SFS-based approaches suffer from decreased accuracy on localizing the true target of selection in bottleneck scenarios. Furthermore, we present an extensive analysis of the effects of gene flow on selective sweep detection, a problem that has been understudied in selective sweep literature.

Relationships between Iraqi Rice Varieties at the Nuclear and Plastid Genome Levels

Due to the importance of the rice crop in Iraq, this study was conducted to determine the origin of the major varieties and understand the evolutionary relationships between Iraqi rice varieties and other Asian rice accessions that could be significant in the improvement of this crop. Five varieties of Oryza sativa were obtained from Baghdad/Iraq, and the whole genomic DNA was sequenced, among these varieties, Amber33, Furat, Yasmin, Buhooth1 and Amber al-Baraka. Raw sequence reads of 33 domesticated Asian rice accessions were obtained from the Sequence Read Archive (SRA-NCBI). The sequence of the whole chloroplast-genome was assembled while only the sequence of 916 concatenated nuclear-genes was assembled. The phylogenetic analysis of both chloroplast and nuclear genomes showed that two main clusters, Indica and Japonica, and further five sub-clusters based upon their ecotype, indica, aus, tropical-japonica, temperate-japonica and basmati were created; moreover, Amber33, Furat, Yasmin and Buhooth1 belonged to the basmati, indica and japonica ecotypes, respectively, where Amber33 was placed in the basmati group as a sister of cultivars from Pakistan and India. This confirms the traditional story that Amber was transferred by a group of people who had migrated from India and settled in southern Iraq a long time ago.

The rice genome revolution: from an ancient grain to Green Super Rice

Rice is a staple crop for half the world's population, which is expected to grow by 3 billion over the next 30 years. It is also a key model for studying the genomics of agroecosystems. This dual role places rice at the centre of an enormous challenge facing agriculture: how to leverage genomics to produce enough food to feed an expanding global population. Scientists worldwide are investigating the genetic variation among domesticated rice species and their wild relatives with the aim of identifying loci that can be exploited to breed a new generation of sustainable crops known as Green Super Rice.

Genomic signatures of adaptation to Sahelian and Soudanian climates in sorghum landraces of Senegal

Uncovering the genomic basis of climate adaptation in traditional crop varieties can provide insight into plant evolution and facilitate breeding for climate resilience. In the African cereal sorghum (Sorghum bicolor L. [Moench]), the genomic basis of adaptation to the semiarid Sahelian zone versus the subhumid Soudanian zone is largely unknown. To address this issue, we characterized a large panel of 421 georeferenced sorghum landrace accessions from Senegal and adjacent locations at 213,916 single-nucleotide polymorphisms (SNPs) using genotyping-by-sequencing. Seven subpopulations distributed along the north-south precipitation gradient were identified. Redundancy analysis found that climate variables explained up to 8% of SNP variation, with climate collinear with space explaining most of this variation (6%). Genome scans of nucleotide diversity suggest positive selection on chromosome 2, 4, 5, 7, and 10 in durra sorghums, with successive adaptation during diffusion along the Sahel. Putative selective sweeps were identified, several of which colocalize with stay-green drought tolerance (Stg) loci, and a priori candidate genes for photoperiodic flowering and inflorescence morphology. Genome-wide association studies of photoperiod sensitivity and panicle compactness identified 35 and 13 associations that colocalize with a priori candidate genes, respectively. Climate-associated SNPs colocalize with Stg3a, Stg1, Stg2, and Ma6 and have allelic distribution consistent with adaptation across Sahelian and Soudanian zones. Taken together, the findings suggest an oligogenic basis of adaptation to Sahelian versus Soudanian climates, underpinned by variation in conserved floral regulatory pathways and other systems that are less understood in cereals.

Crop Biodiversity: An Unfinished Magnum Opus of Nature

Crop biodiversity is one of the major inventions of humanity through the process of domestication. It is also an essential resource for crop improvement to adapt agriculture to ever-changing conditions like global climate change and consumer preferences. Domestication and the subsequent evolution under cultivation have profoundly shaped the genetic architecture of this biodiversity. In this review, we highlight recent advances in our understanding of crop biodiversity. Topics include the reduction of genetic diversity during domestication and counteracting factors, a discussion of the relationship between parallel phenotypic and genotypic evolution, the role of plasticity in genotype × environment interactions, and the important role subsistence farmers play in actively maintaining crop biodiversity and in participatory breeding. Linking genotype and phenotype remains the holy grail of crop biodiversity studies.

The Genomics of Oryza Species Provides Insights into Rice Domestication and Heterosis

Here, we review recent progress in genetic and genomic studies of the diversity of Oryza species. In recent years, unlocking the genetic diversity of Oryza species has provided insights into the genomics of rice domestication, heterosis, and complex traits. Genome sequencing and analysis of numerous wild rice (Oryza rufipogon) and Asian cultivated rice (Oryza sativa) accessions have enabled the identification of genome-wide signatures of rice domestication and the unlocking of the origin of Asian cultivated rice. Moreover, similar studies on genome variations of African rice (Oryza glaberrima) cultivars and their closely related wild progenitor Oryza barthii accessions have provided strong evidence to support a theory of independent domestication in African rice. Integrated genomic approaches have efficiently investigated many heterotic loci in hybrid rice underlying yield heterosis advantages and revealed the genomic architecture of rice heterosis. We conclude that in-depth unlocking of genetic variations among Oryza species will further enhance rice breeding.

Genomic analyses of an extensive collection of wild and cultivated accessions provide new insights into peach breeding history

BACKGROUND

Human selection has a long history of transforming crop genomes. Peach (Prunus persica) has undergone more than 5000 years of domestication that led to remarkable changes in a series of agronomically important traits, but genetic bases underlying these changes and the effects of artificial selection on genomic diversity are not well understood.

RESULTS

Here, we report a comprehensive analysis of peach evolution based on genome sequences of 480 wild and cultivated accessions. By focusing on a set of quantitative trait loci (QTLs), we provide evidence supporting that distinct phases of domestication and improvement have led to an increase in fruit size and taste and extended its geographic distribution. Fruit size was predominantly selected during domestication, and selection for large fruits has led to the loss of genetic diversity in several fruit weight QTLs. In contrast, fruit taste-related QTLs were successively selected for by domestication and improvement, with more QTLs selected for during improvement. Genome-wide association studies of 11 agronomic traits suggest a set of candidate genes controlling these traits and potential markers for molecular breeding. Candidate loci for genes that contributed to the adaption to low-chill regions were identified. Furthermore, the genomic bases of divergent selection for fruit texture and local breeding for different flavors between Asian and European/North American cultivars were also determined.

CONCLUSIONS

Our results elucidate the genetic basis of peach evolution and provide new resources for future genomics-guided peach breeding.

Evolutionary Epigenomics of Retrotransposon-Mediated Methylation Spreading in Rice

Plant genomes contain numerous transposable elements (TEs), and many hypotheses on the evolutionary drivers that restrict TE activity have been postulated. Few models, however, have focused on the evolutionary epigenomic interaction between the plant host and its TE. The host genome recruits epigenetic factors, such as methylation, to silence TEs but methylation can spread beyond the TE sequence and influence the expression of nearby host genes. In this study, we investigated this epigenetic trade-off between TE and proximal host gene silencing by studying the epigenomic regulation of repressing long terminal repeat (LTR) retrotransposons (RTs) in Oryza sativa. Results showed significant evidence of methylation spreading originating from the LTR-RT sequences, and the extent of spreading was dependent on five factors: 1) LTR-RT family, 2) time since the LTR-RT insertion, 3) recombination rate of the LTR-RT region, 4) level of LTR-RT sequence methylation, and 5) chromosomal location. Methylation spreading had negative effects by reducing host gene expression, but only on host genes with LTR-RT inserted in its introns. Our results also suggested high levels of LTR-RT methylation might have a role in suppressing TE-mediated deleterious ectopic recombination. In the end, despite the methylation spreading, no strong epigenetic trade-off was detected and majority of LTR-RT may have only minor epigenetic effects on nearby host genes.

Inferring the distribution of selective effects from a time inhomogeneous model

We have developed a Poisson random field model for estimating the distribution of selective effects of newly arisen nonsynonymous mutations that could be observed as polymorphism or divergence in samples of two related species under the assumption that the two species populations are not at mutation-selection-drift equilibrium. The model is applied to 91Drosophila genes by comparing levels of polymorphism in an African population of D. melanogaster with divergence to a reference strain of D. simulans. Based on the difference of gene expression level between testes and ovaries, the 91 genes were classified as 33 male-biased, 28 female-biased, and 30 sex-unbiased genes. Under a Bayesian framework, Markov chain Monte Carlo simulations are implemented to the model in which the distribution of selective effects is assumed to be Gaussian with a mean that may differ from one gene to the other to sample key parameters. Based on our estimates, the majority of newly-arisen nonsynonymous mutations that could contribute to polymorphism or divergence in Drosophila species are mildly deleterious with a mean scaled selection coefficient of -2.81, while almost 86% of the fixed differences between species are driven by positive selection. There are only 16.6% of the nonsynonymous mutations observed in sex-unbiased genes that are under positive selection in comparison to 30% of male-biased and 46% of female-biased genes that are beneficial. We also estimated that D. melanogaster and D. simulans may have diverged 1.72 million years ago.

Decrease of gene expression diversity during domestication of animals and plants

BACKGROUND

The genetic mechanisms underlying the domestication of animals and plants have been of great interest to biologists since Darwin. To date, little is known about the global pattern of gene expression changes during domestication.

RESULTS

We generated and collected transcriptome data for seven pairs of domestic animals and plants including dog, silkworm, chicken, rice, cotton, soybean and maize and their wild progenitors and compared the expression profiles between the domestic and wild species. Intriguingly, although the number of expressed genes varied little, the domestic species generally exhibited lower gene expression diversity than did the wild species, and this lower diversity was observed for both domestic plants and different kinds of domestic animals including insect, bird and mammal in the whole-genome gene set (WGGS), candidate selected gene set (CSGS) and non-CSGS, with CSGS exhibiting a higher degree of decreased expression diversity. Moreover, different from previous reports which found 2 to 4% of genes were selected by human, we identified 6892 candidate selected genes accounting for 7.57% of the whole-genome genes in rice and revealed that fewer than 8% of the whole-genome genes had been affected by domestication.

CONCLUSIONS

Our results showed that domestication affected the pattern of variation in gene expression throughout the genome and generally decreased the expression diversity across species, and this decrease may have been associated with decreased genetic diversity. This pattern might have profound effects on the phenotypic and physiological changes of domestic animals and plants and provide insights into the genetic mechanisms at the transcriptome level other than decreased genetic diversity and increased linkage disequilibrium underpinning artificial selection.