Minimally destructive hDNA extraction method for retrospective genetics of pinned historical Lepidoptera specimens.

The millions of specimens stored in entomological collections provide a unique opportunity to study historical insect diversity. Current technologies allow to sequence entire genomes of historical specimens and estimate past genetic diversity of present-day endangered species, advancing our understanding of anthropogenic impact on genetic diversity and enabling the implementation of conservation strategies. A limiting challenge is the extraction of historical DNA (hDNA) of adequate quality for sequencing platforms. We tested four hDNA extraction protocols on five body parts of pinned false heath fritillary butterflies, Melitaea diamina, aiming to minimise specimen damage, preserve their scientific value to the collections, and maximise DNA quality and yield for whole-genome re-sequencing. We developed a very effective approach that successfully recovers hDNA appropriate for short-read sequencing from a single leg of pinned specimens using silica-based DNA extraction columns and an extraction buffer that includes SDS, Tris, Proteinase K, EDTA, NaCl, PTB, and DTT. We observed substantial variation in the ratio of nuclear to mitochondrial DNA in extractions from different tissues, indicating that optimal tissue choice depends on project aims and anticipated downstream analyses. We found that sufficient DNA for whole genome re-sequencing can reliably be extracted from a single leg, opening the possibility to monitor changes in genetic diversity maintaining the scientific value of specimens while supporting current and future conservation strategies.

Genome-wide analysis reveals genomic diversity and signatures of selection in Qinchuan beef cattle.

BACKGROUND: Indigenous Chinese cattle have abundant genetic diversity and a long history of artificial selection, giving local breeds advantages in adaptability, forage tolerance and resistance. The detection of selective sweeps and comparative genome analysis of selected breeds and ancestral populations provide a basis for understanding differences among breeds and for the identification and utilization of candidate genes. We investigated genetic diversity, population structure, and signatures of selection using genome-wide sequencing data for a new breed of Qinchuan cattle (QNC, n = 21), ancestral Qinchuan cattle (QCC, n = 20), and Zaosheng cattle (ZSC, n = 19). RESULTS: A population structure analysis showed that the ancestry components of QNC and ZSC were similar. In addition, the QNC and ZSC groups showed higher proportions of European taurine ancestry than that of QCC, and this may explain the larger body size of QNC, approaching that of European cattle under long-term domestication and selection. A neighbor-joining tree revealed that QCC individuals were closely related, whereas QNC formed a distinct group. To search for signatures of selection in the QNC genome, we evaluated nucleotide diversity (θπ), the fixation index (FST) and Tajima's D. Overlapping selective sweeps were enriched for one KEGG pathway, the apelin signaling pathway, and included five candidate genes (MEF2A, SMAD2, CAMK4, RPS6, and PIK3CG). We performed a comprehensive review of genomic variants in QNC, QCC, and ZSC using whole-genome sequencing data. QCC was rich in novel genetic diversity, while diversity in QNC and ZSC cattle was reduced due to strong artificial selection, with divergence from the original cattle. CONCLUSIONS: We identified candidate genes associated with production traits. These results support the success of selective breeding and can guide further breeding and resource conservation of Qinchuan cattle.

Discovering a novel genetic variant in 11 family members who had isolated pheochromocytoma linked to von Hippel-Lindau (VHL) syndrome, aligning with the type 2c phenotype.

INTRODUCTION: Von Hippel-Lindau disease (e.g. VHL) is an autosomal dominant multi-organ cancer syndrome caused by a mutation in the VHL tumour suppressor gene. In this study, we introduce a novel genetic variant found in 11 family members diagnosed initially with isolated Pheochromocytoma. Subsequent findings revealed its association with VHL syndrome and corresponds to the Type 2 C phenotype. METHODS: The VHL gene was amplified through the utilisation of the polymerase chain reaction (PCR). PCR fragments were sequenced using bidirectional Sanger sequencing, using BigDye™ Terminator v3.1 Cycle Sequencing Kit, running on the 3500 genetic analyser. Results were assembled and analysed Using Software SeqA and chromas pro. RESULTS: A heterozygous in-frame duplication of three nucleotides, specifically ATG, c.377_379dup; p.Asp126dup in exon 2, was identified in all the patients tested within the pedigree. CONCLUSION: In this study, we disclose the identification of a novel genetic variant in a Jordanian family, affecting eleven family members with pheochromocytoma associated with VHL disease. This finding underscores the importance of screening family members and contemplating genetic testing for individuals newly diagnosed with pheochromocytoma and could enhance our comprehension of the potential adverse consequences associated with VHL germline mutations.

Goal: To study a novel gene change in a family with Von Hippel-Lindau (e.g. VHL) syndrome, which increases cancer chances.Participants: 11 family members with Pheochromocytoma, a tumour linked to VHL.Methods:Used PCR to copy the VHL gene.Analysed the gene using Sanger sequencing.Findings:Found a novel gene change in all family members. This change, called an in-frame duplication, affects a protein.It's in a specific part of the gene.Conclusion:Stressing the importance of genetic testing for Pheochromocytoma patients to grasp VHL mutation risks.

Future carbon sequestration potential in a widespread transcontinental boreal tree species: Standing genetic variation matters!

Climate change (CC) necessitates reforestation/afforestation programs to mitigate its impacts and maximize carbon sequestration. But comprehending how tree growth, a proxy for fitness and resilience, responds to CC is critical to maximize these programs' effectiveness. Variability in tree response to CC across populations can notably be influenced by the standing genetic variation encompassing both neutral and adaptive genetic diversity. Here, a framework is proposed to assess tree growth potential at the population scale while accounting for standing genetic variation. We applied this framework to black spruce (BS, Picea mariana [Mill] B.S.P.), with the objectives to (1) determine the key climate variables having impacted BS growth response from 1974 to 2019, (2) examine the relative roles of local adaptation and the phylogeographic structure in this response, and (3) project BS growth under two Shared Socioeconomic Pathways while taking standing genetic variation into account. We modeled growth using a machine learning algorithm trained with dendroecological and genetic data obtained from over 2600 trees (62 populations divided in three genetic clusters) in four 48-year-old common gardens, and simulated growth until year 2100 at the common garden locations. Our study revealed that high summer and autumn temperatures negatively impacted BS growth. As a consequence of warming, this species is projected to experience a decline in growth by the end of the century, suggesting maladaptation to anticipated CC and a potential threat to its carbon sequestration capacity. This being said, we observed a clear difference in response to CC within and among genetic clusters, with the western cluster being more impacted than the central and eastern clusters. Our results show that intraspecific genetic variation, notably associated with the phylogeographic structure, must be considered when estimating the response of widespread species to CC.

Genome Sequencing for Diagnosing Rare Diseases.

BACKGROUND: Genetic variants that cause rare disorders may remain elusive even after expansive testing, such as exome sequencing. The diagnostic yield of genome sequencing, particularly after a negative evaluation, remains poorly defined. METHODS: We sequenced and analyzed the genomes of families with diverse phenotypes who were suspected to have a rare monogenic disease and for whom genetic testing had not revealed a diagnosis, as well as the genomes of a replication cohort at an independent clinical center. RESULTS: We sequenced the genomes of 822 families (744 in the initial cohort and 78 in the replication cohort) and made a molecular diagnosis in 218 of 744 families (29.3%). Of the 218 families, 61 (28.0%) - 8.2% of families in the initial cohort - had variants that required genome sequencing for identification, including coding variants, intronic variants, small structural variants, copy-neutral inversions, complex rearrangements, and tandem repeat expansions. Most families in which a molecular diagnosis was made after previous nondiagnostic exome sequencing (63.5%) had variants that could be detected by reanalysis of the exome-sequence data (53.4%) or by additional analytic methods, such as copy-number variant calling, to exome-sequence data (10.8%). We obtained similar results in the replication cohort: in 33% of the families in which a molecular diagnosis was made, or 8% of the cohort, genome sequencing was required, which showed the applicability of these findings to both research and clinical environments. CONCLUSIONS: The diagnostic yield of genome sequencing in a large, diverse research cohort and in a small clinical cohort of persons who had previously undergone genetic testing was approximately 8% and included several types of pathogenic variation that had not previously been detected by means of exome sequencing or other techniques. (Funded by the National Human Genome Research Institute and others.).

Novel phages of Pseudomonas syringae unveil numerous potential auxiliary metabolic genes.

Relatively few phages that infect plant pathogens have been isolated and investigated. The Pseudomonas syringae species complex is present in various environments, including plants. It can cause major crop diseases, such as bacterial canker on apricot trees. This study presents a collection of 25 unique phage genomes that infect P. syringae. These phages were isolated from apricot orchards with bacterial canker symptoms after enrichment with 21 strains of P. syringae. This collection comprises mostly virulent phages, with only three being temperate. They belong to 14 genera, 11 of which are newly discovered, and 18 new species, revealing great genetic diversity within this collection. Novel DNA packaging systems have been identified bioinformatically in one of the new phage species, but experimental confirmation is required to define the precise mechanism. Additionally, many phage genomes contain numerous potential auxiliary metabolic genes with diversified putative functions. At least three phages encode genes involved in bacterial tellurite resistance, a toxic metalloid. This suggests that viruses could play a role in bacterial stress tolerance. This research emphasizes the significance of continuing the search for new phages in the agricultural ecosystem to unravel novel ecological diversity and new gene functions. This work contributes to the foundation for future fundamental and applied research on phages infecting phytopathogenic bacteria.

Genetic variation and evolutionary characteristics of Echovirus 11: new variant within genotype D5 associated with neonatal death found in China.

Echovirus 11 (E11) has gained attention owing to its association with severe neonatal infections. From 2018 to 2023, a surge in severe neonatal cases and fatalities linked to a novel variant of genotype D5 was documented in China, France, and Italy. However, the prevention and control of E11 variants have been hampered by limited background data on the virus circulation and genetic variance. Therefore, the present study investigated the circulating dynamics of E11 and the genetic variation and molecular evolution of genotype D5 through the collection of strains from the national acute flaccid paralysis (AFP) and hand, foot, and mouth disease (HFMD) surveillance system in China during 2000-2022 and genetic sequences published in the GenBank database. The results of this study revealed a prevalent dynamic of E11 circulation, with D5 being the predominant genotype worldwide. Further phylogenetic analysis of genotype D5 indicated that it could be subdivided into three important geographic clusters (D5-CHN1: 2014-2019, D5-CHN2: 2016-2022, and D5-EUR: 2022-2023). Additionally, variant-specific (144) amino acid mutation sites and positive-selection pressure sites (132, 262) were identified in the VP1 region. Cluster-specific recombination patterns were also identified, with CVB5, E6, and CVB4 as the major recombinant viruses. These findings provide a preliminary landscape of E11 circulation worldwide and basic scientific data for further study of the pathogenicity of E11 variants.

Whole-genome sequencing reveals genomic diversity and selection signatures in Xia'nan cattle.

BACKGROUND: The crossbreeding of specialized beef cattle breeds with Chinese indigenous cattle is a common method of genetic improvement. Xia'nan cattle, a crossbreed of Charolais and Nanyang cattle, is China's first specialized beef cattle breed with independent intellectual property rights. After more than two decades of selective breeding, Xia'nan cattle exhibit a robust physique, good environmental adaptability, good tolerance to coarse feed, and high meat production rates. This study analyzed the population genetic structure, genetic diversity, and genomic variations of Xia'nan cattle using whole-genome sequencing data from 30 Xia'nan cattle and 178 published cattle genomic data. RESULT: The ancestry estimating composition analysis showed that the ancestry proportions for Xia'nan cattle were mainly Charolais with a small amount of Nanyang cattle. Through the genetic diversity studies (nucleotide diversity and linkage disequilibrium decay), we found that the genomic diversity of Xia'nan cattle is higher than that of specialized beef cattle breeds in Europe but lower than that of Chinese native cattle. Then, we used four methods to detect genome candidate regions influencing the excellent traits of Xia'nan cattle. Among the detected results, 42 genes (θπ and CLR) and 131 genes (FST and XP-EHH) were detected by two different detection strategies. In addition, we found a region in BTA8 with strong selection signals. Finally, we conducted functional annotation on the detected genes and found that these genes may influence body development (NR6A1), meat quality traits (MCCC1), growth traits (WSCD1, TMEM68, MFN1, NCKAP5), and immunity (IL11RA, CNTFR, CCL27, SLAMF1, SLAMF7, NAA35, and GOLM1). CONCLUSION: We elucidated the genomic features and population structure of Xia'nan cattle and detected some selection signals in genomic regions potentially associated with crucial economic traits in Xia'nan cattle. This research provided a basis for further breeding improvements in Xia'nan cattle and served as a reference for genetic enhancements in other crossbreed cattle.

Diversity study of Beauveria bassiana species for finding the most virulent strain to manage Bemisia tabaci in cotton.

Beauveria bassiana (Bal.-Criv.) is an important entomopathogenic fungus being used for the management of various agricultural pests worldwide. However, all strains of B. bassiana may not be effective against whitefly, Bemisia tabaci, or other pests, and strains show diversity in their growth, sporulation, virulence features, and overall bioefficacy. Thus, to select the most effective strain, a comprehensive way needs to be devised. We studied the diversity among the 102 strains of B. bassiana isolated from 19 insect species based on their physiological features, virulence, and molecular phylogeny, to identify promising ones for the management of B. tabaci. Strains showed diversity in mycelial growth, conidial production, and their virulence against B. tabaci nymphs. The highest nymphal mortality (2nd and 3rd instar) was recorded with MTCC-4511 (95.1%), MTCC-6289 (93.8%), and MTCC-4565 (89.9%) at a concentration of 1 × 106 conidia ml-1 under polyhouse conditions. The highest bioefficacy index (BI) was in MTCC-4511 (78.3%), MTCC-4565 (68.2%), and MTCC-4543 (62.1%). MTCC-4511, MTCC-4565, and MTCC-4543 clustered with positive loading of eigenvalues for the first two principal components and the cluster analysis also corresponded well with PCA (principal component analysis) (nymphal mortality and BI). The molecular phylogeny could not draw any distinct relationship between physiological features, the virulence of B. bassiana strains with the host and location. The BI, PCA, and square Euclidean distance cluster were found the most useful tools for selecting potential entomopathogenic strains. The selected strains could be utilized for the management of the B. tabaci nymphal population in the field through the development of effective formulations. KEY POINTS: • 102 B. bassiana strains showed diversity in growth and virulence against B. tabaci. • Bioefficacy index, PCA, and SED group are efficient tools for selecting potential strains. • MTCC-4511, 4565, and 4543 chosen as the most virulent strains to kill whitefly nymphs.

Model-driven characterization of functional diversity of Pseudomonas aeruginosa clinical isolates with broadly representative phenotypes.

Pseudomonas aeruginosa is a leading cause of infections in immunocompromised individuals and in healthcare settings. This study aims to understand the relationships between phenotypic diversity and the functional metabolic landscape of P. aeruginosa clinical isolates. To better understand the metabolic repertoire of P. aeruginosa in infection, we deeply profiled a representative set from a library of 971 clinical P. aeruginosa isolates with corresponding patient metadata and bacterial phenotypes. The genotypic clustering based on whole-genome sequencing of the isolates, multilocus sequence types, and the phenotypic clustering generated from a multi-parametric analysis were compared to each other to assess the genotype-phenotype correlation. Genome-scale metabolic network reconstructions were developed for each isolate through amendments to an existing PA14 network reconstruction. These network reconstructions show diverse metabolic functionalities and enhance the collective P. aeruginosa pangenome metabolic repertoire. Characterizing this rich set of clinical P. aeruginosa isolates allows for a deeper understanding of the genotypic and metabolic diversity of the pathogen in a clinical setting and lays a foundation for further investigation of the metabolic landscape of this pathogen and host-associated metabolic differences during infection.

[Prevalence and genetic diversity of the Alongshan virus (Flaviviridae) circulating in ticks in the south of Eastern Siberia].

INTRODUCTION: Tick-borne infections are of great importance for many regions of Russia, including Eastern Siberia. This unfavorable epidemiological situation can be characterized not only by the circulation of well-known tick-borne infections, but also by the identification of new pathogens, the role of which remains little or generally unexplored. Multicomponent flavi-like viruses can cause infectious diseases in humans and pose a threat to public health. The purpose of the study was the identification and molecular genetic characterization of the Alongshan virus (Flaviviridae, ALSV) isolates, transmitted by ticks in the south of Eastern Siberia. MATERIALS AND METHODS: Total 1060 ticks were collected and analyzed from the territory of the Republics of Khakassia, Tuva, Buryatia, Irkutsk Region and Transbaikal Territory (Zabaykalsky Krai) in the spring-summer period 2023. ALSV RNA was detected by RT-PCR followed by nucleotide sequence determination and phylogenetic analysis for each segment of the genome. RESULTS: The ALSV infection rate in Ixodes persulcatus ticks collected in the Republic of Khakassia was 3.3% (95% CI: 1.4-7.5); in Irkutsk Oblast - 1.0% (95% CI: 0.3-3.7); in the Republic of Tuva - 0.9% (95% CI: 0.3-3.4) and in Transbaikal Krai - 0.7% (95% CI: 0.2-3.6). Sequences of all four segments of ALSV genetic variants circulating in I. persulcatus ticks in the south of Eastern Siberia are grouped with sequences found in China and clustered into the Asian subgroup transmitted by taiga ticks. The level of difference in the nucleotide sequences of genome fragments among the identified genetic variants of ALSV ranged from 2 to 3%. CONCLUSION: The article shows the widespread distribution of ALSV in I. persulcatus ticks in the Republics of Khakassia and Tyva, Irkutsk Oblast and Transbaikal Territory. The obtained data actualize monitoring of changes in the area of distribution of potentially dangerous for humans flavi-like viruses and their vectors.

Phylogenetic analysis of variants of the Puumala virus (Hantaviridae: Orthohantavirus) circulating in the Saratov region.

The objective is to determine the complete nucleotide sequence and conduct a phylogenetic analysis of genome variants of the Puumala virus isolated in the Saratov region. MATERIALS AND METHODS: The samples for the study were field material collected in the Gagarinsky (formerly Saratovsky), Engelssky, Novoburassky and Khvalynsky districts of the Saratov region in the period from 2019 to 2022. To specifically enrich the Puumala virus genome in the samples, were used PCR and developed a specific primer panel. Next, the resulting PCR products were sequenced and the fragments were assembled into one sequence for each segment of the virus genome. To construct phylogenetic trees, the maximum parsimony algorithm was used. RESULTS: Genetic variants of the Puumala virus isolated in the Saratov region have a high degree of genome similarity to each other, which indicates their unity of origin. According to phylogenetic analysis, they all form a separate branch in the cluster formed by hantaviruses from other subjects of the Volga Federal District. The virus variants from the Republics of Udmurtia and Tatarstan, as well as from the Samara and Ulyanovsk regions, are closest to the samples from the Saratov region. CONCLUSION: The data obtained show the presence of a pronounced territorial confinement of strains to certain regions or areas that are the natural biotopes of their carriers. This makes it possible to fairly accurately determine the territory of possible infection of patients and/or the circulation of carriers of these virus variants based on the sequence of individual segments of their genome.

Optimization of genetic distance threshold for inferring the CRF01_AE molecular network based on next-generation sequencing.

Introduction: HIV molecular network based on genetic distance (GD) has been extensively utilized. However, the GD threshold for the non-B subtype differs from that of subtype B. This study aimed to optimize the GD threshold for inferring the CRF01_AE molecular network. Methods: Next-generation sequencing data of partial CRF01_AE pol sequences were obtained for 59 samples from 12 transmission pairs enrolled from a high-risk cohort during 2009 and 2014. The paired GD was calculated using the Tamura-Nei 93 model to infer a GD threshold range for HIV molecular networks. Results: 2,019 CRF01_AE pol sequences and information on recent HIV infection (RHI) from newly diagnosed individuals in Shenyang from 2016 to 2019 were collected to construct molecular networks to assess the ability of the inferred GD thresholds to predict recent transmission events. When HIV transmission occurs within a span of 1-4 years, the mean paired GD between the sequences of the donor and recipient within the same transmission pair were as follow: 0.008, 0.011, 0.013, and 0.023 substitutions/site. Using these four GD thresholds, it was found that 98.9%, 96.0%, 88.2%, and 40.4% of all randomly paired GD values from 12 transmission pairs were correctly identified as originating from the same transmission pairs. In the real world, as the GD threshold increased from 0.001 to 0.02 substitutions/site, the proportion of RHI within the molecular network gradually increased from 16.6% to 92.3%. Meanwhile, the proportion of links with RHI gradually decreased from 87.0% to 48.2%. The two curves intersected at a GD of 0.008 substitutions/site. Discussion: A suitable range of GD thresholds, 0.008-0.013 substitutions/site, was identified to infer the CRF01_AE molecular transmission network and identify HIV transmission events that occurred within the past three years. This finding provides valuable data for selecting an appropriate GD thresholds in constructing molecular networks for non-B subtypes.

Unveiling the genetic structure of pig population in a Himalayan state Uttarakhand through microsatellite and mitochondrial DNA analyses.

This study traced the maternal lineage of the domestic swine populations using mitochondrial DNA control region markers and genetic diversity using microsatellite markers in Uttarakhand, an Indian state situated at the foothills of the world's youngest (geo-dynamically sensitive) mountain system, "the Himalayas". Analysis of 68 maternally unrelated individuals revealed 20 haplotypes. The maternal signature of the Pacific, Southeast Asian, European, and ubiquitously distributed Chinese haplotypes was present in Uttarakhand's domestic pig population. The D3 haplotype reported in wild pigs from North India was also identified in 47 domestic samples. A unique gene pool, UKD (Uttarakhand Domestic), as another lineage specific to this region has been proposed. Genotypes were analyzed, using 13 sets of microsatellite markers. The observed (Ho) and expected (He) heterozygosities were 0.83 ± 0.02 and 0.84 ± 0.01, respectively. The average polymorphic information content value of 0.83 ± 0.01 indicated the high informativeness of the marker. The overall mean FIS value for all the microsatellite markers was low (F = 0.04, P < 0.01). Seven loci deviated from Hardy-Weinberg equilibrium (HWE) at a significant level (p < 0.05). Two clusters were identified, indicating overlapping populations. These results suggested that though belonging to different maternal lineages, the traditional management practices in Uttarakhand have allowed for genetic mixing and the sharing of genetic material among pig populations. It could contribute to increased genetic diversity but might also result in the loss of distinct genetic characteristics or breed purity of the local breeds if not carefully managed.

Prevalence and genetic diversity of Anaplasma and Ehrlichia in ticks and domesticated animals in Suizhou County, Hubei Province, China.

Anaplasma and Ehrlichia are tick-borne bacterial pathogens that cause anaplasmoses and ehrlichioses in humans and animals. In this study, we examined the prevalence of Anaplasma and Ehrlichia species in ticks and domesticated animals in Suizhou County, Hubei Province in the central China. We used PCR amplification and DNA sequencing of the 16S rRNA, groEL, and gltA genes to analyze. We collected 1900 ticks, including 1981 Haemaphysalis longicornis and 9 Rhipicephalus microplus, 159 blood samples of goats (n = 152), cattle (n = 4), and dogs (n = 3) from May to August of 2023. PCR products demonstrated that Anaplasma bovis, Anaplasma capra, and an Ehrlichia species were detected in the H. longicornis with the minimum infection rates (MIR) of 1.11%, 1.32%, and 0.05%, respectively; A. bovis, A. capra, and unnamed Anaplasma sp. were detected in goats with an infection rate of 26.31%, 1.31% and 1.97%, respectively. Anaplasma and Ehrlichia species were not detected from cattle, dogs and R. microplus ticks. The genetic differences in the groEL gene sequences of the Anaplasma in the current study were large, whereas the 16S rRNA and gltA gene sequences were less disparate. This study shows that ticks and goats in Suizhou County, Hubei Province carry multiple Anaplasma species and an Ehrlichia species, with relatively higher infection rate of A. bovis in goats. Our study indicates that multiple Anaplasma and Ehrlichia species exist in ticks and goats in the central China with potential to cause human infection.

Variability of Chalcone Synthase in Chamomile Accessions (Matricaria chamomilla).

Chamomile (Matricaria chamomilla) is an important medicinal plant whose beneficial activities partly rely on certain flavonoids. The first dedicated step in flavonoid biosynthesis is chalcone synthase (CHS, EC 2.3.1.74). The genomic DNA of CHS was studied in six chamomile specimens from different genotypes to describe interspecimen, as well as interspecific, variability. One specimen of M. discoidea was included as an outgroup. The two exons of CHS of M. chamomilla (McCHS) and M. discoidea (MdCHS) were 188 bp and 1,011 bp long, separated by an intron of variable length between 192 and 199 bp in McCHS and 201 bp in MdCHS, respectively. The two exons with 5.3 and 6.2 mutations per 100 bp, respectively, were more conserved than the intron with 11.5 mutations per 100 bp. In total, 96 SNPs were detected in both species, of which 12 SNPs were only present in MdCHS and 80 SNPs only in McCHS. Overall, 70 haplotypes (multilocus genotypes, MLGs) were detected. The samples could be classified into two groups, a 'compact' group of a low number and diversity of haplotypes and a 'variable' group of a high number and diversity of haplotypes. Of the 74 SNPs in McCHS, only six SNPs were non-synonymous. However, the amino acid changes did not affect critical areas of the enzyme. The combination of the six SNPs resulted in nine translated amino acid MLGs. The CHS network located MdCHS, due to the crossing barrier, quite distant from chamomile. MdCHS docked to McCHS at a position from where McCHS divergently evolved into two directions.

dgfr: an R package to assess sequence diversity of gene families.

BACKGROUND: Gene families are groups of homologous genes that often have similar biological functions. These families are formed by gene duplication events throughout evolution, resulting in multiple copies of an ancestral gene. Over time, these copies can acquire mutations and structural variations, resulting in members that may vary in size, motif ordering and sequence. Multigene families have been described in a broad range of organisms, from single-celled bacteria to complex multicellular organisms, and have been linked to an array of phenomena, such as host-pathogen interactions, immune evasion and embryonic development. Despite the importance of gene families, few approaches have been developed for estimating and graphically visualizing their diversity patterns and expression profiles in genome-wide studies. RESULTS: Here, we introduce an R package named dgfr, which estimates and enables the visualization of sequence divergence within gene families, as well as the visualization of secondary data such as gene expression. The package takes as input a multi-fasta file containing the coding sequences (CDS) or amino acid sequences from a multigene family, performs a pairwise alignment among all sequences, and estimates their distance, which is subjected to dimension reduction, optimal cluster determination, and gene assignment to each cluster. The result is a dataset that allows for the visualization of sequence divergence and expression within the gene family, an approximation of the number of clusters present in the family. CONCLUSIONS: dgfr provides a way to estimate and study the diversity of gene families, as well as visualize the dispersion and secondary profile of the sequences. The dgfr package is available at https://github.com/lailaviana/dgfr under the GPL-3 license.

Genome sequencing of captive white tigers from Bangladesh.

OBJECTIVES: The Bengal tiger Panthera tigris tigris, is an emblematic animal for Bangladesh. Despite being the apex predator in the wild, their number is decreasing due to anthropogenic activities such as hunting, urbanization, expansion of agriculture and deforestation. By contrast, captive tigers are flourishing due to practical conservation efforts. Breeding within the small captive population can produce inbreeding depression and genetic bottlenecks, which may limit the success of conservation efforts. Despite past decades of research, a comprehensive database on genetic variation in the captive and wild Bengal tigers in Bangladesh still needs to be included. Therefore, this research aimed to investigate the White Bengal tiger genome to create a resource for future studies to understand variation underlying important functional traits. DATA DESCRIPTION: Blood samples from Chattogram Zoo were collected for three white Bengal tigers. Genomic DNA for all collected samples were extracted using a commercial DNA extraction kit. Whole genome sequencing was performed using a DNBseq platform. We generated 77 Gb of whole-genome sequencing (WGS) data for three white Bengal tigers (Average 11X coverage/sample). The data we generated will establish a paradigm for tiger research in Bangladesh by providing a genomic resource for future functional studies on the Bengal white tiger.

Phylogenetic analysis and divergence time estimation of Lycium species in China based on the chloroplast genomes.

BACKGROUND: Lycium is an economically and ecologically important genus of shrubs, consisting of approximately 70 species distributed worldwide, 15 of which are located in China. Despite the economic and ecological importance of Lycium, its phylogeny, interspecific relationships, and evolutionary history remain relatively unknown. In this study, we constructed a phylogeny and estimated divergence time based on the chloroplast genomes (CPGs) of 15 species, including subspecies, of the genus Lycium from China. RESULTS: We sequenced and annotated 15 CPGs in this study. Comparative analysis of these genomes from these Lycium species revealed a typical quadripartite structure, with a total sequence length ranging from 154,890 to 155,677 base pairs (bp). The CPGs was highly conserved and moderately differentiated. Through annotation, we identified a total of 128-132 genes. Analysis of the boundaries of inverted repeat (IR) regions showed consistent positioning: the junctions of the IRb/LSC region were located in rps19 in all Lycium species, IRb/SSC between the ycf1 and ndhF genes, and SSC/IRa within the ycf1 gene. Sequence variation in the SSC region exceeded that in the IR region. We did not detect major expansions or contractions in the IR region or rearrangements or insertions in the CPGs of the 15 Lycium species. Comparative analyses revealed five hotspot regions in the CPG: trnR(UCU), atpF-atpH, ycf3-trnS(GGA), trnS(GGA), and trnL-UAG, which could potentially serve as molecular markers. In addition, phylogenetic tree construction based on the CPG indicated that the 15 Lycium species formed a monophyletic group and were divided into two typical subbranches and three minor branches. Molecular dating suggested that Lycium diverged from its sister genus approximately 17.7 million years ago (Mya) and species diversification within the Lycium species of China primarily occurred during the recent Pliocene epoch. CONCLUSION: The divergence time estimation presented in this study will facilitate future research on Lycium, aid in species differentiation, and facilitate diverse investigations into this economically and ecologically important genus.

First detection of Tetraparvovirus ungulate 1 in diseased cattle (Chinese Simmental) from Hunan province, China.

Tetraparvovirus is an emerging parvovirus infecting a variety of mammals and humans, and associated with human diseases including severe acute respiratory infection and acute encephalitis syndrome. In the present study, a Tetraparvovirus ungulate 1 (formerly known as bovine hokovirus) strain HNU-CBY-2023 was identified and characterized from diseased Chinese Simmental from Hunan province, China. The nearly complete genome of HNU-CBY-2023 is 5346 nt in size and showed genomic identities of 85-95.5% to the known Tetraparvovirus ungulate 1 strains from GenBank, indicating a rather genetic variation. Phylogenetic and genetic divergence analyses indicated that Tetraparvovirus ungulate 1 could be divided into two genotypes (I and II), and HNU-CBY-2023 was clustered into genotype II. This study, for the first time, identified Tetraparvovirus ungulate 1 from domestic cattle from mainland China, which will be helpful to understand the prevalence and genetic diversity of Tetraparvovirus ungulate 1.