Genomic Regions Research Articles

Regulatory transposable elements in the encyclopedia of DNA elements

Transposable elements (TEs) comprise ~50% of our genome, but knowledge of how TEs affect genome evolution remains incomplete. Leveraging ENCODE4 data, we provide the most comprehensive study to date of TE contributions to the regulatory genome. We find 236,181 (~25%) human candidate cis-regulatory elements (cCREs) are TE-derived, with over 90% lineage-specific since the human-mouse split, accounting for 8–36% of lineage-specific cCREs. Except for SINEs, cCRE-associated transcription factor (TF) motifs in TEs are derived from ancestral TE sequence more than expected by chance. We show that TEs may adopt similar regulatory activities of elements near their integration site. Since human-mouse divergence, TEs have contributed 3–56% of TF binding site turnover events across 30 examined TFs. Finally, TE-derived cCREs are similar to non-TE cCREs in terms of MPRA activity and GWAS variant enrichment. Overall, our results substantiate the notion that TEs have played an important role in shaping the human regulatory genome.

Resistance of the CRISPR-Cas13a Gene-Editing System to Potato Spindle Tuber Viroid Infection in Tomato and Nicotiana benthamiana.

Gene-editing technology, specifically the CRISPR-Cas13a system, has shown promise in breeding plants resistant to RNA viruses. This system targets RNA and, theoretically, can also combat RNA-based viroids. To test this, the CRISPR-Cas13a system was introduced into tomato plants via transient expression and into Nicotiana benthamiana through transgenic methods, using CRISPR RNAs (crRNAs) targeting the conserved regions of both sense and antisense genomes of potato spindle tuber viroid (PSTVd). In tomato plants, the expression of CRISPR-Cas13a and crRNAs substantially reduced PSTVd accumulation and alleviated disease symptoms. In transgenic N. benthamiana plants, the PSTVd levels were lower as compared to wild-type plants. Several effective crRNAs targeting the PSTVd genomic RNA were also identified. These results demonstrate that the CRISPR-Cas13a system can effectively target and combat viroid RNAs, despite their compact structures.

Identification of a major QTL region and candidate genes responsible for harvest date combining QTL mapping and QTL-seq approaches in walnut (Juglans regia L.)

Walnut (Juglans regia L.) takes place within the Juglandaceae family and is one of the most popular nut species. In walnut breeding, several limiting factors need to be revealed regarding the genetic mechanism underlying important agronomical traits. In the current report, the monitoring of the ‘Chandler x Kaplan-86’ F1 walnut population was carried out four successive years (2015, 2016, 2017, 2018) and the distribution of the plants for harvest date from four years was ranged from August 27th to October 22nd. Then, we dissected the genomic regions associated with harvest date using classical Quantitative Trait Loci (QTLs) mapping and a rapid and cost-effective QTL-seq approach in ‘Chandler x Kaplan-86’ F1 walnut population for the first time. In classical QTL mapping, a major QTL was recognized in LG4 corresponding to Chr1 in Chandler female map with 27.6 (%) phenotyping variance. In QTL-seq, we bulked the DNAs of the genotypes with extreme harvest date separately, the whole genome resequencing was performed for a parent Chandler and both bulks (early and late) to calculate Delta SNP indexes for each Single Nucleotide Polymorphisms (SNPs) and Insertion/Deletion (InDels). Approximately 48.04 million paired reads for Chandler, 149.72, and 150.14 million paired reads were obtained from early and late harvest bulks, respectively. A total of 229,307 SNPs and 45,223 InDels were mined from the merged sequences data covered whole chromosomes. The responsible genomic fragments were identified in Chr2 and Chr11 in QTL-seq. The highest linked InDel locus detected in QTL-seq in Chr11 around 17 million base pairs distinguished the plants according to harvest date. InDel Sequence Characterized Amplified Region (SCAR) marker is a candidate for marker-assisted selection for future walnut breeding studies.

Adaptation in landlocked Atlantic salmon links genetics in wild and farmed salmon to smoltification

Increased knowledge of heritable traits in Atlantic Salmon (Salmo salar) is important to overcome bottlenecks in salmonid aquaculture. Atlantic salmonid populations, both landlocked and anadromous, represent an interesting model to gain insight into anadromy related traits, most notably, the probability to smoltify. While a previous study has identified several genomic regions diverging between anadromous and landlocked populations across the species range, the present study explores these data further with the aim to uncover if some of these genomic regions are linked to beneficial genetic traits associated with smoltification. In this study 17 of these loci were monitored in 669 anadromous salmon originating from 36 full-sibling families that had been reared under common garden conditions. The Smolt Index was calculated, using multiple visual markers, and provided a means of assessing smoltification stage. One SNP, located in Ssa04, showed a significant association with probability to smoltify, where individuals homozygous for the landlocked variant (LL) displayed a decrease in probability of smoltifying after one winter when compared with the homozygous for the anadromous variant (AA). This effect was independent of individual fish size. A separate common garden study comprising 200 individuals from either anadromous or landlocked strains showed that expression levels of ncor1, a thyroid mediator hormone located on the same chromosomal region (Ssa04), were significantly reduced in landlocked individuals post smoltification but remained constant in their anadromous counterparts. This study therefore suggests that while size is still the most important trigger for the induction of smoltification, there may also be an additional genetic component or trigger that has been ‘lost’ during the years deprived of SW transfer. In conclusion, the LL genotype identified here could potentially be used by the industry to delay smoltification and may also represent one of the first clues to the genetic regulation of smoltification in Atlantic salmon.

Genomic insights into the modifications of spike morphology traits during wheat breeding.

Over the past century, environmental changes have significantly impacted wheat spike morphology, crucial for adaptation and grain yield. However, the changes in wheat spike modifications during this periodremain largely unknown. This study examines 16 spike morphology traits in 830 accessions released from 1900 to 2020. It finds that spike weight, grain number per spike (GN), and thousand kernel weight have significantly increased, while spike length has no significant change. The increase in fertile spikelets is due to fewer degenerated spikelets, resulting in a higher GN. Genome-wide association studies identified 49,994 significant SNPs, grouped into 293 genomic regions. The accumulation of favorable alleles in these genomic regions indicates the genetic basis for modification in spike morphology traits. Genetic network analysis of these genomic regions reveals the genetic basis for phenotypic correlations among spike morphology traits. The haplotypes of the identified genomic regions display obvious geographical differentiation in global accessions and environmental adaptation over the past 120 years. In summary, we reveal the genetic basis of adaptive evolution and the interactions of spike morphology, offering valuable resources for the genetic improvement of spike morphology to enhance environmental adaptation.

Methylation-Associated Nucleosomal Patterns of Cell-Free DNA in Cancer Patients and Pregnant Women.

Cell-free DNA (cfDNA) analysis offers an attractive noninvasive means of detecting and monitoring diseases. cfDNA cleavage patterns within a short range (e.g., 11 nucleotides) have been reported to correlate with cytosine-phosphate-guanine (CpG) methylation, allowing fragmentomics-based methylation analysis (FRAGMA). Here, we adopted FRAGMA to the extended region harboring multiple nucleosomes, termed FRAGMAXR. We profiled cfDNA nucleosomal patterns over the genomic regions from -800 to 800 bp surrounding differentially methylated CpG sites, harboring approximately 8 nucleosomes, referred to as CpG-associated cfDNA nucleosomal patterns. Such nucleosomal patterns were analyzed by FRAGMAXR in cancer patients and pregnant women. We identified distinct cfDNA nucleosomal patterns around differentially methylated CpG sites. Compared with subjects without cancer, patients with hepatocellular carcinoma (HCC) showed reduced amplitude of nucleosomal patterns, with a gradual decrease over tumor stages. Nucleosomal patterns associated with differentially methylated CpG sites could be used to train a machine learning model, resulting in the detection of HCC patients with an area under the receiver operating characteristic curve of 0.93. We further demonstrated the feasibility of multicancer detection using a dataset comprising lung, breast, and ovarian cancers. The tissue-of-origin analysis of plasma cfDNA from pregnant women and cancer patients revealed that the placental DNA and tumoral DNA contributions deduced by FRAGMAXR correlated well with values measured using genetic variants (Pearson r: 0.85 and 0.94, respectively). CpG-associated cfDNA nucleosomal patterns of cfDNA molecules are influenced by DNA methylation and might be useful for biomarker developments for cancer liquid biopsy and noninvasive prenatal testing.

Temperature affects recombination rate plasticity and meiotic success between thermotolerant and cold tolerant yeast species.

Meiosis is required for the formation of gametes in all sexually reproducing species and the process is well conserved across the tree of life. However, meiosis is sensitive to a variety of external factors, which can impact chromosome pairing, recombination, and fertility. For example, the optimal temperature for successful meiosis varies between species of plants and animals. This suggests that meiosis is temperature sensitive, and that natural selection may act on variation in meiotic success as organisms adapt to different environmental conditions. To understand how temperature alters the successful completion of meiosis, we utilized two species of the budding yeast Saccharomyces with different temperature preferences: thermotolerant Saccharomyces cerevisiae and cold tolerant Saccharomyces uvarum. We surveyed three metrics of meiosis: sporulation efficiency, spore viability, and recombination rate in multiple strains of each species. As per our predictions, the proportion of cells that complete meiosis and form spores is temperature sensitive, with thermotolerant S. cerevisiae having a higher temperature threshold for successful meiosis than cold tolerant S. uvarum. We confirmed previous observations that S. cerevisiae recombination rate varies between strains and across genomic regions, and add new results that S. uvarum has higher recombination rates than S. cerevisiae. We find that temperature significantly influences recombination rate plasticity in S. cerevisiae and S. uvarum, in agreement with studies in animals and plants. Overall, these results suggest that meiotic thermal sensitivity is associated with organismal thermal tolerance, and may even result in temporal reproductive isolation as populations diverge in thermal profiles.

Digest: Large-effect loci involved in maturation timing affect the early survival of Atlantic salmon.

How might variation in genomic regions that impact many traits modulate fitness across different life stages? Aykanat et al. (2024) show that two large-effect loci associated with age at maturity, six6 and vgll3, impact the survival of wild Atlantic salmon (Salmo salar) through nutrient-dependent, indirect genetic effects. Specifically, the late maturation allele in parental vgll3, and the early maturation allele in maternal six6, increase survival in early life under high nutrient conditions.

Combining RNA Interference and RIG-I Activation to Inhibit Hepatitis E Virus Replication.

Hepatitis E virus (HEV) poses a significant global health threat, with an estimated 20 million infections occurring annually. Despite being a self-limiting illness, in most cases, HEV infection can lead to severe outcomes, particularly in pregnant women and individuals with pre-existing liver disease. In the absence of specific antiviral treatments, the exploration of RNAi interference (RNAi) as a targeted strategy provides valuable insights for urgently needed therapeutic interventions against Hepatitis E. We designed small interfering RNAs (siRNAs) against HEV, which target the helicase domain and the open reading frame 3 (ORF3). These target regions will reduce the risk of viral escape through mutations, as they belong to the most conserved regions in the HEV genome. The siRNAs targeting the ORF3 efficiently inhibited viral replication in A549 cells after HEV infection. Importantly, the siRNA was also highly effective at inhibiting HEV in the persistently infected A549 cell line, which provides a suitable model for chronic infection in patients. Furthermore, we showed that a 5' triphosphate modification on the siRNA sense strand activates the RIG-I receptor, a cytoplasmic pattern recognition receptor that recognizes viral RNA. Upon activation, RIG-I triggers a signaling cascade, effectively suppressing HEV replication. This dual-action strategy, combining the activation of the adaptive immune response and the inherent RNAi pathway, inhibits HEV replication successfully and may lead to the development of new therapies.

Genetic Diversity, Selection Signatures, and Genome-Wide Association Study Identify Candidate Genes Related to Litter Size in Hu Sheep.

Hu sheep is a renowned prolific local sheep breed in China, widely distributed across the country due to its excellent reproductive performance. Deciphering the molecular mechanisms underlying the high fecundity of Hu sheep is crucial for improving the litter size of ewes. In this study, we genotyped 830 female Hu sheep using the Illumina OvineSNP50 BeadChip and performed genetic diversity analysis, selection signature detection, and a genome-wide association study (GWAS) for litter size. Our results revealed that the Hu sheep population exhibits relatively high genetic diversity. A total of 4927 runs of homozygosity (ROH) segments were detected, with the majority (74.73%) being short in length. Different genomic inbreeding coefficients (FROH, FHOM, FGRM, and FUNI) ranged from -0.0060 to 0.0126, showing low levels of inbreeding in this population. Additionally, we identified 91 candidate genomic regions through three complementary selection signature methods, including ROH, composite likelihood ratio (CLR), and integrated haplotype score (iHS), and annotated 189 protein-coding genes. Moreover, we observed two significant SNPs related to the litter size of Hu sheep using GWAS analysis based on a repeatability model. Integrating the selection signatures and the GWAS results, we identified 15 candidate genes associated with litter size, among which BMPR1B and UNC5C were particularly noteworthy. These findings provide valuable insights for improving the reproductive performance and breeding of high-fecundity lines of Hu sheep.

Genomic Regions Associated with Growth and Reproduction Traits in Pink-Eyed White Mink.

In mink breeding, balanced selection for growth and reproductive features is essential because these traits are contradictory. The variables of total number born (TNB), number born alive (NBA), and body weight (BW) are highly valuable in terms of their importance in mink production. A comprehensive understanding of the molecular mechanisms that drive these features could offer vital insights into their genetic compositions. In the present study, the single-nucleotide polymorphism (SNP) genotypes of 219 minks were obtained via double digest restriction-site associated DNA sequencing (ddRAD-seq). Following several rounds of screening, about 2,415,121 high-quality SNPs were selected for a genome-wide association study (GWAS). The GWAS was used to determine BW and reproductive traits in pink-eyed white mink. It was suggested that SLC26A36, STXBP5L, and RPS 29 serve as potential genes for the total number of kits born (TNB), while FSCB, PDPN, NKX 2-1, NFKB 1, NFKBIA, and GABBR1 are key genes for the number born alive (NBA). Moreover, RTTN, PRPF31, MACROD1, and KYAT1 are possible BW genes based on association results and available functional data from gene and mammalian phenotype databases. These results offer essential information about the variety of mink and theoretical principles for applying mink breeds.

Peering into the gaps: Long-read sequencing illuminates structural variants and genomic evolution in the Australasian snapper

Even before genome sequencing, genetic resources have supported species management and breeding programs. Current technologies, such as long-read sequencing, resolve complex genomic regions, like those rich in repeats or high in GC content. Improved genome contiguity enhances accuracy in identifying structural variants (SVs) and transposable elements (TEs). We present an improved genome assembly and SV catalogue for the Australasian snapper (Chrysophrys auratus). The new assembly is more contiguous, allowing for putative identification of 14 centromeres and transfer of 26,115 gene annotations from yellowfin seabream. Compared to the previous assembly, 35,000 additional SVs, including larger and more complex rearrangements, were annotated. SVs and TEs exhibit a distribution pattern skewed towards chromosome ends, likely influenced by recombination. Some SVs overlap with growth-related genes, underscoring their significance. This upgraded genome serves as a foundation for studying natural and artificial selection, offers a reference for related species, and sheds light on genome dynamics shaped by evolution.

Human genomic DNA is widely interspersed with i-motif structures.

DNA i-motif structures are formed in the nuclei of human cells and are believed to provide critical genomic regulation. While the existence, abundance, and distribution of i-motif structures in human cells has been demonstrated and studied by immunofluorescent staining, and more recently NMR and CUT&Tag, the abundance and distribution of such structures in human genomic DNA have remained unclear. Here we utilise high-affinity i-motif immunoprecipitation followed by sequencing to map i-motifs in the purified genomic DNA of human MCF7, U2OS and HEK293T cells. Validated by biolayer interferometry and circular dichroism spectroscopy, our approach aimed to identify DNA sequences capable of i-motif formation on a genome-wide scale, revealing that such sequences are widely distributed throughout the human genome and are common in genes upregulated in G0/G1 cell cycle phases. Our findings provide experimental evidence for the widespread formation of i-motif structures in human genomic DNA and a foundational resource for future studies of their genomic, structural, and molecular roles.

Genome wide association study of antioxidant activity in pigeonpea germplasm

Pigeonpea (Cajanus cajan L.) is a nutritious crop with antioxidant potential, making it valuable for health supplements and medicinal products. Pigeonpea seed coat possess a rich source of antioxidant properties, particularly in dark colour seed coat varieties. Recently, antioxidant-based drugs for the prevention and treatment of complex diseases have evoked a great deal of research interest in natural antioxidants. We resolved genomic regions by conducted a genome-wide association study (GWAS) using a high-throughput 62 K SNP “CcSNPnks” genic chip to identify SNP markers associated with antioxidant properties in 287 diverse panel of pigeonpea genotypes. Evaluated their antioxidant activity using a spectrophotometric technique based on ferric ion-reducing antioxidant potential (FRAP) and calculated the gallic acid equivalents (GAE) per 100 g of extract. Significant differences in antioxidant activities were found among genotypes, ranging from 3.27 mg in ICPL-13271 to 16.85 mg in UP 9-I for ferric iron-reducing activity. The GWAS identified four distinct marker trait association (MTA) that exhibit a strong correlation with antioxidant properties. We identified significant marker trait associations with chr 2 (AX-165372487), chr 6 (AX-165348296) code for RGA1 protein, chr 10 (AX-165367685) codes for alkaline ceramidase (transcript variant), and chr 11 (AX-165411944) codes for (DL)-glycerol-3-phosphatase 2. The identified protein, RGA1, is involved in reducing ROS accumulation in response to stress by up-regulating the transcription of superoxide dismutases. It suppresses GA-induced vegetative growth and floral initiation. Rapidly degraded in response to GA and involved in fruit and flower development. Furthermore, the population structure analysis identified four distinct sub-populations with admixture types, effectively grouping all genotypes based on their ancestry. Our findings advance the understanding of the genetic underpinning antioxidant concentrations, holding potential for applications in crop breeding and aligning with consumer preferences. This study provides a cornerstone for future research into exploiting marker-assisted breeding to produce pigeonpea cultivars.

Genome‐wide association analysis of resistance to anthracnose in the Middle American Diversity Panel of common bean (Phaseolus vulgaris L.)

AbstractAnthracnose (ANTH) caused by Colletotrichum lindemuthianum is a major disease of common bean (Phaseolus vulgaris L.). The genetic basis of ANTH resistance in the Middle American Diversity Panel (MDP) is unknown. The objectives of this study were to identify (1) Middle American accessions resistant to races 7, 19, 51, 63, 167, and 1085 of C. lindemuthianum and (ii) genomic regions and positional candidate genes associated with resistance to these races. The MDP composed of 240 Middle American accessions was evaluated for resistance to races 7, 19, 51, 63, 167, and 1085. The MDP was genotyped with 211,763 single nucleotide polymorphisms (SNPs), and mixed linear model analysis was conducted to identify genomic regions associated with resistance to the six races. Seven accessions were highly resistant to all six races, and these can be used as sources of resistance to improve specific market classes in the Middle American gene pool. The genomic region (385,894 bp) on chromosome Pv04 was significantly associated with resistance to race 167. Genomic regions on Pv02 (41,570,325 bp), Pv07 (24,122,343 bp), and Pv11 (51,707,917 bp) were significantly associated with resistance to race 19. Disease resistance (R) genes with the nucleotide binding‐APAF resistance protein and CED‐4 domain were identified as positional candidate genes on Pv04 and Pv11. There were no SNPs significantly associated with resistance to races 7, 51, 63, and 1085. Pyramiding the identified genomic regions on Pv04, Pv07, and Pv11 could provide durable ANTH resistance in Middle American varieties for races 19 and 167.

Association between mitochondrial DNA copy number and production traits in pigs.

Mitochondria are essential organelles in the regulation of cellular energetic metabolism. Mitochondrial DNA copy number (mtDNA_CN) can be used as a proxy for mitochondria number, size, and activity. The aims of our study are to evaluate the effect of mtDNA_CN and mitochondrial haploblocks on production traits in pigs, and to identify the genetic background of this cellular phenotype. We collected performance data of 234 pigs and extracted DNA from skeletal muscle. Whole-genome sequencing data was used to determine mtDNA_CN. We found positive correlations of muscle mtDNA_CN with backfat thickness at 207 d (+0.14; p-value = 0.07) and negative correlations with carcase loin thickness (-0.14; p-value = 0.03). Pigs with mtDNA_CN values below the lower quartile had greater loin thickness (+4.1 mm; p-value = 0.01) and lower backfat thickness (-1.1 mm; p-value = 0.08), which resulted in greater carcase lean percentage (+2.4%; p-value = 0.04), than pigs with mtDNA_CN values above the upper quartile. These results support the hypothesis that a reduction of mitochondrial activity is associated with greater feed efficiency. Higher mtDNA_CN was also positively correlated with higher meat ultimate pH (+0.19; p-value <0.01) but we did not observe significant difference for meat ultimate pH between the two groups with extreme mtDNA_CN. We found no association of the most frequent mitochondrial haploblocks with mtDNA_CN or the production traits, but several genomic regions that harbour potential candidate genes with functions related to mitochondrial biogenesis and homeostasis were associated with mtDNA_CN. These regions provide new insights into the genetic background of this cellular phenotype but it is still uncertain if such associations translate into noticeable effects on the production traits.

Rare germline genetic variation in PAX8 transcription factor binding sites and susceptibility to epithelial ovarian cancer.

Common genetic variation throughout the genome together with rare coding variants identified to date explain about a half of the inherited genetic component of epithelial ovarian cancer risk. It is likely that rare variation in the non-coding genome will explain some of the unexplained heritability, but identifying such variants is challenging. The primary problem is lack of statistical power to identifying individual risk variants by association as power is a function of sample size, effect size and allele frequency. Power can be increased by using burden tests which test for association of carriers of any variant in a specified genomic region. This has the effect of increasing the putative effect allele frequency. PAX8 is a transcription factor that plays a critical role in tumour progression, migration and invasion. Furthermore, regulatory elements proximal to target genes of PAX8 are enriched for common ovarian cancer risk variants. We hypothesised that rare variation in PAX8 binding sites are also associated with ovarian cancer risk, but unlikely to be associated with risk of breast, colorectal or endometrial cancer. We have used publicly available, whole-genome sequencing data from the UK 100,000 Genomes Project to evaluate the burden of rare variation in PAX8 binding sites across the genome. Data were available for 522 ovarian cancers, 2,984 breast cancers, 2,696 colorectal cancers, 836 endometrial cancers and 2253 non-cancer controls. Active binding sites were defined using data from multiple PAX8 and H3K27 ChIPseq experiments. We found no association between the burden of rare variation in PAX8 binding sites (defined in several ways) and risk of ovarian, breast or endometrial cancer. An apparent association with colorectal cancer was likely to be a technical artefact as a similar association was also detected for rare variation in random regions of the genome. Despite the null result this study provides a proof-of -principle for using burden testing to identify rare, non-coding germline genetic variation associated with disease. Larger sample sizes available from large-scale sequencing projects together with improved understanding of the function of the non-coding genome will increase the potential of similar studies in the future.

Multiomics of a rice population identifies genes and genomic regions that bestow low glycemic index and high protein content

To counter the rising incidence of diabetes and to meet the daily protein needs, we created low glycemic index (GI) rice varieties with protein content (PC) surpassing 14%. In the development of recombinant inbred lines using Samba Mahsuri and IR36 amylose extender (IR36ae) as parental lines, we identified quantitative trait loci and genes associated with low GI, high amylose content (AC), and high PC. By integrating genetic techniques with classification models, this comprehensive approach identified candidate genes on chromosome 2 (qGI2.1/qAC2.1 spanning the region from 18.62 Mb to 19.95 Mb), exerting influence on low GI and high amylose. Notably, the phenotypic variant with high value was associated with the recessive allele of the starch branching enzyme 2b (sbeIIb). The genome-edited sbeIIb line confirmed low GI phenotype in milled rice grains. Further, combinations of alleles created by the highly significant SNPs from the targeted associations and epistatically interacting genes showed ultralow GI phenotypes with high amylose and high protein. Metabolomics analysis of rice with varying AC, PC, and GI revealed that the superior lines of high AC and PC, and low GI were preferentially enriched in glycolytic and amino acid metabolisms, whereas the inferior lines of low AC and PC and high GI were enriched with fatty acid metabolism. The high amylose high protein recombinant inbred line (HAHP_101) was enriched in essential amino acids like lysine. Such lines may be highly relevant for food product development to address diabetes and malnutrition.

Chlomito: a novel tool for precise elimination of organelle genome contamination from nuclear genome assembly.

Accurate reference genomes are fundamental to understanding biological evolution, biodiversity, hereditary phenomena and diseases. However, many assembled nuclear chromosomes are often contaminated by organelle genomes, which will mislead bioinformatic analysis, and genomic and transcriptomic data interpretation. To address this issue, we developed a tool named Chlomito, aiming at precise identification and elimination of organelle genome contamination from nuclear genome assembly. Compared to conventional approaches, Chlomito utilized new metrics, alignment length coverage ratio (ALCR) and sequencing depth ratio (SDR), thereby effectively distinguishing true organelle genome sequences from those transferred into nuclear genomes via horizontal gene transfer (HGT). The accuracy of Chlomito was tested using sequencing data from Plum, Mango and Arabidopsis. The results confirmed that Chlomito can accurately detect contigs originating from the organelle genomes, and the identified contigs covered most regions of the organelle reference genomes, demonstrating efficiency and precision of Chlomito. Considering user convenience, we further packaged this method into a Docker image, simplified the data processing workflow. Overall, Chlomito provides an efficient, accurate and convenient method for identifying and removing contigs derived from organelle genomes in genomic assembly data, contributing to the improvement of genome assembly quality.

Contribution of homoeologous exchange to domestication of polyploid Brassica

BackgroundPolyploidy is widely recognized as a significant evolutionary force in the plant kingdom, contributing to the diversification of plants. One of the notable features of allopolyploidy is the occurrence of homoeologous exchange (HE) events between the subgenomes, causing changes in genomic composition, gene expression, and phenotypic variations. However, the role of HE in plant adaptation and domestication remains unclear.ResultsHere we analyze the whole-genome resequencing data from Brassica napus accessions representing the different morphotypes and ecotypes, to investigate the role of HE in domestication. Our findings demonstrate frequent occurrence of HEs in Brassica napus, with substantial HE patterns shared across populations, indicating their potential role in promoting crop domestication. HE events are asymmetric, with the A genome more frequently replacing C genome segments. These events show a preference for specific genomic regions and vary among populations. We also identify candidate genes in HE regions specific to certain populations, which likely contribute to flowering-time diversification across diverse morphotypes and ecotypes. In addition, we assemble a new genome of a swede accession, confirming the HE signals on the genome and their potential involvement in root tuber development. By analyzing HE in another allopolyploid species, Brassica juncea, we characterize a potential broader role of HE in allopolyploid crop domestication.ConclusionsOur results provide novel insights into the domestication of polyploid Brassica species and highlight homoeologous exchange as a crucial mechanism for generating variations that are selected for crop improvement in polyploid species.