BatchPrimer3: a high throughput web application for PCR and sequencing primer design

Multiplex profiling of developmental cis-regulatory elements with quantitative single-cell expression reporters

The inability to scalably and precisely measure the activity of developmental cis-regulatory elements (CREs) in multicellular systems is a bottleneck in genomics. Here we develop a dual RNA cassette that decouples the detection and quantification tasks inherent to multiplex single-cell reporter assays. The resulting measurement of reporter expression is accurate over multiple orders of magnitude, with a precision approaching the limit set by Poisson counting noise. Together with RNA barcode stabilization via circularization, these scalable single-cell quantitative expression reporters provide high-contrast readouts, analogous to classic in situ assays but entirely from sequencing. Screening >200 regions of accessible chromatin in a multicellular in vitro model of early mammalian development, we identify 13 (8 previously uncharacterized) autonomous and cell-type-specific developmental CREs. We further demonstrate that chimeric CRE pairs generate cognate two-cell-type activity profiles and assess gain- and loss-of-function multicellular expression phenotypes from CRE variants with perturbed transcription factor binding sites. Single-cell quantitative expression reporters can be applied in developmental and multicellular systems to quantitatively characterize native, perturbed and synthetic CREs at scale, with high sensitivity and at single-cell resolution.

Comparative chloroplast genomics, phylogenetic relationships and molecular markers development of Aglaonema commutatum and seven green cultivars of Aglaonema

Aglaonema commutatum is a famous species in the Aglaonema genus, which has important ornamental and economic value. However, its chloroplast genome information and phylogenetic relationships among popular green cultivars of Aglaonema in southern China have not been reported. Herein, chloroplast genomes of one variety of A. commutatum and seven green cultivars of Aglaonema, namely, A. commutatum 'San Remo', 'Kai Sa', 'Pattaya Beauty', 'Sapphire', 'Silver Queen', 'Snow White', 'White Gem', and 'White Horse Prince', were sequenced and assembled for comparative analysis and phylogeny. These eight genomes possessed a typical quadripartite structure that consisted of a LSC region (90,799-91,486 bp), an SSC region (20,508-21,137 bp) and a pair of IR regions (26,661-26,750 bp). Each genome contained 112 different genes, comprising 79 protein-coding genes, 29 tRNA genes and 4 rRNA genes. The gene orders, GC contents, codon usage frequency, and IR/SC boundaries were highly conserved among these eight genomes. Long repeats, SSRs, SNPs and indels were analyzed among these eight genomes. Comparative analysis of 15 Aglaonema chloroplast genomes identified 7 highly variable regions, including trnH-GUG-exon1-psbA, trnS-GCU-trnG-UCC-exon1, trnY-GUA-trnE-UUC, psbC-trnS-UGA, trnF-GAA-ndhJ, ccsA-ndhD, and rps15-ycf1-D2. Reconstruction of the phylogenetic trees based on chloroplast genomes, strongly supported that Aglaonema was a sister to Anchomanes, and that the Aglaonema genus was classified into two sister clades including clade I and clade II, which corresponded to two sections, Aglaonema and Chamaecaulon, respectively. One variety and five cultivars, including A. commutatum 'San Remo', 'Kai Sa', 'Pattaya Beauty', 'Silver Queen', 'Snow White', and 'White Horse Prince', were classified into clade I; and the rest of the two cultivars, including 'Sapphire' and 'White Gem', were classified into clade II. Positive selection was observed in 34 protein-coding genes at the level of the amino acid sites among 77 chloroplast genomes of the Araceae family. Based on the highly variable regions and SSRs, 4 DNA markers were developed to differentiate the clade I and clade II in Aglaonema. In conclusion, this study provided chloroplast genomic resources for Aglaonema, which were useful for its classification and phylogeny.

Most Tibetan weedy barleys originated via recombination between Btr1 and Btr2 in domesticated barley

Tibetan weedy barleys reside at the edges of qingke (hulless barley) fields in Tibet (Xizang). The spikes of these weedy barleys contain or lack a brittle rachis, with either two- or six-rowed spikes and either hulled or hulless grains at maturity. Although the brittle rachis trait of Tibetan weedy barleys is similar to that of wild barley (Hordeum vulgare ssp. spontaneum Thell.), these plants share genetic similarity with domesticated barley. The origin of Tibetan weedy barleys continues to be debated. Here, we show that most Tibetan weedy barleys originated from cross-pollinated hybridization of domesticated barleys, followed by hybrid self-pollination and recombination between Non-brittle rachis 1 (btr1) and 2 (btr2). We discovered the specific genetic ancestry of these weedy barleys in South Asian accessions. Tibetan weedy barleys exhibit lower genetic diversity than wild and Chinese landraces/cultivars and share a close relationship with qingke, genetically differing from typical eastern and western barley populations. We classified Tibetan weedy barleys into two groups, brittle rachis (BR) and non-brittle rachis (NBR); these traits align with the haplotypes of the btr1 and btr2 genes. Whereas wild barleys carry haplotype combinations of Btr1 and Btr2, each showing lower proportions in a population, the recombinant haplotype BTR2H8+BTR1H24 is predominant in the BR group. Haplotype block analysis based on whole-genome sequencing revealed two recombination breakpoints, which are present in 80.6% and 16.8% of BR accessions according to marker-assisted analysis. Hybridization events between wild and domesticated barley were rarely detected. These findings support the notion that Tibetan weedy barleys originated via recombination between Btr1 and Btr2 in domesticated barley.

A cluster of putative resistance genes is associated with a dominant resistance to sunflower broomrape

KEY MESSAGE

The HaOr5 resistance gene is located in a large genomic insertion containing putative resistance genes and provides resistance to O. cumana, preventing successful connection to the sunflower root vascular system. Orobanche cumana (sunflower broomrape) is a parasitic plant that is part of the Orobanchaceae family and specifically infests sunflower crops. This weed is an obligate parasitic plant that does not carry out photosynthetic activity or develop roots and is fully dependent on its host for its development. It produces thousands of dust-like seeds per plant. It possesses a high spreading ability and has been shown to quickly overcome resistance genes successively introduced by selection in cultivated sunflower varieties. The first part of its life cycle occurs underground. The connection to the sunflower vascular system is essential for parasitic plant survival and development. The HaOr5 gene provides resistance to sunflower broomrape race E by preventing the connection of O. cumana to the root vascular system. We mapped a single position of the HaOr5 gene by quantitative trait locus mapping using two segregating populations. The same location of the HaOr5 gene was identified by genome-wide association. Using a large population of thousands of F2 plants, we restricted the location of the HaOr5 gene to a genomic region of 193 kb. By sequencing the whole genome of the resistant line harboring the major resistance gene HaOr5, we identified a large insertion of a complex genomic region containing a cluster of putative resistance genes.

Development of a rapid detection method for the macrolide resistance gene in Mycobacterium avium using the amplification refractory mutation system-loop-mediated isothermal amplification method

Macrolide antibiotics such as clarithromycin (CLR) and azithromycin are the key drugs used in multidrug therapy for Mycobacterium avium complex (MAC) diseases. For these antibacterial drugs, drug susceptibility has been correlated with clinical response in MAC diseases. We have previously demonstrated the correlation between drug susceptibility and mutations in the 23S rRNA gene, which confers resistance to macrolides. Herein, we developed a rapid detection method using the amplification refractory mutation system (ARMS)-loop-mediated isothermal amplification (LAMP) technique to identify mutations in the 23S rRNA gene of M. avium. We examined the applicability of the ARMS-LAMP method to genomic DNA extracted from six genotypes of M. avium clinical isolates. The M. avium isolates were classified into 21 CLR-resistant and 9 CLR-susceptible strains based on the results of drug susceptibility tests; the 23S rRNA genes of these strains were sequenced and analyzed using the ARMS-LAMP method. Sequence analysis revealed that the 9 CLR-sensitive strains were wild-type strains, whereas the 21 CLR-resistant strains comprised 20 mutant-type strains and one wild-type strain. Using ARMS-LAMP, no amplification from genomic DNAs of the 10 wild-type strains was observed using the mutant-type mismatch primer sets (MTPSs); however, amplification from the 20 mutant-type strain DNAs was observed using the MTPSs. The rapid detection method developed by us integrates ARMS-LAMP with a real-time turbidimeter, which can help determine drug resistance in a few hours. In conclusion, ARMS-LAMP might be a new clinically beneficial technology for rapid detection of mutations.IMPORTANCEMultidrug therapy for pulmonary Mycobacterium avium complex disease is centered on the macrolide antibiotics clarithromycin and azithromycin, and resistance to macrolides is an important prognosticator for clinical aggravation. Therefore, it is important to develop a quick and easy method for detecting resistance to macrolides. Drug resistance is known to be correlated with mutations in macrolide resistance genes. We developed a rapid detection method using amplification refractory mutation system (ARMS)-loop-mediated isothermal amplification (LAMP) to identify a mutation in the 23S rRNA gene, which is a macrolide resistance gene. Furthermore, we examined the applicability of this method using M. avium clinical isolates. The rapid method developed by us for detection of the macrolide resistance gene by integrating ARMS-LAMP and a real-time turbidimeter can help in detection of drug resistance within a few hours. Since this method does not require expensive equipment or special techniques and shows high analytical speed, it would be very useful in clinical practice.

Low cycle number multiplex PCR: A novel strategy for the construction of amplicon libraries for next-generation sequencing

Multiplex PCR is a critical step when preparing amplicon library for next-generation sequencing. However, there are several challenges related to multiplex PCR including poor uniformity, nonspecific amplification, and primer-dimers. To address these issues, we propose a novel solution strategy that involves using a low cycle number (<10 cycles) in multiplex PCR and then employing carrier DNAs and magnetic beads for the selection of targeted products. This technique improves the amplicon uniformity while also reducing primer-dimers and PCR artifacts. To evaluate our technique, we initially utilized 120 DNA fragments from mouse genome containing single nucleotide polymorphism (SNP) sites. Sequencing results demonstrated that with only 7 cycles of multiplex PCR, 95.8% of the targeted SNP sites were mapped, with a coverage of at least 1×. The average sequencing depth of all amplicons was 1705.79 ± 1205.30×; 87% of them reached a coverage depth that exceeded 0.2-fold of the average sequencing depth. Our method had a greater uniformity (87%) when compared to Hi-Plex PCR (53.3%). Furthermore, we validated our strategy by randomly selecting 90 primer pairs twice from the initial set of 120 primer-pairs. Next, we used the same protocol to prepare amplicon libraries. The two groups had an average sequencing depth of 1013.30 ± 585.57× and 219.10 ± 158.27×, respectively; over 84% of the amplicons had a sequencing depth that exceeded 0.2-fold of average depth. These results suggest that the use of a low cycle number in multiplex PCR is a cost-effective and efficient approach for the preparation of amplicon libraries.

The Hessian fly resistance gene HvRHF1 is localized in an NBS-LRR gene cluster in barley

Hessian fly (Mayetiola destructor Say) is a significant pest in cereal crops, causing substantial yield losses worldwide. While host resistance is the most efficient method for pest control, research on genetic characterization of Hessian fly resistance in barley (Hordeum vulgare L.) has been limited, and the underlying resistance mechanism remains largely unknown. In this study, we conducted fine mapping of a crucial Hessian fly resistance locus, known as HvRHF1, using a biparental population. Assisted with genetic markers and robust phenotyping assay, we pinpointed the HvRHF1 gene to an ~ 82 kb region on chromosome 4H. Gene prediction and annotation revealed that the HvRHF1 locus comprises three complete NBS-LRR genes, which are characteristic of disease resistance genes. As a result, our study not only provides valuable resources for resistance in barley and genetic tools for breeding, but also identifies candidate genes that lay the foundation for cloning HvRHF1. This endeavor will significantly contribute to our understanding of the molecular mechanisms underlying cereal resistance to Hessian fly.

FBPP: software to design PCR primers and probes for nucleic acid base detection of foodborne pathogens

Foodborne pathogens can be found in various foods, and it is important to detect foodborne pathogens to provide a safe food supply and to prevent foodborne diseases. The nucleic acid base detection method is one of the most rapid and widely used methods in the detection of foodborne pathogens; it depends on hybridizing the target nucleic acid sequence to a synthetic oligonucleotide (probes or primers) that is complementary to the target sequence. Designing primers and probes for this method is a preliminary and critical step. However, new bioinformatics tools are needed to automate, specific and improve the design sets to be used in the nucleic acid‒base method. Thus, we developed foodborne pathogen primer probe design (FBPP), an open-source, user-friendly graphical interface Python-based application supported by the SQL database for foodborne pathogen virulence factors, for (i) designing primers/probes for detection purposes, (ii) PCR and gel electrophoresis photo simulation, and (iii) checking the specificity of primers/probes.

New Set of EST-STR Markers for Discrimination of Related Papaver somniferum L. Varieties

Papaver somniferum L. is cultivated for its edible seeds and for the production of alkaloids. A serious problem in seed trade and processing is the intentional mixing of excellent food-quality seeds with non-food-grade-quality seeds. Tracking the correct or illegitimate handling of seeds requires an efficient method for discrimination and individualization of poppy varieties. As in human and animal forensics, DNA variable regions containing short tandem repeats (STRs) located either in non-coding DNA or in gene sequences (EST-STRs) are preferred markers for discrimination between genotypes. Primers designed for 10 poppy EST-STR loci not analyzed so far were tested for their discriminatory ability on a set of 23 related P. somniferum L. genotypes. Thirty-three EST-STR alleles were identified together. Their polymorphic information content (PIC) values were in the range of 0.175-0.649. The PI value varied in the range of 0.140-0.669, and the cumulative PI was 1.2 × 10-5. PIsibs values varied between 0.436 and 0.820 and the cumulative value was lower (5.0 × 10-3). All analyzed genotypes were distinguished mutually, each with its own unique EST-STR profile. These newly developed EST-STR markers more effectively discriminated P. somniferum L. genotypes, even those genotypes whose DNA profiles were previously identical.

Genome-resolved metatranscriptomics reveals conserved root colonization determinants in a synthetic microbiota

The identification of processes activated by specific microbes during microbiota colonization of plant roots has been hampered by technical constraints in metatranscriptomics. These include lack of reference genomes, high representation of host or microbial rRNA sequences in datasets, or difficulty to experimentally validate gene functions. Here, we recolonized germ-free Arabidopsis thaliana with a synthetic, yet representative root microbiota comprising 106 genome-sequenced bacterial and fungal isolates. We used multi-kingdom rRNA depletion, deep RNA-sequencing and read mapping against reference microbial genomes to analyse the in planta metatranscriptome of abundant colonizers. We identified over 3,000 microbial genes that were differentially regulated at the soil-root interface. Translation and energy production processes were consistently activated in planta, and their induction correlated with bacterial strains' abundance in roots. Finally, we used targeted mutagenesis to show that several genes consistently induced by multiple bacteria are required for root colonization in one of the abundant bacterial strains (a genetically tractable Rhodanobacter). Our results indicate that microbiota members activate strain-specific processes but also common gene sets to colonize plant roots.

Transcriptome analysis of Haemaphysalis flava female using Illumina HiSeq 4000 sequencing: de novo assembly, functional annotation and discovery of SSR markers

BACKGROUND

Ticks are ectoparasites capable of directly damaging their hosts and transmitting vector-borne diseases. The ixodid tick Haemaphysalis flava has a broad distribution that extends from East to South Asia. This tick is a reservoir of severe fever with thrombocytopenia syndrome virus (SFTSV) that causes severe hemorrhagic disease, with cases reported from China, Japan and South Korea. Recently, the distribution of H. flava in South Korea was found to overlap with the occurrence of SFTSV.

METHODS

This study was undertaken to discover the molecular resources of H. flava female ticks using the Illumina HiSeq 4000 system, the Trinity de novo sequence assembler and annotation against public databases. The locally curated Protostome database (PANM-DB) was used to screen the putative adaptation-related transcripts classified to gene families, such as angiotensin-converting enzyme, aquaporin, adenylate cyclase, AMP-activated protein kinase, glutamate receptors, heat shock proteins, molecular chaperones, insulin receptor, mitogen-activated protein kinase and solute carrier family proteins. Also, the repeats and simple sequence repeats (SSRs) were screened from the unigenes using RepeatMasker (v4.0.6) and MISA (v1.0) software tools, followed by the designing of SSRs flanking primers using BatchPrimer 3 (v1.0) software.

RESULTS

The transcriptome produced a total of 69,822 unigenes, of which 46,175 annotated to the homologous proteins in the PANM-DB. The unigenes were also mapped to the EuKaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) specializations. Promiscuous presence of protein kinase, zinc finger (C2H2-type), reverse transcriptase, and RNA recognition motif domains was observed in the unigenes. A total of 3480 SSRs were screened, of which 1907 and 1274 were found as tri- and dinucleotide repeats, respectively. A list of primer sequences flanking the SSR motifs was detailed for validation of polymorphism in H. flava and the related tick species.

CONCLUSIONS

The reference transcriptome information on H. flava female ticks will be useful for an enriched understanding of tick biology, its competency to act as a vector and the study of species diversity related to disease transmission.

Tetraose steroidal glycoalkaloids from potato provide resistance against Alternaria solani and Colorado potato beetle

Plants with innate disease and pest resistance can contribute to more sustainable agriculture. Natural defence compounds produced by plants have the potential to provide a general protective effect against pathogens and pests, but they are not a primary target in resistance breeding. Here, we identified a wild relative of potato, Solanum commersonii, that provides us with unique insight in the role of glycoalkaloids in plant immunity. We cloned two atypical resistance genes that provide resistance to Alternaria solani and Colorado potato beetle through the production of tetraose steroidal glycoalkaloids (SGA). Moreover, we provide in vitro evidence to show that these compounds have potential against a range of different (potato pathogenic) fungi. This research links structural variation in SGAs to resistance against potato diseases and pests. Further research on the biosynthesis of plant defence compounds in different tissues, their toxicity, and the mechanisms for detoxification, can aid the effective use of such compounds to improve sustainability of our food production.

Fine mapping of two recessive powdery mildew resistance genes from Aegilops tauschii accession CIae8

KEY MESSAGE

Two recessive powdery mildew resistance loci pmAeCIae8_2DS and pmAeCIae8_7DS from Aegilops tauschii were mapped and two synthesized hexaploid wheat lines were developed by distant hybridization. Wheat powdery mildew (Pm), one of the worldwide destructive fungal diseases, causes significant yield loss up to 30%. The identification of new Pm resistance genes will enrich the genetic diversity of wheat breeding for Pm resistance. Aegilops tauschii is the ancestor donor of sub-genome D of hexaploid wheat. It provides beneficial genes that can be easily transferred into wheat by producing synthetic hexaploid wheat followed by genetic recombination. We assessed the Pm resistance level of 35 Ae. tauschii accessions from different origins. Accession CIae8 exhibited high Pm resistance. Inheritance analysis and gene mapping were performed using F2 and F2:3 populations derived from the cross between CIae8 and a Pm susceptible accession PI574467. The Pm resistance of CIae8 was controlled by two independent recessive genes. Bulked segregate analysis using a 55 K SNP array revealed the SNPs were mainly enriched into genome regions, i.e. 2DS (13.5-20 Mb) and 7DS (4.0-15.5 Mb). The Pm resistance loci were named as pmAeCIae8_2DS and pmAeCIae8_7DS, respectively. By recombinant screening, we narrowed the pmAeCIae8_2DS into a 370-kb interval flanked by markers CINAU-AE7800 (14.89 Mb) and CINAU-AE20 (15.26 Mb), and narrowed the pmAeCIae8_7DS into a 260-kb interval flanked by markers CINAU-AE58 (4.72 Mb) and CINAU-AE25 (4.98 Mb). The molecular markers closely linked with the resistance loci were developed, and two synthesized hexaploid wheat (SHW) lines were produced. These laid the foundation for cloning of the two resistance loci and for transferring the resistance into common wheat.

Assessing the recovery of Y chromosome microsatellites with population genomic data using Papio and Theropithecus genomes

Y chromosome markers can shed light on male-specific population dynamics but for many species no such markers have been discovered and are available yet, despite the potential for recovering Y-linked loci from available genome sequences. Here, we investigated how effective available bioinformatic tools are in recovering informative Y chromosome microsatellites from whole genome sequence data. In order to do so, we initially explored a large dataset of whole genome sequences comprising individuals at various coverages belonging to different species of baboons (genus: Papio) using Y chromosome references belonging to the same genus and more distantly related species (Macaca mulatta). We then further tested this approach by recovering Y-STRs from available Theropithecus gelada genomes using Papio and Macaca Y chromosome as reference sequences. Identified loci were validated in silico by a) comparing within-species relationships of Y chromosome lineages and b) genotyping male individuals in available pedigrees. Each STR was selected not to extend in its variable region beyond 100 base pairs, so that loci can be developed for PCR-based genotyping of non-invasive DNA samples. In addition to assembling a first set of Papio and Theropithecus Y-specific microsatellite markers, we released TYpeSTeR, an easy-to-use script to identify and genotype Y chromosome STRs using population genomic data which can be modulated according to available male reference genomes and genomic data, making it widely applicable across taxa.

EasySSR: a user-friendly web application with full command-line features for large-scale batch microsatellite mining and samples comparison

Microsatellites, also known as SSRs or STRs, are polymorphic DNA regions with tandem repetitions of a nucleotide motif of size 1-6 base pairs with a broad range of applications in many fields, such as comparative genomics, molecular biology, and forensics. However, the majority of researchers do not have computational training and struggle while running command-line tools or very limited web tools for their SSR research, spending a considerable amount of time learning how to execute the software and conducting the post-processing data tabulation in other tools or manually-time that could be used directly in data analysis. We present EasySSR, a user-friendly web tool with command-line full functionality, designed for practical use in batch identifying and comparing SSRs in sequences, draft, or complete genomes, not requiring previous bioinformatic skills to run. EasySSR requires only a FASTA and an optional GENBANK file of one or more genomes to identify and compare STRs. The tool can automatically analyze and compare SSRs in whole genomes, convert GenBank to PTT files, identify perfect and imperfect SSRs and coding and non-coding regions, compare their frequencies, abundancy, motifs, flanking sequences, and iterations, producing many outputs ready for download such as PTT files, interactive charts, and Excel tables, giving the user the data ready for further analysis in minutes. EasySSR was implemented as a web application, which can be executed from any browser and is available for free at https://computationalbiology.ufpa.br/easyssr/. Tutorials, usage notes, and download links to the source code can be found at https://github.com/engbiopct/EasySSR.

Genome-wide identification, characterization and primer designing of simple sequence repeats across Leguminosae family

Legumes are important clade of commercially important family Leguminosae that mainly include medicinal, flowering and edible plants. Although the genomic sequence of legumes is accessible, only the limited number of effective simple sequence repeat markers has been identified by prior research. Additional polymorphic simple sequence repeats marker discovery will aid in the genetics and breeding of legumes. In this study, 13 complete genome sequences were screened for the identification of chromosome-wise simple sequence repeats (SSRs) and 1,866,861 SSRs were identified. Based on the study, it was observed that the number of SSRs in non-coding region was more as compared to coding region and frequency of mononucleotides was highest followed by di-nucleotides while penta- and hexa-nucleotide repeats were least frequent one. The identified genome-wide SSRs and newly developed SSR markers, primers and their mapping will provide a powerful means for genetic researches across Leguminosae plants, including genetic diversity and evolutionary origin analysis, fingerprinting, QTL mapping and marker-assisted selection for breeding as well as comparative genomic analysis studies.

Transcriptome analysis provides insights into the stress response in cultivated peanut (Arachis hypogaea L.) subjected to drought-stress

BACKGROUND

Peanut (Arachis hypogaea L.) is one of the valuable oilseed crops grown in drought-prone areas worldwide. Drought severely limits peanut production and productivity significantly.

METHOD AND RESULTS

In order to decipher the drought tolerance mechanism in peanut under drought stress, RNA sequencing was performed in TAG - 24 (drought tolerant genotype) and JL-24 (drought susceptible genotype). Approximately 51 million raw reads were generated from four different libraries of two genotypes subjected to drought stress exerted by 20% PEG 6000 stress and control conditions, of which ~ 41 million (80.87%) filtered reads were mapped to the Arachis hypogaea L. reference genome. The transcriptome analysis detected 1,629 differentially expressed genes (DEGs), 186 genes encoding transcription factors (TFs) and 30,199 SSR among the identified DEGs. Among the differentially expressed TF encoding genes, the highest number of genes were WRKY followed by bZIP, C2H2, and MYB during drought stress. The comparative analysis between the two genotypes revealed that TAG-24 exhibits activation of certain key genes and transcriptional factors that are involved in essential biological processes. Specifically, TAG-24 showed activation of genes involved in the plant hormone signaling pathway such as PYL9, Auxin response receptor gene, and ABA. Additionally, genes related to water deprivation such as LEA protein and those involved in combating oxidative damage such as Glutathione reductase were also found to be activated in TAG-24.

CONCLUSION

This genome-wide transcription map, therefore, provides a valuable tool for future transcript profiling under drought stress and enriches the genetic resources available for this important oilseed crop.

Next-Generation-Sequencing-Based Simple Sequence Repeat (SSR) Marker Development and Linkage Mapping in Lentil (Lens culinaris L.)

Simple sequence repeats (SSRs) are highly versatile markers in genetic diversity analysis and plant breeding, making them widely applicable. They hold potential in lentil (Lens culinaris) breeding for genetic diversity analysis, marker-assisted selection (MAS), and linkage mapping. However, the availability and diversity of SSR markers in lentil is limited. We used next-generation sequencing (NGS) technology to develop SSR markers in lentil. NGS allowed us to identify regions of the lentil genome that contained SSRs. Illumina Hiseq-2000 sequencing of the lentil genotype "Karacadağ" resulted in 1,727,734 sequence reads comprising more than 48,390 Mb, and contigs were mined for SSRs, resulting in the identification of a total of 8697 SSR motifs. Among these, dinucleotide repeats were the most abundant (53.38%), followed by trinucleotides (30.38%), hexanucleotides (6.96%), tetranucleotides (6.59%), and pentanucleotides (3.19%). The most frequent repeat in dinucleotides was the TC (21.80%), followed by the GA (17.60%). A total of 2000 primer pairs were designed from these motifs, and 458 SSR markers were validated following their amplified PCR products. A linkage map was constructed using these new SSRs with high linkage disequilibrium (209) and previously known SSRs (11). The highest number of SSR markers (43) was obtained in LG2, while the lowest number of SSR markers (19) was obtained in LG7. The longest linkage group (LG) was LG2 (86.84 cM), whereas the shortest linkage group was LG7 (53.46 cM). The average length between markers ranged from 1.86 cM in LG1 to 2.81 cM in LG7, and the map density was 2.16 cM. The developed SSRs and created linkage map may provide useful information and offer a new library for genetic diversity analyses, linkage mapping studies, and lentil breeding programs.

De Novo Mining and Validating Novel Microsatellite Markers to Assess Genetic Diversity in Maruca vitrata (F.), a Legume Pod Borer

Maruca vitrata (Fabricius) is an invasive insect pest capable of causing enormous economic losses to a broad spectrum of leguminous crops. Microsatellites are valuable molecular markers for population genetic studies; however, an inadequate number of M. vitrata microsatellite loci are available to carry out population association studies. Thus, we utilized this insect's public domain databases for mining expressed sequence tags (EST)-derived microsatellite markers. In total, 234 microsatellite markers were identified from 10053 unigenes. We discovered that trinucleotide repeats were the most predominant microsatellite motifs (61.53%), followed by dinucleotide repeats (23.50%) and tetranucleotide repeats (14.95%). Based on the analysis, twenty-five markers were selected for validation in M. vitrata populations collected from various regions of India. The number of alleles (Na), observed heterozygosity (Ho), and expected heterozygosity (He) ranged from 2 to 5; 0.00 to 0.80; and 0.10 to 0.69, respectively. The polymorphic loci showed polymorphism information content (PIC), ranging from 0.09 to 0.72. Based on the genetic distance matrix, the unrooted neighbor-joining dendrogram differentiated the selected populations into two discrete groups. The SSR markers developed and validated in this study will be helpful in population-level investigations of M. vitrata to understand the gene flow, demography, dispersal patterns, biotype differentiation, and host dynamics.

Transcriptomic classes of BCR-ABL1 lymphoblastic leukemia

In BCR-ABL1 lymphoblastic leukemia, treatment heterogeneity to tyrosine kinase inhibitors (TKIs), especially in the absence of kinase domain mutations in BCR-ABL1, is poorly understood. Through deep molecular profiling, we uncovered three transcriptomic subtypes of BCR-ABL1 lymphoblastic leukemia, each representing a maturation arrest at a stage of B-cell progenitor differentiation. An earlier arrest was associated with lineage promiscuity, treatment refractoriness and poor patient outcomes. A later arrest was associated with lineage fidelity, durable leukemia remissions and improved patient outcomes. Each maturation arrest was marked by specific genomic events that control different transition points in B-cell development. Interestingly, these events were absent in BCR-ABL1⁺ preleukemic stem cells isolated from patients regardless of subtype, which supports that transcriptomic phenotypes are determined downstream of the leukemia-initialing event. Overall, our data indicate that treatment response and TKI efficacy are unexpected outcomes of the differentiation stage at which this leukemia transforms.

Genome-wide identification and expression profile analysis of metal tolerance protein gene family in Eucalyptus grandis under metal stresses

Metal tolerance proteins (MTPs) as Me²⁺/H⁺(K⁺) antiporters participate in the transport of divalent cations, leading to heavy metal stress resistance and mineral utilization in plants. In the present study, to obtain better knowledge of the biological functions of the MTPs family, 20 potential EgMTPs genes were identified in Eucalyptus grandis and classified into seven groups belonging to three cation diffusion facilitator groups (Mn-CDFs, Zn/Fe-CDFs, and Zn-CDFs) and seven groups. EgMTP-encoded amino acids ranged from 315 to 884, and most of them contained 4-6 recognized transmembrane domains and were clearly prognosticated to localize into the cell vacuole. Almost all EgMTP genes experienced gene duplication events, in which some might be uniformly distributed in the genome. The numbers of cation efflux and the zinc transporter dimerization domain were highest in EgMTP proteins. The promoter regions of EgMTP genes have different cis-regulatory elements, indicating that the transcription rate of EgMTP genes can be a controlled response to different stimuli in multiple pathways. Our findings provide accurate perception on the role of the predicted miRNAs and the presence of SSR marker in the Eucalyptus genome and clarify their functions in metal tolerance regulation and marker-assisted selection, respectively. Gene expression profiling based on previous RNA-seq data indicates a probable function for EgMTP genes during development and responses to biotic stress. Additionally, the upregulation of EgMTP6, EgMTP5, and EgMTP11.1 to excess Cd²⁺ and Cu²⁺ exposure might be responsible for metal translocation from roots to leaves.

Identification of qBK2.1, a novel QTL controlling rice resistance against Fusarium fujikuroi

BACKGROUND

Bakanae disease caused by Fusarium fujikuroi is an increasing threat to rice production. The infected plants show symptoms such as elongation, slenderness, chlorosis, a large leaf angle, and even death. Bakanae disease is traditionally managed by seed treatment. However, fungicide-resistant F. fujikuroi isolates have emerged in several Asian areas, including Taiwan. This study aimed to identify new bakanae resistance quantitative trait loci (QTLs) and provide molecular markers to assist future breeding.

RESULTS

A population of F_2:9 recombinant inbred lines (RILs) was derived from the cross between an elite japonica Taiwanese cultivar 'Taikeng 16 (TK16)' and an indica variety 'Budda'. 'Budda' was found highly resistant to all 24 representative isolates of the F. fujikuroi population in Taiwan. For the RIL population, 6,492 polymorphic single nucleotide polymorphisms (SNPs) spanning the rice genome were obtained by genotyping-by-sequencing (GBS) technique, and the disease severity index (DSI) was evaluated by inoculation with a highly virulent F. fujikuroi isolate Ff266. Trait-marker association analysis of 166 RILs identified two QTLs in 'Budda'. qBK2.1 (21.97-30.15 Mb) is a novel and first bakanae resistance QTL identified on chromosome 2. qBK1.8 (5.24-8.66 Mb) partially overlaps with the previously reported qBK1.3 (4.65-8.41 Mb) on chromosome 1. The log of odds (LOD) scores of qBK1.8 and qBK2.1 were 4.75 and 6.13, accounting for 4.9% and 8.1% of the total phenotypic variation, respectively. 64 RILs carrying both qBK1.8 and qBK2.1 showed lower DSI (7%) than the lines carrying only qBK1.8 (15%), only qBK2.1 (13%), or none of the two QTLs (21%). For the future application of identified QTLs, 11 KBioscience competitive allele-specific PCR (KASP) markers and 3 insertion-deletion (InDel) markers were developed.

CONCLUSIONS

Compared to other important rice diseases, knowledge of bakanae resistance has been insufficient, which limited the development and deployment of resistant cultivars. The discovery of qBK2.1 has provided a new source of bakanae resistance. The resistant RILs inheriting good plant type, good taste, and high yield characteristics from 'TK16' can be used as good resistance donors. Our newly developed markers targeting qBK2.1 and qBK1.8 can also serve as an important basis for future fine-mapping and resistance breeding.

Population Genetic Structure of a Rare Butterfly in a Fragmented South Florida Ecosystem

We investigated the genetic structure and diversity between populations of a rare butterfly, the Florida duskywing (Ephyriades brunnea floridensis E. Bell and W. Comstock, 1948) (Lepidoptera: Hesperiidae) across a network of South Florida pine rockland habitat fragments. Based on 81 individuals from seven populations and using multiple polymorphic microsatellite loci, our analyses support the presence of mainland Florida (peninsular) and Florida Keys (island) population groupings, with a moderate, asymmetrical gene flow connecting them, and the presence of private alleles providing unique identities to each. We additionally found that despite a prevalence in many Lepidoptera, the presence of Wolbachia was not identified in any of the samples screened. Our findings can be used to inform conservation and recovery decisions, including population monitoring, organism translocation, and priority areas for management, restoration or stepping-stone creation to help maintain the complex genetic structure of separate populations.

Periodontal disease is associated with increased gut colonization of pathogenic Haemophilus parainfluenzae in patients with Crohn's disease

Intestinal colonization of the oral bacterium Haemophilus parainfluenzae has been associated with Crohn's disease (CD) severity and progression. This study examines the role of periodontal disease (PD) as a modifier for colonization of H. parainfluenzae in patients with CD and explores the mechanisms behind H. parainfluenzae-mediated intestinal inflammation. Fifty subjects with and without CD were evaluated for the presence of PD, and their oral and fecal microbiomes were characterized. PD is associated with increased levels of H. parainfluenzae strains in subjects with CD. Oral inoculation of H. parainfluenzae elicits strain-dependent intestinal inflammation in murine models of inflammatory bowel disease, which is associated with increased intestinal interferon-γ (IFN-γ)+ CD4+ T cells and disruption of the host hypusination pathway. In summary, this study establishes a strain-specific pathogenic role of H. parainfluenzae in intestinal inflammation and highlights the potential effect of PD on intestinal colonization by pathogenic H. parainfluenzae strains in patients with CD.

Tiller Number1 encodes an ankyrin repeat protein that controls tillering in bread wheat

Wheat (Triticum aestivum L.) is a major staple food for more than one-third of the world's population. Tiller number is an important agronomic trait in wheat, but only few related genes have been cloned. Here, we isolate a wheat mutant, tiller number1 (tn1), with much fewer tillers. We clone the TN1 gene via map-based cloning: TN1 encodes an ankyrin repeat protein with a transmembrane domain (ANK-TM). We show that a single amino acid substitution in the third conserved ankyrin repeat domain causes the decreased tiller number of tn1 mutant plants. Resequencing and haplotype analysis indicate that TN1 is conserved in wheat landraces and modern cultivars. Further, we reveal that the expression level of the abscisic acid (ABA) biosynthetic gene TaNCED3 and ABA content are significantly increased in the shoot base and tiller bud of the tn1 mutants; TN1 but not tn1 could inhibit the binding of TaPYL to TaPP2C via direct interaction with TaPYL. Taken together, we clone a key wheat tiller number regulatory gene TN1, which promotes tiller bud outgrowth probably through inhibiting ABA biosynthesis and signaling.

Can allele-specific loop-mediated isothermal amplification be used for rapid detection of target-site herbicide resistance in Lolium spp.?

BACKGROUND

Herbicide resistance is one of the threats to modern agriculture and its early detection is one of the most effective components for sustainable resistance management strategies. Many techniques have been used for target-site-resistance detection. Allele-Specific Loop-Mediated Isothermal Amplification (AS-LAMP) was evaluated as a possible rapid diagnostic method for acetyl-CoA carboxylase (ACCase) and acetolactate synthase (ALS) inhibiting herbicides resistance in Lolium spp.

RESULTS

AS-LAMP protocols were set up for the most frequent mutations responsible for herbicide resistance to ALS (positions 197, 376 and 574) and ACCase (positions 1781, 2041 and 2078) inhibitors in previously characterized and genotyped Lolium spp.

POPULATIONS

A validation step on new putative resistant populations gave the overview of a possible use of this tool for herbicide resistance diagnosis in Lolium spp. Regarding the ACCase inhibitor pinoxaden, in more than 65% of the analysed plants, the LAMP assay and genotyping were in keeping, whereas the results were not consistent when ALS inhibitors resistance was considered. Limitations on the use of this technique for herbicide resistance detection in the allogamous Lolium spp. are discussed.

CONCLUSIONS

The LAMP method used for the detection of target-site resistance in weed species could be applicable with target genes that do not have high genetic variability, such as ACCase gene in Lolium spp.

Initiation of genomics-assisted breeding in Virginia-type peanuts through the generation of a de novo reference genome and informative markers

INTRODUCTION

Virginia-type peanut, Arachis hypogaea subsp. hypogaea, is the second largest market class of peanut cultivated in the United States. It is mainly used for large-seeded, in-shell products. Historically, Virginia-type peanut cultivars were developed through long-term recurrent phenotypic selection and wild species introgression projects. Contemporary genomic technologies represent a unique opportunity to revolutionize the traditional breeding pipeline. While there are genomic tools available for wild and cultivated peanuts, none are tailored specifically to applied Virginia-type cultivar development programs.

METHODS AND RESPECTIVE RESULTS

Here, the first Virginia-type peanut reference genome, "Bailey II", was assembled. It has improved contiguity and reduced instances of manual curation in chromosome arms. Whole-genome sequencing and marker discovery was conducted on 66 peanut lines which resulted in 1.15 million markers. The high marker resolution achieved allowed 34 unique wild species introgression blocks to be cataloged in the A. hypogaea genome, some of which are known to confer resistance to one or more pathogens. To enable marker-assisted selection of the blocks, 111 PCR Allele Competitive Extension assays were designed. Forty thousand high quality markers were selected from the full set that are suitable for mid-density genotyping for genomic selection. Genomic data from representative advanced Virginia-type peanut lines suggests this is an appropriate base population for genomic selection.

DISCUSSION

The findings and tools produced in this research will allow for rapid genetic gain in the Virginia-type peanut population. Genomics-assisted breeding will allow swift response to changing biotic and abiotic threats, and ultimately the development of superior cultivars for public use and consumption.

Whole genome resequencing identifies candidate genes and allelic diagnostic markers for resistance to Ralstonia solanacearum infection in cultivated peanut (Arachis hypogaea L.)

Bacterial wilt disease (BWD), caused by Ralstonia solanacearum is a major challenge for peanut production in China and significantly affects global peanut field productivity. It is imperative to identify genetic loci and putative genes controlling resistance to R. solanacearum (RRS). Therefore, a sequencing-based trait mapping approach termed "QTL-seq" was applied to a recombination inbred line population of 581 individuals from the cross of Yueyou 92 (resistant) and Xinhuixiaoli (susceptible). A total of 381,642 homozygous single nucleotide polymorphisms (SNPs) and 98,918 InDels were identified through whole genome resequencing of resistant and susceptible parents for RRS. Using QTL-seq analysis, a candidate genomic region comprising of 7.2 Mb (1.8-9.0 Mb) was identified on chromosome 12 which was found to be significantly associated with RRS based on combined Euclidean Distance (ED) and SNP-index methods. This candidate genomic region had 180 nonsynonymous SNPs and 14 InDels that affected 75 and 11 putative candidate genes, respectively. Finally, eight nucleotide binding site leucine rich repeat (NBS-LRR) putative resistant genes were identified as the important candidate genes with high confidence. Two diagnostic SNP markers were validated and revealed high phenotypic variation in the different resistant and susceptible RIL lines. These findings advocate the expediency of the QTL-seq approach for precise and rapid identification of candidate genomic regions, and the development of diagnostic markers that are applicable in breeding disease-resistant peanut varieties.

Stepwise Optimization of Real-Time RT-PCR Analysis

Quantitative real-time reverse transcription PCR (qRT-PCR) analysis has been used routinely to quantify gene expression levels. Primer design and the optimization of qRT-PCR parameters are critical for the accuracy and reproducibility of qRT-PCR analysis. Computational tool-assisted primer design often overlooks the presence of homologous sequences of the gene of interest and the sequence similarities between homologous genes in a plant genome. This sometimes results in skipping the optimization of qRT-PCR parameters due to the false confidence in the quality of the designed primers. Here we present a stepwise optimization protocol for single nucleotide polymorphisms (SNPs)-based sequence-specific primer design and sequential optimization of primer sequences, annealing temperatures, primer concentrations, and cDNA concentration range for each reference and target gene. The goal of this optimization protocol is to achieve a standard cDNA concentration curve with an R² ≥ 0.9999 and efficiency (E) = 100 ± 5% for the best primer pair of each gene, which serves as the prerequisite for using the 2^-ΔΔCT method for data analysis.

Resistance of QYm.nau-2D to wheat yellow mosaic virus was derived from an alien introgression into common wheat

The QYm.nau-2D locus conferring wheat yellow mosaic virus resistance is an exotic introgression and we developed 11 diagnostic markers tightly linked to QYm.nau-2D. Wheat yellow mosaic virus (WYMV) is a serious disease of winter wheat in China. Breeding resistant varieties is the most effective strategy for WYMV control. A WYMV resistant locus QYm.nau-2D on the chromosome arm 2DL has been repeatedly reported but the mapped region is large. In the present study, we screened recombinants using a biparental population and mapped QYm.nau-2D into an 18.8 Mb physical interval. By genome-wide association studies of 372 wheat varieties for WYMV resistance in four environments, we narrowed down QYm.nau-2D into a 16.4 Mb interval. Haplotype analysis indicated QYm.nau-2D were present as six different states due to recombination during hybridization breeding. QYm.nau-2D was finally mapped into a linkage block of 11.2 Mb. Chromosome painting using 2D specific probes and collinearity analysis among the published sequences corresponding to QYm.nau-2D region indicated the block was an exotic introgression. The Illumina-sequenced reads of four diploid Aegilops species were mapped to the sequence of Fielder, a variety having the introgression. The mapping reads were significantly increased at the putative introgression regions of Fielder. Ae. uniaristata (NN) had the highest mapping reads, suggesting that QYm.nau-2D was possibly an introgression from genome N. We investigated the agronomic performances of different haplotypes and observed no linkage drag of the alien introgression for the 15 tested traits. For marker-assisted selection of QYm.nau-2D, we developed 11 diagnostic markers tightly linked to the locus. This research provided a case study of an exotic introgression, which has been utilized in wheat improvement for WYMV resistance.

MicroRNA Identification, Target Prediction, and Validation for Crop Improvement

Micro-RNAs (mi-RNAs) are regulatory elements that play a vital role in the growth, development, and metabolic regulation of plants. In current research, the isolation of miRNAs is a tedious and difficult task using in vitro methods. However, recent exploration into the remarkably highly conserved nature of nucleotide sequences of miRNAs assists in the identification of miRNAs in plant species through homologous approaches. Here, we describe the in silico-based method for identification of miRNAs from the EST database which is emerging as a faster and more reliable approach along with the development of miRNA-SSR markers. This approach has the potential to accelerate research into the regulation of gene expression in various plant species such as tea, potato, tomato, tobacco, and orphan crops like cluster bean.

Alkaloid production and response to natural adverse conditions in Peganum harmala: in silico transcriptome analyses

Peganum harmala is a valuable wild plant that grows and survives under adverse conditions and produces pharmaceutical alkaloid metabolites. Using different assemblers to develop a transcriptome improves the quality of assembled transcriptome. In this study, a concrete and accurate method for detecting stress-responsive transcripts by comparing stress-related gene ontology (GO) terms and public domains was designed. An integrated transcriptome for P. harmala including 42 656 coding sequences was created by merging de novo assembled transcriptomes. Around 35 000 transcripts were annotated with more than 90% resemblance to three closely related species of Citrus, which confirmed the robustness of the assembled transcriptome; 4853 stress-responsive transcripts were identified. CYP82 involved in alkaloid biosynthesis showed a higher number of transcripts in P. harmala than in other plants, indicating its diverse alkaloid biosynthesis attributes. Transcription factors (TFs) and regulatory elements with 3887 transcripts comprised 9% of the transcriptome. Among the TFs of the integrated transcriptome, cystein2/histidine2 (C2H2) and WD40 repeat families were the most abundant. The Kyoto Encyclopedia of Genes and Genomes (KEGG) MAPK (mitogen-activated protein kinase) signaling map and the plant hormone signal transduction map showed the highest assigned genes to these pathways, suggesting their potential stress resistance. The P. harmala whole-transcriptome survey provides important resources and paves the way for functional and comparative genomic studies on this plant to discover stress-tolerance-related markers and response mechanisms in stress physiology, phytochemistry, ecology, biodiversity, and evolution. P. harmala can be a potential model for studying adverse environmental cues and metabolite biosynthesis and a major source for the production of various alkaloids.

High-throughput primer design by scoring in piecewise logistic model for multiple polymerase chain reaction variants

Polymerase chain reaction (PCR) variants requiring specific primer types are widely used in various PCR experiments, including generic PCR, inverse PCR, anchored PCR, and ARMS PCR. Few tools can be adapted for multiple PCR variants, and many tools select primers by filtration based on the given parameters, which result in frequent design failures. Here we introduce PrimerScore2, a robust high-throughput primer design tool that can design primers in one click for multiple PCR variants. It scores primers using a piecewise logistic model and the highest-scored primers are selected avoiding the issue of design failure and the necessity to loosen parameters to redesign, and it creatively evaluates specificity by predicting the efficiencies of all target/non-target products. To assess the prediction accuracy of the scores and efficiencies, two next generation sequencing (NGS) libraries were constructed-a 12-plex and a 57-plex-and the results showed that 17 out of 19 (89.5%) low-scoring pairs had a poor depth, 18 out of 19 (94.7%) high-scoring pairs had a high depth, and the depth ratios of the products were linearly correlated with the predicted efficiencies with a slope of 1.025 and a coefficient of determination (R²) 0.935. 116-plex and 114-plex anchored PCR panels designed by PrimerScore2 were applied to 26 maternal plasma samples with male fetuses, the results showed that the predicted fetal DNA fractions were concordant with fractions measured in gold standard method (Y fractions). PrimerScore2 was also used to design 77 monoplex Sanger sequencing primers, the sequencing results indicated that all the primers were effective.

Maximizing genetic gain through unlocking genetic variation in different ecotypes of kalmegh (Andrographis paniculata (Burm. f.) Nee)

Andrographis paniculata, commonly known as kalmegh is among the most popular medicinal herbs in Southeast Asia. It is widely cultivated for medicinal purposes. The bioactive molecule, Andrographolide accumulated in herb leaves has immense therapeutic and economic potential. However, comprehensive information regarding genetic diversity is very limited in this species. The present study assessed genetic diversity between and within the six populations (ecotypes) of twenty-four kalmegh accessions using multiple datasets (agro-morphological traits, phytochemical traits, and genic markers). This is the established report where EST-SSR (Expressed sequence tags-Simple Sequence Repeat) markers have been used to unlock genetic variation in kalmegh. Here, we identified and developed ninety-one metabolic pathway-specific EST-SSR markers. Finally, 32 random EST-SSR primer pairs were selected for genetic diversity assessment. Multivariate analysis to unveil the agro-morphological, phytochemical and genotypic variability was helpful in discriminating various germplasms studied in the present study. Among all the morphological discriptors used in present study, days to fifty percent flowering and dry herb yield were found as potential selection index for AP genetic improvement. Hierarchical cluster analysis built with agro-morphological data identified three major groups. However, corresponding analysis with phytochemical and molecular data generated two clear-cut groups among the studied individuals. Moreover, the grouping of individuals into different clusters using multiple datasets was geographically independent, and also showed inconsistency in grouping among agromorphological, phytochemical and molecular dataset based clusters. However, joint analysis using agro-morphological, phytochemical and genotypic information generated two genetic groups, which could be a valuable resource for identifying complementary crossing panels in the kalmegh breeding program. The accessions AP7, AP13, AP5, AP3 belong to cluster I and accessions AP17, AP18 belong to cluster II could be utilized as potential donors for high dry herb yield and andrographolide content, respectively in different selective breeding programs of AP. Thus, our results provided useful information about the overall genetic diversity and variation in economic traits useful for initiating selective breeding programs for contrasting traits of interest and maximizing genetic gain in kalmegh.

Genome-Wide Analysis of SAUR Gene Family Identifies a Candidate Associated with Fruit Size in Loquat (Eriobotrya japonica Lindl.)

Fruit size is an important fruit quality trait that influences the production and commodity values of loquats (Eriobotrya japonica Lindl.). The Small Auxin Upregulated RNA (SAUR) gene family has proven to play a vital role in the fruit development of many plant species. However, it has not been comprehensively studied in a genome-wide manner in loquats, and its role in regulating fruit size remains unknown. In this study, we identified 95 EjSAUR genes in the loquat genome. Tandem duplication and segmental duplication contributed to the expansion of this gene family in loquats. Phylogenetic analysis grouped the SAURs from Arabidopsis, rice, and loquat into nine clusters. By analyzing the transcriptome profiles in different tissues and at different fruit developmental stages and comparing two sister lines with contrasting fruit sizes, as well as by functional predictions, a candidate gene (EjSAUR22) highly expressed in expanding fruits was selected for further functional investigation. A combination of Indoleacetic acid (IAA) treatment and virus-induced gene silencing revealed that EjSAUR22 was not only responsive to auxin, but also played a role in regulating cell size and fruit expansion. The findings from our study provide a solid foundation for understanding the molecular mechanisms controlling fruit size in loquats, and also provide potential targets for manipulation of fruit size to accelerate loquat breeding.

Risk variants of obesity associated genes demonstrate BMI raising effect in a large cohort

Obesity is highly polygenic disease where several genetic variants have been reportedly associated with obesity in different ethnicities of the world. In the current study, we identified the obesity risk or protective association and BMI raising effect of the minor allele of adiponectin, C1Q and collagen domain containing (ADIPOQ), cholesteryl ester transfer protein (CEPT), FTO alpha-ketoglutarate dependent dioxygenase (FTO), leptin (LEP), and leptin receptor (LEPR) genes in a large cohort stratified into four BMI-based body weight categories i.e., normal weight, lean, over-weight, and obese. Based on selected candidate genetic markers, the genotyping of all study subjects was performed by PCR assays, and genotypes and allele frequencies were calculated. The minor allele frequencies (MAFs) of all genetic markers were computed for total and BMI-based body weight categories and compared with MAFs of global and South Asian (SAS) populations. Genetic associations of variants with obesity risk were calculated and BMI raising effect per copy of the minor allele were estimated. The genetic variants with higher MAFs in obese BMI group were; rs2241766 (G = 0.43), rs17817449 (G = 0.54), rs9939609 (A = 0.51), rs1421085 (C = 0.53), rs1558902 (A = 0.63), and rs1137101 (G = 0.64) respectively. All these variants were significantly associated with obesity (OR = 1.03–4.42) and showed a high BMI raising effect (β = 0.239–0.31 Kg/m²) per copy of the risk allele. In contrast, the MAFs of three variants were higher in lean-normal BMI groups; rs3764261 A = 0.38, rs9941349 T = 0.43, and rs7799039 G = 0.40–0.43). These variants showed obesity protective associations (OR = 0.68–0.76), and a BMI lowering effect per copy of the protective allele (β = -0.103–0.155 Kg/m²). The rs3764261 variant also showed significant and positive association with lean body mass (OR = 2.38, CI = 1.30–4.34). Overall, we report six genetic variants of ADIPOQ, FTO and LEPR genes as obesity-risk markers and a CETP gene variant as lean mass/obesity protective marker in studied Pakistani cohort.

Genetic Variants of ABC and SLC Transporter Genes and Chronic Myeloid Leukaemia: Impact on Susceptibility and Prognosis

Solute carrier (SLC) and ATP-binding cassette (ABC) transporters comprise a variety of proteins expressed on cell membranes responsible for intrusion or extrusion of substrates, respectively, including nutrients, xenobiotics, and chemotherapeutic agents. These transporters mediate the cellular disposition of tyrosine kinase inhibitors (TKIs), and their genetic variants could affect its function, potentially predisposing patients to chronic myeloid leukaemia (CML) and modulating treatment response. We explored the impact of genetic variability (single nucleotide variants—SNVs) of drug transporter genes (ABCB1, ABCG2, SLC22A1, and SLC22A5) on CML susceptibility, drug response, and BCR-ABL1 mutation status. We genotyped 10 SNVs by tetra-primers-AMRS-PCR in 198 CML patients and 404 controls, and assessed their role in CML susceptibility and prognosis. We identified five SNVs associated with CML predisposition, with some variants increasing disease risk, including TT genotype ABCB1 (rs1045642), and others showing a protective effect (GG genotype SLC22A5 rs274558). We also observed different haplotypes and genotypic profiles associated with CML predisposition. Relating to drug response impact, we found that CML patients with the CC genotype (rs2231142 ABCG2) had an increased risk of TKI resistance (six-fold). Additionally, CML patients carrying the CG genotype (rs683369 SLC22A1) presented a 4.54-fold higher risk of BCR-ABL1 mutations. Our results suggest that drug transporters’ SNVs might be involved in CML susceptibility and TKI response, and predict the risk of BCR-ABL1 mutations, highlighting the impact that SNVs could have in therapeutic selection.

Fine mapping of Rf2, a minor Restorer-of-fertility (Rf) gene for cytoplasmic male sterility in chili pepper G164 (Capsicum annuum L.)

Genome re-sequencing and recombination analyses identified Capana06g000193 as a strong candidate for the minor male fertility restoration locus Rf2 in chili pepper G164 harboring two dominant male fertility restoration genes. Male fertility restoration genes of chili pepper restorer line G164 (Capsicum annuum L.) were studied using molecular marker genotypes of an F₂ population (7G) of G164 crossed with the cytoplasmic male sterility line 77013A. The ratio of sterile to fertile single plants in the F₂ population was 1:15. This result indicates that chili pepper G164 has two dominant restoration genes, which we designated as Rf1 and Rf2. An individual plant recessive for Rf1 and heterozygous for Rf2, 7G-112 (rf1rf1Rf2rf2), was identified by molecular marker selection and genetic analysis, and a single Rf2 gene-segregating population with a 3:1 ratio of fertile to sterile plants was developed from the self-pollination of male fertile individuals of 77013A and 7G-112 hybrid progeny. Bulk segregant analysis of fertile and sterile pools from the segregating populations was used to genetically map Rf2 to a 3.1-Mb region on chromosome 6. Rf2 was further narrowed to a 179.3-kb interval through recombination analysis of molecular markers and obtained the most likely candidate gene, Capana06g000193.

High-Resolution Mapping of Barley mild mosaic virus Resistance Gene rym15

Barley yellow mosaic virus (BaYMV) and Barley mild mosaic virus (BaMMV), which are transmitted by the soil-borne plasmodiophorid Polymyxa graminis, cause high yield losses in barley. In previous studies, the recessive BaMMV resistance gene rym15, derived from the Japanese landrace Chikurin Ibaraki 1, was mapped on chromosome 6HS of Hordeum vulgare. In this study, 423 F₄ segmental recombinant inbred lines (RILs) were developed from crosses of Chikurin Ibaraki 1 with two BaMMV-susceptible cultivars, Igri (139 RILs) and Uschi (284 RILs). A set of 32 competitive allele-specific PCR (KASP) assays, designed using single nucleotide polymorphisms (SNPs) from the barley 50 K Illumina Infinium iSelect SNP chip, genotyping by sequencing (GBS) and whole-genome sequencing (WGS), was used as a backbone for construction of two high-resolution maps. Using this approach, the target locus was narrowed down to 0.161 cM and 0.036 cM in the Igri × Chikurin Ibaraki 1 (I × C) and Chikurin Ibaraki 1 × Uschi (C × U) populations, respectively. Corresponding physical intervals of 11.3 Mbp and 0.281 Mbp were calculated for I × C and C × U, respectively, according to the Morex v3 genome sequence. In the 0.281 Mbp target region, six high confidence (HC) and two low confidence (LC) genes were identified. Genome assemblies of BaMMV-susceptible cultivars Igri and Golden Promise from the barley pan-genome, and a HiFi assembly of Chikurin Ibaraki 1 together with re-sequencing data for the six HC and two LC genes in susceptible parental cultivar Uschi revealed functional SNPs between resistant and susceptible genotypes only in two of the HC genes. These SNPs are the most promising candidates for the development of functional markers and the two genes represent promising candidates for functional analysis.

Development and characterization of nitrogen and phosphorus use efficiency responsive genic and miRNA derived SSR markers in wheat

Among all the nutrients, nitrogen (N) and phosphorous (P) are the most limiting factors reducing wheat production and productivity world-wide. These macronutrients are directly applied to soil in the form of fertilizers. However, only 30-40% of these applied fertilizers are utilized by crop plants, while the rest is lost through volatilization, leaching, and surface run off. Therefore, to overcome the deficiency of N and P, it becomes necessary to improve their use efficiency. Marker-assisted selection (MAS) combined with traditional plant breeding approaches is considered best to improve the N and P use efficiency (N/PUE) of wheat varieties. In this study, we developed and evaluated a total of 98 simple sequence repeat (SSR) markers including 66 microRNAs and 32 gene-specific SSRs on a panel of 10 (N and P efficient/deficient) wheat genotypes. Out of these, 35 SSRs were found polymorphic and have been used for the study of genetic diversity and population differentiation. A set of two SSRs, namely miR171a and miR167a were found candidate markers able to discriminate contrasting genotypes for N/PUE, respectively. Therefore, these two markers could be used as functional markers for characterization of wheat germplasm for N and P use efficiency. Target genes of these miRNAs were found to be highly associated with biological processes (24 GO terms) as compared to molecular function and cellular component and shows differential expression under various P starving conditions and abiotic stresses.

GrainGenes: a data-rich repository for small grains genetics and genomics

As one of the US Department of Agriculture-Agricultural Research Service flagship databases, GrainGenes (https://wheat.pw.usda.gov) serves the data and community needs of globally distributed small grains researchers for the genetic improvement of the Triticeae family and Avena species that include wheat, barley, rye and oat. GrainGenes accomplishes its mission by continually enriching its cross-linked data content following the findable, accessible, interoperable and reusable principles, enhancing and maintaining an intuitive web interface, creating tools to enable easy data access and establishing data connections within and between GrainGenes and other biological databases to facilitate knowledge discovery. GrainGenes operates within the biological database community, collaborates with curators and genome sequencing groups and contributes to the AgBioData Consortium and the International Wheat Initiative through the Wheat Information System (WheatIS). Interactive and linked content is paramount for successful biological databases and GrainGenes now has 2917 manually curated gene records, including 289 genes and 254 alleles from the Wheat Gene Catalogue (WGC). There are >4.8 million gene models in 51 genome browser assemblies, 6273 quantitative trait loci and >1.4 million genetic loci on 4756 genetic and physical maps contained within 443 mapping sets, complete with standardized metadata. Most notably, 50 new genome browsers that include outputs from the Wheat and Barley PanGenome projects have been created. We provide an example of an expression quantitative trait loci track on the International Wheat Genome Sequencing Consortium Chinese Spring wheat browser to demonstrate how genome browser tracks can be adapted for different data types. To help users benefit more from its data, GrainGenes created four tutorials available on YouTube. GrainGenes is executing its vision of service by continuously responding to the needs of the global small grains community by creating a centralized, long-term, interconnected data repository. Database URL:https://wheat.pw.usda.gov.

Genomics-Assisted Improvement of Super High-Yield Hybrid Rice Variety "Super 1000" for Resistance to Bacterial Blight and Blast Diseases

The two-line rice hybrid "Super 1000" (GX24S × R900) represents a major landmark achievement of breeding for super-hybrid rice in China. However, both male parent R900 and hybrid "Super 1000" have an obvious defect of high susceptibility to rice bacterial blight (BB) and blast. Thus, improving disease resistance and maintaining the original high-yield capacity are essential for the sustainable application of "Super 1000." In this study, the application of closely linked single-nucleotide polymorphism (SNP) markers for foreground selection of dominant resistance gene loci together with genome-wide SNP markers for the background selection rapidly improved the disease resistance of R900 without disturbing its high-yield capacity. A series of improved R900 lines (iR900, in BC2Fn and BC3Fn generations) were developed to stack resistance genes (Xa23+Pi9, Xa23+Pi1+Pi2/9) by marker-assisted backcrossing and field selection for phenotypes, and further crossed with the female line GX24S to obtain improved hybrid variety Super 1000 (iS1000). The genetic backgrounds of iS1000 and "Super 1000" were profiled by using a 56 K SNP-Chip, and results showed that they shared 98.76% of similarity. Meanwhile, evaluation of the field disease resistance showed that the iR900 lines and iS1000 hybrids possess significantly enhanced resistance to both BB and rice blast. Resistance spectrum assays revealed that the iR900 lines and their derived hybrids exhibited high-level resistance to 28 Xoo strains tested, and enhanced resistance to leaf blast at the seedling stage when infected with 38 Magnaporthe oryzae isolates. Between 2019 and 2020, the multi-location field trials across the middle and lower reaches of the Yangtze River were launched and showed that the iS1000 slightly out-yielded than the original variety. In a large-scale demonstration site (6.73 ha, Yunnan, China), the iS1000 achieved 17.06 t/hm2 of yield in 2019. Moreover, the high similarity was observed in main agronomic traits and grain quality when comparing the improved lines/hybrids to original ones (iR900 vs. R900, iS1000 vs. S1000). This work presented a typical genomics-assisted breeding strategy and practice, which involves in directional introgression and rapid stack of multiple disease resistance genes, endowing the super-high-yield hybrid rice variety with holistic disease resistance but without yield penalty.

Determining Mitochondrial 3243A>G Heteroplasmy Using an ARMS-ddPCR Strategy

OBJECTIVES

Determining mitochondrial DNA (mtDNA) A-to-G substitution at nucleotide 3243 (m.3243A>G) heteroplasmy is essential for both precision diagnosis of m.3243A>G-associated mitochondrial disease and genetic counseling. Precise determination of m.3243A>G heteroplasmy is challenging, however, without appropriate strategies to accommodate heteroplasmic levels ranging from 1% to 100% in samples carrying thousands to millions of mtDNA copies.

METHODS

We used a combined strategy of amplification-refractory mutation system-quantitative polymerase chain reaction (ARMS-qPCR) and droplet digital PCR (ddPCR) to determine m.3243A>G heteroplasmy. Primers were specifically designed and screened for both ARMS-qPCR and ddPCR to determine m.3243A>G heteroplasmy. An optimized ARMS-qPCR-ddPCR-based strategy was established using artificial standards, with different mixtures of m.3243A-containing and m.3243G-containing plasmids and further tested using clinical samples containing the m.3243A>G mutation.

RESULTS

One of 20 primer pairs designed in the study was omitted for ARMS-qPCR-ddPCR strategy application according to criteria of 85% to 110%, R2> 0.98 amplification efficiency, melt curve with a single clear peak, and specificity for m.3243A and m.3243G artificial standards (|CtWt-CtMut|max). Using plasmid standards with various m.3243A>G heteroplasmy (1%-100%) at low, mid, and high copy numbers (3,000, 104, and 105-107, respectively) and DNA from the blood of 20 patients carrying m.3243A>G with mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes, we found that ARMS-qPCR was reliable for determining m.3243A>G at 3% to 100% for low copy number and 1% to 100% for mid to high copy number samples. Meanwhile, ddPCR was reliable for determining m.3243A>G at 1% to 100% at low to mid copy number samples.

CONCLUSIONS

An ARMS-qPCR-ddPCR-based strategy was successfully established for precise determination of m.3243A>G heteroplasmy in complex clinical samples.

SNP Discovery Using BSR-Seq Approach for Spot Blotch Resistance in Wheat (Triticum aestivum L.), an Essential Crop for Food Security

The pathogenic fungus, Bipolaris sorokiniana, that causes spot blotch (SB) disease of wheat, is a major production constraint in the Eastern Gangetic Plains of South Asia and other warm, humid regions of the world. A recombinant inbred line population was developed and phenotyped at three SB-prone locations in India. The single nucleotide polymorphism (SNP) for SB resistance was identified using a bulked segregant RNA-Seq-based approach, referred to as “BSR-Seq.” Transcriptome sequencing of the resistant parent (YS#24), the susceptible parent (YS#58), and their resistant and susceptible bulks yielded a total of 429.67 million raw reads. The bulk frequency ratio (BFR) of SNPs between the resistant and susceptible bulks was estimated, and selection of SNPs linked to resistance was done using sixfold enrichments in the corresponding bulks (BFR >6). With additional filtering criteria, the number of transcripts was further reduced to 506 with 1055 putative polymorphic SNPs distributed on 21 chromosomes of wheat. Based on SNP enrichment on chromosomal loci, five transcripts were found to be associated with SB resistance. Among the five SB resistance-associated transcripts, four were distributed on the 5B chromosome with putative 52 SNPs, whereas one transcript with eight SNPs was present on chromosome 3B. The SNPs linked to the trait were exposed to a tetra-primer ARMS-PCR assay, and an SNP-based allele-specific marker was identified for SB resistance. The in silico study of these five transcripts showed homology with pathogenesis-related genes; the metabolic pathway also exhibits similar results, suggesting their role in the plant defense mechanism.

SSR individual identification system construction and population genetics analysis for Chamaecyparis formosensis

Chamaecyparis formosensis is an endemic species of Taiwan, threatened from intensive use and illegal felling. An individual identification system for C. formosensis is required to provide scientific evidence for court use and deter illegal felling. In this study, 36 polymorphic simple sequence repeat markers were developed. By applying up to 28 non-linked of the developed markers, it is calculated that the cumulative random probability of identity (CP_I) is as low as 1.652 × 10^–12, and the identifiable population size is up to 60 million, which is greater than the known C. formosensis population size in Taiwan. Biogeographical analysis data show that C. formosensis from four geographic areas belong to the same genetic population, which can be further divided into three clusters: SY (Eastern Taiwan), HV and GW (Northwestern Taiwan), and MM (Southwestern Taiwan). The developed system was applied to assess the provenance of samples with 88.44% accuracy rate and therefore can serve as a prescreening tool to reduce the range required for comparison. The system developed in this study is a potential crime-fighting tool against illegal felling.

R3-MYB transcription factor LcMYBx from Litchi chinensis negatively regulates anthocyanin biosynthesis by ectopic expression in tobacco

Anthocyanin accumulation is one of the remarkable physiological changes during fruit ripening. In plants, anthocyanin synthesis is regulated by MYB activators, but the MYB repressors has been recognized recently. Here, we isolated a repressor of anthocyanin synthesis, LcMYBx, from Litchi chinensis Sonn. LcMYBx encoded a typical R3-MYB protein and contained a conserved [D/E]Lx2[R/K]x3Lx6Lx3R motif for interacting with bHLH proteins. Overexpression of LcMYBx in tobacco suppressed anthocyanin accumulation resulting in faded petals from pale-pink to almost white. Gene expression analysis showed the strong down-regulation of endogenous anthocyanin structural and regulatory genes by LcMYBx overexpression. Yeast two-hybrid and bimolecular fluorescence complementation assays indicated that LcMYBx could interact with the transcription factors LcbHLH1 and LcbHLH3. Transient promoter activation assays showed that LcMYBx could inhibit the activation capacity of LcMYB1-LcbHLH3 complex for LcDFR gene. These results suggest that LcMYBx competed with LcMYB1 to LcbHLHs, thus preventing the activation of LcDFR by LcMYB1-LcbHLHs complex and negatively controlling anthocyanin biosynthesis.