Genetic architecture of root and shoot ionomes in rice (Oryza sativa L.).

KEY MESSAGE: Association analysis for ionomic concentrations of 20 elements identified independent genetic factors underlying the root and shoot ionomes of rice, providing a platform for selecting and dissecting causal genetic variants. Understanding the genetic basis of mineral nutrient acquisition is key to fully describing how terrestrial organisms interact with the non-living environment. Rice (Oryza sativa L.) serves both as a model organism for genetic studies and as an important component of the global food system. Studies in rice ionomics have primarily focused on above ground tissues evaluated from field-grown plants. Here, we describe a comprehensive study of the genetic basis of the rice ionome in both roots and shoots of 6-week-old rice plants for 20 elements using a controlled hydroponics growth system. Building on the wealth of publicly available rice genomic resources, including a panel of 373 diverse rice lines, 4.8 M genome-wide single-nucleotide polymorphisms, single- and multi-marker analysis pipelines, an extensive tome of 321 candidate genes and legacy QTLs from across 15 years of rice genetics literature, we used genome-wide association analysis and biparental QTL analysis to identify 114 genomic regions associated with ionomic variation. The genetic basis for root and shoot ionomes was highly distinct; 78 loci were associated with roots and 36 loci with shoots, with no overlapping genomic regions for the same element across tissues. We further describe the distribution of phenotypic variation across haplotypes and identify candidate genes within highly significant regions associated with sulfur, manganese, cadmium, and molybdenum. Our analysis provides critical insight into the genetic basis of natural phenotypic variation for both root and shoot ionomes in rice and provides a comprehensive resource for dissecting and testing causal genetic variants.

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.

Loss of function at RAE2, a previously unidentified EPFL, is required for awnlessness in cultivated Asian rice.

Domestication of crops based on artificial selection has contributed numerous beneficial traits for agriculture. Wild characteristics such as red pericarp and seed shattering were lost in both Asian (Oryza sativa) and African (Oryza glaberrima) cultivated rice species as a result of human selection on common genes. Awnedness, in contrast, is a trait that has been lost in both cultivated species due to selection on different sets of genes. In a previous report, we revealed that at least three loci regulate awn development in rice; however, the molecular mechanism underlying awnlessness remains unknown. Here we isolate and characterize a previously unidentified EPIDERMAL PATTERNING FACTOR-LIKE (EPFL) family member named REGULATOR OF AWN ELONGATION 2 (RAE2) and identify one of its requisite processing enzymes, SUBTILISIN-LIKE PROTEASE 1 (SLP1). The RAE2 precursor is specifically cleaved by SLP1 in the rice spikelet, where the mature RAE2 peptide subsequently induces awn elongation. Analysis of RAE2 sequence diversity identified a highly variable GC-rich region harboring multiple independent mutations underlying protein-length variation that disrupt the function of the RAE2 protein and condition the awnless phenotype in Asian rice. Cultivated African rice, on the other hand, retained the functional RAE2 allele despite its awnless phenotype. Our findings illuminate the molecular function of RAE2 in awn development and shed light on the independent domestication histories of Asian and African cultivated rice.

Genome-Wide Convergence during Evolution of Mangroves from Woody Plants.

When living organisms independently invade a new environment, the evolution of similar phenotypic traits is often observed. An interesting but contentious issue is whether the underlying molecular biology also converges in the new habitat. Independent invasions of tropical intertidal zones by woody plants, collectively referred to as mangrove trees, represent some dramatic examples. The high salinity, hypoxia, and other stressors in the new habitat might have affected both genomic features and protein structures. Here, we developed a new method for detecting convergence at conservative Sites (CCS) and applied it to the genomic sequences of mangroves. In simulations, the CCS method drastically reduces random convergence at rapidly evolving sites as well as falsely inferred convergence caused by the misinferences of the ancestral character. In mangrove genomes, we estimated ∼400 genes that have experienced convergence over the background level of convergence in the nonmangrove relatives. The convergent genes are enriched in pathways related to stress response and embryo development, which could be important for mangroves' adaptation to the new habitat.

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.

Next-generation phenotyping: requirements and strategies for enhancing our understanding of genotype-phenotype relationships and its relevance to crop improvement.

More accurate and precise phenotyping strategies are necessary to empower high-resolution linkage mapping and genome-wide association studies and for training genomic selection models in plant improvement. Within this framework, the objective of modern phenotyping is to increase the accuracy, precision and throughput of phenotypic estimation at all levels of biological organization while reducing costs and minimizing labor through automation, remote sensing, improved data integration and experimental design. Much like the efforts to optimize genotyping during the 1980s and 1990s, designing effective phenotyping initiatives today requires multi-faceted collaborations between biologists, computer scientists, statisticians and engineers. Robust phenotyping systems are needed to characterize the full suite of genetic factors that contribute to quantitative phenotypic variation across cells, organs and tissues, developmental stages, years, environments, species and research programs. Next-generation phenotyping generates significantly more data than previously and requires novel data management, access and storage systems, increased use of ontologies to facilitate data integration, and new statistical tools for enhancing experimental design and extracting biologically meaningful signal from environmental and experimental noise. To ensure relevance, the implementation of efficient and informative phenotyping experiments also requires familiarity with diverse germplasm resources, population structures, and target populations of environments. Today, phenotyping is quickly emerging as the major operational bottleneck limiting the power of genetic analysis and genomic prediction. The challenge for the next generation of quantitative geneticists and plant breeders is not only to understand the genetic basis of complex trait variation, but also to use that knowledge to efficiently synthesize twenty-first century crop varieties.

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.

Rapid Evolution of Autosomal Binding Sites of the Dosage Compensation Complex in Drosophila melanogaster and Its Association With Transcription Divergence.

How pleiotropy influences evolution of protein sequence remains unclear. The male-specific lethal (MSL) complex in Drosophila mediates dosage compensation by 2-fold upregulation of the X chromosome in males. Nevertheless, several MSL proteins also bind autosomes and likely perform functions not related to dosage compensation. Here, we study the evolution of MOF, MSL1, and MSL2 biding sites in Drosophila melanogaster and its close relative Drosophila simulans. We found pervasive expansion of the MSL binding sites in D. melanogaster, particularly on autosomes. The majority of these newly-bound regions are unlikely to function in dosage compensation and associated with an increase in expression divergence between D. melanogaster and D. simulans. While dosage-compensation related sites show clear signatures of adaptive evolution, these signatures are even more marked among autosomal regions. Our study points to an intriguing avenue of investigation of pleiotropy as a mechanism promoting rapid protein sequence evolution.

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.