Search

Abstract Several studies have documented the significant impact of methodological choices in microbiome analyses. The myriad of methodological options available complicate the replication of results and generally limit the comparability of findings between independent studies that use differing techniques and measurement pipelines. Here we describe the Mosaic Standards Challenge (MSC), an international interlaboratory study designed to assess the impact of methodological variables on the results. The MSC did not prescribe methods but rather asked participating labs to analyze 7 shared reference samples (5 × human stool samples and 2 × mock communities) using their standard laboratory methods. To capture the array of methodological variables, each participating lab completed a metadata reporting sheet that included 100 different questions regarding the details of their protocol. The goal of this study was to survey the methodological landscape for microbiome metagenomic sequencing (MGS) analyses and the impact of methodological decisions on metagenomic sequencing results. A total of 44 labs participated in the MSC by submitting results (16S or WGS) along with accompanying metadata; thirty 16S rRNA gene amplicon datasets and 14 WGS datasets were collected. The inclusion of two types of reference materials (human stool and mock communities) enabled analysis of both MGS measurement variability between different protocols using the biologically-relevant stool samples, and MGS bias with respect to ground truth values using the DNA mixtures. Owing to the compositional nature of MGS measurements, analyses were conducted on the ratio of Firmicutes: Bacteroidetes allowing us to directly apply common statistical methods. The resulting analysis demonstrated that protocol choices have significant effects, including both bias of the MGS measurement associated with a particular methodological choices, as well as effects on measurement robustness as observed through the spread of results between labs making similar methodological choices. In the analysis of the DNA mock communities, MGS measurement bias was observed even when there was general consensus among the participating laboratories. This study was the result of a collaborative effort that included academic, commercial, and government labs. In addition to highlighting the impact of different methodological decisions on MGS result comparability, this work also provides insights for consideration in future microbiome measurement study design.

Common bean has two distinct domestication centers in Mesoamerica and in the Andes. The authors show that the Andean is the first gene pool successfully introduced in Europe and identify signature of pervasive introgression among gene pools and of selection for flowering underlying adaptation.

The microbial cells colonizing the human body form an ecosystem that is integral to the regulation and maintenance of human health. Elucidation of specific associations between the human microbiome and health outcomes is facilitating the development of microbiome-targeted recommendations and treatments (e.g., fecal microbiota transplant; pre-, pro-, and post-biotics) to help prevent and treat disease. However, the potential of such recommendations and treatments to improve human health has yet to be fully realized. Technological advances have led to the development and proliferation of a wide range of tools and methods to collect, store, sequence, and analyze microbiome samples. However, differences in methodology at each step in these analytic processes can lead to variability in results due to the unique biases and limitations of each component. This technical variability hampers the detection and validation of associations with small to medium effect sizes. Therefore, the American Society for Nutrition (ASN) Nutritional Microbiology Group Engaging Members (GEM), sponsored by the Institute for the Advancement of Food and Nutrition Sciences (IAFNS), hosted a satellite session on methods in nutrition and gut microbiome research to review currently available methods for microbiome research, best practices, as well as tools and standards to aid in comparability of methods and results. This manuscript summarizes the topics and research discussed at the session. Consideration of the guidelines and principles reviewed in this session will increase the accuracy, precision, and comparability of microbiome research and ultimately the understanding of the associations between the human microbiome and health.

Crop wild species are increasingly important for crop improvement. Peanut (Arachis hypogaea L.) wild relatives comprise a diverse genetic pool that is being used to broaden its narrow genetic base. Peanut is an allotetraploid species extremely susceptible to peanut root-knot nematode (PRKN) Meloidogyne arenaria. Current resistant cultivars rely on a single introgression for PRKN resistance incorporated from the wild relative Arachis cardenasii, which could be overcome as a result of the emergence of virulent nematode populations. Therefore, new sources of resistance may be needed. Near-immunity has been found in the peanut wild relative Arachis stenosperma. The two loci controlling the resistance, present on chromosomes A02 and A09, have been validated in tetraploid lines and have been shown to reduce nematode reproduction by up to 98%. To incorporate these new resistance QTL into cultivated peanut, we used a marker-assisted backcrossing approach, using PRKN A. stenosperma-derived resistant lines as donor parents. Four cycles of backcrossing were completed, and SNP assays linked to the QTL were used for foreground selection. In each backcross generation seed weight, length, and width were measured, and based on a statistical analysis we observed that only one generation of backcrossing was required to recover the elite peanut’s seed size. A populating of 271 BC3F1 lines was genome-wide genotyped to characterize the introgressions across the genome. Phenotypic information for leaf spot incidence and domestication traits (seed size, fertility, plant architecture, and flower color) were recorded. Correlations between the wild introgressions in different chromosomes and the phenotypic data allowed us to identify candidate regions controlling these domestication traits. Finally, PRKN resistance was validated in BC3F3 lines. We observed that the QTL in A02 and/or large introgression in A09 are needed for resistance. This present work represents an important step toward the development of new high-yielding and nematode-resistant peanut cultivars.

Although seed and pod traits are important for peanut breeding, little is known about the inheritance of these traits. A recombinant inbred line (RIL) population of 156 lines from a cross of Tifrunner x NC 3033 was genotyped with the Axiom_Arachis1 SNP array and SSRs to generate a genetic map composed of 1524 markers in 29 linkage groups (LG). The genetic positions of markers were compared with their physical positions on the peanut genome to confirm the validity of the linkage map and explore the distribution of recombination and potential chromosomal rearrangements. This linkage map was then used to identify Quantitative Trait Loci (QTL) for seed and pod traits that were phenotyped over three consecutive years for the purpose of developing trait-associated markers for breeding. Forty-nine QTL were identified in 14 LG for seed size index, kernel percentage, seed weight, pod weight, single-kernel, double-kernel, pod area and pod density. Twenty QTL demonstrated phenotypic variance explained (PVE) greater than 10% and eight more than 20%. Of note, seven of the eight major QTL for pod area, pod weight and seed weight (PVE >20% variance) were attributed to NC 3033 and located in a single linkage group, LG B06_1. In contrast, the most consistent QTL for kernel percentage were located on A07/B07 and derived from Tifrunner.

Abstract Soybean [Glycinemax (L.) Merr.] maturity determines the growing region of a given soybean variety and is a primary factor in yield and other agronomic traits. The objectives of this research were to identify the quantitative trait loci (QTL) associated with maturity groups (MGs) and determine the genetic control of soybean maturity in each MG. Using data from 16,879 soybean accessions, genome‐wide association (GWA) analyses were conducted for each paired MG and across MGs 000 through IX. Genome‐wide association analyses were also performed using 184 genotypes (MGs V–IX) with days to flowering (DTF) and maturity (DTM) collected in the field. A total of 58 QTL were identified to be significantly associated with MGs in individual GWAs, which included 12 reported maturity loci and two stem termination genes. Genome‐wide associations across MGs 000–IX detected a total of 103 QTL and confirmed 54 QTL identified in the individual GWAs. Of significant loci identified, qMG‐5.2 had effects on the highest number (9) of MGs, followed by E2, E3, Dt2, qMG‐15.5, E1, qMG‐13.1, qMG‐7.1, and qMG‐16.1, which affected five to seven MGs. A high number of genetic loci (8–25) that affected MGs 0–V were observed. Stem termination genes Dt1 and Dt2 mainly had significant allele variation in MGs II–V. Genome‐wide associations for DTF, DTM, and reproductive period (RP) in the diversity panel confirmed 15 QTL, of which seven were observed in MGs V–IX. The results generated can help soybean breeders manipulate the maturity loci for genetic improvement of soybean yield.

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed. Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Abstract Background Plant domestication provides a unique model to study genome evolution. Many studies have been conducted to examine genes, genetic diversity, genome structure, and epigenome changes associated with domestication. Interestingly, domesticated accessions have significantly higher [A] and [T] values across genome-wide polymorphic sites than accessions sampled from the corresponding progenitor species. However, the relative contributions of different genomic regions to this genome divergence pattern and underlying mechanisms have not been well characterized. Results Here, we investigate the genome-wide base-composition patterns by analyzing millions of SNPs segregating among 100 accessions from a teosinte-maize comparison set and among 302 accessions from a wild-domesticated soybean comparison set. We show that non-genic part of the genome has a greater contribution than genic SNPs to the [AT]-increase observed between wild and domesticated accessions in maize and soybean. The separation between wild and domesticated accessions in [AT] values is significantly enlarged in non-genic and pericentromeric regions. Motif frequency and sequence context analyses show the motifs (PyCG) related to solar-UV signature are enriched in these regions, particularly when they are methylated. Additional analysis using population-private SNPs also implicates the role of these motifs in relatively recent mutations. With base-composition across polymorphic sites as a genome phenotype, genome scans identify a set of putative candidate genes involved in UV damage repair pathways. Conclusions The [AT]-increase is more pronounced in genomic regions that are non-genic, pericentromeric, transposable elements; methylated; and with low recombination. Our findings establish important links among UV radiation, mutation, DNA repair, methylation, and genome evolution.

Abstract Machine learning has demonstrated potential in analyzing large, complex biological data. In practice, however, biological information is required in addition to machine learning for successful application.

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed. Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Regulatory changes include divergence in both cis-elements and trans-factors, which play roles in organismal evolution. Whole genome duplications (WGD) followed by diploidization are a recurrent feature in the evolutionary history of angiosperms. Prior studies have shown that duplicated genes have different evolutionary fates due to variable selection constraints and results in genomic compositions with hallmarks of paleopolyploidy. The recent sequential WGDs and post-WGD evolution in the common ancestor of cultivated soybean (Glycine max) and wild soybean (Glycine soja), together with other models of gene duplication, have resulted in a highly duplicated genome. In this study, we investigated the transcriptional changes in G. soja and G. max. We identified a sizable proportion of interspecific differentially expressed genes (DEGs) and found parental expression level dominance of G. max in their F1 hybrids. By classifying genes into different regulatory divergence types, we found the trans-regulatory changes played a predominant role in transcriptional divergence between wild and cultivated soybean. The same gene ontology (GO) and protein family (Pfam) terms were found to be over-represented in DEGs and genes of cis-only between JY47 and GS, suggesting the substantial contribution of cis-regulatory divergences to the evolution of wild and cultivated soybeans. By further dissecting genes into five different duplication modes, we found genes in different duplication modes tend to accumulate different types of regulatory differences. A relatively higher proportion of cis-only regulatory divergences was detected in singleton, dispersed, proximal, and tandem duplicates than WGD duplicates and genome-wide level, which is in line with the prediction of gene balance hypothesis for the differential fates of duplicated genes post-WGD. The numbers of cis-only and trans-only regulated genes were similar for singletons, whereas there were more genes of trans-only than cis-only in the rest duplication types, especially in WGD in which there were two times more trans-only genes than that in cis-only type. Tandem duplicates showed the highest proportion of trans-only genes probably due to some special features of this class. In summary, our results demonstrate that genes in different duplication modes have different fates in transcriptional evolution underpinned by cis- or trans-regulatory divergences in soybean and likely in other paleopolyploid higher organisms.

Cronobacter species are opportunistic pathogens capable of causing life-threatening infections in humans, with serious complications arising in neonates, infants, immuno-compromised individuals, and elderly adults. The genus is comprised of seven species: Cronobacter sakazakii, Cronobacter malonaticus, Cronobacter turicensis, Cronobacter muytjensii, Cronobacter dublinensis, Cronobacter universalis, and Cronobacter condimenti. Despite a multiplicity of genomic data for the genus, little is known about likely transmission vectors. Using DNA microarray analysis, in parallel with whole genome sequencing, and targeted PCR analyses, the total gene content of two C. malonaticus, three C. turicensis, and 14 C. sakazaki isolated from various filth flies was assessed. Phylogenetic relatedness among these and other strains obtained during surveillance and outbreak investigations were comparatively assessed. Specifically, microarray analysis (MA) demonstrated its utility to cluster strains according to species-specific and sequence type (ST) phylogenetic relatedness, and that the fly strains clustered among strains obtained from clinical, food and environmental sources from United States, Europe, and Southeast Asia. This combinatorial approach was useful in data mining for virulence factor genes, and phage genes and gene clusters. In addition, results of plasmidotyping were in agreement with the species identity for each strain as determined by species-specific PCR assays, MA, and whole genome sequencing. Microarray and BLAST analyses of Cronobacter fly sequence datasets were corroborative and showed that the presence and absence of virulence factors followed species and ST evolutionary lines even though such genes were orthologous. Additionally, zebrafish infectivity studies showed that these pathotypes were as virulent to zebrafish embryos as other clinical strains. In summary, these findings support a striking phylogeny amongst fly, clinical, and surveillance strains isolated during 2010–2015, suggesting that flies are capable vectors for transmission of virulent Cronobacter spp.; they continue to circulate among United States and European populations, environments, and that this “pattern of circulation” has continued over decades.

In a rapidly growing global probiotic market, end-users have difficulty distinguishing between high quality and poor quality products. This ambiguity threatens the trust consumers and healthcare providers have in probiotic products. To address this problem, we recommend that companies undergo third-party evaluations to certify probiotic quality and label accuracy. In order to communicate about product quality to end-users, indication of certification on product labels is helpful, although not all manufacturers choose to use this approach. Herein we discuss: third-party certification, the process of setting standards for identity, purity, and quantification of probiotics; some emerging methodologies useful for quality assessment; and some technical challenges unique to managing quality of live microbial products. This review provides insights of an Expert Panel engaged in this process and aims to update the reader on relevant current scientific methodologies. Establishing validated methodologies for all aspects of quality assessment is an essential component of this process and can be facilitated by established organizations, such as United States Pharmacopeia. Emerging methodologies including whole genome sequencing and flow cytometry are poised to play important roles in these processes.

Polyploidization has played a prominent role in the evolutionary history of plants. Two recent and sequential allopolyploidization events have resulted in the formation of wheat species with different ploidies, and which provide a model to study the effects of polyploidization on the evolution of gene expression. In this study, we identified differentially expressed genes (DEGs) between four BBAA tetraploid wheats of three different ploidy backgrounds. DEGs were found to be unevenly distributed among functional categories and duplication modes. We observed more DEGs in the extracted tetraploid wheat (ETW) than in natural tetraploid wheats (TD and TTR13) as compared to a synthetic tetraploid (AT2). Furthermore, DEGs showed higher Ka/Ks ratios than those that did not show expression changes (non-DEGs) between genotypes, indicating DEGs and non-DEGs experienced different selection pressures. For A-B homeolog pairs with DEGs, most of them had only one differentially expressed copy, however, when both copies of a homeolog pair were DEGs, the A and B copies were more likely to be regulated to the same direction. Our results suggest that both cis- and inter-subgenome trans-regulatory changes are important drivers in the evolution of homeologous gene expression in polyploid wheat, with ploidy playing a significant role in the process.

Transposons are ubiquitous genomic components that play pivotal roles in plant gene and genome evolution. We analyzed two genome sequences of common bean (Phaseolus vulgaris) and identified a new CACTA transposon family named pvCACTA1. The family is extremely abundant, as more than 12,000 pvCACTA1 elements were found. To our knowledge, this is the most abundant CACTA family reported thus far. The computational and fluorescence in situ hybridization (FISH) analyses indicated that the pvCACTA1 elements were concentrated in terminal regions of chromosomes and frequently generated AT-rich 3 bp target site duplications (TSD, WWW, W is A or T). Comparative analysis of the common bean genomes from two domesticated genetic pools revealed that new insertions or excisions of pvCACTA1 elements occurred after the divergence of the two common beans, and some of the polymorphic elements likely resulted in variation in gene sequences. pvCACTA1 elements were detected in related species but not outside the Phaseolus genus. We calculated the molecular evolutionary rate of pvCACTA1 transposons using orthologous elements that indicated that most transposition events likely occurred before the divergence of the two gene pools. These results reveal unique features and evolution of this new transposon family in the common bean genome.

Fluorescence in situ hybridization (FISH)-based karyotyping is a powerful cytogenetics tool to study chromosome organization, behavior, and chromosome evolution. Here, we developed a FISH-based karyotyping system using a probe mixture comprised of centromeric and subtelomeric satellite repeats, 5S rDNA, and chromosome-specific BAC clones in common bean, which enables one to unambiguously distinguish all 11 chromosome pairs. Furthermore, we applied the karyotyping system to several wild relatives and landraces of common bean from two distinct gene pools, as well as other related Phaseolus species, to investigate repeat evolution in the genus Phaseolus. Comparison of karyotype maps within common bean indicates that chromosomal distribution of the centromeric and subtelomeric satellite repeats is stable, whereas the copy number of the repeats was variable, indicating rapid amplification/reduction of the repeats in specific genomic regions. In Phaseolus species that diverged approximately 2–4 million yr ago, copy numbers of centromeric repeats were largely reduced or diverged, and chromosomal distributions have changed, suggesting rapid evolution of centromeric repeats. We also detected variation in the distribution pattern of subtelomeric repeats in Phaseolus species. The FISH-based karyotyping system revealed that satellite repeats are actively and rapidly evolving, forming genomic features unique to individual common bean accessions and Phaseolus species.

Rice bean (Vigna umbellata Thunb.), a warm-season annual legume, is grown in Asia mainly for dried grain or fodder and plays an important role in human and animal nutrition because the grains are rich in protein and some essential fatty acids and minerals. With the aim of expediting the genetic improvement of rice bean, we initiated a project to develop genomic resources and tools for molecular breeding in this little-known but important crop. Here we report the construction of an SSR-enriched genomic library from DNA extracted from pooled young leaf tissues of 22 rice bean genotypes and developing SSR markers. In 433,562 reads generated by a Roche 454 GS-FLX sequencer, we identified 261,458 SSRs, of which 48.8% were of compound form. Dinucleotide repeats were predominant with an absolute proportion of 81.6%, followed by trinucleotides (17.8%). Other types together accounted for 0.6%. The motif AC/GT accounted for 77.7% of the total, followed by AAG/CTT (14.3%), and all others accounted for 12.0%. Among the flanking sequences, 2928 matched putative genes or gene models in the protein database of Arabidopsis thaliana, corresponding with 608 non-redundant Gene Ontology terms. Of these sequences, 11.2% were involved in cellular components, 24.2% were involved molecular functions, and 64.6% were associated with biological processes. Based on homolog analysis, 1595 flanking sequences were similar to mung bean and 500 to common bean genomic sequences. Comparative mapping was conducted using 350 sequences homologous to both mung bean and common bean sequences. Finally, a set of primer pairs were designed, and a validation test showed that 58 of 220 new primers can be used in rice bean and 53 can be transferred to mung bean. However, only 11 were polymorphic when tested on 32 rice bean varieties. We propose that this study lays the groundwork for developing novel SSR markers and will enhance the mapping of qualitative and quantitative traits and marker-assisted selection in rice bean and other Vigna species.

Abstract Epigenetic profiling in diploid, allopolyploid, and domesticated cotton shows that despite most DNA methylation being conserved and stably inherited, alterations likely due to hybridization and domestication affect gene expression.

Plant genome evolution can be very complex and challenging to describe, even within a genus. Mechanisms that underlie genome variation are complex and can include whole-genome duplications, gene duplication and/or loss, and, importantly, multiple chromosomal rearrangements. Lupins (Lupinus) diverged from other legumes approximately 60 mya. In contrast to New World lupins, Old World lupins show high variability not only for chromosome numbers (2n = 32–52), but also for the basic chromosome number (x = 5–9, 13) and genome size. The evolutionary basis that underlies the karyotype evolution in lupins remains unknown, as it has so far been impossible to identify individual chromosomes. To shed light on chromosome changes and evolution, we used comparative chromosome mapping among 11 Old World lupins, with Lupinus angustifolius as the reference species. We applied set of L. angustifolius-derived bacterial artificial chromosome clones for fluorescence in situ hybridization. We demonstrate that chromosome variations in the species analyzed might have arisen from multiple changes in chromosome structure and number. We hypothesize about lupin karyotype evolution through polyploidy and subsequent aneuploidy. Additionally, we have established a cytogenomic map of L. angustifolius along with chromosome markers that can be used for related species to further improve comparative studies of crops and wild lupins.

Common bean (Phaseolus vulgaris L.) is an important legume crop worldwide. However, abiotic and biotic stress limits bean yields to

Cytogenetics played an essential role in studies of chromosome structure, behavior, and evolution in numerous plant species. The advent of molecular cytogenetics combined with recent development of genomic resources has ushered in a new era of chromosome studies that have greatly advanced our knowledge of karyotypic diversity, genome and chromosome organization, and chromosomal evolution in legumes. This review summarizes some of the achievements of cytogenetic studies in legumes in the Phaseoloid clade, which includes several important legume crops such as common bean ( L.), cowpea [ (L.) Walp.], soybean [ (L.) Merr.], and pigeonpea [ (L.) Huth]. In the Phaseoloid clade, karyotypes are mostly stable. There are, however, several species with extensive chromosomal changes. Fluorescence in situ hybridization has been useful to reveal chromosomal structure by physically mapping transposons, satellite repeats, ribosomal DNA genes, and bacterial artificial chromosome clones onto chromosomes. Polytene chromosomes, which are much longer than the mitotic chromosomes, have been successfully found and used in cytogenetic studies in some and species. Molecular cytogenetics will continue to be an important tool in legume genetics and genomics, and we discuss future applications of molecular cytogenetics to better understand chromosome and genome structure and evolution in legumes.

A small fast neutron mutant population has been established from Phaseolus vulgaris cv. Red Hawk. We leveraged the available P. vulgaris genome sequence and high throughput next generation DNA sequencing to examine the genomic structure of five Phaseolus vulgaris cv. Red Hawk fast neutron mutants with striking visual phenotypes. Analysis of these genomes identified three classes of structural variation; between cultivar variation, natural variation within the fast neutron mutant population, and fast neutron induced mutagenesis. Our analyses focused on the latter two classes. We identified 23 large deletions (>40 bp) common to multiple individuals, illustrating residual heterogeneity and regions of structural variation within the common bean cv. Red Hawk. An additional 18 large deletions were identified in individual mutant plants. These deletions, ranging in size from 40 bp to 43,000 bp, are potentially the result of fast neutron mutagenesis. Six of the 18 deletions lie near or within gene coding regions, identifying potential candidate genes causing the mutant phenotype.