Bridging old and new: diversity and evaluation of high iron-associated stress response of rice cultivated in West Africa.

Adoption of rice varieties that perform well under high iron-associated (HIA) stress environments can enhance rice production in West Africa. This study reports the genetic characterization of 323 rice accessions and breeding lines cultivated in West Africa using genotyping-by-sequencing and their phenotypic response to HIA treatments in hydroponic solution (1500 mg l-1 FeSO4·7H2O) and hot-spot fields. The germplasm consisted of four genetic subpopulations: Oryza glaberrima (14%), O. sativa-japonica (7%), O. sativa-indica Group 1 (45%), and O. sativa-indica Group 2 (25%). Severe versus mild stress in the field was associated with a reduced SPAD value (12%), biomass (56%), and grain yield (57%), with leaf bronzing explaining 30% and 21% of the variation for biomass and grain yield, respectively. Association mapping using 175 indica genotypes identified 23 significant single nucleotide polymorphism (SNP) markers that mapped to 14 genomic regions. Genome-wide association study (GWAS) signals associated with leaf bronzing, a routinely used indicator of HIA stress, differed in hydroponic compared with field conditions. Contrastingly, six significant SNPs on chromosomes 8 and 9 were associated with the SPAD value under HIA stress in both field and hydroponic experiments, and a candidate potassium transporter gene mapped under the peak on chromosome 8. This study helps define criteria for assessing rice performance under HIA environments.

Prevalence of cryptic species in morphologically uniform taxa - Fast speciation and evolutionary radiation in Asian frogs.

Diversity and distributions of cryptic species have long been a vexing issue. Identification of species boundaries is made difficult by the lack of obvious morphological differences. Here, we investigate the cryptic diversity and evolutionary history of an underappreciated group of Asian frog species (Megophrys) to explore the pattern and dynamic of amphibian cryptic species. We sequenced four mitochondrial genes and five nuclear genes and delineated species using multiple approaches, combining DNA and mating-call data. A Bayesian species tree was generated to estimate divergence times and to reconstruct ancestral ranges. Macroevolutionary analyses and hybridization tests were conducted to explore the evolutionary dynamics of this cryptic group. Our phylogenies support the current subgenera. We revealed 43 cryptic species, 158% higher than previously thought. The species-delimitation results were further confirmed by mating-call data and morphological divergence. We found that these Asian frogs entered China from the Sunda Shelf 48 Mya, followed by an ancient radiation event during middle Miocene. We confirmed the efficiency of the multispecies coalescent model for delimitation of species with low morphological diversity. Species diversity of Megophrys is severely underappreciated, and species distributions have been misestimated as a result.

Extremely low genetic diversity across mangrove taxa reflects past sea level changes and hints at poor future responses.

The projected increases in sea levels are expected to affect coastal ecosystems. Tropical communities, anchored by mangrove trees and having experienced frequent past sea level changes, appear to be vibrant at present. However, any optimism about the resilience of these ecosystems is premature because the impact of past climate events may not be reflected in the current abundance. To assess the impact of historical sea level changes, we conducted an extensive genetic diversity survey on the Indo-Malayan coast, a hotspot with a large global mangrove distribution. A survey of 26 populations in six species reveals extremely low genome-wide nucleotide diversity and hence very small effective population sizes (Ne ) in all populations. Whole-genome sequencing of three mangrove species further shows the decline in Ne to be strongly associated with the speed of past changes in sea level. We also used a recent series of flooding events in Yalong Bay, southern China, to test the robustness of mangroves to sea level changes in relation to their genetic diversity. The events resulted in the death of half of the mangrove trees in this area. Significantly, less genetically diverse mangrove species suffered much greater destruction. The dieback was accompanied by a drastic reduction in local invertebrate biodiversity. We thus predict that tropical coastal communities will be seriously endangered as the global sea level rises. Well-planned coastal development near mangrove forests will be essential to avert this crisis.

The emergence of the hyperinvasive vine, Mikania micrantha (Asteraceae), via admixture and founder events inferred from population transcriptomics.

Biological invasions that involve well-documented rapid adaptations to new environments provide unequalled opportunities for testing evolutionary hypotheses. Mikania micrantha Kunth (Asteraceae), a perennial herbaceous vine native to tropical Central and South America, successfully invaded tropical Asia in the early 20th century. It is regarded as one of the most aggressive weeds in the world. To elucidate the molecular and evolutionary processes underlying this invasion, we extensively sampled this weed throughout its invaded range in South-East and South Asia and surveyed its genetic structure using variants detected from population transcriptomics. Clustering results suggest that more than one source population contributed to this invasion. Computer simulations using genomewide genetic variation support a scenario of admixture and founder events during invasion. The genes differentially expressed between native and invasive populations were found to be involved in oxidative and high light intensity stress responses, pointing to a possible ecological mechanism of adaptation. Our results provide a foundation for further detailed mechanistic and population studies of this ecologically and economically important invasion. This line of research promises to provide new mitigation strategies for invasive species as well as insights into mechanisms of adaptation.

Robust phenotyping strategies for evaluation of stem non-structural carbohydrates (NSC) in rice.

Rice plants (Oryza sativa) accumulate excess photoassimilates in the form of non-structural carbohydrates (NSCs) in their stems prior to heading that can later be mobilized to supplement photosynthate production during grain-filling. Despite longstanding interest in stem NSC for rice improvement, the dynamics of NSC accumulation, remobilization, and re-accumulation that have genetic potential for optimization have not been systematically investigated. Here we conducted three pilot experiments to lay the groundwork for large-scale diversity studies on rice stem NSC. We assessed the relationship of stem NSC components with 21 agronomic traits in large-scale, tropical yield trials using 33 breeder-nominated lines, established an appropriate experimental design for future genetic studies using a Bayesian framework to sample sub-datasets from highly replicated greenhouse data using 36 genetically diverse genotypes, and used 434 phenotypically divergent rice stem samples to develop two partial least-squares (PLS) models using near-infrared (NIR) spectra for accurate, rapid prediction of rice stem starch, sucrose, and total non-structural carbohydrates. We find evidence that stem reserves are most critical for short-duration varieties and suggest that pre-heading stem NSC is worthy of further experimentation for breeding early maturing rice.

Age-dependent transition from cell-level to population-level control in murine intestinal homeostasis revealed by coalescence analysis.

In multi-cellular organisms, tissue homeostasis is maintained by an exquisite balance between stem cell proliferation and differentiation. This equilibrium can be achieved either at the single cell level (a.k.a. cell asymmetry), where stem cells follow strict asymmetric divisions, or the population level (a.k.a. population asymmetry), where gains and losses in individual stem cell lineages are randomly distributed, but the net effect is homeostasis. In the mature mouse intestinal crypt, previous evidence has revealed a pattern of population asymmetry through predominantly symmetric divisions of stem cells. In this work, using population genetic theory together with previously published crypt single-cell data obtained at different mouse life stages, we reveal a strikingly dynamic pattern of stem cell homeostatic control. We find that single-cell asymmetric divisions are gradually replaced by stochastic population-level asymmetry as the mouse matures to adulthood. This lifelong process has important developmental and evolutionary implications in understanding how adult tissues maintain their homeostasis integrating the trade-off between intrinsic and extrinsic regulations.

High-Resolution Inflorescence Phenotyping Using a Novel Image-Analysis Pipeline, PANorama.

Variation in inflorescence development is an important target of selection for numerous crop species, including many members of the Poaceae (grasses). In Asian rice (Oryza sativa), inflorescence (panicle) architecture is correlated with yield and grain-quality traits. However, many rice breeders continue to use composite phenotypes in selection pipelines, because measuring complex, branched panicles requires a significant investment of resources. We developed an open-source phenotyping platform, PANorama, which measures multiple architectural and branching phenotypes from images simultaneously. PANorama automatically extracts skeletons from images, allows users to subdivide axes into individual internodes, and thresholds away structures, such as awns, that normally interfere with accurate panicle phenotyping. PANorama represents an improvement in both efficiency and accuracy over existing panicle imaging platforms, and flexible implementation makes PANorama capable of measuring a range of organs from other plant species. Using high-resolution phenotypes, a mapping population of recombinant inbred lines, and a dense single-nucleotide polymorphism data set, we identify, to our knowledge, the largest number of quantitative trait loci (QTLs) for panicle traits ever reported in a single study. Several areas of the genome show pleiotropic clusters of panicle QTLs, including a region near the rice Green Revolution gene SEMIDWARF1. We also confirm that multiple panicle phenotypes are distinctly different among a small collection of diverse rice varieties. Taken together, these results suggest that clusters of small-effect QTLs may be responsible for varietal or subpopulation-specific panicle traits, representing a significant opportunity for rice breeders selecting for yield performance across different genetic backgrounds.

Two evolutionary histories in the genome of rice: the roles of domestication genes.

Genealogical patterns in different genomic regions may be different due to the joint influence of gene flow and selection. The existence of two subspecies of cultivated rice provides a unique opportunity for analyzing these effects during domestication. We chose 66 accessions from the three rice taxa (about 22 each from Oryza sativa indica, O. sativa japonica, and O. rufipogon) for whole-genome sequencing. In the search for the signature of selection, we focus on low diversity regions (LDRs) shared by both cultivars. We found that the genealogical histories of these overlapping LDRs are distinct from the genomic background. While indica and japonica genomes generally appear to be of independent origin, many overlapping LDRs may have originated only once, as a result of selection and subsequent introgression. Interestingly, many such LDRs contain only one candidate gene of rice domestication, and several known domestication genes have indeed been "rediscovered" by this approach. In summary, we identified 13 additional candidate genes of domestication.

Genome-wide misexpression of X-linked versus autosomal genes associated with hybrid male sterility.

Postmating reproductive isolation is often manifested as hybrid male sterility, for which X-linked genes are overrepresented (the so-called large X effect). In contrast, X-linked genes are significantly under-represented among testis-expressing genes. This seeming contradiction may be germane to the X:autosome imbalance hypothesis on hybrid sterility, in which the X-linked effect is mediated mainly through the misexpression of autosomal genes. In this study, we compared gene expression in fertile and sterile males in the hybrids between two Drosophila species. These hybrid males differ only in a small region of the X chromosome containing the Ods-site homeobox (OdsH) (also known as Odysseus) locus of hybrid sterility. Of genes expressed in the testis, autosomal genes were, indeed, more likely to be misexpressed than X-linked genes under the sterilizing action of OdsH. Since this mechanism of X:autosome interaction is only associated with spermatogenesis, a connection between X:autosome imbalance and the high rate of hybrid male sterility seems plausible.

Odorant Binding Causes Cytoskeletal Rearrangement, Leading to Detectable Changes in Endothelial and Epithelial Barrier Function and Micromotion.

Non-olfactory cells have excellent biosensor potential because they express functional olfactory receptors (ORs) and are non-neuronal cells that are easy to culture. ORs are G-protein coupled receptors (GPCRs), and there is a well-established link between different classes of G-proteins and cytoskeletal structure changes affecting cellular morphology that has been unexplored for odorant sensing. Thus, the present study was conducted to determine if odorant binding in non-olfactory cells causes cytoskeletal changes that will lead to cell changes detectable by electric cell-substrate impedance sensing (ECIS). To this end, we used the human umbilical vein endothelial cells (HUVECs), which express OR10J5, and the human keratinocyte (HaCaT) cells, which express OR2AT4. Using these two different cell barriers, we showed that odorant addition, lyral and Sandalore, respectively, caused an increase in cAMP, changes in the organization of the cytoskeleton, and a decrease in the integrity of the junctions between the cells, causing a decrease in cellular electrical resistance. In addition, the random cellular movement of the monolayers (micromotion) was significantly decreased after odorant exposure. Collectively, these data demonstrate a new physiological role of olfactory receptor signaling in endothelial and epithelial cell barriers and represent a new label-free method to detect odorant binding.

Environmental and genetic perturbations reveal different networks of metabolic regulation.

Progress in systems biology depends on accurate descriptions of biological networks. Connections in a regulatory network are identified as correlations of gene expression across a set of environmental or genetic perturbations. To use this information to predict system behavior, we must test how the nature of perturbations affects topologies of networks they reveal. To probe this question, we focused on metabolism of Drosophila melanogaster. Our source of perturbations is a set of crosses among 92 wild-derived lines from five populations, replicated in a manner permitting separate assessment of the effects of genetic variation and environmental fluctuation. We directly assayed activities of enzymes and levels of metabolites. Using a multivariate Bayesian model, we estimated covariance among metabolic parameters and built fine-grained probabilistic models of network topology. The environmental and genetic co-regulation networks are substantially the same among five populations. However, genetic and environmental perturbations reveal qualitative differences in metabolic regulation, suggesting that environmental shifts, such as diet modifications, produce different systemic effects than genetic changes, even if the primary targets are the same.

Eradication of avian leukosis virus subgroups J and K in broiler cross chickens by selection against infected birds using multilocus PCR.

The avian leukosis virus (ALV) is a serious threat to sustainable and economically viable commercial poultry management world-wide. Active infections can result in more than 20% flock loss, resulting in significant economic damage. ALV detection and elimination from flocks and breeding programs is complicated by high sequence variability and the presence of endogenous virus copies which show up as false positives in assays. Previously-developed approaches to virus detection are either too labor-intensive to implement on an industrial scale or suffer from high false negative or positive rates. We developed a novel multi-locus multiplex quantitative real-time PCR system to detect viruses belonging to the J and K genetic subgroups that are particularly prevalent in our region. We used this system to eradicate ALV from our broiler breeding program comprising thousands of individuals. Our approach can be generalized to other ALV subgroups and other highly genetically diverse pathogens.

Phenotypic Variation and the Impact of Admixture in the Oryza rufipogon Species Complex (ORSC).

Crop wild relatives represent valuable reservoirs of variation for breeding, but their populations are threatened in natural habitats, are sparsely represented in genebanks, and most are poorly characterized. The focus of this study is the Oryza rufipogon species complex (ORSC), wild progenitor of Asian rice (Oryza sativa L.). The ORSC comprises perennial, annual and intermediate forms which were historically designated as O. rufipogon, O. nivara, and O. sativa f. spontanea (or Oryza spp., an annual form of mixed O. rufipogon/O. nivara and O. sativa ancestry), respectively, based on non-standardized morphological, geographical, and/or ecologically-based species definitions and boundaries. Here, a collection of 240 diverse ORSC accessions, characterized by genotyping-by-sequencing (113,739 SNPs), was phenotyped for 44 traits associated with plant, panicle, and seed morphology in the screenhouse at the International Rice Research Institute, Philippines. These traits included heritable phenotypes often recorded as characterization data by genebanks. Over 100 of these ORSC accessions were also phenotyped in the greenhouse for 18 traits in Stuttgart, Arkansas, and 16 traits in Ithaca, New York, United States. We implemented a Bayesian Gaussian mixture model to infer accession groups from a subset of these phenotypic data and ascertained three phenotype-based group assignments. We used concordance between the genotypic subpopulations and these phenotype-based groups to identify a suite of phenotypic traits that could reliably differentiate the ORSC populations, whether measured in tropical or temperate regions. The traits provide insight into plant morphology, life history (perenniality versus annuality) and mating habit (self- versus cross-pollinated), and are largely consistent with genebank species designations. One phenotypic group contains predominantly O. rufipogon accessions characterized as perennial and largely out-crossing and one contains predominantly O. nivara accessions characterized as annual and largely inbreeding. From these groups, 42 "core" O. rufipogon and 25 "core" O. nivara accessions were identified for domestication studies. The third group, comprising 20% of our collection, has the most accessions identified as Oryza spp. (51.2%) and levels of O. sativa admixture accounting for more than 50% of the genome. This third group is potentially useful as a "pre-breeding" pool for breeders attempting to incorporate novel variation into elite breeding lines.

Evolution of protein-coding genes in Drosophila.

Several contributing factors have been implicated in evolutionary rate heterogeneity among proteins, but their evolutionary mechanisms remain poorly characterized. The recently sequenced 12 Drosophila genomes provide a unique opportunity to shed light on these unresolved issues. Here, we focus on the role of natural selection in shaping evolutionary rates. We use the Drosophila genomic data to distinguish between factors that increase the strength of purifying selection on proteins and factors that affect the amount of positive selection experienced by proteins. We confirm the importance of translational selection in shaping protein evolution in Drosophila and show that factors such as tissue bias in expression, gene essentiality, intron number, and recombination rate also contribute to evolutionary rate variation among proteins.

Low Additive Genetic Variation in a Trait Under Selection in Domesticated Rice.

Quantitative traits are important targets of both natural and artificial selection. The genetic architecture of these traits and its change during the adaptive process is thus of fundamental interest. The fate of the additive effects of variants underlying a trait receives particular attention because they constitute the genetic variation component that is transferred from parents to offspring and thus governs the response to selection. While estimation of this component of phenotypic variation is challenging, the increasing availability of dense molecular markers puts it within reach. Inbred plant species offer an additional advantage because phenotypes of genetically identical individuals can be measured in replicate. This makes it possible to estimate marker effects separately from the contribution of the genetic background not captured by genotyped loci. We focused on root growth in domesticated rice, Oryza sativa, under normal and aluminum (Al) stress conditions, a trait under recent selection because it correlates with survival under drought. A dense single nucleotide polymorphism (SNP) map is available for all accessions studied. Taking advantage of this map and a set of Bayesian models, we assessed additive marker effects. While total genetic variation accounted for a large proportion of phenotypic variance, marker effects contributed little information, particularly in the Al-tolerant tropical japonica population of rice. We were unable to identify any loci associated with root growth in this population. Models estimating the aggregate effects of all measured genotypes likewise produced low estimates of marker heritability and were unable to predict total genetic values accurately. Our results support the long-standing conjecture that additive genetic variation is depleted in traits under selection. We further provide evidence that this depletion is due to the prevalence of low-frequency alleles that underlie the trait.

Multiple Small-Effect Alleles of Indica Origin Enhance High Iron-Associated Stress Tolerance in Rice Under Field Conditions in West Africa.

Understanding the genetics of field-based tolerance to high iron-associated (HIA) stress in rice can accelerate the development of new varieties with enhanced yield performance in West African lowland ecosystems. To date, few field-based studies have been undertaken to rigorously evaluate rice yield performance under HIA stress conditions. In this study, two NERICA × O. sativa bi-parental rice populations and one O.sativa diversity panel consisting of 296 rice accessions were evaluated for grain yield and leaf bronzing symptoms over multiple years in four West African HIA stress and control sites. Mapping of these traits identified a large number of QTLs and single nucleotide polymorphisms (SNPs) associated with stress tolerance in the field. Favorable alleles associated with tolerance to high levels of iron in anaerobic rice soils were rare and almost exclusively derived from the indica subpopulation, including the most favorable alleles identified in NERICA varieties. These findings highlight the complex genetic architecture underlying rice response to HIA stress and suggest that a recurrent selection program focusing on an expanded indica genepool could be productively used in combination with genomic selection to increase the efficiency of selection in breeding programs designed to enhance tolerance to this prevalent abiotic stress in West Africa.

Evolutionary constraint and adaptation in the metabolic network of Drosophila.

Organisms must carefully control their metabolism in order to survive. On the other hand, enzymes must adapt in response to evolutionary pressures on the pathways in which they are imbedded. Taking advantage of the newly available whole-genome sequences of 12 Drosophila species, we examined how protein function and metabolic network architecture influence rates of enzyme evolution. We found that despite high overall constraint, there were significant differences in rates of amino acid substitution among functional classes of enzymes. This heterogeneity arises because proteins involved in the metabolism of foreign compounds evolve relatively rapidly, whereas enzymes that act in "core" metabolism exhibit much slower rates of amino acid replacement, suggesting strong selective constraint. Network architecture also influences enzymes' rates of amino acid replacement. In particular, enzymes that share metabolites with many other enzymes are relatively constrained, although apparently not because they are more likely to be essential. Our analyses suggest that this pattern is driven by strong constraint of enzymes acting at branch points in metabolic pathways. We conclude that metabolic network architecture and enzyme function separately affect enzyme evolution rates.

Fast Ordered Sampling of DNA Sequence Variants.

Explosive growth in the amount of genomic data is matched by increasing power of consumer-grade computers. Even applications that require powerful servers can be quickly tested on desktop or laptop machines if we can generate representative samples from large data sets. I describe a fast and memory-efficient implementation of an on-line sampling method developed for tape drives 30 years ago. Focusing on genotype files, I test the performance of this technique on modern solid-state and spinning hard drives, and show that it performs well compared to a simple sampling scheme. I illustrate its utility by developing a method to quickly estimate genome-wide patterns of linkage disequilibrium (LD) decay with distance. I provide open-source software that samples loci from several variant format files, a separate program that performs LD decay estimates, and a C++ library that lets developers incorporate these methods into their own projects.

Proper control of genetic background with precise allele substitution: a comment on Coyne and Elwyn.

Coyne and Elwyn report that, using Drosophila stocks provided by us, they were unable to replicate our experiments measuring the effects of desaturase-2 on stress tolerance. In this note, we provide evidence that these authors did not properly control for the differences in genetic background between the lines. Their experiments are thus not meaningful replications of our previous ones. We discuss ways of circumventing this problem in studies involving induced mutations in single genes.

Open access resources for genome-wide association mapping in rice.

Increasing food production is essential to meet the demands of a growing human population, with its rising income levels and nutritional expectations. To address the demand, plant breeders seek new sources of genetic variation to enhance the productivity, sustainability and resilience of crop varieties. Here we launch a high-resolution, open-access research platform to facilitate genome-wide association mapping in rice, a staple food crop. The platform provides an immortal collection of diverse germplasm, a high-density single-nucleotide polymorphism data set tailored for gene discovery, well-documented analytical strategies, and a suite of bioinformatics resources to facilitate biological interpretation. Using grain length, we demonstrate the power and resolution of our new high-density rice array, the accompanying genotypic data set, and an expanded diversity panel for detecting major and minor effect QTLs and subpopulation-specific alleles, with immediate implications for rice improvement.