Comparative Genome Analysis of Scutellaria baicalensis and Scutellaria barbata Reveals the Evolution of Active Flavonoid Biosynthesis

Scutellaria baicalensis (S. baicalensis) and Scutellaria barbata (S. barbata) are common medicinal plants of the Lamiaceae family. Both produce specific flavonoid compounds, including baicalein, scutellarein, norwogonin, and wogonin, as well as their glycosides, which exhibit antioxidant and antitumor activities. Here, we report chromosome-level genome assemblies of S. baicalensis and S. barbata with quantitative chromosomal variation (2n = 18 and 2n = 26, respectively). The divergence of S. baicalensis and S. barbata occurred far earlier than previously reported, and a whole-genome duplication (WGD) event was identified. The insertion of long terminal repeat elements after speciation might be responsible for the observed chromosomal expansion and rearrangement. Comparative genome analysis of the congeneric species revealed the species-specific evolution of chrysin and apigenin biosynthetic genes, such as the S. baicalensis-specific tandem duplication of genes encoding phenylalanine ammonia lyase and chalcone synthase, and the S. barbata-specific duplication of genes encoding 4-CoA ligase. In addition, the paralogous duplication, colinearity, and expression diversity of CYP82D subfamily members revealed the functional divergence of genes encoding flavone hydroxylase between S. baicalensis and S. barbata. Analyzing these Scutellaria genomes reveals the common and species-specific evolution of flavone biosynthetic genes. Thus, these findings would facilitate the development of molecular breeding and studies of biosynthesis and regulation of bioactive compounds.     

Population genomic footprints of selection and associations with climate in natural populations of Arabidopsis halleri from the Alps

Natural genetic variation is essential for the adaptation of organisms to their local environment and to changing environmental conditions. Here, we examine genomewide patterns of nucleotide variation in natural populations of the outcrossing herb Arabidopsis halleri and associations with climatic variation among populations in the Alps. Using a pooled population sequencing (Pool‐Seq) approach, we discovered more than two million SNPs in five natural populations and identified highly differentiated genomic regions and SNPs using F_ST‐based analyses. We tested only the most strongly differentiated SNPs for associations with a nonredundant set of environmental factors using partial Mantel tests to identify topo‐climatic factors that may underlie the observed footprints of selection. Possible functions of genes showing signatures of selection were identified by Gene Ontology analysis. We found 175 genes to be highly associated with one or more of the five tested topo‐climatic factors. Of these, 23.4% had unknown functions. Genetic variation in four candidate genes was strongly associated with site water balance and solar radiation, and functional annotations were congruent with these environmental factors. Our results provide a genomewide perspective on the distribution of adaptive genetic variation in natural plant populations from a highly diverse and heterogeneous alpine environment.     

Pedigree‐based analysis of derivation of genome segments of an elite rice reveals key regions during its breeding

Analyses of genome variations with high‐throughput assays have improved our understanding of genetic basis of crop domestication and identified the selected genome regions, but little is known about that of modern breeding, which has limited the usefulness of massive elite cultivars in further breeding. Here we deploy pedigree‐based analysis of an elite rice, Huanghuazhan, to exploit key genome regions during its breeding. The cultivars in the pedigree were resequenced with 7.6× depth on average, and 2.1 million high‐quality single nucleotide polymorphisms (SNPs) were obtained. Tracing the derivation of genome blocks with pedigree and information on SNPs revealed the chromosomal recombination during breeding, which showed that 26.22% of Huanghuazhan genome are strictly conserved key regions. These major effect regions were further supported by a QTL mapping of 260 recombinant inbred lines derived from the cross of Huanghuazhan and a very dissimilar cultivar, Shuanggui 36, and by the genome profile of eight cultivars and 36 elite lines derived from Huanghuazhan. Hitting these regions with the cloned genes revealed they include numbers of key genes, which were then applied to demonstrate how Huanghuazhan were bred after 30 years of effort and to dissect the deficiency of artificial selection. We concluded the regions are helpful to the further breeding based on this pedigree and performing breeding by design. Our study provides genetic dissection of modern rice breeding and sheds new light on how to perform genomewide breeding by design.     

Resequencing core accessions of a pedigree identifies derivation of genomic segments and key agronomic trait loci during cotton improvement

Upland cotton (Gossypium hirsutum) is the world's largest source of natural fibre and dominates the global textile industry. Hybrid cotton varieties exhibit strong heterosis that confers high fibre yields, yet the genome‐wide effects of artificial selection that have influenced Upland cotton during its breeding history are poorly understood. Here, we resequenced Upland cotton genomes and constructed a variation map of an intact breeding pedigree comprising seven elite and 19 backbone parents. Compared to wild accessions, the 26 pedigree accessions underwent strong artificial selection during domestication that has resulted in reduced genetic diversity but stronger linkage disequilibrium and higher extents of selective sweeps. In contrast to the backbone parents, the elite parents have acquired significantly improved agronomic traits, with an especially pronounced increase in the lint percentage. Notably, identify by descent (IBD) tracking revealed that the elite parents inherited abundant beneficial trait segments and loci from the backbone parents and our combined analyses led to the identification of a core genomic segment which was inherited in the elite lines from the parents Zhong 7263 and Ejing 1 and that was strongly associated with lint percentage. Additionally, SNP correlation analysis of this core segment showed that a non‐synonymous SNP (A‐to‐G) site in a gene encoding the cell wall‐associated receptor‐like kinase 3 (GhWAKL3) protein was highly correlated with increased lint percentage. Our results substantially increase the valuable genomics resources available for future genetic and functional genomics studies of cotton and reveal insights that will facilitate yield increases in the molecular breeding of cotton.     

Temporal dynamics of three culturable γ-Proteobacteria taxa in salt marsh sediments

Temporal dynamics of three marine -Proteobacteria taxa withculturable members were followed in sediment communities in southeastern USsaltmarshes using a whole-genome hybridization approach. Labeled DNA from threebacterial isolates was used to probe sediment community DNA from two saltmarshes on Sapelo Island, GA at multiple time points over a three-year period.The relative abundance of the three isolates (and their close relatives)accounted for up to 28% of the total sediment community DNA. Hybridizationsignals for the taxon represented by isolate SIGA28M (with a 16S rRNA sequencesimilarity of 97% to Vibrio proteolyticus) showedsignificant temporal variation in both marshes, varying from 1% to 28% over thethree-year period. The taxa represented by isolates SIGA172a and SIGA198 (16SrRNA sequence similarities of 90% and 92% to Shewanellafrigidimarina and Vibrio nigripulchritudo,respectively) accounted for less than 6% of the sediment community at any timepoint. The variation in relative abundance of these three groups did not appearto follow clear seasonal trends, and could not be readily correlated withenvironmental variables.     

Pipeliner: software to evaluate the performance of bioinformatics pipelines for next‐generation resequencing

The choice of technology and bioinformatics approach is critical in obtaining accurate and reliable information from next‐generation sequencing (NGS) experiments. An increasing number of software and methodological guidelines are being published, but deciding upon which approach and experimental design to use can depend on the particularities of the species and on the aims of the study. This leaves researchers unable to produce informed decisions on these central questions. To address these issues, we developed pipeliner – a tool to evaluate, by simulation, the performance of NGS pipelines in resequencing studies. Pipeliner provides a graphical interface allowing the users to write and test their own bioinformatics pipelines with publicly available or custom software. It computes a number of statistics summarizing the performance in SNP calling, including the recovery, sensitivity and false discovery rate for heterozygous and homozygous SNP genotypes. Pipeliner can be used to answer many practical questions, for example, for a limited amount of NGS effort, how many more reliable SNPs can be detected by doubling coverage and halving sample size or what is the false discovery rate provided by different SNP calling algorithms and options. Pipeliner thus allows researchers to carefully plan their study's sampling design and compare the suitability of alternative bioinformatics approaches for their specific study systems. Pipeliner is written in C++ and is freely available from http://github.com/brunonevado/Pipeliner .     

Co‐evolutionary dynamics between public good producers and cheats in the bacterium Pseudomonas aeruginosa

The production of beneficial public goods is common in the microbial world, and so is cheating – the exploitation of public goods by nonproducing mutants. Here, we examine co‐evolutionary dynamics between cooperators and cheats and ask whether cooperators can evolve strategies to reduce the burden of exploitation, and whether cheats in turn can improve their exploitation abilities. We evolved cooperators of the bacterium Pseudomonas aeruginosa, producing the shareable iron‐scavenging siderophore pyoverdine, together with cheats, defective in pyoverdine production but proficient in uptake. We found that cooperators managed to co‐exist with cheats in 56% of all replicates over approximately 150 generations of experimental evolution. Growth and competition assays revealed that co‐existence was fostered by a combination of general adaptions to the media and specific adaptions to the co‐evolving opponent. Phenotypic screening and whole‐genome resequencing of evolved clones confirmed this pattern, and suggest that cooperators became less exploitable by cheats because they significantly reduced their pyoverdine investment. Cheats, meanwhile, improved exploitation efficiency through mutations blocking the costly pyoverdine‐signalling pathway. Moreover, cooperators and cheats evolved reduced motility, a pattern that likely represents adaptation to laboratory conditions, but at the same time also affects social interactions by reducing strain mixing and pyoverdine sharing. Overall, we observed parallel evolution, where co‐existence of cooperators and cheats was enabled by a combination of adaptations to the abiotic and social environment and their interactions.     

Sequence Capture and Next Generation Resequencing of the MHC Region Highlights Potential Transplantation Determinants in HLA Identical Haematopoietic Stem Cell Transplantation

How cells coordinate the immune system activities is important for potentially life-saving organ or stem cell transplantations. Polymorphic immunoregulatory genes, many of them located in the human major histocompatibility complex, impact the process and assure the proper execution of tolerance-versus-activity mechanisms. In haematopoietic stem cell transplantation, on the basis of fully human leukocyte antigen (HLA)-matched donor–recipient pairs, adverse effects like graft versus leukaemia and graft versus host are observed and difficult to handle. So far, high-resolution HLA typing was performed with Sanger sequencing, but for methodological reasons information on additional immunocompetent major histocompatibility complex loci has not been revealed. Now, we have used microarray sequence capture and targeted enrichment combined with next generation pyrosequencing for 3.5 million base pair human major histocompatibility complex resequencing in a clinical transplant setting and describe 3025 variant single nucleotide polymorphisms, insertions and deletions among recipient and donor in a single sequencing experiment. Taken together, the presented data show that sequence capture and massively parallel pyrosequencing can be used as a new tool for risk assessment in the setting of allogeneic stem cell transplantation.