Genomic Structural Variation Research Articles

Modular networks and genomic variation during progression from stable angina pectoris through ischemic cardiomyopathy to chronic heart failure

BackgroundAnalyzing disease–disease relationships plays an important role for understanding etiology, disease classification, and drug repositioning. However, as cardiovascular diseases with causative links, the molecular relationship among stable angina pectoris (SAP), ischemic cardiomyopathy (ICM) and chronic heart failure (CHF) is not clear.MethodsIn this study, by integrating the multi-database data, we constructed paired disease progression modules (PDPMs) to identified relationship among SAP, ICM and CHF based on module reconstruction pairs (MRPs) of K-value calculation (a Euclidean distance optimization by integrating module topology parameters and their weights) methods. Finally, enrichment analysis, literature validation and structural variation (SV) were performed to verify the relationship between the three diseases in PDPMs.ResultsTotal 16 PDPMs were found with K > 0.3777 among SAP, ICM and CHF, in which 6 pairs in SAP–ICM, 5 pairs for both ICM–CHF and SAP–CHF. SAP–ICM was the most closely related by having the smallest average K-value (K = 0.3899) while the maximum is SAP–CHF (K = 0.4006). According to the function of the validation gene, inflammatory response were through each stage of SAP–ICM–CHF, while SAP–ICM was uniquely involved in fibrosis, and genes were related in affecting the upstream of PI3K–Akt signaling pathway. 4 of the 11 genes (FLT1, KDR, ANGPT2 and PGF) in SAP–ICM–CHF related to angiogenesis in HIF-1 signaling pathway. Furthermore, we identified 62.96% SVs were protein deletion in SAP–ICM–CHF, and 53.85% SVs were defined as protein replication in SAP–ICM, while ICM–CHF genes were mainly affected by protein deletion.ConclusionThe PDPMs analysis approach combined with genomic structural variation provides a new avenue for determining target associations contributing to disease progression and reveals that inflammation and angiogenesis may be important links among SAP, ICM and CHF progression.

Pan-genome and multi-parental framework for high-resolution trait dissection in melon (Cucumis melo).

Linking genotype with phenotype is a fundamental goal in biology and requires robust data for both. Recent advances in plant-genome sequencing have expedited comparisons among multiple-related individuals. The abundance of structural genomic within-species variation that has been discovered indicates that a single reference genome cannot represent the complete sequence diversity of a species, leading to the expansion of the pan-genome concept. For high-resolution forward genetics, this unprecedented access to genomic variation should be paralleled and integrated with phenotypic characterization of genetic diversity. We developed a multi-parental framework for trait dissection in melon (Cucumis melo), leveraging a novel pan-genome constructed for this highly variable cucurbit crop. A core subset of 25 diverse founders (MelonCore25), consisting of 24 accessions from the two widely cultivated subspecies of C. melo, encompassing 12 horticultural groups, and 1 feral accession was sequenced using a combination of short- and long-read technologies, and their genomes were assembled de novo. The construction of this melon pan-genome exposed substantial variation in genome size and structure, including detection of ~300 000 structural variants and ~9 million SNPs. A half-diallel derived set of 300 F2 populations, representing all possible MelonCore25 parental combinations, was constructed as a framework for trait dissection through integration with the pan-genome. We demonstrate the potential of this unified framework for genetic analysis of various melon traits, including rind color intensity and pattern, fruit sugar content, and resistance to fungal diseases. We anticipate that utilization of this integrated resource will enhance genetic dissection of important traits and accelerate melon breeding.

Copy Number Variation among Resistance Genes Analogues in Brassica napus.

Copy number variations (CNVs) are defined as deletions, duplications and insertions among individuals of a species. There is growing evidence that CNV is a major factor underlining various autoimmune disorders and diseases in humans; however, in plants, especially oilseed crops, the role of CNVs in disease resistance is not well studied. Here, we investigate the genome-wide diversity and genetic properties of CNVs in resistance gene analogues (RGAs) across eight Brassica napus lines. A total of 1137 CNV events (704 deletions and 433 duplications) were detected across 563 RGAs. The results show CNVs are more likely to occur across clustered RGAs compared to singletons. In addition, 112 RGAs were linked to a blackleg resistance QTL, of which 25 were affected by CNV. Overall, we show that the presence and abundance of CNVs differ between lines, suggesting that in B. napus, the distribution of CNVs depends on genetic background. Our findings advance the understanding of CNV as an important type of genomic structural variation in B. napus and provide a resource to support breeding of advanced canola lines.

26. Identification of novel genomic structural variations in Angiosarcoma by Optical Genome Mapping

Angiosarcomas are malignant vascular tumors with a number of known molecular aberrations identified by standard technologies, such as FISH and next generation sequencing. However, potentially clinically significant large genomic aberrations, such as gene fusion, large amplification and deletions, may not be detectable by these technologies. To gain further insights into the genomic landscape of angiosarcoma, we performed optical genome mapping (OGM) of 3 primary human angiosarcoma cell lines. OGM analysis of these cells revealed multiple unique structural aberrations, including large insertions/deletions, inversions, duplications, inter-chromosomal and intra-chromosomal translocations. Further analysis showed a 354-Kb deletion in Plexin D1 gene locus that regulates the migration of endothelial cells and is required for normal development of the heart and vasculature; reciprocal t(4;9) translocation in UGT8 (UDP Glycosyltransferase 8) gene locus that is involved in biosynthesis of galactocerebrosidase; reciprocal t(6;22) translocation involving HIRA (Histone Cell Cycle Regulator) and EYS (Eyes Shut Homolog) gene loci; and 14 Mb inversion interrupting RYR2 (Ryanodine Receptor 2) and TAF1A (TATA-Box Binding Protein Associated Factor) gene loci. These findings demonstrate the utility of optical genome mapping as a new tool to uncover large genomic structural variations to better understand the molecular pathogenesis of malignant vascular tumors. Future studies are aimed to determine the biological significance of these structural variations in angiosarcoma.

Time-ordering japonica/geng genomes analysis indicates the importance of large structural variants in rice breeding.

Temperate japonica/geng (GJ) rice yield has significantly improved due to intensive breeding efforts, dramatically enhancing global food security. However, little is known about the underlying genomic structural variations (SVs) responsible for this improvement. We compared 58 long-read assemblies comprising cultivated and wild rice species in the present study, revealing 156 319 SVs. The phylogenomic analysis based on the SV dataset detected the putatively selected region of GJ sub-populations. A significant portion of the detected SVs overlapped with genic regions were found to influence the expression of involved genes inside GJ assemblies. Integrating the SVs and causal genetic variants underlying agronomic traits into the analysis enables the precise identification of breeding signatures resulting from complex breeding histories aimed at stress tolerance, yield potential and quality improvement. Further, the results demonstrated genomic and genetic evidence that the SV in the promoter of LTG1 is accounting for chilling sensitivity, and the increased copy numbers of GNP1 were associated with positive effects on grain number. In summary, the current study provides genomic resources for retracing the properties of SVs-shaped agronomic traits during previous breeding procedures, which will assist future genetic, genomic and breeding research on rice.

Optical genome mapping identifies clinically relevant genomic rearrangements in prostate cancer biopsy sample

BackgroundProstate cancer (PCa) is characterized by complex genomic rearrangements such as the ETS oncogene family fusions, yet the clinical relevance is not well established. While paneled genetic tests of DNA repair genes are recommended in advanced PCa, conventional genomic or cytogenetic tools are not ideal for genome-wide screening of structural variations (SVs) such as balanced translocation due to cost and/or resolution issues.MethodsIn this study, we tested the feasibility of whole-genome optical genomic mapping (OGM), a newly developed platform for genome-wide SV analysis to detect complex genomic rearrangements in consecutive unselected PCa samples from MRI/US-fusion targeted biopsy.ResultsWe tested ten samples, and nine (90%) passed quality check. Average mapping rate and coverage depth were 58.1 ± 23.7% and 157.3 ± 97.7×, respectively (mean ± SD). OGM detected copy number alterations such as chr6q13 loss and chr8q12-24 gain. Two adjacent tumor samples were distinguished by inter/intra-chromosomal translocations, revealing that they’re from the same ancestor. Furthermore, OGM detected large deletion of chr13q13.1 accompanied by inter-chromosomal translocation t(13;20)(q13.1;p13) occurring within BRCA2 gene, suggesting complete loss of function.ConclusionIn conclusion, clinically relevant genomic SVs were successfully detected in PCa samples by OGM. We suggest that OGM can complement panel sequencing of DNA repair genes BRCA1/2 or ATM in high-risk PCa.

Influences of rare copy-number variation on human complex traits

The human genome contains hundreds of thousands of regions harboring copy-number variants (CNV). However, the phenotypic effects of most such polymorphisms are unknown because only larger CNVs have been ascertainable from SNP-array data generated by large biobanks. We developed a computational approach leveraging haplotype sharing in biobank cohorts to more sensitively detect CNVs. Applied to UKBiobank, this approach accounted for approximately half of all rare gene inactivation events produced by genomic structural variation. This CNV call set enabled a detailed analysis of associations between CNVs and 56 quantitative traits, identifying 269 independent associations (p<5×10-8) likely to be causally driven by CNVs. Putative target genes were identifiable for nearly half of the loci, enabling insights intodosage sensitivity of these genes and uncovering several gene-trait relationships. These results demonstrate the ability of haplotype-informed analysis to provide insights into the genetic basis of human complex traits.

In-depth analysis of genomes and functional genomics of orchid using cutting-edge high-throughput sequencing.

High-throughput sequencing technology has been facilitated the development of new methodologies and approaches for studying the origin and evolution of plant genomes and subgenomes, population domestication, and functional genomics. Orchids have tens of thousands of members in nature. Many of them have promising application potential in the extension and conservation of the ecological chain, the horticultural use of ornamental blossoms, and the utilization of botanical medicines. However, a large-scale gene knockout mutant library and a sophisticated genetic transformation system are still lacking in the improvement of orchid germplasm resources. New gene editing tools, such as the favored CRISPR-Cas9 or some base editors, have not yet been widely applied in orchids. In addition to a large variety of orchid cultivars, the high-precision, high-throughput genome sequencing technology is also required for the mining of trait-related functional genes. Nowadays, the focus of orchid genomics research has been directed to the origin and classification of species, genome evolution and deletion, gene duplication and chromosomal polyploidy, and flower morphogenesis-related regulation. Here, the progressing achieved in orchid molecular biology and genomics over the past few decades have been discussed, including the evolution of genome size and polyploidization. The frequent incorporation of LTR retrotransposons play important role in the expansion and structural variation of the orchid genome. The large-scale gene duplication event of the nuclear genome generated plenty of recently tandem duplicated genes, which drove the evolution and functional divergency of new genes. The evolution and loss of the plastid genome, which mostly affected genes related to photosynthesis and autotrophy, demonstrated that orchids have experienced more separate transitions to heterotrophy than any other terrestrial plant. Moreover, large-scale resequencing provide useful SNP markers for constructing genetic maps, which will facilitate the breeding of novel orchid varieties. The significance of high-throughput sequencing and gene editing technologies in the identification and molecular breeding of the trait-related genes in orchids provides us with a representative trait-improving gene as well as some mechanisms worthy of further investigation. In addition, gene editing has promise for the improvement of orchid genetic transformation and the investigation of gene function. This knowledge may provide a scientific reference and theoretical basis for orchid genome studies.

Two independent allohexaploidizations and genomic fractionation in Solanales

Solanales, an order of flowering plants, contains the most economically important vegetables among all plant orders. To date, many Solanales genomes have been sequenced. However, the evolutionary processes of polyploidization events in Solanales and the impact of polyploidy on species diversity remain poorly understood. We compared two representative Solanales genomes (Solanum lycopersicum L. and Ipomoea triloba L.) and the Vitis vinifera L. genome and confirmed two independent polyploidization events. Solanaceae common hexaploidization (SCH) and Convolvulaceae common hexaploidization (CCH) occurred ∼43–49 and ∼40–46 million years ago (Mya), respectively. Moreover, we identified homologous genes related to polyploidization and speciation and constructed multiple genomic alignments with V. vinifera genome, providing a genomic homology framework for future Solanales research. Notably, the three polyploidization-produced subgenomes in both S. lycopersicum and I. triloba showed significant genomic fractionation bias, suggesting the allohexaploid nature of the SCH and CCH events. However, we found that the higher genomic fractionation bias of polyploidization-produced subgenomes in Solanaceae was likely responsible for their more abundant species diversity than that in Convolvulaceae. Furthermore, through genomic fractionation and chromosomal structural variation comparisons, we revealed the allohexaploid natures of SCH and CCH, both of which were formed by two-step duplications. In addition, we found that the second step of two paleohexaploidization events promoted the expansion and diversity of β-amylase (BMY) genes in Solanales. These current efforts provide a solid foundation for future genomic and functional exploration of Solanales.

Evolutionary divergence of duplicated genomes in newly described allotetraploid cottons

Allotetraploid cotton (Gossypium) species represents a model system for the study of plant polyploidy, molecular evolution, and domestication. Here, chromosome-scale genome sequences were obtained and assembled for two recently described wild species of tetraploid cotton, Gossypium ekmanianum [(AD)6, Ge] and Gossypium stephensii [(AD)7, Gs], and one early form of domesticated Gossypium hirsutum, race punctatum [(AD)1, Ghp]. Based on phylogenomic analysis, we provide a dated whole-genome level perspective for the evolution of the tetraploid Gossypium clade and resolved the evolutionary relationships of Gs, Ge, and domesticated G. hirsutum. We describe genomic structural variation that arose during Gossypium evolution and describe its correlates-including phenotypic differentiation, genetic isolation, and genetic convergence-that contributed to cotton biodiversity and cotton domestication. Presence/absence variation is prominent in causing cotton genomic structural variations. A presence/absence variation-derived gene encoding a phosphopeptide-binding protein is implicated in increasing fiber length during cotton domestication. The relatively unimproved Ghp offers the potential for gene discovery related to adaptation to environmental challenges. Expanded gene families enoyl-CoA δ isomerase 3 and RAP2-7 may have contributed to abiotic stress tolerance, possibly by targeting plant hormone-associated biochemical pathways. Our results generate a genomic context for a better understanding of cotton evolution and for agriculture.

Direct Comparison of HPV16 Viral Genomic Integration, Copy Loss, and Structural Variants in Oropharyngeal and Uterine Cervical Cancers Reveal Distinct Relationships to E2 Disruption and Somatic Alteration.

Simple SummaryHPV16 causes approximately 60% of uterine cervical cancer and 95% of HPV-driven oropharynx cancers. Despite both being HPV-associated, these tumors are very different. We directly compare integration of the HPV16 genome into the human genome in these two diseases, finding that the viral gene E2 is frequently lost in cervical cancer, but usually maintained in oropharyngeal cancer. We also found that oropharyngeal cancers with integration have many more integration sites per tumor and these more frequently occur in genomic regions with a high density of genes.Squamous cell carcinoma of the oropharynx caused by HPV type 16 (HPV16+ OPSCC) is the most common HPV-associated malignancy in the USA and has many molecular differences from uterine cervical squamous cell carcinoma (UCSCC). Our understanding of HPV oncogenesis relied on studies of UCSCC revealing a consensus model reliant on HPV integration with a loss of E2. Here, we compare patterns of HPV integration in UCSCC and OPSCC by analysis of affinity capture sequencing of the HPV16 genome in 104 OPSCC and 44 UCSCC tumors. These cohorts were contemporaneously sequenced using an identical strategy. Integration was identified using discordant read pair clustering and assembly-based approaches. Viral integration sites, structural variants, and copy losses were examined. While large-scale deep losses of HPV16 genes were common in UCSCC and were associated with E2 loss, deep copy losses of the HPV16 genome were infrequent in HPV16+ OPSCC. Similarly, structural variants within HPV16 favored E2 loss in UCSCC but not OPSCC. HPV16 integration sites were non-random, with recurrent integration hot-spots identified. OPSCC tumors had many more integration sites per tumor when compared to UCSCC and had more integration sites in genomic regions with high gene density. These data show that viral integration and E2 disruption are distinct in UCSCC and OPSCC. Our findings also add to growing literature suggesting that HPV tumorigenesis in OPSCC does not follow the model developed based on UCSCC.

Genome-wide detection of copy number variation in American mink using whole-genome sequencing

BackgroundCopy number variations (CNVs) represent a major source of genetic diversity and contribute to the phenotypic variation of economically important traits in livestock species. In this study, we report the first genome-wide CNV analysis of American mink using whole-genome sequence data from 100 individuals. The analyses were performed by three complementary software programs including CNVpytor, DELLY and Manta.ResultsA total of 164,733 CNVs (144,517 deletions and 20,216 duplications) were identified representing 5378 CNV regions (CNVR) after merging overlapping CNVs, covering 47.3 Mb (1.9%) of the mink autosomal genome. Gene Ontology and KEGG pathway enrichment analyses of 1391 genes that overlapped CNVR revealed potential role of CNVs in a wide range of biological, molecular and cellular functions, e.g., pathways related to growth (regulation of actin cytoskeleton, and cAMP signaling pathways), behavior (axon guidance, circadian entrainment, and glutamatergic synapse), lipid metabolism (phospholipid binding, sphingolipid metabolism and regulation of lipolysis in adipocytes), and immune response (Wnt signaling, Fc receptor signaling, and GTPase regulator activity pathways). Furthermore, several CNVR-harbored genes associated with fur characteristics and development (MYO5A, RAB27B, FGF12, SLC7A11, EXOC2), and immune system processes (SWAP70, FYN, ORAI1, TRPM2, and FOXO3).ConclusionsThis study presents the first genome-wide CNV map of American mink. We identified 5378 CNVR in the mink genome and investigated genes that overlapped with CNVR. The results suggest potential links with mink behaviour as well as their possible impact on fur quality and immune response. Overall, the results provide new resources for mink genome analysis, serving as a guideline for future investigations in which genomic structural variations are present.

Genome-wide analysis of CNVs in three populations of Tibetan sheep using whole-genome resequencing.

Copy number variation (CNV), an important source of genomic structural variation, can disturb genetic structure, dosage, regulation and expression, and is associated with phenotypic diversity and adaptation to local environments in mammals. In the present study, 24 resequencing datasets were used to characterize CNVs in three ecotypic populations of Tibetan sheep and assess CNVs related to domestication and adaptation in Qinghai-Tibetan Plateau. A total of 87,832 CNV events accounting for 0.3% of the sheep genome were detected. After merging the overlapping CNVs, 2777 CNV regions (CNVRs) were obtained, among which 1098 CNVRs were shared by the three populations. The average length of these CNVRs was more than 3 kb, and duplication events were more frequent than deletions. Functional analysis showed that the shared CNVRs were significantly enriched in 56 GO terms and 18 KEGG pathways that were mainly concerned with ABC transporters, olfactory transduction and oxygen transport. Moreover, 188 CNVRs overlapped with 97 quantitative trait loci (QTLs), such as growth and carcass QTLs, immunoglobulin QTLs, milk yield QTLs and fecal egg counts QTLs. PCDH15, APP and GRID2 overlapped with body weight QTLs. Furthermore, Vst analysis showed that RUNX1, LOC101104348, LOC105604082 and PAG11 were highly divergent between Highland-type Tibetan Sheep (HTS) and Valley-type Tibetan sheep (VTS), and RUNX1 and LOC101111988 were significantly differentiated between VTS and Oura-type Tibetan sheep (OTS). The duplication of RUNX1 may facilitate the hypoxia adaptation of OTS and HTS in Qinghai-Tibetan Plateau, which deserves further research in detail. In conclusion, for the first time, we represented the genome-wide distribution characteristics of CNVs in Tibetan sheep by resequencing, and provided a valuable genetic variation resource, which will facilitate the elucidation of the genetic basis underlying the distinct phenotypic traits and local adaptation of Tibetan sheep.

Copy Number Variants in Two Northernmost Cattle Breeds Are Related to Their Adaptive Phenotypes.

Copy number variations (CNVs) are genomic structural variants with potential functional and evolutionary effects on phenotypes. In this study, we report the identification and characterization of CNVs from the whole-genome resequencing data of two northernmost cattle breeds from Russia: the Yakut and Kholmogory cattle and their phylogenetically most related breeds, Hanwoo and Holstein, respectively. Comparisons of the CNV regions (CNVRs) among the breeds led to the identification of breed-specific CNVRs shared by cold-adapted Kholmogory and Yakut cattle. An investigation of their overlap with genes, regulatory domains, conserved non-coding elements (CNEs), enhancers, and quantitative trait loci (QTLs) was performed to further explore breed-specific biology and adaptations. We found CNVRs enriched for gene ontology terms related to adaptation to environments in both the Kholmogory and Yakut breeds and related to thermoregulation specifically in Yakut cattle. Interestingly, the latter has also been supported when exploring the enrichment of breed-specific CNVRs in the regulatory domains and enhancers, CNEs, and QTLs implying the potential contribution of CNVR to the Yakut and Kholmogory cattle breeds’ adaptation to a harsh environment.

Exploring Genes and Phenotypes Within Chromosomal Regions Using OMIM's GeneScout.

Structural variation in genomes, such as copy number variants (CNVs), is under scrutiny for its contribution to phenotypic expression and evolution. Regions of homozygosity (ROH) are ripe for phenotype-gene discovery. Determining the genes and related phenotypes within genomic regions is key to studying potential functional and phenotypic consequences. Because individuals have multiple CNVs and ROHs in their genome, identifying genomic regions that are phenotypically significant is challenging. GeneScout is a web-based tool that can be used to search genomic regions to display and filter the genes and their associated phenotypes within regions of interest. Phenotypes and their associated gene(s) can then be filtered to show only the genes with phenotypes that have a particular inheritance pattern and/or specific clinical feature(s). Phenotypes can then be selected to compare the clinical synopses side-by-side in Online Mendelian Inheritance in Man (OMIM® ).Additionally, two coordinate sets can be compared to determine either the regions of overlap or the unique regions (subtraction). The resulting coordinate ranges are displayed on the results page, and the results table displays only the genes and phenotypes present within the coordinate ranges. The interactive table includes gene-specific links to external resources such as ClinVar, ClinGen validity, ClinGen dosage, and gnomAD, and a diamond symbol next to the gene name indicates a gene that spans the start or end of a coordinate range. Searches and comparisons may be performed for coordinates in assemblies GRCh37 (hg19) and GRCh38 (hg38). The results page offers the option to liftover coordinates entered in GRCh37 to GRCh38 and updates the results table to display the gene content based on assembly GRCh38. The search coordinates and results table can be downloaded in a tab-delimited or Excel file. © 2022 Wiley Periodicals LLC. Basic Protocol: Searching GeneScout.

Large Fragment InDels Reshape Genome Structure of Porcine Alveolar Macrophage 3D4/21 Cells.

The porcine monomyeloid cell line, or 3D4/21 cells, is an effective tool to study the immune characteristics and virus infection mechanism of pigs. Due to the introduction of the neomycin resistance gene and the SV40 large T antigen gene, its genome has undergone essential changes, which are still unknown. Studying the variation in genome structure, especially the large fragments of insertions and deletions (InDels), is one of the proper ways to reveal these issues. In this study, an All-seq method was established by combining Mate-pair and Shotgun sequencing methods, and the detection and verification of large fragments of InDels were performed on 3D4/21 cells. The results showed that there were 844 InDels with a length of more than 1 kb, of which 12 regions were deletions of more than 100 kb in the 3D4/21 cell genome. In addition, compared with porcine primary alveolar macrophages, 82 genes including the CD163 had lost transcription in 3D4/21 cells, and 72 genes gained transcription as well. Further referring to the Hi-C structure, it was found that the fusion of the topologically associated domains (TADs) caused by the deletion may lead to abnormal gene function. The results of this study provide a basis for elaborating the genome structure and functional variation in 3D4/21 cells, provide a method for rapid and convenient detection of large-scale InDels, and provide useful clues for the study of the porcine immune function genome and the molecular mechanism of virus infection.

EVOLUTIONARY PATTERNS OF GENOME SIZE AND CHROMOSOME NUMBER VARIATION IN BEGONIACEAE

Cytological data resources are crucial to the study and understanding of the evolution of complex taxa. Recent research on the genus Begonia L. has provided a robust phylogenetic background for the analysis of evolutionary patterns across the group and has established that Begonia is variable in genome size and chromosome number. This paper provides an overview of the genome structural variation present in Begonia and an updated chromosome number and genome size dataset for the genus. Chromosome numbers of more than 400 species are presented and discussed within their current taxonomic and phylogenetic context. A more complete chromosome number dataset is available for Neotropical and Asian Begonia sections than for those from Africa. The distribution of chromosome numbers across phylogenetic trees supports the idea of Begonia sections as natural groups, because most variation is found between sections rather than within them. Some larger Begonia clades were found to have larger chromosome number variation. Moreover, groups with the most variable chromosome numbers belong to some of the taxonomically complex or unresolved Begonia sections. Genome size variation was found not to correlate with changes in chromosome number. It suggests that Begonia genome dynamics are caused not only by large-scale duplications, rearrangements, and changes in ploidy levels but also by changes in the repetitive fraction of the genome, which probably cause changes in chromosome size. This could potentially play an important role in species radiations.

Comparative Genomic and Phylogenetic Analysis of Forty Gentiana Chloroplast Genomes.

Gentiana plants, which have great medicinal and ornamental value, are widely distributed in diverse habitats and have complex taxonomy. Here 40 Gentiana chloroplast genomes were used for comparative genomic analysis and divergence time estimation. The complete chloroplast genome of G. rhodantha was sequenced, assembled, and annotated. Comparative genomic and phylogenetic analysis were provided for variation analysis of Gentiana. Gentiana species satisfy the characteristics of intra-Sect conservation and inter-Sect variation in chloroplast genome structure and IR boundaries. All Gentiana Sects can be clustered into a single one and separated from each other; however, Ser. Apteroideae and Ser. Confertifoliae in Sect. Monopodiae are more closely related to Sect. Frigida and Sect. Cruciata, respectively. Gentiana has experienced two large gene loss events; the first, the collective loss of the rps16 gene at genus formation and the second, the collective loss of the ndh gene when Ser. Ornatae and Ser. Verticillatae completed their differentiation. Comparative genomic analysis support that Sect. Stenogyne and Sect. Otophora became the independent genera Metagentiana and Kuepferia. Seven divergence hotspot regions were screened based on Pi values, and could serve as DNA-specific barcodes for Gentiana. This study provides a further theoretical basis for taxonomic analysis, genetic diversity, evolutionary mechanism and molecular identification in Gentiana.

Evolutionary insights into porcine genomic structural variations based on a novel‐constructed dataset from 24 worldwide diverse populations

AbstractStructural variations (SVs) are important DNA polymorphisms that contribute to genetic diversity and evolution in humans, animals, and plants. In this study, we present a novel swine SV dataset of 79,919 deletions, 23,638 duplications, and 9333 inversions with average sequence depths of 24.1× from 24 varieties of worldwide pig populations, encompassing 305 individuals. Genotypes of SVs, particularly deletions, can accurately group individuals based on their population identity. We showed that exon‐covering deletions were subject to negative selection. Fixation index and differential allele frequency analysis identified highly differentiated SVs between European and Asian indigenous pigs, including deletions in NR6A1 and PLAG1, which are significantly associated with vertebrate numbers and growth performances, respectively. The growth‐enhancing allele at the deletion in PLAG1 was shared by European commercial pigs and Northern Chinese indigenous pigs including Laiwu and Min pigs, suggesting potential introgression from European commercial breeds into Chinese indigenous breeds. Moreover, we uncovered highly differentiated SVs in 139 genes between domesticated pigs and wild boars in Asia, temperature and altitude adaptation‐associated SVs in 41 genes, and population‐specific SVs in 718 genes. This study provides novel insights into the role of porcine SVs in domestication, environmental adaptation, and breed formation.

Influence of genomic structural variations and nutritional conditions on the emergence of quorum sensing-dependent gene regulation defects in Burkholderia glumae.

Bacteria often change their genetic and physiological traits to survive in harsh environments. To determine whether, in various strains of Burkholderia glumae, genomic diversity is associated with the ability to adapt to ever-changing environments, whole genomes of 44 isolates from different hosts and regions were analyzed. Whole-genome phylogenetic analysis of the 44 isolates revealed six clusters and two divisions. While all isolates possessed chromosomes 1 and 2, strains BGR80S and BGR81S had one chromosome resulting from the merging of the two chromosomes. Upon comparison of genomic structures to the prototype BGR1, inversions, deletions, and rearrangements were found within or between chromosomes 1 and/or 2 in the other isolates. When three isolates—BGR80S, BGR15S, and BGR21S, representing clusters III, IV, and VI, respectively—were grown in Luria-Bertani medium, spontaneous null mutations were identified in qsmR encoding a quorum-sensing master regulator. Six days after subculture, qsmR mutants were found at detectable frequencies in BGR15S and BGR21S, and reached approximately 40% at 8 days after subculture. However, the qsmR mutants appeared 2 days after subculture in BGR80S and dominated the population, reaching almost 80%. No qsmR mutant was detected at detectable frequency in BGR1 or BGR13S. The spontaneous qsmR mutants outcompeted their parental strains in the co-culture. Daily addition of glucose or casamino acids to the batch cultures of BGR80S delayed emergence of qsmR mutants and significantly reduced their incidence. These results indicate that spontaneous qsmR mutations are correlated with genomic structures and nutritional conditions.