A viral histone H4 (CpBV-H4) is encoded in a polydnavirus, Cotesia plutellae bracovirus (CpBV). It plays a crucial role in parasitism of an endoparasitoid wasp, C. plutellae, against diamondback moth, Plutella xylostella, by altering host gene expression in an epigenetic mode by its N-terminal tail after joining host nucleosomes. Comparative transcriptomic analysis between parasitized and nonparasitized P. xylostella by RNA-Seq indicated that 1,858 genes were altered at more than two folds in expression levels at late parasitic stage, including 877 up-regulated genes and 981 down-regulated genes. Among parasitic factors altering host gene expression, CpBV-H4 alone explained 16.3% of these expressional changes. To characterize the joining sites of CpBV-H4 on host chromosomes, ChIP-Seq (chromatin immunoprecipitation followed by deep sequencing) was applied to chromatins extracted from parasitized larvae. It identified specific 538 ChIP targets. Joining sites were rich (60.2%) in AT sequence. Almost 40% of ChIP targets included short nucleotide repeat sequences presumably recognizable by transcriptional factors and chromatin remodeling factors. To further validate these CpBV-H4 targets, CpBV-H4 was transiently expressed in nonparasitized host at late larval stage and subjected to ChIP-Seq. Two kinds of ChIP-Seqs shared 51 core joining sites. Common targets were close (within 1 kb) to genes regulated at expression levels by CpBV-H4. However, other host genes not close to CpBV-H4 joining sites were also regulated by CpBV-H4. These results indicate that CpBV-H4 joins specific chromatin regions of P. xylostella and controls about one sixth of the total host genes that were regulated by C. plutellae parasitism in an epigenetic mode.

A very common protein modification in multicellular organisms is protein glycosylation or the addition of carbohydrate structures to the peptide backbone. Although the Class of the Insecta is the largest animal taxon on Earth, almost all information concerning glycosylation in insects is derived from studies with only one species, namely the fruit fly Drosophila melanogaster.

R2 non-LTR retrotransposons exclusively insert into the 28S rRNA genes of their host, and are expressed by co-transcription with the rDNA unit. The grasshopper Eyprepocnemis plorans contains transcribed rDNA clusters on most of its A chromosomes, as well as non-transcribed rDNA clusters on the parasitic B chromosomes found in many populations. Here the structure of the E. plorans R2 element, its abundance relative to the number of rDNA units and its retrotransposition activity were determined. Animals screened from five populations contained on average over 12,000 rDNA units on their A chromosomes, but surprisingly only about 100 R2 elements. Monitoring the patterns of R2 insertions in individuals from these populations revealed only low levels of retrotransposition. The low rates of R2 insertion observed in E. plorans differ from the high levels of R2 insertion previously observed in insect species that have many fewer rDNA units. It is proposed that high levels of R2 are strongly selected against in E. plorans, because the rDNA transcription machinery in this species is unable to differentiate between R2-inserted and uninserted units. The B chromosomes of E. plorans contain an additional 7,000 to 15,000 rDNA units, but in contrast to the A chromosomes, from 150 to over 1,500 R2 elements. The higher concentration of R2 in the inactive B chromosomes rDNA clusters suggests these chromosomes can act as a sink for R2 insertions thus further reducing the level of insertions on the A chromosomes. These studies suggest an interesting evolutionary relationship between the parasitic B chromosomes and R2 elements.

Tissue culture is performed to maintain isolated portions of multicellular organisms in an artificial milieu that is outside the individual organism and for considerable periods of time; cells derived from cultured explants are, in general, different from cells of the corresponding tissue in a living organism. The changes in cultured tissues that precede and often explain the subsequent cell proliferation of explant-derived cells have been partially studied, but little is known about the molecular and genomic basis of these changes. Comparative transcriptomics of intact and cultured (90 hours in MGM-450 insect medium) Bombyx mori tissues revealed that fewer genes represented a larger portion of the transcriptome of intact fat body tissues than of cultured fat body tissues. This analysis also indicated that expression of genes encoding sugar transporters and immune response proteins increased during culture and that expression of genes encoding lipoproteins and cuticle proteins decreased during culture. These results provide support for hypotheses that cultured tissues respond immunologically to surgery, adapt to the medium by accelerating sugar uptake, and terminate their identity as part of an intact organism by becoming independent of that organism.

Unlocking the vast genomic diversity stored in natural history collections would create unprecedented opportunities for genome-scale evolutionary, phylogenetic, domestication and population genomic studies. Many researchers have been discouraged from using historical specimens in molecular studies because of both generally limited success of DNA extraction and the challenges associated with PCR-amplifying highly degraded DNA. In today's next-generation sequencing (NGS) world, opportunities and prospects for historical DNA have changed dramatically, as most NGS methods are actually designed for taking short fragmented DNA molecules as templates. Here we show that using a standard multiplex and paired-end Illumina sequencing approach, genome-scale sequence data can be generated reliably from dry-preserved plant, fungal and insect specimens collected up to 115 years ago, and with minimal destructive sampling. Using a reference-based assembly approach, we were able to produce the entire nuclear genome of a 43-year-old Arabidopsis thaliana (Brassicaceae) herbarium specimen with high and uniform sequence coverage. Nuclear genome sequences of three fungal specimens of 22-82 years of age (Agaricus bisporus, Laccaria bicolor, Pleurotus ostreatus) were generated with 81.4-97.9% exome coverage. Complete organellar genome sequences were assembled for all specimens. Using de novo assembly we retrieved between 16.2-71.0% of coding sequence regions, and hence remain somewhat cautious about prospects for de novo genome assembly from historical specimens. Non-target sequence contaminations were observed in 2 of our insect museum specimens. We anticipate that future museum genomics projects will perhaps not generate entire genome sequences in all cases (our specimens contained relatively small and low-complexity genomes), but at least generating vital comparative genomic data for testing (phylo)genetic, demographic and genetic hypotheses, that become increasingly more horizontal. Furthermore, NGS of historical DNA enables recovering crucial genetic information from old type specimens that to date have remained mostly unutilized and, thus, opens up a new frontier for taxonomic research as well.

The nuclear genomes of most animal species include NUMTs, segments of the mitogenome incorporated into their chromosomes. Although NUMT counts are known to vary greatly among species, there has been no comprehensive study of their frequency/attributes in the most diverse group of terrestrial organisms, insects. This study examines NUMTs derived from a 658 bp 5' segment of the cytochrome c oxidase I (COI) gene, the barcode region for the animal kingdom. This assessment is important because unrecognized NUMTs can elevate estimates of species richness obtained through DNA barcoding and derived approaches (eDNA, metabarcoding). This investigation detected nearly 10,000 COI NUMTs ≥ 100 bp in the genomes of 1,002 insect species (range = 0-443). Variation in nuclear genome size explained 56% of the mitogenome-wide variation in NUMT counts. Although insect orders with the largest genome sizes possessed the highest NUMT counts, there was considerable variation among their component lineages. Two thirds of COI NUMTs possessed an IPSC (indel and/or premature stop codon) allowing their recognition and exclusion from downstream analyses. The remainder can elevate species richness as they showed 10.1% mean divergence from their mitochondrial homologue. The extent of exposure to "ghost species" is strongly impacted by the target amplicon's length. NUMTs can raise apparent species richness by up to 22% when a 658 bp COI amplicon is examined versus a doubling of apparent richness when 150 bp amplicons are targeted. Given these impacts, metabarcoding and eDNA studies should target the longest possible amplicons while also avoiding use of 12S/16S rDNA as they triple NUMT exposure because IPSC screens cannot be employed.

Fungi are key dietary resources for many animals. Fungi, in consequence, have evolved sophisticated physical and chemical defences for repelling and impairing fungivores. Expression of such defences may entail costs, requiring diversion of energy and nutrients away from fungal growth and reproduction. Inducible resistance that is mounted after attack by fungivores may allow fungi to circumvent the potential costs of defence when not needed. However, no information exists on whether fungi display inducible resistance. We combined organism and fungal gene expression approaches to investigate whether fungivory induces resistance in fungi.

Genome sequencing projects have been completed for several species representing four highly diverged holometabolous insect orders, Diptera, Hymenoptera, Coleoptera, and Lepidoptera. The striking evolutionary diversity of insects argues a need for efficient methods to apply genome information from such models to genetically uncharacterized species. Constructing conserved synteny maps plays a crucial role in this task. Here, we demonstrate the use of fluorescence in situ hybridization with bacterial artificial chromosome probes as a powerful tool for physical mapping of genes and comparative genome analysis in Lepidoptera, which have numerous and morphologically uniform holokinetic chromosomes.

Among insect taxa, ants exhibit one of the most variable chromosome numbers ranging from n = 1 to n = 60. This high karyotype diversity is suggested to be correlated to ants diversification. The karyotype evolution of ants is usually understood in terms of Robertsonian rearrangements towards an increase in chromosome numbers. The ant genus Mycetophylax is a small monogynous basal Attini ant (Formicidae: Myrmicinae), endemic to sand dunes along the Brazilian coastlines. A recent taxonomic revision validates three species, Mycetophylax morschi, M. conformis and M. simplex. In this paper, we cytogenetically characterized all species that belongs to the genus and analyzed the karyotypic evolution of Mycetophylax in the context of a molecular phylogeny and ancestral character state reconstruction. M. morschi showed a polymorphic number of chromosomes, with colonies showing 2n = 26 and 2n = 30 chromosomes. M. conformis presented a diploid chromosome number of 30 chromosomes, while M. simplex showed 36 chromosomes. The probabilistic models suggest that the ancestral haploid chromosome number of Mycetophylax was 17 (Likelihood framework) or 18 (Bayesian framework). The analysis also suggested that fusions were responsible for the evolutionary reduction in chromosome numbers of M. conformis and M. morschi karyotypes whereas fission may determines the M. simplex karyotype. These results obtained show the importance of fusions in chromosome changes towards a chromosome number reduction in Formicidae and how a phylogenetic background can be used to reconstruct hypotheses about chromosomes evolution.

Insects and plants have coexisted for over 350 million years and their interactions have affected ecosystems and agricultural practices worldwide. Variation in herbivorous insects' virulence to circumvent host resistance has been extensively documented. However, despite decades of investigation, the genetic foundations of virulence are currently unknown. The brown planthopper (Nilaparvata lugens) is the most destructive rice (Oryza sativa) pest in the world. The identification of the resistance gene Bph1 and its introduction in commercial rice varieties prompted the emergence of a new virulent brown planthopper biotype that was able to break the resistance conferred by Bph1. In this study, we aimed to construct a high density linkage map for the brown planthopper and identify the loci responsible for its virulence in order to determine their genetic architecture. Based on genotyping data for hundreds of molecular markers in three mapping populations, we constructed the most comprehensive linkage map available for this species, covering 96.6% of its genome. Fifteen chromosomes were anchored with 124 gene-specific markers. Using genome-wide scanning and interval mapping, the Qhp7 locus that governs preference for Bph1 plants was mapped to a 0.1 cM region of chromosome 7. In addition, two major QTLs that govern the rate of insect growth on resistant rice plants were identified on chromosomes 5 (Qgr5) and 14 (Qgr14). This is the first study to successfully locate virulence in the genome of this important agricultural insect by marker-based genetic mapping. Our results show that the virulence which overcomes the resistance conferred by Bph1 is controlled by a few major genes and that the components of virulence originate from independent genetic characters. The isolation of these loci will enable the elucidation of the molecular mechanisms underpinning the rice-brown planthopper interaction and facilitate the development of durable approaches for controlling this most destructive agricultural insect.

F-box proteins play critical roles in plant responses to biotic/abiotic stresses. In the present study, a total of 68 wheat F-box/Kelch (TaFBK) genes, unevenly distributed across 21 chromosomes and encoding 74 proteins, were identified in EnsemblPlants. Protein sequences were compared with those of Arabidopsis and three cereal species by phylogenetic and domain analyses, where the wheat sequences were resolved into 6 clades. In silico analysis of a digital PCR dataset revealed that TaFBKs were expressed at multiple developmental stages and tissues, and in response to drought and/or heat stresses. The TaFBK19 gene, a homolog of the Attenuated Far-Red Response (AFR) genes in other plant species, and hence named TaAFR, was selected for further analysis. Reverse-transcription quantitative real-time PCR (RT-qPCR) was carried out to determine tissue-specific, hormone and stress (abiotic/biotic) responsive expression patterns. Of interest, TaAFR was expressed most abundantly in the leaves, and its expression in response to leaf rust variants suggests a potential role in compatible vs incompatible rust responses. The protein was predicted to localize in cytosol, but it was shown experimentally to localize in both the cytosol and the nucleus of tobacco. A series of protein interaction studies, starting with a yeast-2-hybrid (Y2H) library screen (wheat leaf infected with incompatible leaf rust pathogens), led to the identification of three TaAFR interacting proteins. Skp1/ASK1-like protein (Skp1) was found to interact with the F-box domain of TaAFR, while ADP-ribosylation factor 2-like isoform X1 (ARL2) and phenylalanine ammonia-lyase (PAL) were shown to interact with its Kelch domain. The data presented herein provides a solid foundation from which the function and metabolic network of TaAFR and other wheat FBKs can be further explored.

The sophisticated colony organization of eusocial insects is primarily maintained through the utilization of pheromones. The regulation of these complex social interactions requires intricate chemoreception systems. The recent publication of the genome of Zootermopsis nevadensis opened a new avenue to study molecular basis of termite caste systems. Although there has been a growing interest in the termite chemoreception system that regulates their sophisticated caste system, the relationship between division of labor and expression of chemoreceptor genes remains to be explored. Using high-throughput mRNA sequencing (RNA-seq), we found several chemoreceptors that are differentially expressed among castes and between sexes in a subterranean termite Reticulitermes speratus. In total, 53 chemoreception-related genes were annotated, including 22 odorant receptors, 7 gustatory receptors, 12 ionotropic receptors, 9 odorant-binding proteins, and 3 chemosensory proteins. Most of the chemoreception-related genes had caste-related and sex-related expression patterns; in particular, some chemoreception genes showed king-biased or queen-biased expression patterns. Moreover, more than half of the genes showed significant age-dependent differences in their expression in female and/or male reproductives. These results reveal a strong relationship between the evolution of the division of labor and the regulation of chemoreceptor gene expression, thereby demonstrating the chemical communication and underlining chemoreception mechanism in social insects.

In most insect species, a variety of serine protease inhibitors (SPIs) have been found in multiple tissues, including integument, gonad, salivary gland, and hemolymph, and are required for preventing unwanted proteolysis. These SPIs belong to different families and have distinct inhibitory mechanisms. Herein, we predicted and characterized potential SPI genes based on the genome sequences of silkworm, Bombyx mori. As a result, a total of eighty SPI genes were identified in B. mori. These SPI genes contain 10 kinds of SPI domains, including serpin, Kunitz_BPTI, Kazal, TIL, amfpi, Bowman-Birk, Antistasin, WAP, Pacifastin, and alpha-macroglobulin. Sixty-three SPIs contain single SPI domain while the others have at least two inhibitor units. Some SPIs also contain non-inhibitor domains for protein-protein interactions, including EGF, ADAM_spacer, spondin_N, reeler, TSP_1 and other modules. Microarray analysis showed that fourteen SPI genes from lineage-specific TIL family and Group F of serpin family had enriched expression in the silk gland. The roles of SPIs in resisting pathogens were investigated in silkworms when they were infected by four pathogens. Microarray and qRT-PCR experiments revealed obvious up-regulation of 8, 4, 3 and 3 SPI genes after infection with Escherichia coli, Bacillus bombysepticus, Beauveria bassiana or B. mori nuclear polyhedrosis virus (BmNPV), respectively. On the contrary, 4, 11, 7 and 9 SPI genes were down-regulated after infection with E. coli, B. bombysepticus, B. bassiana or BmNPV, respectively. These results suggested that these SPI genes may be involved in resistance to pathogenic microorganisms. These findings may provide valuable information for further clarifying the roles of SPIs in the development, immune defence, and efficient synthesis of silk gland protein.

Insect saliva plays an important role in modulation of plant-insect interactions. Although this area of research has generated much attention in recent years, mechanisms of how saliva affects plant responses remain poorly understood. To address this void, the present study investigated the impact of the brown planthopper (Nilaparvata lugens, Stål; hereafter BPH) salivary gland extract (SGE) on rice (Oryza sativa) systemic responses at the mRNA level. Differentially expressed rice mRNAs were generated through suppression subtractive hybridization (SSH) and classified into six functional groups. Those with the most representatives were from the primary metabolism (28%), signaling-defense (22%) and transcription-translation-regulation group (16%). To validate SSH library results, six genes were further analyzed by One-Step Real-Time Reverse Transcriptase-PCR. Five of these genes exhibited up-regulation levels of more than 150% of those in the control group in at least one post-application time point. Results of this study allow assignment of at least two putative roles of BPH saliva: First, application of SGE induces immediate systemic responses at the mRNA level, suggesting that altering of the rice transcriptome at sites distant to hoppers feeding locations may play an important role in BPH-rice interactions. Second, 58% of SGE-responsive up-regulated genes have a secondary function associated with senescence, a process characterized by remobilization of nutrients. This suggests that BPH salivary secretions may reprogram the rice transcriptome for nutritional enhancement. When these findings are translated onto 'whole plant' scale, they indicate that BPH saliva may play the 'wise investment' role of 'minimum input today, maximum output tomorrow'.

Parasitic lice demonstrate an unusual array of mitochondrial genome architectures and gene arrangements. We characterized the mitochondrial genome of Geomydoecus aurei, a chewing louse (Phthiraptera: Trichodectidae) found on pocket gophers (Rodentia: Geomyidae) using reads from both Illumina and Oxford Nanopore sequencing coupled with PCR, cloning, and Sanger sequencing to verify structure and arrangement for each chromosome. The genome consisted of 12 circular mitochondrial chromosomes ranging in size from 1,318 to 2,088 nucleotides (nt). Total genome size was 19,015 nt. All 37 genes typical of metazoans (2 rRNA genes, 22 tRNA genes, and 13 protein-coding genes) were present. An average of 26% of each chromosome was composed of non-gene sequences. Within the non-gene region of each chromosome, there was a 79-nt nucleotide sequence that was identical among chromosomes and a conserved sequence with secondary structure that was always followed by a poly-T region. We hypothesize that these regions may be important in the initiation of transcription and DNA replication, respectively. The G. aurei genome shares 8 derived gene clusters with other chewing lice of mammals, but in G. aurei, genes on several chromosomes are not contiguous.

The Adaptor Protein (AP)-3 complex is an evolutionary conserved, molecular sorting device that mediates the intracellular trafficking of proteins to lysosomes and related organelles. Genetic defects in AP-3 subunits lead to impaired biogenesis of lysosome-related organelles (LROs) such as mammalian melanosomes and insect eye pigment granules. In this work, we have performed a forward screening for genetic modifiers of AP-3 function in the fruit fly, Drosophila melanogaster. Specifically, we have tested collections of large multi-gene deletions--which together covered most of the autosomal chromosomes-to identify chromosomal regions that, when deleted in single copy, enhanced or ameliorated the eye pigmentation phenotype of two independent AP-3 subunit mutants. Fine-mapping led us to define two non-overlapping, relatively small critical regions within fly chromosome 3. The first critical region included the Atg2 gene, which encodes a conserved protein involved in autophagy. Loss of one functional copy of Atg2 ameliorated the pigmentation defects of mutants in AP-3 subunits as well as in two other genes previously implicated in LRO biogenesis, namely Blos1 and lightoid, and even increased the eye pigment content of wild-type flies. The second critical region included the ArfGAP1 gene, which encodes a conserved GTPase-activating protein with specificity towards GTPases of the Arf family. Loss of a single functional copy of the ArfGAP1 gene ameliorated the pigmentation phenotype of AP-3 mutants but did not to modify the eye pigmentation of wild-type flies or mutants in Blos1 or lightoid. Strikingly, loss of the second functional copy of the gene did not modify the phenotype of AP-3 mutants any further but elicited early lethality in males and abnormal eye morphology when combined with mutations in Blos1 and lightoid, respectively. These results provide genetic evidence for new functional links connecting the machinery for biogenesis of LROs with molecules implicated in autophagy and small GTPase regulation.

MicroRNAs (miRNAs) are a group of small RNAs involved in various biological processes through negative regulation of mRNAs at the post-transcriptional level. Although miRNA profiles have been documented in over two dozen insect species, few are agricultural pests. In this study, both conserved and novel miRNAs in the diamondback moth, Plutella xylostella L., a devastating insect pest of cruciferous crops worldwide, were documented. High-throughput sequencing of a small RNA library constructed from a mixed life stages of P. xylostella, including eggs, 1st to 4th (last) instar larvae, pupae and adults, identified 384 miRNAs, of which 174 were P. xylostella specific. In addition, temporal expressions of 234 miRNAs at various developmental stages were investigated using a customized microarray analysis. Among the 91 differentially expressed miRNAs, qRT-PCR analysis was used to validate highly expressed miRNAs at each stage. The combined results not only systematically document miRNA profiles in an agriculturally important insect pest, but also provide molecular targets for future functional analysis and, ultimately, genetic-based pest control practice.

Satellite repetitive sequences that accumulate in the heterochromatin consist a large fraction of a genome and due to their properties are suggested to be implicated in centromere function. Current knowledge of heterochromatic regions of Bactrocera oleae genome, the major pest of the olive tree, is practically nonexistent. In our effort to explore the repetitive DNA portion of B. oleae genome, a novel satellite sequence designated BoR300 was isolated and cloned. The present study describes the genomic organization, abundance and chromosomal distribution of BoR300 which is organized in tandem, forming arrays of 298 bp-long monomers. Sequence analysis showed an AT content of 60.4%, a CENP-B like-motif and a high curvature value based on predictive models. Comparative analysis among randomly selected monomers demonstrated a high degree of sequence homogeneity (88%-97%) of BoR300 repeats, which are present at approximately 3,000 copies per haploid genome accounting for about 0.28% of the total genomic DNA, based on two independent qPCR approaches. In addition, expression of the repeat was also confirmed through RT-PCR, by which BoR300 transcripts were detected in both sexes. Fluorescence in situ hybridization (FISH) of BoR300 on mitotic metaphases and polytene chromosomes revealed signals to the centromeres of two out of the six chromosomes which indicated a chromosome-specific centromeric localization. Moreover, BoR300 is not conserved in the closely related Bactrocera species tested and it is also absent in other dipterans, but it's rather restricted to the B. oleae genome. This feature of species-specificity attributed to BoR300 satellite makes it a good candidate as an identification probe of the insect among its relatives at early development stages.

Prophages are highly dynamic components in the bacterial genome and play an important role in intraspecies variations. There are at least two prophages in the chromosomes of Candidatus Liberibacter asiaticus' (Las) Floridian isolates. Las is both unculturable and the most prevalent species of Liberibacter pathogens that cause huanglongbing (HLB), a worldwide destructive disease of citrus. In this study, seven new prophage variants resulting from two hyper-variable regions were identified by screening clone libraries of infected citrus, periwinkle and psyllids. Among them, Types A and B share highly conserved sequences and localize within the two prophages, FP1 and FP2, respectively. Although Types B and C were abundant in all three libraries, Type A was much more abundant in the libraries from the Las-infected psyllids than from the Las-infected plants, and Type D was only identified in libraries from the infected host plants but not from the infected psyllids. Sequence analysis of these variants revealed that the variations may result from recombination and rearrangement events. Conventional PCR results using type-specific molecular markers indicated that A, B, C and D are the four most abundant types in Las-infected citrus and periwinkle. However, only three types, A, B and C are abundant in Las-infected psyllids. Typing results for Las-infected citrus field samples indicated that mixed populations of Las bacteria present in Floridian isolates, but only the Type D population was correlated with the blotchy mottle symptom. Extended cloning and sequencing of the Type D region revealed a third prophage/phage in the Las genome, which may derive from the recombination of FP1 and FP2. Dramatic variations in these prophage regions were also found among the global Las isolates. These results are the first to demonstrate the prophage/phage-mediated dynamics of Las populations in plant and insect hosts, and their correlation with insect transmission and disease development.

The common cutworm (CCW; Spodoptera litura Fabricius) is a serious herbivorous insect pest of soybean (Glycine max) in Asia and Oceania. Previously, we identified quantitative trait loci (QTLs) for CCW-antibiosis-resistance, CCW-1 and CCW-2, and antixenosis-resistance, qRslx1 and qRslx2, in the cultivar 'Himeshirazu'. The effects of these QTLs are useful in the breeding of CCW-resistant cultivars. In this study, we conducted an antixenosis bioassay on CCW using recombinant inbred lines derived from a cross between a wild soybean (Glycine soja) and the leading cultivar 'Fukuyutaka' to identify CCW-resistance genes in G. soja. The QTL analysis revealed six and four novel antixenosis-resistance QTLs in 2012 and 2013, respectively. Among them, the QTLs on chromosomes 2 and 7, designated qRslx4 and qRslx3, respectively, were stably detected in both years. qRslx3 exhibited the largest effect in both years, suggesting that qRslx3 can be exploited in the breeding of CCW-resistant soybean. Furthermore, qRslx3 and qRslx4 can be used, along with previously reported QTLs from 'Himeshirazu', to enhance the CCW-resistance of soybean cultivars because their chromosomal positions are unique. These new CCW-resistance QTLs from G. soja should play important roles in the breeding of CCW-resistant soybean cultivars.