The barley pan-genome reveals the hidden legacy of mutation breeding.

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

HvGST plays a key role in anthocyanin accumulation in colored barley.

Blue aleurone of barley is caused by the accumulation of delphinidin-based derivatives. Although these compounds are ideal nutrients for human health, they are undesirable contaminants in malt brewing. Therefore, the ability to add and remove this trait easily would facilitate breeding barley for different purposes. Here we identified a glutathione S-transferase gene (HvGST) that was responsible for the blue aleurone trait in Tibetan qingke barley by performing a genome-wide association study and RNA-sequencing analysis. Gene variation and expression analysis indicated that HvGST also participates in the transport and accumulation of anthocyanin in purple barley. Haplotype and the geographic distribution analyses of HvGST alleles revealed two independent natural variants responsible for the emergence of white aleurone: a 203-bp deletion causing premature termination of translation in qingke barley and two key single nucleotide polymorphisms in the promoter resulting in low transcription in Western barley. This study contributes to a better understanding of mechanisms of colored barley formation, and provides a comprehensive reference for marker-assisted barley breeding.

TGFBR3 supports anoikis through suppressing ATF4 signaling.

Epithelial morphogenesis and oncogenic transformation can cause loss of cell adhesion, and detached cells are eliminated by anoikis. Here, we reveal that transforming growth factor ß receptor 3 (TGFBR3) acts as an anoikis mediator through the coordination of activating transcription factor 4 (ATF4). In breast cancer tissues, TGFBR3 is progressively lost, but elevated TGFBR3 is associated with a histologic subtype characterized by cellular adhesion defects. Dissecting the impact of extracellular matrix (ECM) deprivation, we demonstrate that ECM loss promotes TGFBR3 expression, which in turn causes differentiation of cell aggregates, conferring a low-adhesion phenotype, and drives the intrinsic apoptotic pathway. We demonstrate that inhibition of TGFBR3 impairs epithelial anoikis by activating ATF4 signaling. These preclinical findings provide a rationale for therapeutic inhibition of ATF4 in the subgroup of breast cancer patients with low TGFBR3 expression.

Genomic variation in 3,010 diverse accessions of Asian cultivated rice.

Here we analyse genetic variation, population structure and diversity among 3,010 diverse Asian cultivated rice (Oryza sativa L.) genomes from the 3,000 Rice Genomes Project. Our results are consistent with the five major groups previously recognized, but also suggest several unreported subpopulations that correlate with geographic location. We identified 29 million single nucleotide polymorphisms, 2.4 million small indels and over 90,000 structural variations that contribute to within- and between-population variation. Using pan-genome analyses, we identified more than 10,000 novel full-length protein-coding genes and a high number of presence-absence variations. The complex patterns of introgression observed in domestication genes are consistent with multiple independent rice domestication events. The public availability of data from the 3,000 Rice Genomes Project provides a resource for rice genomics research and breeding.

CoreSNP: an efficient pipeline for core marker profile selection from genome-wide SNP datasets in crops.

BACKGROUND: DNA marker profiles play a crucial role in the identification and registration of germplasm, as well as in the distinctness, uniformity, and stability (DUS) testing of new plant variety protection. However, selecting minimal marker sets from large-scale SNP dataset can be challenging to distinguish a maximum number of samples. RESULTS: Here, we developed the CoreSNP pipeline using a "divide and conquer" strategy and a "greedy" algorithm. The pipeline offers adjustable parameters to guarantee the distinction of each sample pair with at least two markers. Additionally, it allows datasets with missing loci as input. The pipeline was tested in barley, soybean, wheat, rice and maize. A few dozen of core SNPs were efficiently selected in different crops with SNP array, GBS, and WGS dataset, which can differentiate thousands of individual samples. The core SNPs were distributed across all chromosomes, exhibiting lower pairwise linkage disequilibrium (LD) and higher polymorphism information content (PIC) and minor allele frequencies (MAF). It was shown that both the genetic diversity of the population and the characteristics of the original dataset can significantly influence the number of core markers. In addition, the core SNPs capture a certain level of the original population structure. CONCLUSIONS: CoreSNP is an efficiency way of core marker sets selection based on Genome-wide SNP datasets of crops. Combined with low-density SNP chip or genotyping technologies, it can be a cost-effective way to simplify and expedite the evaluation of genetic resources and differentiate different crop varieties. This tool is expected to have great application prospects in the rapid comparison of germplasm and intellectual property protection of new varieties.

Positional cloning identified HvTUBULIN8 as the candidate gene for round lateral spikelet (RLS) in barley (Hordeum vulgare L.).

KEY MESSAGE: Map-based cloning, subcellular localization, virus-induced-gene-silencing and transcriptomic analysis reveal HvTUB8 as a candidate gene with pleiotropic effects on barley spike and leaf development via ethylene and chlorophyll metabolism. Barley lateral spikelet morphology and grain shape play key roles in grain physical quality and yield. Several genes and QTLs for these traits have been cloned or fine mapped previously. Here, we report the phenotypic and genotypic analysis of a barley mutant with round lateral spikelet (rls) from cv. Edamai 934. rls had round lateral spikelet, short but round grain, shortened awn, thick glume and dark green leaves. Histocytologic and ultrastructural analysis revealed that the difference of grain shape of rls was caused by change of cell arrangement in glume, and the dark leaf color resulted from enlarged chloroplast. HvTUBULIN8 (HvTUB8) was identified as the candidate gene for rls by combination of RNA-Seq, map-based-cloning, virus-induced-gene-silencing (VIGS) and protein subcellular location. A single G-A substitution at the third exon of HvTUB8 resulted in change of Cysteine 354 to tyrosine. Furthermore, the mutant isoform Hvtub8 could be detected in both nucleus and cytoplasm, whereas the wild-type protein was only in cytoplasm and granular organelles of wheat protoplasts. Being consistent with the rare phenotype, the "A" allele of HvTUB8 was only detected in rls, but not in a worldwide barley germplasm panel with 400 accessions. VIGS confirmed that HvTUB8 was essential to maintain spike integrity. RNA-Seq results suggested that HvTUB8 may control spike morphogenesis via ethylene homeostasis and signaling, and control leaf color through chlorophyll metabolism. Collectively, our results support HvTUB8 as a candidate gene for barley spike and leaf morphology and provide insight of a novel mechanism of it in barley development.

Genome-Wide Association Study Identifies Resistance Loci for Bacterial Blight in a Collection of Asian Temperate Japonica Rice Germplasm.

Growing resistant rice cultivars is the most effective strategy to control bacterial blight (BB), a devastating disease caused by Xanthomonas oryzae pv. oryzae (Xoo). Screening resistant germplasm and identifying resistance (R) genes are prerequisites for breeding resistant rice cultivars. We conducted a genome-wide association study (GWAS) to detect quantitative trait loci (QTL) associated with BB resistance using 359 East Asian temperate Japonica accessions inoculated with two Chinese Xoo strains (KS6-6 and GV) and one Philippine Xoo strain (PXO99A). Based on the 55K SNPs Array dataset of the 359 Japonica accessions, eight QTL were identified on rice chromosomes 1, 2, 4, 10, and 11. Four of the QTL coincided with previously reported QTL, and four were novel loci. Six R genes were localized in the qBBV-11.1, qBBV-11.2, and qBBV-11.3 loci on chromosome 11 in this Japonica collection. Haplotype analysis revealed candidate genes associated with BB resistance in each QTL. Notably, LOC_Os11g47290 in qBBV-11.3, encoding a leucine-rich repeat receptor-like kinase, was a candidate gene associated with resistance to the virulent strain GV. Knockout mutants of Nipponbare with the susceptible haplotype of LOC_Os11g47290 exhibited significantly improved BB resistance. These results will be useful for cloning BB resistance genes and breeding resistant rice cultivars.

Identification and allele mining of new candidate genes underlying rice grain weight and grain shape by genome-wide association study.

BACKGROUND: Grain weight and grain shape are important agronomic traits that affect the grain yield potential and grain quality of rice. Both grain weight and grain shape are controlled by multiple genes. The 3,000 Rice Genomes Project (3 K RGP) greatly facilitates the discovery of agriculturally important genetic variants and germplasm resources for grain weight and grain shape. RESULTS: Abundant natural variations and distinct phenotic differentiation among the subgroups in grain weight and grain shape were observed in a large population of 2,453 accessions from the 3 K RGP. A total of 21 stable quantitative trait nucleotides (QTNs) for the four traits were consistently identified in at least two of 3-year trials by genome-wide association study (GWAS), including six new QTNs (qTGW3.1, qTGW9, qTGW11, qGL4/qRLW4, qGL10, and qRLW1) for grain weight and grain shape. We further predicted seven candidate genes (Os03g0186600, Os09g0544400, Os11g0163600, Os04g0580700, Os10g0399700, Os10g0400100 and Os01g0171000) for the six new QTNs by high-density association and gene-based haplotype analyses. The favorable haplotypes of the seven candidate genes and five previously cloned genes in elite accessions with high TGW and RLW are also provided. CONCLUSIONS: Our results deepen the understanding of the genetic basis of grain weight and grain shape in rice and provide valuable information for improving rice grain yield and grain quality through molecular breeding.

Aggregation behavior of aqu/nC60 produced via extended mixing: Influence of sunlight and agitation intensity.

Aggregation of C60, as an important process governing its mobility and toxicity, has been quantitatively investigated. However, effects of sunlight and agitation intensity on the aggregation behavior of aqu/nC60 produced via extended mixing, have not been clarified. Therefore, in the present study, the aggregation behavior of aqu/nC60 produced at 500 and 800 rpm in the absence and presence of sunlight was investigated. Aggregation with increasing concentrations could be accelerated, while changes of Zave and zeta potential were not obvious. Critical coagulation concentrations (CCCs) of aqu/nC60 obtained at 800 rpm in the absence/presence of sunlight and that at 500 rpm under sunlight were 330, 205 and 170 mM NaCl, and 10.0, 2.6 and 3.1 mM CaCl2, respectively. These CCCs indicated that the aqu/nC60 prepared by the extended mixing were more stable than those produced by other methods. Salt-induced aggregation occurred more easily for aqu/nC60 formed under sunlight than that formed in the dark. Extra surface oxidation induced by high agitation intensity remarkably increased the stability of aqu/nC60 in NaCl solutions. In contrast, in CaCl2 solutions, aqu/nC60 formed at high agitation intensity had similar stability or even inadequate stability to that obtained at low agitation intensity due to the charge neutralization and cross-link bridging.

RPAN: rice pan-genome browser for ∼3000 rice genomes.

A pan-genome is the union of the gene sets of all the individuals of a clade or a species and it provides a new dimension of genome complexity with the presence/absence variations (PAVs) of genes among these genomes. With the progress of sequencing technologies, pan-genome study is becoming affordable for eukaryotes with large-sized genomes. The Asian cultivated rice, Oryza sativa L., is one of the major food sources for the world and a model organism in plant biology. Recently, the 3000 Rice Genome Project (3K RGP) sequenced more than 3000 rice genomes with a mean sequencing depth of 14.3×, which provided a tremendous resource for rice research. In this paper, we present a genome browser, Rice Pan-genome Browser (RPAN), as a tool to search and visualize the rice pan-genome derived from 3K RGP. RPAN contains a database of the basic information of 3010 rice accessions, including genomic sequences, gene annotations, PAV information and gene expression data of the rice pan-genome. At least 12 000 novel genes absent in the reference genome were included. RPAN also provides multiple search and visualization functions. RPAN can be a rich resource for rice biology and rice breeding. It is available at http://cgm.sjtu.edu.cn/3kricedb/ or http://www.rmbreeding.cn/pan3k.

ccNET: Database of co-expression networks with functional modules for diploid and polyploid Gossypium.

Plant genera with both diploid and polyploid species are a common evolutionary occurrence. Polyploids, especially allopolyploids such as cotton and wheat, are a great model system for heterosis research. Here, we have integrated genome sequences and transcriptome data of Gossypium species to construct co-expression networks and identified functional modules from different cotton species, including 1155 and 1884 modules in G. arboreum and G. hirsutum, respectively. We overlayed the gene expression results onto the co-expression network. We further provided network comparison analysis for orthologous genes across the diploid and allotetraploid Gossypium We also constructed miRNA-target networks and predicted PPI networks for both cotton species. Furthermore, we integrated in-house ChIP-seq data of histone modification (H3K4me3) together with cis-element analysis and gene sets enrichment analysis tools for studying possible gene regulatory mechanism in Gossypium species. Finally, we have constructed an online ccNET database (http://structuralbiology.cau.edu.cn/gossypium) for comparative gene functional analyses at a multi-dimensional network and epigenomic level across diploid and polyploid Gossypium species. The ccNET database will be beneficial for community to yield novel insights into gene/module functions during cotton development and stress response, and might be useful for studying conservation and diversity in other polyploid plants, such as T. aestivum and Brassica napus.

Differential deposition of H2A.Z in combination with histone modifications within related genes in Oryza sativa callus and seedling.

As a histone variant, H2A.Z is highly conserved among species and plays a significant role in diverse cellular processes. Here, we generated genome-wide maps of H2A.Z in Oryza sativa (rice) callus and seedling by combining chromatin immunoprecipitation using H2A.Z antibody and high-throughput sequencing. We found a significantly high peak and a small peak of H2A.Z distributed at the 5' and 3' ends of highly expressed genes, respectively. H2A.Z was also associated with inactive genes in both tissues. H3 lysine 4 trimethylation was associated with H2A.Z deposition at the 5' end of expressed genes, and H3 lysine 27 trimethylation peaks were partially associated with H2A.Z. In summary, our study provides global analysis data for the distribution of H2A.Z in the rice genome. Our results demonstrate that the differential deposition of H2A.Z might play important roles in gene transcription during rice development.

Novel tests for evaluating two ROC curves under paired samples.

Disease prevention is important and can be accomplished by developing diagnostic tests. The receiver operating characteristic (ROC) curve and the area under the ROC curve (AUC) are used to assess the accuracy of diagnostic tests. The assessment for the superiority between evaluating two diagnostic tests is needed when comparing two diagnostic tests. Existing tests are constructed by comparing two AUCs under the paired samples. Nevertheless, it is problematic when two ROC curves are crossing. This article proposes a test that takes into account the possible correlation between pairs. Simulations are conducted to evaluate the feasibility of the test.

Down-regulation of OsSPX1 caused semi-male sterility, resulting in reduction of grain yield in rice.

OsSPX1, a rice SPX domain gene, involved in the phosphate (Pi)-sensing mechanism plays an essential role in the Pi-signalling network through interaction with OsPHR2. In this study, we focused on the potential function of OsSPX1 during rice reproductive phase. Based on investigation of OsSPX1 antisense and sense transgenic rice lines in the paddy fields, we discovered that the down-regulation of OsSPX1 caused reduction of seed-setting rate and filled grain number. Through examination of anthers and pollens of the transgenic and wild-type plants by microscopy, we found that the antisense of OsSPX1 gene led to semi-male sterility, with lacking of mature pollen grains and phenotypes with a disordered surface of anthers and pollens. We further conducted rice whole-genome GeneChip analysis to elucidate the possible molecular mechanism underlying why the down-regulation of OsSPX1 caused deficiencies in anthers and pollens and lower seed-setting rate in rice. The down-regulation of OsSPX1 significantly affected expression of genes involved in carbohydrate metabolism and sugar transport, anther development, cell cycle, etc. These genes may be related to pollen fertility and male gametophyte development. Our study demonstrated that down-regulation of OsSPX1 disrupted rice normal anther and pollen development by affecting carbohydrate metabolism and sugar transport, leading to semi-male sterility, and ultimately resulted in low seed-setting rate and grain yield.

JAZ7 negatively regulates dark-induced leaf senescence in Arabidopsis.

JASMONATE ZIM-domain (JAZ) proteins play important roles in plant defence and growth by regulating jasmonate signalling. Through data mining, we discovered that the JAZ7 gene was up-regulated in darkness. In the dark, the jaz7 mutant displayed more severe leaf yellowing, quicker chlorophyll degradation, and higher hydrogen peroxide accumulation compared with wild-type (WT) plants. The mutant phenotype of dark-induced leaf senescence could be rescued in the JAZ7-complemented and -overexpression lines. Moreover, the double mutants of jaz7 myc2 and jaz7 coi1 exhibited delayed leaf senescence. We further employed GeneChip analysis to study the molecular mechanism. Some key genes down-regulated in the triple mutant myc2 myc3 myc4 were up-regulated in the jaz7 mutant under darkness. The Gene Ontology terms 'leaf senescence' and 'cell death' were significantly enriched in the differentially expressed genes. Combining the genetic and transcriptomic analyses together, we proposed a model whereby darkness can induce JAZ7, which might further block MYC2 to suppress dark-induced leaf senescence. In darkness, the mutation of JAZ7 might partially liberate MYC2/MYC3/MYC4 from suppression, leading the MYC proteins to bind to the G-box/G-box-like motifs in the promoters, resulting in the up-regulation of the downstream genes related to indole-glucosinolate biosynthesis, sulphate metabolism, callose deposition, and JA-mediated signalling pathways. In summary, our genetic and transcriptomic studies established the JAZ7 protein as an important regulator in dark-induced leaf senescence.

Joint model of multiple longitudinal measures and a binary outcome: An application to predict orthostatic hypertension for subacute stroke patients.

Stroke patients with orthostatic hypertensive responses that are one of the blood pressure regulation problems can easily fall down while doing rehabilitation, which may result in prolonged hospitalization and delayed treatment and recovery. This may result in increasing the medical cost and burden. In turn, developing a diagnostic test for the orthostatic hypertension (OH) is clinically important for patients who are suffering from stroke. Clinically, a patient needs to have a tilt testing that requires measuring the change of blood pressures and heart rate at all angles to determine whether a stroke patient has OH. It takes lots of time and effort to perform the test. Assuming there exist measurement errors when obtaining the blood pressures and heart rate at all angles, this paper proposes using multiple mixed-effect models to obtain the true trajectories of these measurements, which take into account the measurement error and the possible correlation among multiple measurements, and a logistic regression uses these true trajectories at a given time and other fixed-effect covariates as predictors to predict the status of OH. The joint likelihood function is derived to estimate parameters and the area under the receiver operating characteristics curve is used to estimate the predictive power of the model. Monte Carlo simulations are performed to evaluate the feasibility of the proposed methods. Also, the proposed model is implemented in the real data and provides an acceptable predictive power.

Natural variations of HvSRN1 modulate the spike rachis node number in barley.

Grain number, one of the major determinants of yield in Triticeae crops, is largely determined by spikelet number and spike rachis node number (SRN). Here, we identified three quantitative trait loci (QTLs) for SRN using 145 recombinant inbred lines derived from a barley R90/1815D cross. qSRN1, the major-effect QTL, was mapped to chromosome 2H and explained up to 38.77% of SRN variation. Map-based cloning revealed that qSRN1 encodes the RAWUL domain-containing protein HvSRN1. Further analysis revealed that two key SNPs in the HvSRN1 promoter region (∼2 kb upstream of the transcription start site) affect the transcript level of HvSRN1 and contribute to variation in SRN. Similar to its orthologous proteins OsLAX2 and ZmBA2, HvSRN1 showed protein-protein interactions with HvLAX1, suggesting that the LAX2-LAX1 model for spike morphology regulation may be conserved in Poaceae crops. CRISPR-Cas9-induced HvSRN1 mutants showed reduced SRN but increased grain size and weight, demonstrating a trade-off effect. Our results shed light on the role of HvSRN1 variation in regulating the balance between grain number and weight in barley.

Identifying single-cell molecular programs by stochastic profiling.

Cells in tissues can be morphologically indistinguishable yet show molecular expression patterns that are remarkably heterogeneous. Here we describe an approach to comprehensively identify co-regulated, heterogeneously expressed genes among cells that otherwise appear identical. The technique, called stochastic profiling, involves repeated, random selection of very small cell populations via laser-capture microdissection followed by a customized single-cell amplification procedure and transcriptional profiling. Fluctuations in the resulting gene-expression measurements are then analyzed statistically to identify transcripts that are heterogeneously coexpressed. We stochastically profiled matrix-attached human epithelial cells in a three-dimensional culture model of mammary-acinar morphogenesis. Of 4,557 transcripts, we identified 547 genes with strong cell-to-cell expression differences. Clustering of this heterogeneous subset revealed several molecular 'programs' implicated in protein biosynthesis, oxidative-stress responses and NF-kappaB signaling, which we independently confirmed by RNA fluorescence in situ hybridization. Thus, stochastic profiling can reveal single-cell heterogeneities without the need to measure expression in individual cells.