Search

Abstract Background RNA-seq is a fundamental technique in genomics, yet reference bias, where transcripts derived from non-reference alleles are quantified less accurately, can undermine the accuracy of RNA-seq quantification and thus the conclusions made downstream. Reference bias in RNA-seq analysis has yet to be explored in complex polyploid genomes despite evidence that they are often a complex mosaic of wild relative introgressions, which introduce blocks of highly divergent genes. Results Here we use hexaploid wheat as a model complex polyploid, using both simulated and experimental data to show that RNA-seq alignment in wheat suffers from widespread reference bias which is largely driven by divergent introgressed genes. This leads to underestimation of gene expression and incorrect assessment of homoeologue expression balance. By incorporating gene models from ten wheat genome assemblies into a pantranscriptome reference, we present a novel method to reduce reference bias, which can be readily scaled to capture more variation as new genome and transcriptome data becomes available. Conclusions This study shows that the presence of introgressions can lead to reference bias in wheat RNA-seq analysis. Caution should be exercised by researchers using non-sample reference genomes for RNA-seq alignment and novel methods, such as the one presented here, should be considered.

Abstract Date palm (Phoenix dactylifera L.) is able to grow and complete its life cycle while being rooted in highly saline soils. Which of the many well‐known salt‐tolerance strategies are combined to fine‐tune this remarkable resilience is unknown. The precise location, whether in the shoot or the root, where these strategies are employed remains uncertain, leaving us unaware of how the various known salt‐tolerance mechanisms are integrated to fine‐tune this remarkable resilience. To address this shortcoming, we exposed date palm to a salt stress dose equivalent to seawater for up to 4 weeks and applied integrative multi‐omics analyses followed by targeted metabolomics, hormone, and ion analyses. Integration of proteomic into transcriptomic data allowed a view beyond simple correlation, revealing a remarkably high degree of convergence between gene expression and protein abundance. This sheds a clear light on the acclimatization mechanisms employed, which depend on reprogramming of protein biosynthesis. For growth in highly saline habitats, date palm effectively combines various salt‐tolerance mechanisms found in both halophytes and glycophytes: “avoidance” by efficient sodium and chloride exclusion at the roots, and “acclimation” by osmotic adjustment, reactive oxygen species scavenging in leaves, and remodeling of the ribosome‐associated proteome in salt‐exposed root cells. Combined efficiently as in P. dactylifera L., these sets of mechanisms seem to explain the palm's excellent salt stress tolerance.

Abstract Finger millet is a key food security crop widely grown in eastern Africa, India and Nepal. Long considered a ‘poor man’s crop’, finger millet has regained attention over the past decade for its climate resilience and the nutritional qualities of its grain. To bring finger millet breeding into the 21st century, here we present the assembly and annotation of a chromosome-scale reference genome. We show that this ~1.3 million years old allotetraploid has a high level of homoeologous gene retention and lacks subgenome dominance. Population structure is mainly driven by the differential presence of large wild segments in the pericentromeric regions of several chromosomes. Trait mapping, followed by variant analysis of gene candidates, reveals that loss of purple coloration of anthers and stigma is associated with loss-of-function mutations in the finger millet orthologs of the maize R1/B1 and Arabidopsis GL3/EGL3 anthocyanin regulatory genes. Proanthocyanidin production in seed is not affected by these gene knockouts.

Tinker et al. identified the position and effects of major QTLs relative to a new fully annotated reference genome in five recombinant inbred line populations representing nine diverse oat (Avena sativa) varieties. They also characterized two major chromosome rearrangements that may impact breeding targets affected by QTL that are located in these regions.

Abstract Aegilops tauschiiAe. tauschiiAe. tauschii

Abstract Chromosome-scale genome sequence assemblies underpin pan-genomic studies. Recent genome assembly efforts in the large-genome Triticeae crops wheat and barley have relied on the commercial closed-source assembly algorithm DeNovoMagic. We present TRITEX, an open-source computational workflow that combines paired-end, mate-pair, 10X Genomics linked-read with chromosome conformation capture sequencing data to construct sequence scaffolds with megabase-scale contiguity ordered into chromosomal pseudomolecules. We evaluate the performance of TRITEX on publicly available sequence data of tetraploid wild emmer and hexaploid bread wheat, and construct an improved annotated reference genome sequence assembly of the barley cultivar Morex as a community resource.

Abstract Background Recent improvements in DNA sequencing and genome scaffolding have paved the way to generate high-quality de novo assemblies of pseudomolecules representing complete chromosomes of wheat and its wild relatives. These assemblies form the basis to compare the dynamics of wheat genomes on a megabase scale. Results Here, we provide a comparative sequence analysis of the 700-megabase chromosome 2D between two bread wheat genotypes—the old landrace Chinese Spring and the elite Swiss spring wheat line ‘CH Campala Lr22a’. Both chromosomes were assembled into megabase-sized scaffolds. There is a high degree of sequence conservation between the two chromosomes. Analysis of large structural variations reveals four large indels of more than 100 kb. Based on the molecular signatures at the breakpoints, unequal crossing over and double-strand break repair were identified as the molecular mechanisms that caused these indels. Three of the large indels affect copy number of NLRs, a gene family involved in plant immunity. Analysis of SNP density reveals four haploblocks of 4, 8, 9 and 48 Mb with a 35-fold increased SNP density compared to the rest of the chromosome. Gene content across the two chromosomes was highly conserved. Ninety-nine percent of the genic sequences were present in both genotypes and the fraction of unique genes ranged from 0.4 to 0.7%. Conclusions This comparative analysis of two high-quality chromosome assemblies enabled a comprehensive assessment of large structural variations and gene content. The insight obtained from this analysis will form the basis of future wheat pan-genome studies.

Vibrio vulnificus is a bacterium responsible for severe gastroenteritis, sepsis and wound infections. Gastroenteritis and sepsis are commonly associated with the consumption of raw oysters, whereas wound infection is often associated with the handling of contaminated fish. Although classical virulence factors of this emerging pathogen are well characterised, there remains a paucity of knowledge regarding the general biology of this species. To investigate the presence of previously unreported virulence factors, we applied whole genome sequencing to a panel of ten V. vulnificus strains with varying virulence potentials. This identified two novel type 6 secretion systems (T6SSs), systems that are known to have a role in bacterial virulence and population dynamics. By utilising a range of molecular techniques and assays we have demonstrated the functionality of one of these T6SSs. Furthermore, we have shown that this system is subject to thermoregulation and is negatively regulated by increasing salinity concentrations. This secretion system was also shown to be involved in the killing of V. vulnificus strains that did not possess this system and a model is proposed as to how this interaction may contribute to population dynamics within V. vulnificus strains. In addition to this intra-species killing, this system also contributes to the killing of inter bacterial species and may have a role in the general composition of Vibrio species in the environment.

Most eukaryotic oleaginous species are yeasts and among them the basidiomycete red yeast, Rhodotorula (Rhodosporidium) toruloides (Pucciniomycotina) is known to produce high quantities of lipids when grown in nitrogen-limiting media, and has potential for biodiesel production. The genome of the CGMCC 2.1609 strain of this oleaginous red yeast was sequenced using a hybrid of Roche 454 and Illumina technology generating 13× coverage. The de novo assembly was carried out using MIRA and scaffolded using MAQ and BAMBUS. The sequencing and assembly resulted in 365 scaffolds with total genome size of 33.4 Mb. The complete genome sequence of this strain was deposited in GenBank and the accession number is LKER00000000. The annotation is available on Figshare (doi:10.6084/m9.figshare.4754251).

We study the structural correlations in the Italian overnight money market over the period 1999–2010. We show that the structural correlations vary across different versions of the network. Moreover, we employ different configuration models and examine whether higher-level characteristics of the observed network can be statistically reconstructed by maximizing the entropy of a randomized ensemble of networks restricted only by the lower-order features of the observed network. We find that often many of the high order correlations in the observed network can be considered emergent from the information embedded in the degree sequence in the binary version and in both the degree and strength sequences in the weighted version. However, this information is not enough to allow the models to account for all the patterns in the observed higher order structural correlations. In particular, one of the main features of the observed network that remains unexplained is the abnormally high level of weighted clustering in the years preceding the crisis, i.e., the huge increase in various indirect exposures generated via more intensive interbank credit links.

BACKGROUND: In the current era of strong worldwide market couplings the global financial village became highly prone to systemic collapses, events that can rapidly sweep throughout the entire village. METHODOLOGY/PRINCIPAL FINDINGS: We present a new methodology to assess and quantify inter-market relations. The approach is based on the correlations between the market index, the index volatility, the market Index Cohesive Force and the meta-correlations (correlations between the intra-correlations.) We investigated the relations between six important world markets--U.S., U.K., Germany, Japan, China and India--from January 2000 until December 2010. We found that while the developed "western" markets (U.S., U.K., Germany) are highly correlated, the interdependencies between these markets and the developing "eastern" markets (India and China) are volatile and with noticeable maxima at times of global world events. The Japanese market switches "identity"--it switches between periods of high meta-correlations with the "western" markets and periods when it behaves more similarly to the "eastern" markets. CONCLUSIONS/SIGNIFICANCE: The methodological framework presented here provides a way to quantify the evolvement of interdependencies in the global market, evaluate a world financial network and quantify changes in the world inter market relations. Such changes can be used as precursors to the agitation of the global financial village. Hence, the new approach can help to develop a sensitive "financial seismograph" to detect early signs of global financial crises so they can be treated before they develop into worldwide events.