An online database for einkorn wheat to aid in gene discovery and functional genomics studies

Diploid A-genome wheat (einkorn wheat) presents a nutrition-rich option as an ancient grain crop and a resource for the improvement of bread wheat against abiotic and biotic stresses. Realizing the importance of this wheat species, reference-level assemblies of two einkorn wheat accessions were generated (wild and domesticated). This work reports an einkorn genome database that provides an interface to the cereals research community to perform comparative genomics, applied genetics and breeding research. It features queries for annotated genes, the use of a recent genome browser release, and the ability to search for sequence alignments using a modern BLAST interface. Other features include a comparison of reference einkorn assemblies with other wheat cultivars through genomic synteny visualization and an alignment visualization tool for BLAST results. Altogether, this resource will help wheat research and breeding. Database URL  https://wheat.pw.usda.gov/GG3/pangenome.

[Serum immunoglobulin concentration and the risk of type 2 diabetes mellitus in adults in Tianjin City: a prospective cohort study]

Objective: To evaluate the association between immunoglobulin concentration and the risk of type 2 diabetes mellitus (T2DM) in adults in Tianjin City. Methods: Based on the Tianjin Chronic Low-grade Systemic Inflammation and Health (TCLSIHealth) cohort from January 2010 to December 2018, subjects who had completed the measurement of baseline immunoglobulin concentration and blood glucose concentration and not been diagnosed with any type of diabetes at baseline were selected in this study. The collected data included the concentration of serum immunoglobulin (IgG, IgM, IgA and IgE), fasting blood glucose and other potential confounders. The subjects were divided into four groups from Q1 to Q4 according to the quartiles of baseline immunoglobulin concentration. The multivariable Cox regression model was used to assess the association between the baseline immunoglobulin concentration and T2DM. Results: A total of 6 315 subjects aged (50.1±10.0) years were included. About 390 subjects were newly diagnosed with T2DM during the follow-up period. The incidence rate was 16.8/1 000 person-years. After adjusting for age, sex, waist circumference, smoking status, drinking status, eosinophil ratio, metabolic syndrome, first-or second-degree family history, and reciprocal adjusting for other immunoglobulin concentrations, compared to the lowest quartile concentration group Q1, subjects in group Q4 with the highest quartile of IgG concentration showed a lower risk of T2DM (HR=0.71, 95%CI: 0.52-0.97), and subjects in group Q4 with the highest quartile of IgM concentration also had a decreased risk of T2DM (HR=0.66, 95%CI: 0.47-0.91). Subjects in group Q4 with the highest quartile of IgA concentration had an increased risk of T2DM (HR=1.56, 95%CI: 1.18-2.07). The risk of T2DM decreased with the increase of serum IgG and IgM concentrations (Ptrend=0.018, Ptrend=0.010) and increased with the increase of serum IgA concentrations (Ptrend<0.001). No association was found between the concentration of IgE and T2DM risk (HR=0.99, 95%CI: 0.74-1.31, Ptrend=0.891). Conclusion: The concentration of IgG and IgM is negatively associated with the risk of T2DM, and the concentration of IgA is positively associated with the risk of T2DM in Tianjin City. The concentrations of IgG, IgM and IgA could be a predictor of hyperglycemia and T2DM.

Comparative genomic analysis of 5Mg chromosome of Aegilops geniculata and 5Uu chromosome of Aegilops umbellulata reveal genic diversity in the tertiary gene pool

Wheat is one of the most important cereal crops for the global food security. Due to its narrow genetic base, modern bread wheat cultivars face challenges from increasing abiotic and biotic stresses. Since genetic improvement is the most sustainable approach, finding novel genes and alleles is critical for enhancing the genetic diversity of wheat. The tertiary gene pool of wheat is considered a gold mine for genetic diversity as novel genes and alleles can be identified and transferred to wheat cultivars. Aegilops geniculata and Ae. umbellulata are the key members of the tertiary gene pool of wheat and harbor important genes against abiotic and biotic stresses. Homoeologous-group five chromosomes (5U^u and 5M^g) have been extensively studied from Ae. geniculata and Ae. umbellulata as they harbor several important genes including Lr57, Lr76, Yr40, Yr70, Sr53 and chromosomal pairing loci. In the present study, using chromosome DNA sequencing and RNAseq datasets, we performed comparative analysis to study homoeologous gene evolution in 5M^g, 5U^u, and group 5 wheat chromosomes. Our findings highlight the diversity of transcription factors and resistance genes, resulting from the differential expansion of the gene families. Both the chromosomes were found to be enriched with the "response to stimulus" category of genes providing resistance against biotic and abiotic stress. Phylogenetic study positioned the M genome closer to the D genome, with higher proximity to the A genome than the B genome. Over 4000 genes were impacted by SNPs on 5D, with 4-5% of those genes displaying non-disruptive variations that affect gene function.

Toward reducing the immunogenic potential of wheat flour: identification and characterization of wheat lines missing omega-5 gliadins encoded by the 1D chromosome

Eleven wheat lines that are missing genes for the 1D-encoded omega-5 gliadins will facilitate breeding efforts to reduce the immunogenic potential of wheat flour for patients susceptible to wheat allergy. Efforts to reduce the levels of allergens in wheat flour that cause wheat-dependent exercise-induced anaphylaxis are complicated by the presence of genes encoding omega-5 gliadins on both chromosomes 1B and 1D of hexaploid wheat. In this study, we screened 665 wheat germplasm samples using gene specific DNA markers for omega-5 gliadins encoded by the genes on 1D chromosome that were obtained from the reference wheat Chinese Spring. Eleven wheat lines missing the PCR product corresponding to 1D omega-5 gliadin gene sequences were identified. Two of the lines contained the 1BL·1RS translocation. Relative quantification of gene copy numbers by qPCR revealed that copy numbers of 1D omega-5 gliadins in the other nine lines were comparable to those in 1D null lines of Chinese Spring, while copy numbers of 1B omega-5 gliadins were like those of Chinese Spring. 2-D immunoblot analysis of total flour proteins from the selected lines using a specific monoclonal antibody against the N-terminal sequence of omega-5 gliadin showed no reactivity in regions of the blots containing previously identified 1D omega-5 gliadins. Interestingly, RP-UPLC analysis of the gliadin fractions of the selected lines indicated that the expression of omega-1,2 gliadins was also significantly reduced in seven of the lines, implying that 1D omega-5 gliadin and 1D omega-1,2 gliadin genes are tightly linked on the Gli-D1 loci of chromosome 1D. Wheat lines missing the omega-5 gliadins encoded by the genes on 1D chromosome should be useful in future breeding efforts to reduce the immunogenic potential of wheat flour.

The CC-NB-LRR protein BSR1 from Brachypodium confers resistance to Barley stripe mosaic virus in gramineous plants by recognising TGB1 movement protein

Although some nucleotide binding, leucine-rich repeat immune receptor (NLR) proteins conferring resistance to specific viruses have been identified in dicot plants, NLR proteins involved in viral resistance have not been described in monocots. We have used map-based cloning to isolate the CC-NB-LRR (CNL) Barley stripe mosaic virus (BSMV) resistance gene barley stripe resistance 1 (BSR1) from Brachypodium distachyon Bd3-1 inbred line. Stable BSR1 transgenic Brachypodium line Bd21-3, barley (Golden Promise) and wheat (Kenong 199) plants developed resistance against BSMV ND18 strain. Allelic variation analyses indicated that BSR1 is present in several Brachypodium accessions collected from countries in the Middle East. Protein domain swaps revealed that the intact LRR domain and the C-terminus of BSR1 are required for resistance. BSR1 interacts with the BSMV ND18 TGB1 protein in planta and shows temperature-sensitive antiviral resistance. The R³⁹⁰ and T³⁹² residues of TGB1_ND (ND18 strain) and the G¹⁹⁶ and K¹⁹⁷ residues within the BSR1 P-loop motif are key amino acids required for immune activation. BSR1 is the first cloned virus resistance gene encoding a typical CNL protein in monocots, highlighting the utility of the Brachypodium model for isolation and analysis of agronomically important genes for crop improvement.

[Associations between internet addiction, screen time and depressive symptoms]

Objective: To understand the associations between internet addiction, screen time (computer/mobile devices use and television watching time) and depressive symptoms in adults. Methods: A total of 6 932 adults aged <60 years from the Tianjin Chronic Low-grade Sgstemic Inflammation and Health (TCLSIH) Cohort of 2013-2019 were surveyed. The information about their computer/mobile devices use and television watching time were collected by using a self-reported questionnaire. The depressive symptoms were assessed using the self-rating depression scale (SDS). The adults surveyed were divided into two groups: non-depressive symptom group (SDS score <45) and depressive symptom group (SDS score ≥45). The associations between internet addiction, screen time and depressive symptoms were estimated using Cox proportional hazard regression models, with adjusting for multiple confounders. Results: After adjusting for confounding factors, the hazard ratios (HRs) of depressive symptom in the adults who had internet addiction before, had light internet addiction and had moderate or severe internet addiction were 0.83 (95%CI: 0.56-1.23) , 1.20 (95%CI: 1.03-1.41) for light and 1.48 (95%CI: 1.16-1.89), respectively, compared with those without internet addiction. The linear trend test results of the association between internet addiction and depressive symptoms was significant (trend P<0.001). Compared with the adults who used computer/mobile devices for <1 hour/day, the HRs of depressive symptoms in those who used computer/mobile devices for >1 hour, >3 hours, >5 hours and >10 hours were 0.59 (95%CI: 0.40-0.88), 0.58 (95%CI: 0.40-0.85), 0.52 (95%CI: 0.36-0.76) and 0.69 (95%CI: 0.45-1.05) respectively, a U-shaped association was found between computer/mobile devices use time and depressive symptoms (trend P<0.001). Compared with the adults who never watch TV, the HR of depressive symptoms was 1.36 (95%CI:1.09-1.69) for those watching TV for ≥3 hours/day in crude model and 1.34 (95%CI: 1.07-1.68) for those watching TV for ≥3 hours/day in adjusted model (trend P<0.001). Conclusion: Our findings suggested that internet addiction and television watching time were associated with an increased risk of depressive symptoms, while computer/mobile device use time was associated with a reduced risk of depressive symptoms.

CRISPR-Cas9 Gene Editing of the Sal1 Gene Family in Wheat

The highly conserved Sal1 encodes a bifunctional enzyme with inositol polyphosphate-1-phosphatase and 3′ (2′), 5′-bisphosphate nucleotidase activity and has been shown to alter abiotic stress tolerance in plants when disrupted. Precise gene editing techniques were used to generate Sal1 mutants in hexaploid bread wheat. The CRISPR (Clustered Regulatory Interspaced Short Palindromic Repeats) Cas9 system with three guide RNAs (gRNAs) was used to inactivate six Sal1 homologous genes within the Bobwhite wheat genome. The resulting mutant wheat plants with all their Sal1 genes disabled had slimmer stems, had a modest reduction in biomass and senesced more slowly in water limiting conditions, but did not exhibit improved yield under drought conditions. Our results show that multiplexed gRNAs enabled effective targeted gene editing of the Sal1 gene family in hexaploid wheat. These Sal1 mutant wheat plants will be a resource for further research studying the function of this gene family in wheat.

[Autopsies and placental examinations of perinatal fetal deaths: a clinicopathological analysis of 105 cases]

Objective: To summarize the clinicopathological factors related to perinatal fetal death and to evaluate importance of fetal autopsy and placental pathology. Methods: The clinicopathological data of 105 perinatal fetal deaths in Beijing Haidian Maternal and Child Health Hospital from November 2012 to December 2020 were retrospectively analyzed. Relevant literature was also reviewed. Results: The maternal age of the deceased fetuses ranged from 22 to 43 years with the average (31.35±4.04 years), and the gestational weeks were 28-40⁺⁶ weeks. Among them, 101 were singleton cases and 4 twin cases. 103 fetuses died in uterus and 2 died during delivery. Relevant factors analysis of the 105 perinatal fetal deaths showed that 86 cases (81.9%, 86/105) were related to umbilical cord/placental abnormality, 10 cases (9.5%, 10/105) uterine infection, 6 cases (5.7%, 6/105) fetal factors, 1 case was fetal maternal blood transfusion syndrome, 1 case twin blood transfusion syndrome, and 1 case died of complete uterine rupture. Among the 86 cases related to umbilical cord/placental abnormality, the diagnosis was most often based on the gross examination of placenta. The most common cause of death was umbilical cord torsion with thin root, followed by placental abruption, tight umbilical cord winding, vascular rupture and umbilical cord true knot. The morphology of placenta revealed mainly functional changes. Among the 10 cases related to intrauterine infections, the placenta generally showed lobular placental edema. The morphological characteristics of ascending infection were mainly acute chorioamnionitis, and the morphological characteristics of blood-borne infection were mainly acute or chronic villitis, as well as villous interstitial inflammation. Identification of viral inclusions suggested viral etiology, while the final diagnosis was relied on laboratory testing. Among the 6 cases related to fetal abnormality, the diagnostic value of placenta was limited and the diagnosis could be made with fetal autopsy. Conclusion: The causes of perinatal fetal death are complex, diverse, and often the synergistic result of multiple factors. Fetal autopsy and placental pathology are the key technical means to identify the cause of death and deserve more attention and utilization.

[Singleton placentas with abnormal shape: a clinicopathological analysis of 130 cases]

Objective: To investigate the pathological characteristics of singleton placenta with abnormal shape and its influence on the outcome of maternal-fetal pregnancy. Methods: The clinicopathological data of singleton placentas with abnormal shape from January 2014 to December 2020 in the Department of Pathology, Haidian Maternal and Children Health Hospital were analyzed retrospectively. Results: There were 130 singleton placentas with abnormal shape in this cohort, including 48 succenturiate placentas, 12 bilobed placentas, 50 marginate placentas, 13 circumvallate placentas, 3 annular placentas, 2 membranous placentas and 2 fenestrated placentas. Gestational age ranged from 29⁺⁵ to 40⁺⁴ weeks. There were 51 cases of premature rupture of membranes, 11 cases of placenta previa, 5 cases of placental abruption, 15 cases of placental adhesion/implantation and 27 cases of postpartum hemorrhage. There were 46 preterm fetuses,28 fetuses with fetal growth restriction, 22 fetuses with intrauterine distress, and 1 fetus with intrauterine death. Grossly, the placental lobules of succenturiate placentas had apparent size difference, while two lobules of bilobate placenta were more consistent. The chorionic plate size was smaller than the bottom plate of circumvallate placenta, the folded fetal membrane in the rim of placenta was thickened (termed marginate placenta if there was no thickening). The membranous placenta was characterized by a thin, large membrane-like shape. Annular placenta showed characteristic hollow cylinder, ring or horseshoe-shape. Fenestrated placenta was characterized by tissue defects near central area. Microscopically, functional/morphologic changes were the main manifestations of inadequate maternal-fetal perfusion, including villous infarction, distal villous dysplasia and excessive villous maturation. Conclusions: The abnormal shaped singleton placentas showed variable extent of inadequate maternal-fetal perfusion, which may lead to adverse pregnancy outcomes such as premature delivery, fetal growth restriction, intrauterine distress or fetal death.

Aegilops tauschii genome assembly Aet v5.0 features greater sequence contiguity and improved annotation

Aegilops tauschii is the donor of the D subgenome of hexaploid wheat and an important genetic resource. The reference-quality genome sequence Aet v4.0 for Ae. tauschii acc. AL8/78 was therefore an important milestone for wheat biology and breeding. Further advances in sequencing acc. AL8/78 and release of the Aet v5.0 sequence assembly are reported here. Two new optical maps were constructed and used in the revision of pseudomolecules. Gaps were closed with Pacific Biosciences long-read contigs, decreasing the gap number by 38,899. Transposable elements and protein-coding genes were reannotated. The number of annotated high-confidence genes was reduced from 39,635 in Aet v4.0 to 32,885 in Aet v5.0. A total of 2245 biologically important genes, including those affecting plant phenology, grain quality, and tolerance of abiotic stresses in wheat, was manually annotated and disease-resistance genes were annotated by a dedicated pipeline. Disease-resistance genes encoding nucleotide-binding site domains, receptor-like protein kinases, and receptor-like proteins were preferentially located in distal chromosome regions, whereas those encoding transmembrane coiled-coil proteins were dispersed more evenly along the chromosomes. Discovery, annotation, and expression analyses of microRNA (miRNA) precursors, mature miRNAs, and phasiRNAs are reported, including miRNA target genes. Other small RNAs, such as hc-siRNAs and tRFs, were characterized. These advances enhance the utility of the Ae. tauschii genome sequence for wheat genetics, biotechnology, and breeding.

Proteomic Determination of Low-Molecular-Weight Glutenin Subunit Composition in Aroona Near-Isogenic Lines and Standard Wheat Cultivars

The low-molecular weight glutenin subunit (LMW-GS) composition of wheat (Triticum aestivum) flour has important effects on end-use quality. However, assessing the contributions of each LMW-GS to flour quality remains challenging because of the complex LMW-GS composition and allelic variation among wheat cultivars. Therefore, accurate and reliable determination of LMW-GS alleles in germplasm remains an important challenge for wheat breeding. In this study, we used an optimized reversed-phase HPLC method and proteomics approach comprising 2-D gels coupled with liquid chromatography–tandem mass spectrometry (MS/MS) to discriminate individual LMW-GSs corresponding to alleles encoded by the Glu-A3, Glu-B3, and Glu-D3 loci in the ‘Aroona’ cultivar and 12 ‘Aroona’ near-isogenic lines (ARILs), which contain unique LMW-GS alleles in the same genetic background. The LMW-GS separation patterns for ‘Aroona’ and ARILs on chromatograms and 2-D gels were consistent with those from a set of 10 standard wheat cultivars for Glu-3. Furthermore, 12 previously uncharacterized spots in ‘Aroona’ and ARILs were excised from 2-D gels, digested with chymotrypsin, and subjected to MS/MS. We identified their gene haplotypes and created a 2-D gel map of LMW-GS alleles in the germplasm for breeding and screening for desirable LMW-GS alleles for wheat quality improvement.

Optical maps refine the bread wheat Triticum aestivum cv. Chinese Spring genome assembly

Until recently, achieving a reference-quality genome sequence for bread wheat was long thought beyond the limits of genome sequencing and assembly technology, primarily due to the large genome size and > 80% repetitive sequence content. The release of the chromosome scale 14.5-Gb IWGSC RefSeq v1.0 genome sequence of bread wheat cv. Chinese Spring (CS) was, therefore, a milestone. Here, we used a direct label and stain (DLS) optical map of the CS genome together with a prior nick, label, repair and stain (NLRS) optical map, and sequence contigs assembled with Pacific Biosciences long reads, to refine the v1.0 assembly. Inconsistencies between the sequence and maps were reconciled and gaps were closed. Gap filling and anchoring of 279 unplaced scaffolds increased the total length of pseudomolecules by 168 Mb (excluding Ns). Positions and orientations were corrected for 233 and 354 scaffolds, respectively, representing 10% of the genome sequence. The accuracy of the remaining 90% of the assembly was validated. As a result of the increased contiguity, the numbers of transposable elements (TEs) and intact TEs have increased in IWGSC RefSeq v2.1 compared with v1.0. In total, 98% of the gene models identified in v1.0 were mapped onto this new assembly through development of a dedicated approach implemented in the MAGAAT pipeline. The numbers of high-confidence genes on pseudomolecules have increased from 105 319 to 105 534. The reconciled assembly enhances the utility of the sequence for genetic mapping, comparative genomics, gene annotation and isolation, and more general studies on the biology of wheat.

Characterization of Pm68, a new powdery mildew resistance gene on chromosome 2BS of Greek durum wheat TRI 1796

New powdery mildew resistance gene Pm68 was found in the terminal region of chromosome 2BS of Greek durum wheat TRI 1796. The co-segregated molecular markers could be used for MAS. Durum wheat (Triticum turgidum L. var. durum Desf.) is not only an important cereal crop for pasta making, but also a genetic resource for common wheat improvement. In the present study, a Greek durum wheat TRI 1796 was found to confer high resistance to all 22 tested isolates of Blumeria graminis f. sp. tritici (Bgt). Inheritance study on the F₁ plants and the F₂ population derived from the cross TRI 1796/PI 584832 revealed that the resistance in TRI 1796 was controlled by a single dominant gene, herein designated Pm68. Using the bulked segregant RNA-Seq (BSR-Seq) analysis combined with molecular analysis, Pm68 was mapped to the terminal part of the short arm of chromosome 2B and flanked by markers Xdw04 and Xdw12/Xdw13 with genetic distances of 0.22 cM each. According to the reference genome of durum wheat cv. Svevo, the corresponding physical region spanned the Pm68 locus was about 1.78-Mb, in which a number of disease resistance-related genes were annotated. This study reports the new powdery mildew resistance gene Pm68 that would be a valuable resource for improvement of both common wheat and durum wheat. The co-segregated markers (Xdw05-Xdw11) developed here would be useful tools for marker-assisted selection (MAS) in breeding.

A high-quality reference genome sequence of Salvia miltiorrhiza provides insights into tanshinone synthesis in its red rhizomes

Salvia miltiorrhiza Bunge, also known as red sage or Danshen, is an important traditional Chinese medicine (TCM) that has been used for thousands of years to treat cardiovascular and other diseases. It is also considered an important model TCM plant. Here, a high-quality reference genome of S. miltiorrhiza was generated by combining PacBio long-read sequencing and chromatin interaction mapping (Hi-C) technologies, resulting in the chromosome-scale assembly of a 594.75-Mb genome sequence with a contig N50 of 2.70 Mb. This assembly shows the highest level of continuity for a Danshen genome generated thus far. The S. miltiorrhiza genome contained 32,483 protein-coding genes, with a repetitive DNA content of approximately 64.84%. The high percentage of young LTRs suggests that multiple TE transposition bursts occurred recently in S. miltiorrhiza. Genes unique to secondary metabolism pathways were expanded in the S. miltiorrhiza genome. A new CYP450 gene cluster was identified in the phloem of red roots where active components were synthesized. This reference genome sequence will facilitate future studies aimed at the elucidation of the secondary metabolism synthesis pathway and the genetic improvement of S. miltiorrhiza.

[Method of dietary nutritional status assessment and its application in cohort study of nutritional epidemiology]

Minimizing the burden on study subjects and assessing the general dietary nutritional status as accurately as possible are the basis of a nutritional epidemiological cohort study in the general population. While introducing the main dietary nutrition assessment methods, this paper manly describes the basic contents and principles for the development of food frequency questionnaire, and briefly illustrates the problems and solutions for the development of area specific food frequency questionnaires by taking the example of Tianjin Chronic Low-grade Systemic Inflammation and Health (TCLSIH) cohort study. Finally, discusses preliminarily the necessity and possibility of developing a national food frequency questionnaire.

[Prognosis analysis of 112 cases with Japanese encephalitis in adults]

Objective: To analyze the improvement of clinical symptoms and recovery of neurological function in adult Japanese encephalitis, and study the prognostic factors. Methods: Follow-up was conducted for 112 hospitalized patients with Japanese encephalitis (JE) in adults at the Department of Neurology of three hospitals in Gansu province from July to October 2016, from July to October 2017, 6 months and 1 year after onset, respectively. The neurological functional recovery was evaluated by modified Ranking Scale (mRS).The influencing factors were analyzed by logistic regression model. Results: Among the 112 adult patients with JE after 1year follow-up, 57% (64/112) were completely recovered (mRS score=0), and 14%(16/112) had mild neurological dysfunction (mRS score=1 or 2 points), 20% (22/112) had moderate to severe neurological dysfunction (mRS score 3 to 5), and 9% (10/112) died. In 102 survivors, decreased consciousness were fully recovered (100%), 75% of the mental and behavior disorders, 64% of cognitive/memory impairment, 71% of language function disorder, 61% of paralysis, 73% of extrapyramidal symptoms were fully recovered, and 92% of the seizures were controlled. Comparison of clinical data of initial on-set between good prognosis group (mRS score≤2, 80 cases) and poor prognosis group (mRS score>2, 32 cases) showed that initial clinical manifestation with seizures, consciousness (GCS score), cerebrospinal fluid pressure, and lesion of MRI involved in midbrain had statistically significant differences (all P<0.05) . Multivariate analysis demonstrated that cerebrospinal fluid (CSF) pressure>250 mmH(2)O and lesion of midbrain in MRI were independent risk factors of poor prognosis in adult patients with JE. Conclusion: JE is an acute and infectious viral encephalitis of the central nervous system with high disability and mortality. Most patients were completely recovered, and some had neurological sequelae. CSF pressure>250 mmH(2)O and lesion of midbrain in MRI are independent risk factors for poor prognosis.

GrainGenes: centralized small grain resources and digital platform for geneticists and breeders

GrainGenes (https://wheat.pw.usda.gov or https://graingenes.org) is an international centralized repository for curated, peer-reviewed datasets useful to researchers working on wheat, barley, rye and oat. GrainGenes manages genomic, genetic, germplasm and phenotypic datasets through a dynamically generated web interface for facilitated data discovery. Since 1992, GrainGenes has served geneticists and breeders in both the public and private sectors on six continents. Recently, several new datasets were curated into the database along with new tools for analysis. The GrainGenes homepage was enhanced by making it more visually intuitive and by adding links to commonly used pages. Several genome assemblies and genomic tracks are displayed through the genome browsers at GrainGenes, including the Triticum aestivum (bread wheat) cv. 'Chinese Spring' IWGSC RefSeq v1.0 genome assembly, the Aegilops tauschii (D genome progenitor) Aet v4.0 genome assembly, the Triticum turgidum ssp. dicoccoides (wild emmer wheat) cv. 'Zavitan' WEWSeq v.1.0 genome assembly, a T. aestivum (bread wheat) pangenome, the Hordeum vulgare (barley) cv. 'Morex' IBSC genome assembly, the Secale cereale (rye) select 'Lo7' assembly, a partial hexaploid Avena sativa (oat) assembly and the Triticum durum cv. 'Svevo' (durum wheat) RefSeq Release 1.0 assembly. New genetic maps and markers were added and can be displayed through CMAP. Quantitative trait loci, genetic maps and genes from the Wheat Gene Catalogue are indexed and linked through the Wheat Information System (WheatIS) portal. Training videos were created to help users query and reach the data they need. GSP (Genome Specific Primers) and PIECE2 (Plant Intron Exon Comparison and Evolution) tools were implemented and are available to use. As more small grains reference sequences become available, GrainGenes will play an increasingly vital role in helping researchers improve crops.

Comparison of MALDI-TOF-MS and RP-HPLC as Rapid Screening Methods for Wheat Lines With Altered Gliadin Compositions

The wheat gliadins are a complex group of flour proteins that can trigger celiac disease and serious food allergies. As a result, mutation breeding and biotechnology approaches are being used to develop new wheat lines with reduced immunogenic potential. Key to these efforts is the development of rapid, high-throughput methods that can be used as a first step in selecting lines with altered gliadin contents. In this paper, we optimized matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) and reversed-phase high-performance liquid chromatography (RP-HPLC) methods for the separation of gliadins from Triticum aestivum cv. Chinese Spring (CS). We evaluated the quality of the resulting profiles using the complete set of gliadin gene sequences recently obtained from this cultivar as well as a set of aneuploid lines in CS. The gliadins were resolved into 13 peaks by MALDI-TOF-MS. α- or γ-gliadins that contain abundant celiac disease epitopes and are likely targets for efforts to reduce the immunogenicity of flour were found in several peaks. However, other peaks contained multiple α- and γ-gliadins, including one peak with as many as 12 different gliadins. In comparison, separation of proteins by RP-HPLC yielded 28 gliadin peaks, including 13 peaks containing α-gliadins and eight peaks containing γ-gliadins. While the separation of α- and γ-gliadins gliadins achieved by RP-HPLC was better than that achieved by MALDI-TOF-MS, it was not possible to link peaks with individual protein sequences. Both MALDI-TOF-MS and RP-HPLC provided adequate separation of ω-gliadins. While MALDI-TOF-MS is faster and could prove useful in studies that target specific gliadins, RP-HPLC is an effective method that can be applied more broadly to detect changes in gliadin composition.

OrthoVenn2: a web server for whole-genome comparison and annotation of orthologous clusters across multiple species

OrthoVenn is a powerful web platform for the comparison and analysis of whole-genome orthologous clusters. Here we present an updated version, OrthoVenn2, which provides new features that facilitate the comparative analysis of orthologous clusters among up to 12 species. Additionally, this update offers improvements to data visualization and interpretation, including an occurrence pattern table for interrogating the overlap of each orthologous group for the queried species. Within the occurrence table, the functional annotations and summaries of the disjunctions and intersections of clusters between the chosen species can be displayed through an interactive Venn diagram. To facilitate a broader range of comparisons, a larger number of species, including vertebrates, metazoa, protists, fungi, plants and bacteria, have been added in OrthoVenn2. Finally, a stand-alone version is available to perform large dataset comparisons and to visualize results locally without limitation of species number. In summary, OrthoVenn2 is an efficient and user-friendly web server freely accessible at https://orthovenn2.bioinfotoolkits.net.

Exploiting the reference genome sequence of hexaploid wheat: a proteomic study of flour proteins from the cultivar Chinese Spring

Although the economic value of wheat flour is determined by the complement of gluten proteins, these proteins have been challenging to study because of the complexity of the major protein groups and the tremendous sequence diversity among wheat cultivars. The completion of a high-quality wheat genome sequence from the reference wheat Chinese Spring recently facilitated the assembly and annotation of a complete set of gluten protein genes from a single cultivar, making it possible to link individual proteins in the flour to specific gene sequences. In a proteomic analysis of total wheat flour protein from Chinese Spring using quantitative two-dimensional gel electrophoresis combined with tandem mass spectrometry, gliadins or low-molecular-weight glutenin subunits were identified as the predominant proteins in 72 protein spots. Individual spots were associated with 40 of 56 Chinese Spring gene sequences, including 16 of 26 alpha gliadins, 10 of 11 gamma gliadins, six of seven omega gliadins, one of two delta gliadins, and nine of ten LMW-GS. Most genes that were not associated with protein spots were either expressed at low levels in endosperm or encoded proteins with high similarity to other proteins. A wide range of protein accumulation levels were observed and discrepancies between transcript levels and protein levels were noted. This work together with similar studies using other commercial cultivars should provide new insight into the molecular basis of wheat flour quality and allergenic potential.

Rapid evolution of α-gliadin gene family revealed by analyzing Gli-2 locus regions of wild emmer wheat

α-Gliadins are a major group of gluten proteins in wheat flour that contribute to the end-use properties for food processing and contain major immunogenic epitopes that can cause serious health-related issues including celiac disease (CD). α-Gliadins are also the youngest group of gluten proteins and are encoded by a large gene family. The majority of the gene family members evolved independently in the A, B, and D genomes of different wheat species after their separation from a common ancestral species. To gain insights into the origin and evolution of these complex genes, the genomic regions of the Gli-2 loci encoding α-gliadins were characterized from the tetraploid wild emmer, a progenitor of hexaploid bread wheat that contributed the AABB genomes. Genomic sequences of Gli-2 locus regions for the wild emmer A and B genomes were first reconstructed using the genome sequence scaffolds along with optical genome maps. A total of 24 and 16 α-gliadin genes were identified for the A and B genome regions, respectively. α-Gliadin pseudogene frequencies of 86% for the A genome and 69% for the B genome were primarily caused by C to T substitutions in the highly abundant glutamine codons, resulting in the generation of premature stop codons. Comparison with the homologous regions from the hexaploid wheat cv. Chinese Spring indicated considerable sequence divergence of the two A genomes at the genomic level. In comparison, conserved regions between the two B genomes were identified that included α-gliadin pseudogenes containing shared nested TE insertions. Analyses of the genomic organization and phylogenetic tree reconstruction indicate that although orthologous gene pairs derived from speciation were present, large portions of α-gliadin genes were likely derived from differential gene duplications or deletions after the separation of the homologous wheat genomes ~ 0.5 MYA. The higher number of full-length intact α-gliadin genes in hexaploid wheat than that in wild emmer suggests that human selection through domestication might have an impact on α-gliadin evolution. Our study provides insights into the rapid and dynamic evolution of genomic regions harboring the α-gliadin genes in wheat.

Towards reducing the immunogenic potential of wheat flour: omega gliadins encoded by the D genome of hexaploid wheat may also harbor epitopes for the serious food allergy WDEIA

BACKGROUND

Omega-5 gliadins are a group of highly repetitive gluten proteins in wheat flour encoded on the 1B chromosome of hexaploid wheat. These proteins are the major sensitizing allergens in a severe form of food allergy called wheat-dependent exercise-induced anaphylaxis (WDEIA). The elimination of omega-5 gliadins from wheat flour through biotechnology or breeding approaches could reduce the immunogenic potential and adverse health effects of the flour.

RESULTS

A mutant line missing low-molecular weight glutenin subunits encoded at the Glu-B3 locus was selected previously from a doubled haploid population generated from two Korean wheat cultivars. Analysis of flour from the mutant line by 2-dimensional gel electrophoresis coupled with tandem mass spectrometry revealed that the omega-5 gliadins and several gamma gliadins encoded by the closely linked Gli-B1 locus were also missing as a result of a deletion of at least 5.8 Mb of chromosome 1B. Two-dimensional immunoblot analysis of flour proteins using sera from WDEIA patients showed reduced IgE reactivity in the mutant relative to the parental lines due to the absence of the major omega-5 gliadins. However, two minor proteins showed strong reactivity to patient sera in both the parental and the mutant lines and also reacted with a monoclonal antibody against omega-5 gliadin. Analysis of the two minor reactive proteins by mass spectrometry revealed that both proteins correspond to omega-5 gliadin genes encoded on chromosome 1D that were thought previously to be pseudogenes.

CONCLUSIONS

While breeding approaches can be used to reduce the levels of the highly immunogenic omega-5 gliadins in wheat flour, these approaches are complicated by the genetic linkage of different classes of gluten protein genes and the finding that omega-5 gliadins may be encoded on more than one chromosome. The work illustrates the importance of detailed knowledge about the genomic regions harboring the major gluten protein genes in individual wheat cultivars for future efforts aimed at reducing the immunogenic potential of wheat flour.

Reassessment of the evolution of wheat chromosomes 4A, 5A, and 7B

Comparison of genome sequences of wild emmer wheat and Aegilops tauschii suggests a novel scenario of the evolution of rearranged wheat chromosomes 4A, 5A, and 7B. Past research suggested that wheat chromosome 4A was subjected to a reciprocal translocation T(4AL;5AL)1 that occurred in the diploid progenitor of the wheat A subgenome and to three major rearrangements that occurred in polyploid wheat: pericentric inversion Inv(4AS;4AL)1, paracentric inversion Inv(4AL;4AL)1, and reciprocal translocation T(4AL;7BS)1. Gene collinearity along the pseudomolecules of tetraploid wild emmer wheat (Triticum turgidum ssp. dicoccoides, subgenomes AABB) and diploid Aegilops tauschii (genomes DD) was employed to confirm these rearrangements and to analyze the breakpoints. The exchange of distal regions of chromosome arms 4AS and 4AL due to pericentric inversion Inv(4AS;4AL)1 was detected, and breakpoints were validated with an optical Bionano genome map. Both breakpoints contained satellite DNA. The breakpoints of reciprocal translocation T(4AL;7BS)1 were also found. However, the breakpoints that generated paracentric inversion Inv(4AL;4AL)1 appeared to be collocated with the 4AL breakpoints that had produced Inv(4AS;4AL)1 and T(4AL;7BS)1. Inv(4AS;4AL)1, Inv(4AL;4AL)1, and T(4AL;7BS)1 either originated sequentially, and Inv(4AL;4AL)1 was produced by recurrent chromosome breaks at the same breakpoints that generated Inv(4AS;4AL)1 and T(4AL;7BS)1, or Inv(4AS;4AL)1, Inv(4AL;4AL)1, and T(4AL;7BS)1 originated simultaneously. We prefer the latter hypothesis since it makes fewer assumptions about the sequence of events that produced these chromosome rearrangements.

Towards high-energy, high-resolution computed tomography via a laser driven micro-spot gamma-ray source

Computed Tomography (CT) is a powerful method for non-destructive testing (NDT) and metrology awakes with expanding application fields. To improve the spatial resolution of high energy CT, a micro-spot gamma-ray source based on bremsstrahlung from a laser wakefield accelerator was developed. A high energy CT using the source was performed, which shows that the resolution of reconstruction can reach 100 μm at 10% contrast. Our proof-of-principle demonstration indicates that laser driven micro-spot gamma-ray sources provide a prospective way to increase the spatial resolution and toward to high energy micro CT. Due to the advantage in spatial resolution, laser based high energy CT represents a large potential for many NDT applications.

Analysis of the Gli-D2 locus identifies a genetic target for simultaneously improving the breadmaking and health-related traits of common wheat

Gliadins are a major component of wheat seed proteins. However, the complex homoeologous Gli-2 loci (Gli-A2, -B2 and -D2) that encode the α-gliadins in commercial wheat are still poorly understood. Here we analyzed the Gli-D2 locus of Xiaoyan 81 (Xy81), a winter wheat cultivar. A total of 421.091 kb of the Gli-D2 sequence was assembled from sequencing multiple bacterial artificial clones, and 10 α-gliadin genes were annotated. Comparative genomic analysis showed that Xy81 carried only eight of the α-gliadin genes of the D genome donor Aegilops tauschii, with two of them each experiencing a tandem duplication. A mutant line lacking Gli-D2 (DLGliD2) consistently exhibited better breadmaking quality and dough functionalities than its progenitor Xy81, but without penalties in other agronomic traits. It also had an elevated lysine content in the grains. Transcriptome analysis verified the lack of Gli-D2 α-gliadin gene expression in DLGliD2. Furthermore, the transcript and protein levels of protein disulfide isomerase were both upregulated in DLGliD2 grains. Consistent with this finding, DLGliD2 had increased disulfide content in the flour. Our work sheds light on the structure and function of Gli-D2 in commercial wheat, and suggests that the removal of Gli-D2 and the gliadins specified by it is likely to be useful for simultaneously enhancing the end-use and health-related traits of common wheat. Because gliadins and homologous proteins are widely present in grass species, the strategy and information reported here may be broadly useful for improving the quality traits of diverse cereal crops.

Structural variation and rates of genome evolution in the grass family seen through comparison of sequences of genomes greatly differing in size

Homology was searched with genes annotated in the Aegilops tauschii pseudomolecules against genes annotated in the pseudomolecules of tetraploid wild emmer wheat, Brachypodium distachyon, sorghum and rice. Similar searches were performed with genes annotated in the rice pseudomolecules. Matrices of collinear genes and rearrangements in their order were constructed. Optical BioNano genome maps were constructed and used to validate rearrangements unique to the wild emmer and Ae. tauschii genomes. Most common rearrangements were short paracentric inversions and short intrachromosomal translocations. Intrachromosomal translocations outnumbered segmental intrachromosomal duplications. The densities of paracentric inversion lengths were approximated by exponential distributions in all six genomes. Densities of collinear genes along the Ae. tauschii chromosomes were highly correlated with meiotic recombination rates but those of rearrangements were not, suggesting different causes of the erosion of gene collinearity and evolution of major chromosome rearrangements. Frequent rearrangements sharing breakpoints suggested that chromosomes have been rearranged recurrently at some sites. The distal 4 Mb of the short arms of rice chromosomes Os11 and Os12 and corresponding regions in the sorghum, B. distachyon and Triticeae genomes contain clusters of interstitial translocations including from 1 to 7 collinear genes. The rates of acquisition of major rearrangements were greater in the large wild emmer wheat and Ae. tauschii genomes than in the lineage preceding their divergence or in the B. distachyon, rice and sorghum lineages. It is suggested that synergy between large quantities of dynamic transposable elements and annual growth habit have been the primary causes of the fast evolution of the Triticeae genomes.

Experimental Evidence of Kinetic Effects in Indirect-Drive Inertial Confinement Fusion Hohlraums

We present the first experimental evidence supported by simulations of kinetic effects launched in the interpenetration layer between the laser-driven hohlraum plasma bubbles and the corona plasma of the compressed pellet at the Shenguang-III prototype laser facility. Solid plastic capsules were coated with carbon-deuterium layers; as the implosion neutron yield is quenched, DD fusion yield from the corona plasma provides a direct measure of the kinetic effects inside the hohlraum. An anomalous large energy spread of the DD neutron signal (∼282  keV) and anomalous scaling of the neutron yield with the thickness of the carbon-deuterium layers cannot be explained by the hydrodynamic mechanisms. Instead, these results can be attributed to kinetic shocks that arise in the hohlraum-wall-ablator interpenetration region, which result in efficient acceleration of the deuterons (∼28.8  J, 0.45% of the total input laser energy). These studies provide novel insight into the interactions and dynamics of a vacuum hohlraum and near-vacuum hohlraum.

Modular assembly of transposable element arrays by microsatellite targeting in the guayule and rice genomes

Background

Guayule (Parthenium argentatum A. Gray) is a rubber-producing desert shrub native to Mexico and the United States. Guayule represents an alternative to Hevea brasiliensis as a source for commercial natural rubber. The efficient application of modern molecular/genetic tools to guayule improvement requires characterization of its genome.

Results

The 1.6 Gb guayule genome was sequenced, assembled and annotated. The final 1.5 Gb assembly, while fragmented (N₅₀ = 22 kb), maps > 95% of the shotgun reads and is essentially complete. Approximately 40,000 transcribed, protein encoding genes were annotated on the assembly. Further characterization of this genome revealed 15 families of small, microsatellite-associated, transposable elements (TEs) with unexpected chromosomal distribution profiles. These SaTar (Satellite Targeted) elements, which are non-autonomous Mu-like elements (MULEs), were frequently observed in multimeric linear arrays of unrelated individual elements within which no individual element is interrupted by another. This uniformly non-nested TE multimer architecture has not been previously described in either eukaryotic or prokaryotic genomes. Five families of similarly distributed non-autonomous MULEs (microsatellite associated, modularly assembled) were characterized in the rice genome. Families of TEs with similar structures and distribution profiles were identified in sorghum and citrus.

Conclusion

The sequencing and assembly of the guayule genome provides a foundation for application of current crop improvement technologies to this plant. In addition, characterization of this genome revealed SaTar elements with distribution profiles unique among TEs. Satar targeting appears based on an alternative MULE recombination mechanism with the potential to impact gene evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4653-6) contains supplementary material, which is available to authorized users.

Dynamic Evolution of α-Gliadin Prolamin Gene Family in Homeologous Genomes of Hexaploid Wheat

Wheat Gli-2 loci encode complex groups of α-gliadin prolamins that are important for breadmaking, but also major triggers of celiac disease (CD). Elucidation of α-gliadin evolution provides knowledge to produce wheat with better end-use properties and reduced immunogenic potential. The Gli-2 loci contain a large number of tandemly duplicated genes and highly repetitive DNA, making sequence assembly of their genomic regions challenging. Here, we constructed high-quality sequences spanning the three wheat homeologous α-gliadin loci by aligning PacBio-based sequence contigs with BioNano genome maps. A total of 47 α-gliadin genes were identified with only 26 encoding intact full-length protein products. Analyses of α-gliadin loci and phylogenetic tree reconstruction indicate significant duplications of α-gliadin genes in the last ~2.5 million years after the divergence of the A, B and D genomes, supporting its rapid lineage-independent expansion in different Triticeae genomes. We showed that dramatic divergence in expression of α-gliadin genes could not be attributed to sequence variations in the promoter regions. The study also provided insights into the evolution of CD epitopes and identified a single indel event in the hexaploid wheat D genome that likely resulted in the generation of the highly toxic 33-mer CD epitope.

Gene Duplication and Evolution Dynamics in the Homeologous Regions Harboring Multiple Prolamin and Resistance Gene Families in Hexaploid Wheat

Improving end-use quality and disease resistance are important goals in wheat breeding. The genetic loci controlling these traits are highly complex, consisting of large families of prolamin and resistance genes with members present in all three homeologous A, B, and D genomes in hexaploid bread wheat. Here, orthologous regions harboring both prolamin and resistance gene loci were reconstructed and compared to understand gene duplication and evolution in different wheat genomes. Comparison of the two orthologous D regions from the hexaploid wheat Chinese Spring and the diploid progenitor Aegilops tauschii revealed their considerable difference due to the presence of five large structural variations with sizes ranging from 100 kb to 2 Mb. As a result, 44% of the Ae. tauschii and 71% of the Chinese Spring sequences in the analyzed regions, including 79 genes, are not shared. Gene rearrangement events, including differential gene duplication and deletion in the A, B, and D regions, have resulted in considerable erosion of gene collinearity in the analyzed regions, suggesting rapid evolution of prolamin and resistance gene families after the separation of the three wheat genomes. We hypothesize that this fast evolution is attributed to the co-evolution of the two gene families dispersed within a high recombination region. The identification of a full set of prolamin genes facilitated transcriptome profiling and revealed that the A genome contributes the least to prolamin expression because of its smaller number of expressed intact genes and their low expression levels, while the B and D genomes contribute similarly.

New insights into structural organization and gene duplication in a 1.75-Mb genomic region harboring the α-gliadin gene family in Aegilops tauschii, the source of wheat D genome

Among the wheat prolamins important for its end-use traits, α-gliadins are the most abundant, and are also a major cause of food-related allergies and intolerances. Previous studies of various wheat species estimated that between 25 and 150 α-gliadin genes reside in the Gli-2 locus regions. To better understand the evolution of this complex gene family, the DNA sequence of a 1.75-Mb genomic region spanning the Gli-2 locus was analyzed in the diploid grass, Aegilops tauschii, the ancestral source of D genome in hexaploid bread wheat. Comparison with orthologous regions from rice, sorghum, and Brachypodium revealed rapid and dynamic changes only occurring to the Ae. tauschii Gli-2 region, including insertions of high numbers of non-syntenic genes and a high rate of tandem gene duplications, the latter of which have given rise to 12 copies of α-gliadin genes clustered within a 550-kb region. Among them, five copies have undergone pseudogenization by various mutation events. Insights into the evolutionary relationship of the duplicated α-gliadin genes were obtained from their genomic organization, transcription patterns, transposable element insertions and phylogenetic analyses. An ancestral glutamate-like receptor (GLR) gene encoding putative amino acid sensor in all four grass species has duplicated only in Ae. tauschii and generated three more copies that are interspersed with the α-gliadin genes. Phylogenetic inference and different gene expression patterns support functional divergence of the Ae. tauschii GLR copies after duplication. Our results suggest that the duplicates of α-gliadin and GLR genes have likely taken different evolutionary paths; conservation for the former and neofunctionalization for the latter.

Femtosecond Probing of Plasma Wakefields and Observation of the Plasma Wake Reversal Using a Relativistic Electron Bunch

We show that a high-energy electron bunch can be used to capture the instantaneous longitudinal and transverse field structures of the highly transient, microscopic, laser-excited relativistic wake with femtosecond resolution. The spatiotemporal evolution of wakefields in a plasma density up ramp is measured and the reversal of the plasma wake, where the wake wavelength at a particular point in space increases until the wake disappears completely only to reappear at a later time but propagating in the opposite direction, is observed for the first time by using this new technique.

Brilliant petawatt gamma-ray pulse generation in quantum electrodynamic laser-plasma interaction

We show a new resonance acceleration scheme for generating ultradense relativistic electron bunches in helical motions and hence emitting brilliant vortical γ-ray pulses in the quantum electrodynamic (QED) regime of circularly-polarized (CP) laser-plasma interactions. Here the combined effects of the radiation reaction recoil force and the self-generated magnetic fields result in not only trapping of a great amount of electrons in laser-produced plasma channel, but also significant broadening of the resonance bandwidth between laser frequency and that of electron betatron oscillation in the channel, which eventually leads to formation of the ultradense electron bunch under resonant helical motion in CP laser fields. Three-dimensional PIC simulations show that a brilliant γ-ray pulse with unprecedented power of 6.7 PW and peak brightness of 10²⁵ photons/s/mm²/mrad²/0.1% BW (at 15 MeV) is emitted at laser intensity of 1.9 × 10²³ W/cm².

Large-charge quasimonoenergetic electron beams produced by off-axis colliding laser pulses in underdense plasma

Electrons can be efficiently injected into a plasma wave by colliding two counterpropagating laser pulses in a laser wakefield acceleration. However, the generation of a high-quality electron beam with a large charge is difficult in the traditional on-axis colliding scheme due to the growth of the electron beam duration coming from the increase of the beam charge. To solve this problem, we propose an off-axis colliding scheme, in which the collision point is away from the axis of the driver pulse. We show that the electrons injected from the off-axis region are highly concentered on the tail of the bubble even for a large trapped charge, thus feeling almost the same accelerating field. As a result, quasimonoenergetic electron beams with a large charge can be produced. The validity of this scheme is confirmed by both the particle-in-cell simulations and the Hamiltonian model. Furthermore, it is shown that a Laguerre-Gauss (LG) laser can be adopted as the injection pulse to realize the off-axis colliding injection in three dimensions symmetrically, which may be useful in simplifying the technical layout of the real experiment setup.

Genome resources for climate-resilient cowpea, an essential crop for food security

Cowpea (Vigna unguiculata L. Walp.) is a legume crop that is resilient to hot and drought-prone climates, and a primary source of protein in sub-Saharan Africa and other parts of the developing world. However, genome resources for cowpea have lagged behind most other major crops. Here we describe foundational genome resources and their application to the analysis of germplasm currently in use in West African breeding programs. Resources developed from the African cultivar IT97K-499-35 include a whole-genome shotgun (WGS) assembly, a bacterial artificial chromosome (BAC) physical map, and assembled sequences from 4355 BACs. These resources and WGS sequences of an additional 36 diverse cowpea accessions supported the development of a genotyping assay for 51 128 SNPs, which was then applied to five bi-parental RIL populations to produce a consensus genetic map containing 37 372 SNPs. This genetic map enabled the anchoring of 100 Mb of WGS and 420 Mb of BAC sequences, an exploration of genetic diversity along each linkage group, and clarification of macrosynteny between cowpea and common bean. The SNP assay enabled a diversity analysis of materials from West African breeding programs. Two major subpopulations exist within those materials, one of which has significant parentage from South and East Africa and more diversity. There are genomic regions of high differentiation between subpopulations, one of which coincides with a cluster of nodulin genes. The new resources and knowledge help to define goals and accelerate the breeding of improved varieties to address food security issues related to limited-input small-holder farming and climate stress.

Physical Mechanism of the Transverse Instability in Radiation Pressure Ion Acceleration

The transverse stability of the target is crucial for obtaining high quality ion beams using the laser radiation pressure acceleration (RPA) mechanism. In this Letter, a theoretical model and supporting two-dimensional (2D) particle-in-cell (PIC) simulations are presented to clarify the physical mechanism of the transverse instability observed in the RPA process. It is shown that the density ripples of the target foil are mainly induced by the coupling between the transverse oscillating electrons and the quasistatic ions, a mechanism similar to the oscillating two stream instability in the inertial confinement fusion research. The predictions of the mode structure and the growth rates from the theory agree well with the results obtained from the PIC simulations in various regimes, indicating the model contains the essence of the underlying physics of the transverse breakup of the target.

Improved spatial calibration for the CXRS system on EAST

A Charge eXchange Recombination Spectroscopy (CXRS) diagnostic system has been developed to measure profiles of ion temperature and rotation since 2014 on EAST. Several techniques have been developed to improve the spatial calibration of the CXRS diagnostic. The sightline location was obtained by measuring the coordinates of three points on each sightline using an articulated flexible coordinate measuring arm when the vessel was accessible. After vacuum pumping, the effect of pressure change in the vacuum vessel was evaluated by observing the movement of the light spot from back-illuminated sightlines on the first wall using the newly developed articulated inspection arm. In addition, the rotation of the periscope after vacuum pumping was derived by using the Doppler shift of neutral beam emission spectra without magnetic field. Combining these techniques, improved spatial calibration was implemented to provide a complete and accurate description of the EAST CXRS system. Due to the effects of the change of air pressure, a ∼0.4° periscope rotation, yielding a ∼20 mm movement of the major radius of observation positions to the lower field side, was derived. Results of Zeeman splitting of neutral beam emission spectra with magnetic field also showed good agreement with the calibration results.

PIECE 2.0: an update for the plant gene structure comparison and evolution database

PIECE (Plant Intron Exon Comparison and Evolution) is a web-accessible database that houses intron and exon information of plant genes. PIECE serves as a resource for biologists interested in comparing intron–exon organization and provides valuable insights into the evolution of gene structure in plant genomes. Recently, we updated PIECE to a new version, PIECE 2.0 (http://probes.pw.usda.gov/piece or http://aegilops.wheat.ucdavis.edu/piece). PIECE 2.0 contains annotated genes from 49 sequenced plant species as compared to 25 species in the previous version. In the current version, we also added several new features: (i) a new viewer was developed to show phylogenetic trees displayed along with the structure of individual genes; (ii) genes in the phylogenetic tree can now be also grouped according to KOG (The annotation of Eukaryotic Orthologous Groups) and KO (KEGG Orthology) in addition to Pfam domains; (iii) information on intronless genes are now included in the database; (iv) a statistical summary of global gene structure information for each species and its comparison with other species was added; and (v) an improved GSDraw tool was implemented in the web server to enhance the analysis and display of gene structure. The updated PIECE 2.0 database will be a valuable resource for the plant research community for the study of gene structure and evolution.

Sequencing and comparative analyses of Aegilops tauschii chromosome arm 3DS reveal rapid evolution of Triticeae genomes

Bread wheat (Triticum aestivum, AABBDD) is an allohexaploid species derived from two rounds of interspecific hybridizations. A high-quality genome sequence assembly of diploid Aegilops tauschii, the donor of the wheat D genome, will provide a useful platform to study polyploid wheat evolution. A combined approach of BAC pooling and next-generation sequencing technology was employed to sequence the minimum tiling path (MTP) of 3176 BAC clones from the short arm of Ae. tauschii chromosome 3 (At3DS). The final assembly of 135 super-scaffolds with an N50 of 4.2 Mb was used to build a 247-Mb pseudomolecule with a total of 2222 predicted protein-coding genes. Compared with the orthologous regions of rice, Brachypodium, and sorghum, At3DS contains 38.67% more genes. In comparison to At3DS, the short arm sequence of wheat chromosome 3B (Ta3BS) is 95-Mb large in size, which is primarily due to the expansion of the non-centromeric region, suggesting that transposable element (TE) bursts in Ta3B likely occurred there. Also, the size increase is accompanied by a proportional increase in gene number in Ta3BS. We found that in the sequence of short arm of wheat chromosome 3D (Ta3DS), there was only less than 0.27% gene loss compared to At3DS. Our study reveals divergent evolution of grass genomes and provides new insights into sequence changes in the polyploid wheat genome.

Identification of new single nucleotide polymorphisms affecting total number born and candidate genes related to ovulation rate in Chinese Erhualian pigs

The Chinese Erhualian pig has the highest record for litter size in the world. However, the genetic mechanism of its high prolificacy remains poorly understood. In our study, large phenotypic variations in litter size were found among Erhualian sows. Significant differences in total number born (TNB) and corpora lutea numbers were observed between sows with high and low estimated breeding values (EBVs) for TNB. To identify single nucleotide polymorphisms (SNPs) associated with TNB, a selective genomic scan was conducted on 18 sows representing the top 10% and 18 sows representing the bottom 10% of EBVs of 177 sows using Illumina Porcine SNP60 genotype data. Genome-wide fixation coefficient (F_ST ) values were calculated for each SNP between the high- and low-EBV groups. A total of 154 SNPs were significantly differentiated loci between the two groups. Of the top 10 highest F_ST SNPs, rs81399474, rs81400131 and rs81405013 on SSC8 and rs81434499 and rs81434489 on SSC 12 corresponded to previously reported QTL for litter size. The other five SNPs, rs81367039 on SSC2, rs80891106 on SSC7, rs81477883 on SSC12 and rs80938898 and rs80971725 on SSC14, appeared to be novel QTL for TNB. Significant associations between rs81399474 on SSC8 and TNB were confirmed in 313 Erhualian sows. Forty genes were identified around the top 10 highest F_ST SNPs, of which UCHL1, adjacent to rs81399474, and RPS6KB1 and CLTC, adjacent to rs81434499, have been reported to affect the ovulation rate in pig. The findings can advance understanding of the genetic variations in litter size of pigs.

Rapid evolutionary dynamics in a 2.8-Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi-gene families, the DNA sequence of a 2.8-Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance-like gene families, was analyzed in the diploid grass Aegilops tauschii, the D-genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non-syntenic genes and only 13 ancestral genes are shared among these grass species. These non-syntenic genes, including prolamin and resistance-like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non-syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin-related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance-like genes and subsequent differential expansion of each R gene family. The high frequency of non-syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8-Mb region with nine-fold higher than the genome-wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.

Fast-forward genetics by radiation hybrids to saturate the locus regulating nuclear-cytoplasmic compatibility in Triticum

The nuclear-encoded species cytoplasm specific (scs) genes control nuclear-cytoplasmic compatibility in wheat (genus Triticum). Alloplasmic cells, which have nucleus and cytoplasm derived from different species, produce vigorous and vital organisms only when the correct version of scs is present in their nucleus. In this study, bulks of in vivo radiation hybrids segregating for the scs phenotype have been genotyped by sequencing with over 1.9 million markers. The high marker saturation obtained for a critical region of chromosome 1D allowed identification of 3318 reads that mapped in close proximity of the scs. A novel in silico approach was deployed to extend these short reads to sequences of up to 70 Kb in length and identify candidate open reading frames (ORFs). Markers were developed to anchor the short contigs containing ORFs to a radiation hybrid map of 650 individuals with resolution of 288 Kb. The region containing the scs locus was narrowed to a single Bacterial Artificial Chromosome (BAC) contig of Aegilops tauschii. Its sequencing and assembly by nano-mapping allowed rapid identification of a rhomboid gene as the only ORF existing within the refined scs locus. Resequencing of this gene from multiple germplasm sources identified a single nucleotide mutation, which gives rise to a functional amino acid change. Gene expression characterization revealed that an active copy of this rhomboid exists on all homoeologous chromosomes of wheat, and depending on the specific cytoplasm each copy is preferentially expressed. Therefore, a new methodology was applied to unique genetic stocks to rapidly identify a strong candidate gene for the control of nuclear-cytoplasmic compatibility in wheat.

Glutamine synthetase in Durum Wheat: Genotypic Variation and Relationship with Grain Protein Content

Grain protein content (GPC), is one of the most important trait in wheat and its characterized by a very complex genetic control. The identification of wheat varieties with high GPC (HGPC), as well as the characterization of central enzymes involved in these processes, are important for more sustainable agricultural practices. In this study, we focused on Glutamine synthetase (GS) as a candidate to study GPC in wheat. We analyzed GS expression and its enzymatic activity in different tissues and phenological stages in 10 durum wheat genotypes with different GPC. Although each genotype performed quite differently from the others, both because their genetic variability and their adaptability to specific environmental conditions, the highest GS activity and expression were found in genotypes with HGPC and vice versa the lowest ones in genotypes with low GPC (LGPC). Moreover, in genotypes contrasting in GPC bred at different nitrogen regimes (0, 60, 140 N Unit/ha) GS behaved differently in diverse organs. Nitrogen supplement increased GS expression and activity in roots of all genotypes, highlighting the key role of this enzyme in nitrogen assimilation and ammonium detoxification in roots. Otherwise, nitrogen treatments decreased GS expression and activity in the leaves of HGPC genotypes and did not affect GS in the leaves of LGPC genotypes. Finally, no changes in GS and soluble protein content occurred at the filling stage in the caryopses of all analyzed genotypes.

Capturing relativistic wakefield structures in plasmas using ultrashort high-energy electrons as a probe

A new method capable of capturing coherent electric field structures propagating at nearly the speed of light in plasma with a time resolution as small as a few femtoseconds is proposed. This method uses a few femtoseconds long relativistic electron bunch to probe the wake produced in a plasma by an intense laser pulse or an ultra-short relativistic charged particle beam. As the probe bunch traverses the wake, its momentum is modulated by the electric field of the wake, leading to a density variation of the probe after free-space propagation. This variation of probe density produces a snapshot of the wake that can directly give many useful information of the wake structure and its evolution. Furthermore, this snapshot allows detailed mapping of the longitudinal and transverse components of the wakefield. We develop a theoretical model for field reconstruction and verify it using 3-dimensional particle-in-cell (PIC) simulations. This model can accurately reconstruct the wakefield structure in the linear regime, and it can also qualitatively map the major features of nonlinear wakes. The capturing of the injection in a nonlinear wake is demonstrated through 3D PIC simulations as an example of the application of this new method.

GSP: a web-based platform for designing genome-specific primers in polyploids

MOTIVATION

The sequences among subgenomes in a polyploid species have high similarity, making it difficult to design genome-specific primers for sequence analysis.

RESULTS

We present GSP, a web-based platform to design genome-specific primers that distinguish subgenome sequences in a polyploid genome. GSP uses BLAST to extract homeologous sequences of the subgenomes in existing databases, performs a multiple sequence alignment, and design primers based on sequence variants in the alignment. An interactive primers diagram, a sequence alignment viewer and a virtual electrophoresis are displayed as parts of the primer design result. GSP also designs specific primers from multiple sequences uploaded by users.

AVAILABILITY AND IMPLEMENTATION

GSP is a user-friendly and efficient web platform freely accessible at http://probes.pw.usda.gov/GSP Source code and command-line application are available at https://github.com/bioinfogenome/GSP CONTACTS: yong.gu@ars.usda.gov or devin.coleman-derr@ars.usda.gov

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

WheatExp: an RNA-seq expression database for polyploid wheat

BACKGROUND

For functional genomics studies, it is important to understand the dynamic expression profiles of transcribed genes in different tissues, stages of development and in response to environmental stimuli. The proliferation in the use of next-generation sequencing technologies by the plant research community has led to the accumulation of large volumes of expression data. However, analysis of these datasets is complicated by the frequent occurrence of polyploidy among economically-important crop species. In addition, processing and analyzing such large volumes of sequence data is a technical and time-consuming task, limiting their application in functional genomics studies, particularly for smaller laboratories which lack access to high-powered computing infrastructure. Wheat is a good example of a young polyploid species with three similar genomes (97 % identical among homoeologous genes), rapidly accumulating RNA-seq datasets and a large research community.

DESCRIPTION

We present WheatExp, an expression database and visualization tool to analyze and compare homoeologue-specific transcript profiles across a broad range of tissues from different developmental stages in polyploid wheat. Beginning with publicly-available RNA-seq datasets, we developed a pipeline to distinguish between homoeologous transcripts from annotated genes in tetraploid and hexaploid wheat. Data from multiple studies is processed and compiled into a database which can be queried either by BLAST or by searching for a known gene of interest by name or functional domain. Expression data of multiple genes can be displayed side-by-side across all expression datasets providing immediate access to a comprehensive panel of expression data for specific subsets of wheat genes.

CONCLUSIONS

The development of a publicly accessible expression database hosted on the GrainGenes website - http://wheat.pw.usda.gov/WheatExp/ - coupled with a simple and readily-comparable visualization tool will empower the wheat research community to use RNA-seq data and to perform functional analyses of target genes. The presented expression data is homoeologue-specific allowing for the analysis of relative contributions from each genome to the overall expression of a gene, a critical consideration for breeding applications. Our approach can be expanded to other polyploid species by adjusting sequence mapping parameters according to the specific divergence of their genomes.