A simple and effective method to enrich endogenous DNA from mammalian faeces

Utilization of faeces has long been a popular approach for genetic and ecological studies of wildlife. However, the success of molecular marker genotyping and genome resequencing is often unpredictable due to insufficient enrichment of endogenous DNA in the total faecal DNA that is dominated by bacterial DNA. Here, we report a simple and cheap method named PEERS to predominantly lyse animal cells over bacteria by using sodium dodecyl sulphate so as to discharge endogenous DNA into liquid phase before bacterial DNA. By brief centrifugation, total DNA with enriched endogenous fraction can be extracted from the supernatant using routine methods. Our assessments showed that the endogenous DNA extracted by PEERS was significantly enriched for various types of faeces from different species, preservation time and conditions. It significantly improves the genotyping correctness and efficiency of genome resequencing with the total additional cost of $ 0.1 and a short incubation step to treat a faecal sample. We also provide methods to assess the enrichment efficiency of mitochondrial and nuclear DNA and models to predict the usability of faecal DNA for genotyping of short tandem repeat, single-nucleotide polymorphism and whole-genome resequencing.

Study on the origin of the Baise horse based on whole-genome resequencing

The Baise horse, an indigenous horse breed mainly distributed in the Baise region of Guangxi province in southwest China, has a long history as draft animal. However, there is a lack of research regarding the origin and ancestral composition of the Baise horse. In this study, whole-genome resequencing data from 236 horses of seven Chinese indigenous horse breeds, five foreign horse breeds, and four Przewalski's horses were used to investigate the relationships between the Baise horse and other horse breeds. The results showed that foreign horse breeds had no significant impact on the formation of the Baise horse. The two southwestern horse populations, the Debao pony and the Jinjiang horse, exhibit the closest genetic affinity with the Baise horse. This is consistent with their adjacent geographical distribution. Analysis of the migration route revealed a gene flow from the Chakouyi horse into the Baise horse. In summary, our results confirm that the formation of the Baise horse did not involve participation from foreign breeds. Geographical distance emerges as a crucial factor in determining the genetic relationships with the Baise horse. Gene flows of indigenous horse breeds along ancient routes of trade activities had played a role in the formation of the Baise horse.

Variation in genomic vulnerability to climate change across temperate populations of eelgrass (Zostera marina)

A global decline in seagrass populations has led to renewed calls for their conservation as important providers of biogenic and foraging habitat, shoreline stabilization and carbon storage. Eelgrass (Zostera marina) occupies the largest geographic range among seagrass species spanning a commensurately broad spectrum of environmental conditions. In Canada, eelgrass is managed as a single phylogroup despite occurring across three oceans and a range of ocean temperatures and salinity gradients. Previous research has focused on applying relatively few markers to reveal population structure of eelgrass, whereas a whole-genome approach is warranted to investigate cryptic structure among populations inhabiting different ocean basins and localized environmental conditions. We used a pooled whole-genome re-sequencing approach to characterize population structure, gene flow and environmental associations of 23 eelgrass populations ranging from the Northeast United States to Atlantic, subarctic and Pacific Canada. We identified over 500,000 SNPs, which when mapped to a chromosome-level genome assembly revealed six broad clades of eelgrass across the study area, with pairwise F ST ranging from 0 among neighbouring populations to 0.54 between Pacific and Atlantic coasts. Genetic diversity was highest in the Pacific and lowest in the subarctic, consistent with colonization of the Arctic and Atlantic oceans from the Pacific less than 300 kya. Using redundancy analyses and two climate change projection scenarios, we found that subarctic populations are predicted to be potentially more vulnerable to climate change through genomic offset predictions. Conservation planning in Canada should thus ensure that representative populations from each identified clade are included within a national network so that latent genetic diversity is protected, and gene flow is maintained. Northern populations, in particular, may require additional mitigation measures given their potential susceptibility to a rapidly changing climate.

Whole-genome resequencing facilitates the development of a 50K single nucleotide polymorphism genotyping array for Scots pine (Pinus sylvestris L.) and its transferability to other pine species

Scots pine (Pinus sylvestris L.) is one of the most widespread and economically important conifer species in the world. Applications like genomic selection and association studies, which could help accelerate breeding cycles, are challenging in Scots pine because of its large and repetitive genome. For this reason, genotyping tools for conifer species, and in particular for Scots pine, are commonly based on transcribed regions of the genome. In this article, we present the Axiom Psyl50K array, the first single nucleotide polymorphism (SNP) genotyping array for Scots pine based on whole-genome resequencing, that represents both genic and intergenic regions. This array was designed following a two-step procedure: first, 192 trees were sequenced, and a 430K SNP screening array was constructed. Then, 480 samples, including haploid megagametophytes, full-sib family trios, breeding population, and range-wide individuals from across Eurasia were genotyped with the screening array. The best 50K SNPs were selected based on quality, replicability, distribution across the draft genome assembly, balance between genic and intergenic regions, and genotype-environment and genotype-phenotype associations. Of the final 49 877 probes tiled in the array, 20 372 (40.84%) occur inside gene models, while the rest lie in intergenic regions. We also show that the Psyl50K array can yield enough high-confidence SNPs for genetic studies in pine species from North America and Eurasia. This new genotyping tool will be a valuable resource for high-throughput fundamental and applied research of Scots pine and other pine species.

Gut microbiota facilitate adaptation of invasive moths to new host plants

Gut microbiota are important in the adaptation of phytophagous insects to their plant hosts. However, the interaction between gut microbiomes and pioneering populations of invasive insects during their adaptation to new hosts, particularly in the initial phases of invasion, has been less studied. We studied the contribution of the gut microbiome to host adaptation in the globally recognized invasive pest, Hyphantria cunea, as it expands its range into southern China. The southern population of H. cunea shows effective adaptation to Metasequoia glyptostroboides and exhibits greater larval survival on Metasequoia than the original population. Genome resequencing revealed no significant differences in functions related to host adaptation between the two populations. The compatibility between southern H. cunea populations and M. glyptostroboides revealed a correlation between the abundance of several gut bacteria genera (Bacteroides, Blautia, and Coprococcus) and H. cunea survival. Transplanting the larval gut microbiome from southern to northern populations enhanced the adaptability of the latter to the previously unsuitable plant M. glyptostroboides. This research provides evidence that the gut microbiome of pioneering populations can enhance the compatibility of invasive pests to new hosts and enable more rapid adaptation to new habitats.

Molecular basis of the phenotypic variants arising from a Pseudoalteromonas lipolytica mutator

Bacterial deficiencies in the DNA repair system can produce mutator strains that promote adaptive microevolution. However, the role of mutator strains in marine Pseudoalteromonas, capable of generating various gain-of-function genetic variants within biofilms, remains largely unknown. In this study, inactivation of mutS in Pseudoalteromonas lipolytica conferred an approximately 100-fold increased resistance to various antibiotics, including ciprofloxacin, rifampicin and aminoglycoside. Furthermore, the mutator of P. lipolytica generated variants that displayed enhanced biofilm formation but reduced swimming motility, indicating a high phenotypic diversity within the ΔmutS population. Additionally, we observed a significant production rate of approximately 50 % for the translucent variants, which play important roles in biofilm formation, when the ΔmutS strain was cultured on agar plates or under shaking conditions. Using whole-genome deep-sequencing combined with genetic manipulation, we demonstrated that point mutations in AT00_17115 within the capsular biosynthesis cluster were responsible for the generation of translucent variants in the ΔmutS subpopulation, while mutations in flagellar genes fliI and flgP led to a decrease in swimming motility. Collectively, this study reveals a specific mutator-driven evolution in P. lipolytica, characterized by substantial genetic and phenotypic diversification, thereby offering a reservoir of genetic attributes associated with microbial fitness.

Genome-Wide Scanning Loci and Differentially Expressed Gene Analysis Unveils the Molecular Mechanism of Chlorantraniliprole Resistance in Spodoptera frugiperda

Chlorantraniliprole has been widely used to controlSpodoptera frugiperda, but it has led to the development of chlorantraniliprole resistance. Multiomics analysis of strains with two extreme traits helps to elucidate the complex mechanisms involved. Herein, following genome resequencing and application of the Euclidean distance algorithm, 550 genes within a 16.20-Mb-linked region were identified from chlorantraniliprole-resistant (Ch-R) and chlorantraniliprole-susceptible (Ch-Sus) strains. Using transcriptome sequencing, 2066 differentially expressed genes were identified between Ch-R and Ch-Sus strains. Through association analysis, three glutathione S-transferase family genes and four trehalose transporter genes were selected for functional verification. Notably, SfGSTD1 had the strongest binding ability with chlorantraniliprole and is responsible for chlorantraniliprole tolerance. The Ch-R strain also increased the intracellular trehalose content by upregulating the transcription of SfTret1, thereby contributing to chlorantraniliprole resistance. These findings provide a new perspective to reveal the mechanism of resistance of agricultural pests to insecticides.

A kiss of deep homology: partial convergence in the genomic basis of hypertrophied lips in cichlid fish and human cleft lip

The genomic loci generating both adaptive and maladaptive variation could be surprisingly predictable in deeply homologous vertebrate structures like lips. Variation in highly conserved vertebrate traits such as jaws and teeth in organisms as evolutionarily disparate as teleost fishes and mammals are known to be structured by the same genes. Likewise, hypertrophied lips that have evolved repeatedly in Neotropical and African cichlid fish lineages could share unexpectedly similar genetic bases themselves and even provide surprising insight into the loci underlying human craniofacial anomalies. To isolate the genomic regions underlying adaptive divergence in hypertrophied lips, we first employed genome wide associations (GWA) in several species of African cichlids from Lake Malawi. Then, we tested if these GWA regions were shared through hybridization with another Lake Malawi cichlid lineage that has evolved hypertrophied lips seemingly in parallel. Overall, introgression among hypertrophied lip lineages appeared limited. Among our Malawi GWA regions, one contained the gene kcnj16 that has been implicated in the convergently evolved hypertrophied lips in Central American Midas cichlids that diverged from the Malawi radiation over 50 million years ago. The Malawi hypertrophied lip GWA regions also contained several additional genes that cause human lip-associated birth defects. Cichlid fishes are becoming prominent examples of replicated genomic architecture underlying trait convergence and are increasingly providing insight into human craniofacial anomalies such as cleft lip.

The Genomic Variation and Differentially Expressed Genes on the 6P Chromosomes in Wheat-Agropyron cristatum Addition Lines 5113 and II-30-5 Confer Different Desirable Traits

Wild relatives of wheat are essential gene pools for broadening the genetic basis of wheat. Chromosome rearrangements and genomic variation in alien chromosomes are widespread. Knowledge of the genetic variation between alien homologous chromosomes is valuable for discovering and utilizing alien genes. In this study, we found that 5113 and II-30-5, two wheat-A. cristatum 6P addition lines, exhibited considerable differences in heading date, grain number per spike, and grain weight. Genome resequencing and transcriptome analysis revealed significant differences in the 6P chromosomes of the two addition lines, including 143,511 single-nucleotide polymorphisms, 62,103 insertion/deletion polymorphisms, and 757 differentially expressed genes. Intriguingly, genomic variations were mainly distributed in the middle of the chromosome arms and the proximal centromere region. GO and KEGG analyses of the variant genes and differentially expressed genes showed the enrichment of genes involved in the circadian rhythm, carbon metabolism, carbon fixation, and lipid metabolism, suggesting that the differential genes on the 6P chromosome are closely related to the phenotypic differences. For example, the photosynthesis-related genes PsbA, PsbT, and YCF48 were upregulated in II-30-5 compared with 5113. ACS and FabG are related to carbon fixation and fatty acid biosynthesis, respectively, and both carried modification variations and were upregulated in 5113 relative to II-30-5. Therefore, this study provides important guidance for cloning desirable genes from alien homologous chromosomes and for their effective utilization in wheat improvement.

Screening of candidate genes related to differences in growth and development between Chinese indigenous and Western pig breeds

Neijiang (NJ) and Yacha (YC) are two indigenous pig breeds in the Sichuan basin of China, displaying higher resistance to diseases, lower lean ratio, and slower growth rate than the commercial Western pig breed Yorkshire (YS). The molecular mechanisms underlying the differences in growth and development between these pig breeds are still unknown. In the present study, five pigs from NJ, YC, and YS breeds were subjected to the whole genome resequencing, and then the differential single-nucleotide polymorphisms (SNPs) were screened using a 10-kb window sliding in 1-kb step using the F_st method. Finally, 48,924, 48,543, and 46,228 nonsynonymous single-nucleotide polymorphism loci (nsSNPs) were identified between NJ and YS, NJ and YC, and YC and YS, which highly or moderately affected 2,490, 800, and 444 genes, respectively. Moreover, three nsSNPs were detected in the genes of acetyl-CoA acetyltransferase 1 (ACAT1) insulin-like growth factor 2 receptor (IGF2R), insulin-like growth factor 2 and mRNA-binding protein 3 (IGF2BP3), which potentially affected the transformation of acetyl-CoA to acetoacetyl-CoA and the normal functions of the insulin signaling pathways. Moreover, serous determinations revealed significantly lower acetyl-CoA content in YC than in YS, supporting that ACAT1 might be a reason explaining the differences in growth and development between YC and YS breeds. Contents of phosphatidylcholine (PC) and phosphatidic acid (PA) significantly differed between the pig breeds, suggesting that glycerophospholipid metabolism might be another reason for the differences between Chinese and Western pig breeds. Overall, these results might contribute basic information to understand the genetic differences determining the phenotypical traits in pigs.

Phylogenomic insights into the origin and evolutionary history of evergreen broadleaved forests in East Asia under Cenozoic climate change

The evergreen versus deciduous leaf habit is an important functional trait for adaptation of forest trees and has been hypothesized to be related to the evolutionary processes of the component species under paleoclimatic change, and potentially reflected in the dynamic history of evergreen broadleaved forests (EBLFs) in East Asia. However, knowledge about the shift of evergreen versus deciduous leaf with the impact of paleoclimatic change using genomic data remains rare. Here, we focus on the Litsea complex (Lauraceae), a key lineage with dominant species of EBLFs, to gain insights into how evergreen versus deciduous trait shifted, providing insights into the origin and historical dynamics of EBLFs in East Asia under Cenozoic climate change. We reconstructed a robust phylogeny of the Litsea complex using genome-wide single-nucleotide variants (SNVs) with eight clades resolved. Fossil-calibrated analyses, diversification rate shifts, ancestral habit, ecological niche modelling and climate niche reconstruction were employed to estimate its origin and diversification pattern. Taking into account studies on other plant lineages dominating EBLFs of East Asia, it was revealed that the prototype of EBLFs in East Asia probably emerged in the Early Eocene (55-50 million years ago [Ma]), facilitated by the greenhouse warming. As a response to the cooling and drying climate in the Middle to Late Eocene (48-38 Ma), deciduous habits were evolved in the dominant lineages of the EBLFs in East Asia. Up to the Early Miocene (23 Ma), the prevailing of East Asian monsoon increased the extreme seasonal precipitation and accelerated the emergence of evergreen habits of the dominant lineages, and ultimately shaped the vegetation resembling that of today.

Resequencing of 410 Sesame Accessions Identifies SINST1 as the Major Underlying Gene for Lignans Variation

Sesame is a promising oilseed crop that produces specific lignans of clinical importance. Hence, a molecular description of the regulatory mechanisms of lignan biosynthesis is essential for crop improvement. Here, we resequence 410 sesame accessions and identify 5.38 and 1.16 million SNPs (single nucleotide polymorphisms) and InDels, respectively. Population genomic analyses reveal that sesame has evolved a geographic pattern categorized into northern (NC), middle (MC), and southern (SC) groups, with potential origin in the southern region and subsequent introduction to the other regions. Selective sweeps analysis uncovers 120 and 75 significant selected genomic regions in MC and NC groups, respectively. By screening these genomic regions, we unveiled 184 common genes positively selected in these subpopulations for exploitation in sesame improvement. Genome-wide association study identifies 17 and 72 SNP loci for sesamin and sesamolin variation, respectively, and 11 candidate causative genes. The major pleiotropic SNP^C/A locus for lignans variation is located in the exon of the gene SiNST1. Further analyses revealed that this locus was positively selected in higher lignan content sesame accessions, and the "C" allele is favorable for a higher accumulation of lignans. Overexpression of SiNST1^C in sesame hairy roots significantly up-regulated the expression of SiMYB58, SiMYB209, SiMYB134, SiMYB276, and most of the monolignol biosynthetic genes. Consequently, the lignans content was significantly increased, and the lignin content was slightly increased. Our findings provide insights into lignans and lignin regulation in sesame and will facilitate molecular breeding of elite varieties and marker-traits association studies.

Molecular evidence for adaptive evolution of drought tolerance in wild cereals

The considerable drought tolerance of wild cereal crop progenitors has diminished during domestication in the pursuit of higher productivity. Regaining this trait in cereal crops is essential for global food security but requires novel genetic insight. Here, we assessed the molecular evidence for natural variation of drought tolerance in wild barley (Hordeum spontaneum), wild emmer wheat (Triticum dicoccoides), and Brachypodium species collected from dry and moist habitats at Evolution Canyon, Israel (ECI). We report that prevailing moist vs dry conditions have differentially shaped the stomatal and photosynthetic traits of these wild cereals in their respective habitats. We present the genomic and transcriptomic evidence accounting for differences, including co-expression gene modules, correlated with physiological traits, and selective sweeps, driven by the xeric site conditions on the African Slope (AS) at ECI. Co-expression gene module 'circadian rhythm' was linked to significant drought-induced delay in flowering time in Brachypodium stacei genotypes. African Slope-specific differentially expressed genes are important in barley drought tolerance, verified by silencing Disease-Related Nonspecific Lipid Transfer 1 (DRN1), Nonphotochemical Quenching 4 (NPQ4), and Brassinosteroid-Responsive Ring-H1 (BRH1). Our results provide new genetic information for the breeding of resilient wheat and barley in a changing global climate with increasingly frequent drought events.

Study of genetic variation and its association with tensile strength among bamboo species through whole genome resequencing

Moso bamboo (Phyllostachys edulis) is a versatile plant species that is widely used as a construction material by many low-income countries due to the lack of major construction materials such as steel and reinforced concrete. It is also widely used in China. Bamboo is an economically sustainable material that behaves exceptionally in natural disasters such as earthquakes and it can offer viable solutions for contemporary engineering challenges. Despite bamboo's potential in the engineering sector, biological features such as its long generation time, its large genome size, and its polyploidy are constraining factors for genetic and genomic studies that potentially can assist the breeding efforts. This study re-sequenced 8 Phyllostachys species and 18 natural accessions of Ph. edulis, generating a large set of functionally annotated molecular markers (Single Nucleotide Polymorphisms (SNPs) and InDels) providing key genomic resource information. Moreover, all this genomic information was used to carry out a preliminary genome-wide association analysis and several candidate genes were identified to be correlated with a mechanical property that is of high interest to structural engineers: its tensile strength normal to its fibers (i.e., splitting).

High levels of inbreeding with spatial and host-associated structure in lice of an endangered freshwater seal

Host-specialist parasites of endangered large vertebrates are in many cases more endangered than their hosts. In particular, low host population densities and reduced among-host transmission rates are expected to lead to inbreeding within parasite infrapopulations living on single host individuals. Furthermore, spatial population structures of directly-transmitted parasites should be concordant with those of their hosts. Using population genomic approaches, we investigated inbreeding and population structure in a host-specialist seal louse (Echinophthirius horridus) infesting the Saimaa ringed seal (Phoca hispida saimensis), which is endemic to Lake Saimaa in Finland, and is one of the most endangered pinnipeds in the world. We conducted genome resequencing of pairs of lice collected from 18 individual Saimaa ringed seals throughout the Lake Saimaa complex. Our analyses showed high genetic similarity and inbreeding between lice inhabiting the same individual seal host, indicating low among-host transmission rates. Across the lake, genetic differentiation among individual lice was correlated with their geographic distance, and assignment analyses revealed a marked break in the genetic variation of the lice in the middle of the lake, indicating substantial population structure. These findings indicate that movements of Saimaa ringed seals across the main breeding areas of the fragmented Lake Saimaa complex may in fact be more restricted than suggested by previous population-genetic analyses of the seals themselves.

High-resolution bin-based linkage mapping uncovers the genetic architecture and heterosis-related loci of plant height in indica-japonica derived populations

Plant height (PH) is an important trait affecting the plant architecture, seed yield, and harvest index. However, the molecular mechanisms underlying PH heterosis remain unclear. In addition, useful PH-related genes must be urgently identified to facilitate ideal plant architecture breeding in rice (Oryza sativa L.). In the present study, to explore rice quantitative trait loci (QTLs) and heterosis-related loci of PH in rice, we developed a high-generation (>F₁₅ ) population of 272 recombinant inbred lines (RIL) from a cross of two elite varieties, Luohui 9 (indica/xian) × RPY geng (japonica/geng), and two testcross hybrid populations derived from the crosses of RILs and two cytoplasmic male sterile lines (YTA [indica] and Z7A [japonica]). Using deep resequencing data, a high-density genetic map containing 4758 bin markers was constructed, with a total map distance of 2356.41 cM. Finally, 31 PH-related QTLs for different PH component lengths or tiller numbers across five seasons were identified. Two major environment-specific PH QTLs were stably detected in Hainan (qPH-3.1) or Hubei (qPH-5.1), which have undergone significant functional alterations in rice with changes in geographical environment. Based on comparative genomics, gene function annotation, homolog identification, and existing literature (pioneering studies), candidate genes for multiple QTLs were fine-mapped, and the candidate genes qPH-3.1 and qPH-5.1 for PH were further validated using CRISPR-Cas9 gene editing. Specifically, qPH-3.1 was characterized as a pleiotropic gene, and the qPH-3.1 knockout line showed reduced PH, delayed heading, a decreased seed setting rate, and increased tiller numbers. Importantly, 10 PH heterosis-related QTLs were identified in the testcross populations, and a better-parent heterosis locus (qBPH-5.2) completely covered qPH-5.1. Furthermore, the cross results of fixed-genotype RILs verified the dominant effects of qPH-3.1 and qPH-5.1. Together, these findings further our understanding of the genetic mechanisms of PH and offer multiple highly reliable gene targets for breeding rice varieties with ideal architecture and high yield potential in the immediate future.

Quantitative Trait Locus Mapping and Identification of Candidate Genes for Resistance to Fusarium Wilt Race 7 Using a Resequencing-Based High Density Genetic Bin Map in a Recombinant Inbred Line Population of Gossypium barbadense

Fusarium wilt caused by Fusarium oxysporum f. sp. vasinfectum (FOV) is one of the most destructive diseases in cotton (Gossypium spp.) production, and use of resistant cultivars is the most cost-effective method managing the disease. To understand the genetic basis of cotton resistance to FOV race 7 (FOV7), this study evaluated a recombinant inbred line (RIL) population of 110 lines of G. barbadense from a cross between susceptible Xinhai 14 and resistant 06-146 in eight tests and constructed a high-density genetic linkage map with resequencing-based 933,845 single-nucleotide polymorphism (SNP) markers covering a total genetic distance of 2483.17 cM. Nine quantitative trait loci (QTLs) for FOV7 resistance were identified, including qFOV7-D03-1 on chromosome D03 in two tests. Through a comparative analysis of gene expression and DNA sequence for predicted genes within the QTL region between the two parents and selected lines inoculated with FOV7, GB_D03G0217 encoding for a calmodulin (CaM)-like (CML) protein was identified as a candidate gene. A further analysis confirmed that the expression of GB_D03G0217 was suppressed, leading to increased disease severity in plants of the resistant parent with virus induced gene silencing (VIGS).

When tropical and subtropical congeners met: Multiple ancient hybridization events within Eriobotrya in the Yunnan-Guizhou Plateau, a tropical-subtropical transition area in China

Global climate changes during the Miocene may have created ample opportunities for hybridization between members of tropical and subtropical biomes at the boundary between these zones. Yet, very few studies have explored this possibility. The Yunnan-Guizhou Plateau (YGP) in Southwest China is a biodiversity hotspot for vascular plants, located in a transitional area between the floristic regions of tropical Southeast Asia and subtropical East Asia. The genus Eriobotrya (Rosaceae) comprises both tropical and subtropical taxa, with 12 species recorded in the YGP, making it a suitable basis for testing the hypothesis of between-biome hybridization. Therefore, we surveyed the evolutionary history of Eriobotrya by examining three chloroplast regions and five nuclear genes for 817 individuals (47 populations) of 23 Eriobotrya species (including 19 populations of 12 species in the YGP), plus genome re-sequencing of 33 representative samples. We concluded that: (1) phylogenetic positions for 16 species exhibited strong cytonuclear conflicts, most probably due to ancient hybridization; (2) the YGP is a hotspot for hybridization, with 11 species showing clear evidence of chloroplast capture; and (3) Eriobotrya probably originated in tropical Asia during the Eocene. From the Miocene onwards, the intensification of the Eastern Asia monsoon and global cooling may have shifted the tropical-subtropical boundary and caused secondary contact between species, thus providing ample opportunity for hybridization and diversification of Eriobotrya, especially in the YGP. Our study highlights the significant role that paleoclimate changes probably played in driving hybridization and generating rich species diversity in climate transition zones.

Precise Characterization and Tracking of Stably Inherited Artificial Minichromosomes Made by Telomere-Mediated Chromosome Truncation in Brassica napus

Plant artificial minichromosomes are the next-generation technology for plant genetic engineering and represent an independent platform for expressing foreign genes and the tools for studying the structure and function of chromosomes. Minichromosomes have been successfully produced by telomere-mediated chromosome truncation in several plants. However, previous studies have primarily focused on the construction and rough characterization of minichromosomes, while the development of stably inherited minichromosomes and their precise characterization and tracking over different generations have rarely been demonstrated. In this study, a 0.35-kb direct repeat of the Arabidopsis telomeric sequence was transformed into Brassica napus to produce artificial minichromosomes, which were analyzed by multifluorescence in situ hybridization (multi-FISH), Southern hybridization, and primer extension telomere rapid amplification (PETRA). The stably inherited minichromosomes C2 and C4 were developed by crossing transgenic plants with wild-type plants and then selfing the hybrids. Notably, two truncation sites on chromosomes C2 and C4, respectively, were identified by resequencing; thus, the artificial minichromosomes were tracked over different generations with insertion site-specific PCR. This study provided two stably inherited minichromosomes in oilseed rape and describes approaches to precisely characterize the truncation position and track the minichromosomes in offspring through multi-FISH, genome resequencing, and insertion site-specific PCR.

Karst) by whole genome resequencing and its transferability to other spruce species

Norway spruce (Picea abies L. Karst) is one of the most important forest tree species with significant economic and ecological impact in Europe. For decades, genomic and genetic studies on Norway spruce have been challenging due to the large and repetitive genome (19.6 Gb with more than 70% being repetitive). To accelerate genomic studies, including population genetics, genome‐wide association studies (GWAS) and genomic selection (GS), in Norway spruce and related species, we here report on the design and performance of a 50K single nucleotide polymorphism (SNP) genotyping array for Norway spruce. The array is developed based on whole genome resequencing (WGS), making it the first WGS‐based SNP array in any conifer species so far. After identifying SNPs using genome resequencing data from 29 trees collected in northern Europe, we adopted a two‐step approach to design the array. First, we built a 450K screening array and used this to genotype a population of 480 trees sampled from both natural and breeding populations across the Norway spruce distribution range. These samples were then used to select high‐confidence probes that were put on the final 50K array. The SNPs selected are distributed over 45,552 scaffolds from the P. abies version 1.0 genome assembly and target 19,954 unique gene models with an even coverage of the 12 linkage groups in Norway spruce. We show that the array has a 99.5% probe specificity, >98% Mendelian allelic inheritance concordance, an average sample call rate of 96.30% and an SNP call rate of 98.90% in family trios and haploid tissues. We also observed that 23,797 probes (50%) could be identified with high confidence in three other spruce species (white spruce [Picea glauca], black spruce [P. mariana] and Sitka spruce [P. sitchensis]). The high‐quality genotyping array will be a valuable resource for genetic and genomic studies in Norway spruce as well as in other conifer species of the same genus.

Resequencing of 296 cultivated and wild lotus accessions unravels its evolution and breeding history

Lotus (family: Nelumbonaceae) are perennial aquatic plants that represent one of the most ancient basal dicots. In the present study, we resequenced 296 lotus accessions from various geographical locations and germplasms to explore their genomic diversity and population structure. This germplasm set consisted of four accessions of American wild lotus and 292 accessions of Asian lotus, which were divided into four subgroups: wild, rhizome, flower and seed. Total single nucleotide polymorphisms (SNPs) suggested that the wild lotus had the highest variant number (7 191 010). Population structure and genome diversity analysis indicated that the American wild lotus demonstrated a distant genetic relationship with the Asian lotus. Furthermore, the seed and rhizome lotus groups had not originated from a single source but rather had a more complex multisource origin. Besides that, the seed lotus showed higher genetic diversity, which might have been due to the gene flow from the flower lotus to seed lotus by artificial crossing, and the rhizome lotus showed a much lower genetic diversity than the other groups. The present study provides SNP markers for lotus genomic diversity analysis, which will be useful for guiding lotus breeding.

A complete set of monosomic alien addition lines developed from Gossypium anomalum in a Gossypium hirsutum background: genotypic and phenotypic characterization

Gossypium anomalum (B₁B₁) is a valuable wild resource for the genetic improvement of G. hirsutum (A₁A₁D₁D₁) in terms of fiber quality and disease and pest resistance, but the inherent difficulties in distant hybridization hinder its utilization in breeding programs. Monosomic alien addition lines (MAALs) are powerful tools for interspecific gene transfer. First, to access useful genes from G. anomalum, a fertile hexaploid from G. hirsutum × G. anomalum was obtained and then additional chromosomes were selected using SSR markers in successive backcrosses and self-crossing from BC₂F₁ to BC₄F₄. Finally, a complete set of 13 MAALs were developed. All the MAALs were confirmed by chromosome-specific anchored SSRs and genome-wide resequencing. The MAALs demonstrated abundant variation in morphological, agronomic, yield-related, and fiber quality traits. MAAL_3B had excellent fiber strength and fineness, indicating that the transmitted chromosome may carry desirable genes for the observed phenotypes. This complete set of MAALs will provide important genetic bridge material for the identification and introgression of favorable genes from G. anomalum and lay an important foundation for the genetic improvement of cotton.

Combining QTL Mapping with Genome Resequencing Identifies an Indel in an R Gene that is Associated with Variation in Leaf Rust Disease Resistance in Poplar

Poplar trees (Populus spp.) are important and are widely grown worldwide. However, the extensive occurrence of leaf rust disease caused by Melampsora spp. seriously inhibits their growth and reduces their biomass. In our previous study, a high-quality genetic map was constructed for the poplar F1 population I-69 × XYY by using next-generation sequencing-based genotyping-by-sequencing. Here, we collected phenotypic data on leaf rust disease resistance on three different dates for all 300 progenies of the F1 population. Combining a high-quality genetic map and phenotypic data, we were able to detect 11 major quantitative trait loci (QTLs) for leaf rust disease resistance. Among these 11 QTLs, two pairs were detected on at least two dates. In the corresponding genomic sequence, we found that resistance (R) gene clusters were located in these two QTL regions. By using genome resequencing, PCR confirmation and statistical analysis, a 611-bp deletion within an R gene in one QTL region was found to be associated with variation in leaf rust disease resistance. A PCR-based examination of this 611-bp deletion was performed. This 611-bp deletion was also found to affect mRNA splicing and form a new protein with the loss of some key protein domains. Based on this study, we were able to determine the genetic architecture of variation in poplar leaf rust disease resistance, and the 611-bp deletion in the R gene could be used as a diagnostic marker for future poplar molecular breeding.

Association between blooming time and climatic adaptation in Prunus mume

Prunus mume Sieb. et Zucc. is an important fruit crop of the subtropical region, originating in China. It blooms earlier than other deciduous fruit trees, but different regions have different blooming periods. The time of anthesis is related to the dormancy period, and a certain amount of chilling promotes bud break and blooming. To identify the relationship between blooming time and the climatic adaptation of P. mume cultivars in China, the nuclear and chloroplast genomes of 19 cultivars from the main cultivation areas of P. mume in China were resequenced. The average depth of coverage was 34X-76X, and a total of 388,134 single nucleotide polymorphisms were located within the coding regions of the gene (CDs). Additionally, the 19 cultivar accessions were divided into three groups based on their blooming time: early, mid, and late. Associated with the blooming time groups, 21 selective sweep regions were identified, which could provide evidence supporting the possible model of P. mume domestication originating due to natural selection. Furthermore, we identified a flowering gene, FRIGIDA-LIKE 3 (FRL3), seems to affect the blooming time and the climatic adaptation of P. mume cultivars. This study is a major step toward understanding the climatic adaptation of P. mume cultivars in China.

Whole-Genome Resequencing Identifies KIT New Alleles That Affect Coat Color Phenotypes in Pigs

The Duroc × (Landrace × Large White) hybrid pig (DLY) is the most popular commercial pig used in the Chinese pig industry. DLY pigs are usually white but sometimes show colored phenotypes. Colored DLY pigs are not favored by slaughterhouses and retailers, thus causing certain economic losses to farmers in China. In this study, we first conducted a genome-wide association study and RNA sequencing to demonstrate that KIT variants are responsible for diversifying coat color phenotypes segregating in a DLY population. We then defined the precise sizes and locations of four duplications (DUP1-4), four candidate causative mutations at the KIT locus, in the pig reference genome using the whole-genome sequence data of representative colored individuals. The sequence data also enabled us to identify a list of new KIT alleles. By investigating the association between these new alleles and coat color phenotypes, we provide further evidence that DUP2 is another causative mutation for the solid white coat color in pigs. DUP1 (the KIT gene duplication), DUP2 and the splice mutation are all required for the manifestation of a solid white coat color. DUP4 had a more significant effect on the formation of the belt phenotype compared with DUP3. Given the necessity of DUP2 for the solid white coat color, we detected I^N/I^N homozygotes lacking DUP2 in Large White and Landrace pigs and found that French Landrace pigs had the highest frequency (8.98%) of I^N/I^N individuals. This study not only advances our understanding of the molecular mechanism of the color phenotype in pigs, but also establishes a simple and accurate method for the screening of KIT I^N/I^N homozygotes in Large White and Landrace that would cause colored DLY pigs.

Whole-Genome Resequencing of a Worldwide Collection of Rapeseed Accessions Reveals the Genetic Basis of Ecotype Divergence

Rapeseed (Brassica napus), an important oilseed crop, has adapted to diverse climate zones and latitudes by forming three main ecotype groups, namely winter, semi-winter, and spring types. However, genetic variations underlying the divergence of these ecotypes are largely unknown. Here, we report the global pattern of genetic polymorphisms in rapeseed determined by resequencing a worldwide collection of 991 germplasm accessions. A total of 5.56 and 5.53 million single-nucleotide polymorphisms (SNPs) as well as 1.86 and 1.92 million InDels were identified by mapping reads to the reference genomes of "Darmor-bzh" and "Tapidor," respectively. We generated a map of allelic drift paths that shows splits and mixtures of the main populations, and revealed an asymmetric evolution of the two subgenomes of B. napus by calculating the genetic diversity and linkage disequilibrium parameters. Selective-sweep analysis revealed genetic changes in genes orthologous to those regulating various aspects of plant development and response to stresses. A genome-wide association study identified SNPs in the promoter regions of FLOWERING LOCUS T and FLOWERING LOCUS C orthologs that corresponded to the different rapeseed ecotype groups. Our study provides important insights into the genomic footprints of rapeseed evolution and flowering-time divergence among three ecotype groups, and will facilitate screening of molecular markers for accelerating rapeseed breeding.

Genomic adaptation of flowering-time genes during the expansion of rice cultivation area

The diversification of flowering time in response to natural environments is critical for the spread of crops to diverse geographic regions. In contrast with recent advances in understanding the molecular basis of photoperiodic flowering in rice (Oryza sativa), little is known about how flowering-time diversification is structured within rice subspecies. By analyzing genome sequencing data and a set of 429 chromosome segment substitution lines (CSSLs) originating from 10 diverse rice accessions with wide distributions, we revealed diverse effects of allelic variations for common flowering-time quantitative trait loci in the recipient's background. Although functional variations associated with a few loci corresponded to standing variations among subspecies, the identified functional nucleotide polymorphisms occurred recently after rice subgroup differentiation, indicating that the functional diversity of flowering-time gene sequences was not particularly associated with phylogenetic relationship between rice subspecies. Intensive analysis of the Hd1 genomic region identified the signature of an early introgression of the Hd1 with key mutation(s) in aus and temperate japonica accessions. Our data suggested that, after such key introgressions, new mutations were selected and accelerated the flowering-time diversity within subspecies during the expansion of rice cultivation area. This finding may imply that new genome-wide changes for flowering-time adaptation are one of the critical determinants for establishing genomic architecture of local rice subgroups. In-depth analyses of various rice genomes coupling with the genetically confirmed phenotypic changes in a large set of CSSLs enabled us to demonstrate how rice genome dynamics has coordinated with the adaptation of cultivated rice during the expansion of cultivation area.

Genomic adaptation of flowering-time genes during the expansion of rice cultivation area

The diversification of flowering time in response to natural environments is critical for the spread of crops to diverse geographic regions. In contrast with recent advances in understanding the molecular basis of photoperiodic flowering in rice (Oryza sativa), little is known about how flowering-time diversification is structured within rice subspecies. By analyzing genome sequencing data and a set of 429 chromosome segment substitution lines (CSSLs) originating from 10 diverse rice accessions with wide distributions, we revealed diverse effects of allelic variations for common flowering-time quantitative trait loci in the recipient's background. Although functional variations associated with a few loci corresponded to standing variations among subspecies, the identified functional nucleotide polymorphisms occurred recently after rice subgroup differentiation, indicating that the functional diversity of flowering-time gene sequences was not particularly associated with phylogenetic relationship between rice subspecies. Intensive analysis of the Hd1 genomic region identified the signature of an early introgression of the Hd1 with key mutation(s) in aus and temperate japonica accessions. Our data suggested that, after such key introgressions, new mutations were selected and accelerated the flowering-time diversity within subspecies during the expansion of rice cultivation area. This finding may imply that new genome-wide changes for flowering-time adaptation are one of the critical determinants for establishing genomic architecture of local rice subgroups. In-depth analyses of various rice genomes coupling with the genetically confirmed phenotypic changes in a large set of CSSLs enabled us to demonstrate how rice genome dynamics has coordinated with the adaptation of cultivated rice during the expansion of cultivation area.

In vivo transposon tagging in the nonheterocystous nitrogen-fixing cyanobacterium Leptolyngbya boryana

Nitrogenase is an oxygen-vulnerable metalloenzyme that catalyzes nitrogen fixation. It largely remains unknown how nitrogenase coexists with oxygenic photosynthesis in nonheterocystous cyanobacteria, since there have been no appropriate model cyanobacteria so far. Here, we demonstrate in vivo transposon tagging in the nonheterocystous cyanobacterium Leptolyngbya boryana as a forward genetics approach. By conjugative transfer, a mini-Tn5-derived vector, pKUT-Tn5-Sm/Sp, was transferred from Escherichia coli to L. boryana cells. Of 1839 streptomycin-resistant colonies, we isolated three mutants showing aberrant diazotrophic growth. Genome resequencing identified the insertion sites of the transposon in the mutants. This in vivo transposon tagging mutagenesis of L. boryana provides a promising system to investigate molecular mechanisms to resolve the Oxygen Paradox between nitrogen fixation and oxygenic photosynthesis in cyanobacteria.

Genomic adaptation to drought in wild barley is driven by edaphic natural selection at the Tabigha Evolution Slope

Ecological divergence at a microsite suggests adaptive evolution, and this study examined two abutting wild barley populations, each 100 m across, differentially adapted to drought tolerance on two contrasting soil types, Terra Rossa and basalt at the Tabigha Evolution Slope, Israel. We resequenced the genomes of seven and six wild barley genotypes inhabiting the Terra Rossa and basalt soils, respectively, and identified a total of 69,192,653 single-nucleotide variants (SNVs) and insertions/deletions in comparison with a reference barley genome. Comparative genomic analysis between these abutting wild barley populations involved 19,615,087 high-quality SNVs. The results revealed dramatically different selection sweep regions relevant to drought tolerance driven by edaphic natural selection within 2,577 selected genes in these regions, including key drought-responsive genes associated with ABA synthesis and degradation (such as Cytochrome P450 protein) and ABA receptor complex (such as PYL2, SNF1-related kinase). The genetic diversity of the wild barley population inhabiting Terra Rossa soil is much higher than that from the basalt soil. Additionally, we identified different sets of genes for drought adaptation in the wild barley populations from Terra Rossa soil and from wild barley populations from Evolution Canyon I at Mount Carmel. These genes are associated with abscisic acid signaling, signaling and metabolism of reactive oxygen species, detoxification and antioxidative systems, rapid osmotic adjustment, and deep root morphology. The unique mechanisms for drought adaptation of the wild barley from the Tabigha Evolution Slope may be useful for crop improvement, particularly for breeding of barley cultivars with high drought tolerance.

BIFURCATE FLOWER TRUSS: a novel locus controlling inflorescence branching in tomato contains a defective MAP kinase gene

A mutant line, bifurcate flower truss (bif), was recovered from a tomato genetics programme. Plants from the control line produced a mean of 0.16 branches per truss, whereas the value for bif plants was 4.1. This increase in branching was accompanied by a 3.3-fold increase in flower number and showed a significant interaction with exposure to low temperature during truss development. The control line and bif genomes were resequenced and the bif gene was mapped to a 2.01 Mbp interval on chromosome 12; all coding region polymorphisms in the interval were surveyed, and five candidate genes displaying altered protein sequences were detected. One of these genes, SlMAPK1, encoding a mitogen-activated protein (MAP) kinase, contained a leucine to stop codon mutation predicted to disrupt kinase function. SlMAPK1 is an excellent candidate for bif because knock-out mutations of an Arabidopsis orthologue MPK6 were reported to have increased flower number. An introgression browser was used to demonstrate that the origin of the bif genomic DNA at the BIF locus was Solanum galapagense and that the SlMAPK1 null mutant is a naturally occurring allele widespread only on the Galápagos Islands. This work strongly implicates SlMAPK1 as part of the network of genes controlling inflorescence branching in tomato.

A New Single Nucleotide Polymorphism Database for Rainbow Trout Generated Through Whole Genome Resequencing

Single-nucleotide polymorphisms (SNPs) are highly abundant markers, which are broadly distributed in animal genomes. For rainbow trout (Oncorhynchus mykiss), SNP discovery has been previously done through sequencing of restriction-site associated DNA (RAD) libraries, reduced representation libraries (RRL) and RNA sequencing. Recently we have performed high coverage whole genome resequencing with 61 unrelated samples, representing a wide range of rainbow trout and steelhead populations, with 49 new samples added to 12 aquaculture samples from AquaGen (Norway) that we previously used for SNP discovery. Of the 49 new samples, 11 were double-haploid lines from Washington State University (WSU) and 38 represented wild and hatchery populations from a wide range of geographic distribution and with divergent migratory phenotypes. We then mapped the sequences to the new rainbow trout reference genome assembly (GCA_002163495.1) which is based on the Swanson YY doubled haploid line. Variant calling was conducted with FreeBayes and SAMtools mpileup, followed by filtering of SNPs based on quality score, sequence complexity, read depth on the locus, and number of genotyped samples. Results from the two variant calling programs were compared and genotypes of the double haploid samples were used for detecting and filtering putative paralogous sequence variants (PSVs) and multi-sequence variants (MSVs). Overall, 30,302,087 SNPs were identified on the rainbow trout genome 29 chromosomes and 1,139,018 on unplaced scaffolds, with 4,042,723 SNPs having high minor allele frequency (MAF > 0.25). The average SNP density on the chromosomes was one SNP per 64 bp, or 15.6 SNPs per 1 kb. Results from the phylogenetic analysis that we conducted indicate that the SNP markers contain enough population-specific polymorphisms for recovering population relationships despite the small sample size used. Intra-Population polymorphism assessment revealed high level of polymorphism and heterozygosity within each population. We also provide functional annotation based on the genome position of each SNP and evaluate the use of clonal lines for filtering of PSVs and MSVs. These SNPs form a new database, which provides an important resource for a new high density SNP array design and for other SNP genotyping platforms used for genetic and genomics studies of this iconic salmonid fish species.

Exploring Genomic Variants Related to Residual Feed Intake in Local and Commercial Chickens by Whole Genomic Resequencing

Improving feed efficiency is a major goal in poultry production to reduce production costs and increase profitability. The genomic variants and possible molecular mechanisms responsible for residual feed intake (RFI) in chickens, however, remain poorly understood. In this study, using both local and commercial breeds, genome re-sequencing of low RFI and high RFI chickens was performed to elucidate the genomic variants underlying RFI. Results showed that 8,505,214 and 8,479,041 single nucleotide polymorphisms (SNPs) were detected in low and high RFI Beijing-You chickens, respectively; 8,352,008 and 8,372,769 SNPs were detected in low- and high-RFI Cobb chickens, respectively. Through a series of filtering processes, 3746 candidate SNPs assigned to 1137 genes in Beijing-You chickens and 575 candidate SNPs (448 genes) in Cobb chickens were found. The validation of the selected 191 SNPs showed that 46 SNPs were significantly associated with the RFI in an independent population of 779 Cobb chickens, suggesting that the method of screening associated SNPs with whole genome sequencing (WGS) strategy was reasonable. Functions annotation of RFI-related genes indicated that genes in Beijing-You were enriched in lipid and carbohydrate metabolism, as well as the phosphatase and tensin homolog (PTEN) signaling pathway. In Cobb, however, RFI-related genes were enriched in the feed behavior process and cAMP responsive element binding protein (CREB) signaling pathway. For both breeds, organismal development physiological processes were enriched. Correspondingly, NOS1, PHKG1, NEU3 and PIP5K1B were differentially expressed in Beijing-You, while CDC42, CSK, PIK3R3, CAMK4 and PLCB4 were differentially expressed in Cobb, suggesting that these might be key genes that contribute to RFI. The results of the present study identified numerous novel SNPs for RFI, which provide candidate biomarkers for use in the genetic selection for RFI. The study has improved knowledge of the genomic variants and potential biological pathways underlying RFI in chickens.

Evaluation of Quality Assessment Protocols for High Throughput Genome Resequencing Data

Large-scale initiatives aiming to recover the complete sequence of thousands of human genomes are currently being undertaken worldwide, concurring to the generation of a comprehensive catalog of human genetic variation. The ultimate and most ambitious goal of human population scale genomics is the characterization of the so-called human “variome,” through the identification of causal mutations or haplotypes. Several research institutions worldwide currently use genotyping assays based on Next-Generation Sequencing (NGS) for diagnostics and clinical screenings, and the widespread application of such technologies promises major revolutions in medical science. Bioinformatic analysis of human resequencing data is one of the main factors limiting the effectiveness and general applicability of NGS for clinical studies. The requirement for multiple tools, to be combined in dedicated protocols in order to accommodate different types of data (gene panels, exomes, or whole genomes) and the high variability of the data makes difficult the establishment of a ultimate strategy of general use. While there already exist several studies comparing sensitivity and accuracy of bioinformatic pipelines for the identification of single nucleotide variants from resequencing data, little is known about the impact of quality assessment and reads pre-processing strategies. In this work we discuss major strengths and limitations of the various genome resequencing protocols are currently used in molecular diagnostics and for the discovery of novel disease-causing mutations. By taking advantage of publicly available data we devise and suggest a series of best practices for the pre-processing of the data that consistently improve the outcome of genotyping with minimal impacts on computational costs.

LINE Insertion Polymorphisms are Abundant but at Low Frequencies across Populations of Anolis carolinensis

Vertebrate genomes differ considerably in size and structure. Among the features that show the most variation is the abundance of Long Interspersed Nuclear Elements (LINEs). Mammalian genomes contain 100,000s LINEs that belong to a single clade, L1, and in most species a single family is usually active at a time. In contrast, non-mammalian vertebrates (fish, amphibians and reptiles) contain multiple active families, belonging to several clades, but each of them is represented by a small number of recently inserted copies. It is unclear why vertebrate genomes harbor such drastic differences in LINE composition. To address this issue, we conducted whole genome resequencing to investigate the population genomics of LINEs across 13 genomes of the lizard Anolis carolinensis sampled from two geographically and genetically distinct populations in the Eastern Florida and the Gulf Atlantic regions of the United States. We used the Mobile Element Locator Tool to identify and genotype polymorphic insertions from five major clades of LINEs (CR1, L1, L2, RTE and R4) and the 41 subfamilies that constitute them. Across these groups we found large variation in the frequency of polymorphic insertions and the observed length distributions of these insertions, suggesting these groups vary in their activity and how frequently they successfully generate full-length, potentially active copies. Though we found an abundance of polymorphic insertions (over 45,000) most of these were observed at low frequencies and typically appeared as singletons. Site frequency spectra for most LINEs showed a significant shift toward low frequency alleles compared to the spectra observed for total genomic single nucleotide polymorphisms. Using Tajima's D, FST and the mean number of pairwise differences in LINE insertion polymorphisms, we found evidence that negative selection is acting on LINE families in a length-dependent manner, its effects being stronger in the larger Eastern Florida population. Our results suggest that a large effective population size and negative selection limit the expansion of polymorphic LINE insertions across these populations and that the probability of LINE polymorphisms reaching fixation is extremely low.

Soybean domestication: the origin, genetic architecture and molecular bases

Domestication provides an important model for the study of evolution, and information learned from domestication research aids in the continued improvement of crop species. Recent progress in de novo assembly and whole-genome resequencing of wild and cultivated soybean genomes, in addition to new archeological discoveries, sheds light on the origin of this important crop and provides a clearer view on the modes of artificial selection that drove soybean domestication and diversification. This novel genomic information enables the search for polymorphisms that underlie variation in agronomic traits and highlights genes that exhibit a signature of selection, leading to the identification of a number of candidate genes that may have played important roles in soybean domestication, diversification and improvement. These discoveries provide a novel point of comparison on the evolutionary bases of important agronomic traits among different crop species.

Mechanistic Study of Utilization of Water-Insoluble Saccharomyces cerevisiae Glucans by Bifidobacterium breve Strain JCM1192

Bifidobacteria exert beneficial effects on hosts and are extensively used as probiotics. However, due to the genetic inaccessibility of these bacteria, little is known about their mechanisms of carbohydrate utilization and regulation. Bifidobacterium breve strain JCM1192 can grow on water-insoluble yeast (Saccharomyces cerevisiae) cell wall glucans (YCWG), which were recently considered as potential prebiotics. According to the results of ¹H nuclear magnetic resonance (NMR) spectrometry, the YCWG were composed of highly branched (1→3,1→6)-β-glucans and (1→4,1→6)-α-glucans. Although the YCWG were composed of 78.3% β-glucans and 21.7% α-glucans, only α-glucans were consumed by the B. breve strain. The ABC transporter (malEFG1) and pullulanase (aapA) genes were transcriptionally upregulated in the metabolism of insoluble yeast glucans, suggesting their potential involvement in the process. A nonsense mutation identified in the gene encoding an ABC transporter ATP-binding protein (MalK) led to growth failure of an ethyl methanesulfonate-generated mutant with yeast glucans. Coculture of the wild-type strain and the mutant showed that this protein was responsible for the import of yeast glucans or their breakdown products, rather than the export of α-glucan-catabolizing enzymes. Further characterization of the carbohydrate utilization of the mutant and three of its revertants indicated that this mutation was pleiotropic: the mutant could not grow with maltose, glycogen, dextrin, raffinose, cellobiose, melibiose, or turanose. We propose that insoluble yeast α-glucans are hydrolyzed by extracellular pullulanase into maltose and/or maltooligosaccharides, which are then transported into the cell by the ABC transport system composed of MalEFG1 and MalK. The mechanism elucidated here will facilitate the development of B. breve and water-insoluble yeast glucans as novel synbiotics.IMPORTANCE In general, Bifidobacterium strains are genetically intractable. Coupling classic forward genetics with next-generation sequencing, here we identified an ABC transporter ATP-binding protein (MalK) responsible for the import of insoluble yeast glucan breakdown products by B. breve JCM1192. We demonstrated the pleiotropic effects of the ABC transporter ATP-binding protein in maltose/maltooligosaccharide, raffinose, cellobiose, melibiose, and turanose transport. With the addition of transcriptional analysis, we propose that insoluble yeast glucans are broken down by extracellular pullulanase into maltose and/or maltooligosaccharides, which are then transported into the cell by the ABC transport system composed of MalEFG1 and MalK. The mechanism elucidated here will facilitate the development of B. breve and water-insoluble yeast glucans as novel synbiotics.

Genome-wide DNA polymorphism in the indica rice varieties RGD-7S and Taifeng B as revealed by whole genome re-sequencing

Next-generation sequencing technologies provide opportunities to further understand genetic variation, even within closely related cultivars. We performed whole genome resequencing of two elite indica rice varieties, RGD-7S and Taifeng B, whose F1 progeny showed hybrid weakness and hybrid vigor when grown in the early- and late-cropping seasons, respectively. Approximately 150 million 100-bp pair-end reads were generated, which covered ∼86% of the rice (Oryza sativa L. japonica 'Nipponbare') reference genome. A total of 2,758,740 polymorphic sites including 2,408,845 SNPs and 349,895 InDels were detected in RGD-7S and Taifeng B, respectively. Applying stringent parameters, we identified 961,791 SNPs and 46,640 InDels between RGD-7S and Taifeng B (RGD-7S/Taifeng B). The density of DNA polymorphisms was 256.8 SNPs and 12.5 InDels per 100 kb for RGD-7S/Taifeng B. Copy number variations (CNVs) were also investigated. In RGD-7S, 1989 of 2727 CNVs were overlapped in 218 genes, and 1231 of 2010 CNVs were annotated in 175 genes in Taifeng B. In addition, we verified a subset of InDels in the interval of hybrid weakness genes, Hw3 and Hw4, and obtained some polymorphic InDel markers, which will provide a sound foundation for cloning hybrid weakness genes. Analysis of genomic variations will also contribute to understanding the genetic basis of hybrid weakness and heterosis.

Genomewide ancestry and divergence patterns from low-coverage sequencing data reveal a complex history of admixture in wild baboons

Naturally occurring admixture has now been documented in every major primate lineage, suggesting its key role in primate evolutionary history. Active primate hybrid zones can provide valuable insight into this process. Here, we investigate the history of admixture in one of the best-studied natural primate hybrid zones, between yellow baboons (Papio cynocephalus) and anubis baboons (Papio anubis) in the Amboseli ecosystem of Kenya. We generated a new genome assembly for yellow baboon and low-coverage genomewide resequencing data from yellow baboons, anubis baboons and known hybrids (n = 44). Using a novel composite likelihood method for estimating local ancestry from low-coverage data, we found high levels of genetic diversity and genetic differentiation between the parent taxa, and excellent agreement between genome-scale ancestry estimates and a priori pedigree, life history and morphology-based estimates (r(2)  = 0.899). However, even putatively unadmixed Amboseli yellow individuals carried a substantial proportion of anubis ancestry, presumably due to historical admixture. Further, the distribution of shared vs. fixed differences between a putatively unadmixed Amboseli yellow baboon and an unadmixed anubis baboon, both sequenced at high coverage, is inconsistent with simple isolation-migration or equilibrium migration models. Our findings suggest a complex process of intermittent contact that has occurred multiple times in baboon evolutionary history, despite no obvious fitness costs to hybrids or major geographic or behavioural barriers. In combination with the extensive phenotypic data available for baboon hybrids, our results provide valuable context for understanding the history of admixture in primates, including in our own lineage.

Larva-mediated chalkbrood resistance-associated single nucleotide polymorphism markers in the honey bee Apis mellifera

Chalkbrood is a disease affecting honey bees that seriously impairs brood growth and productivity of diseased colonies. Although honey bees can develop chalkbrood resistance naturally, the details underlying the mechanisms of resistance are not fully understood, and no easy method is currently available for selecting and breeding resistant bees. Finding the genes involved in the development of resistance and identifying single nucleotide polymorphisms (SNPs) that can be used as molecular markers of resistance is therefore a high priority. We conducted genome resequencing to compare resistant (Res) and susceptible (Sus) larvae that were selected following in vitro chalkbrood inoculation. Twelve genomic libraries, including 14.4 Gb of sequence data, were analysed using SNP-finding algorithms. Unique SNPs derived from chromosomes 2 and 11 were analysed in this study. SNPs from resistant individuals were confirmed by PCR and Sanger sequencing using in vitro reared larvae and resistant colonies. We found strong support for an association between the C allele at SNP C2587245T and chalkbrood resistance. SNP C2587245T may be useful as a genetic marker for the selection of chalkbrood resistance and high royal jelly production honey bee lines, thereby helping to minimize the negative effects of chalkbrood on managed honey bees.

Identifying the genome-wide genetic variation between precocious trifoliate orange and its wild type and developing new markers for genetics research

To increase our understanding of the genes involved in flowering in citrus, we performed genome resequencing of an early flowering trifoliate orange mutant (Poncirus trifoliata L. Raf.) and its wild type. At the genome level, 3,932,628 single nucleotide polymorphisms (SNPs), 1,293,383 insertion/deletion polymorphisms (InDels), and 52,135 structural variations were identified between the mutant and its wild type based on the citrus reference genome. Based on integrative analysis of resequencing and transcriptome analysis, 233,998 SNPs and 75,836 InDels were also identified between the mutant and its wild type at the transcriptional level. Also, 272 citrus homologous flowering-time transcripts containing genetic variation were also identified. Gene Ontology and Kyoto Encyclopaedia of Genes and Genomes annotation revealed that the transcripts containing the mutant- and the wild-type-specific InDel were involved in diverse biological processes and molecular function. Among these transcripts, there were 131 transcripts that were expressed differently in the two genotypes. When 268 selected InDels were tested on 32 genotypes of the three genera of Rutaceae for the genetic diversity assessment, these InDel-based markers showed high transferability. This work provides important information that will allow a better understanding of the citrus genome and that will be helpful for dissecting the genetic basis of important traits in citrus.

Efficient Genome-Wide Sequencing and Low-Coverage Pedigree Analysis from Noninvasively Collected Samples

Research on the genetics of natural populations was revolutionized in the 1990s by methods for genotyping noninvasively collected samples. However, these methods have remained largely unchanged for the past 20 years and lag far behind the genomics era. To close this gap, here we report an optimized laboratory protocol for genome-wide capture of endogenous DNA from noninvasively collected samples, coupled with a novel computational approach to reconstruct pedigree links from the resulting low-coverage data. We validated both methods using fecal samples from 62 wild baboons, including 48 from an independently constructed extended pedigree. We enriched fecal-derived DNA samples up to 40-fold for endogenous baboon DNA and reconstructed near-perfect pedigree relationships even with extremely low-coverage sequencing. We anticipate that these methods will be broadly applicable to the many research systems for which only noninvasive samples are available. The lab protocol and software (“WHODAD”) are freely available at www.tung-lab.org/protocols-and-software.html and www.xzlab.org/software.html, respectively.

Whole genome resequencing of Botrytis cinerea isolates identifies high levels of standing diversity

How standing genetic variation within a pathogen contributes to diversity in host/pathogen interactions is poorly understood, partly because most studied pathogens are host-specific, clonally reproducing organisms which complicates genetic analysis. In contrast, Botrytis cinerea is a sexually reproducing, true haploid ascomycete that can infect a wide range of diverse plant hosts. While previous work had shown significant genomic variation between two isolates, we proceeded to assess the level and frequency of standing variation in a population of B. cinerea. To begin measuring standing genetic variation in B. cinerea, we re-sequenced the genomes of 13 different isolates and aligned them to the previously sequenced T4 reference genome. In addition one of these isolates was resequenced from four independently repeated cultures. A high level of genetic diversity was found within the 13 isolates. Within this variation, we could identify clusters of genes with major effect polymorphisms, i.e., polymorphisms that lead to a predicted functional knockout, that surrounded genes involved in controlling vegetative incompatibility. The genotype at these loci was able to partially predict the interaction of these isolates in vegetative fusion assays showing that these loci control vegetative incompatibility. This suggests that the vegetative incompatibility loci within B. cinerea are associated with regions of increased genetic diversity. The genome re-sequencing of four clones from the one isolate (Grape) that had been independently propagated over 10 years showed no detectable spontaneous mutation. This suggests that B. cinerea does not display an elevated spontaneous mutation rate. Future work will allow us to test if, and how, this diversity may be contributing to the pathogen's broad host range.

Detecting rare structural variation in evolving microbial populations from new sequence junctions using breseq

New mutations leading to structural variation (SV) in genomes—in the form of mobile element insertions, large deletions, gene duplications, and other chromosomal rearrangements—can play a key role in microbial evolution. Yet, SV is considerably more difficult to predict from short-read genome resequencing data than single-nucleotide substitutions and indels (SN), so it is not yet routinely identified in studies that profile population-level genetic diversity over time in evolution experiments. We implemented an algorithm for detecting polymorphic SV as part of the breseq computational pipeline. This procedure examines split-read alignments, in which the two ends of a single sequencing read match disjoint locations in the reference genome, in order to detect structural variants and estimate their frequencies within a sample. We tested our algorithm using simulated Escherichia coli data and then applied it to 500- and 1000-generation population samples from the Lenski E. coli long-term evolution experiment (LTEE). Knowledge of genes that are targets of selection in the LTEE and mutations present in previously analyzed clonal isolates allowed us to evaluate the accuracy of our procedure. Overall, SV accounted for ~25% of the genetic diversity found in these samples. By profiling rare SV, we were able to identify many cases where alternative mutations in key genes transiently competed within a single population. We also found, unexpectedly, that mutations in two genes that rose to prominence at these early time points always went extinct in the long term. Because it is not limited by the base-calling error rate of the sequencing technology, our approach for identifying rare SV in whole-population samples may have a lower detection limit than similar predictions of SNs in these data sets. We anticipate that this functionality of breseq will be useful for providing a more complete picture of genome dynamics during evolution experiments with haploid microorganisms.

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 FST -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.

Transgene-free genome editing in Caenorhabditis elegans using CRISPR-Cas. Genetics. 2013;195:1167-1171. [PMID: 23979577 DOI: 10.1534/genetics.113.155879]

CRISPR-Cas is an efficient method for genome editing in organisms from bacteria to human cells. We describe a transgene-free method for CRISPR-Cas-mediated cleavage in nematodes, enabling RNA-homology-targeted deletions that cause loss of gene function; analysis of whole-genome sequencing indicates that the nuclease activity is highly specific.