Genomic characters of Anaplasma bovis and genetic diversity in China

The emergence of Anaplasma bovis or A. bovis-like infection in humans from China and the United States of America has raised concern about the public health importance of this pathogen. Although A. bovis has been detected in a wide range of ticks and mammals in the world, no genome of the pathogen is available up to now, which has prohibited us from better understanding the genetic basis for its pathogenicity. Here we describe an A. bovis genome from metagenomic sequencing of an infected goat in China. Anaplasma bovis had the smallest genome of the genus Anaplasma, and relatively lower GC content. Phylogenetic analysis of single-copy orthologue sequence showed that A. bovis was closely related to A. platys and A. phagocytophilum, but relatively far from intraerythrocytic Anaplasma species. Anaplasma bovis had 116 unique orthogroups and lacked 51 orthogroups in comparison to other Anaplasma species. The virulence factors of A. bovis were significantly less than those of A. phagocytophilum, suggesting less pathogenicity of A. bovis. When tested by specific PCR assays, A. bovis was detected in 23 of 29 goats, with an infection rate up to 79.3% (95% CI: 64.6% ∼94.1%). The phylogenetic analyses based on partial 16S rRNA, gltA and groEL genes indicated that A. bovis had high genetic diversity. The findings of this study lay a foundation for further understanding of the biological characteristics and genetic diversity of A. bovis, and will facilitate the formulation of prevention and control strategies.

Simultaneous determination of subspecies and geographic origins of 110 rice cultivars by microsatellite markers

Rice varieties of different subspecies types (indica rice and japonica rice) across various geographical origins (Hunan, Jiangsu, and Northeast China) were monitored using microsatellite markers (simple sequence repeats, SSR). 110 representative rice cultivars were collected from the main crop areas. Multiple methods including clustering analysis (neighbor-joining (NJ) method, unweighted pair-group method with arithmetic mean (UPGMA) method), principal component analysis (PCA) and model-based grouping were applied. The study revealed that 25 pairs of SSR markers exhibited a broad range of polymorphism information content (PIC) values, ranging from 0.240 to 0.830. Furthermore, our study successfully achieved a higher overall mean correct rate of 99.09% in determining the geographical origin of rice. Simultaneously, it accurately classified indica rice and japonica rice. These findings are significant as they provide an SSR fingerprint of 110 high-quality rice cultivars, serving as a valuable scientific resource for the detection of rice adulteration and traceability of its origin.

Integrative phylogenetic, phylogeographic and morphological characterisation of the Unio crassus species complex reveals cryptic diversity with important conservation implications

The global decline of freshwater mussels and their crucial ecological services highlight the need to understand their phylogeny, phylogeography and patterns of genetic diversity to guide conservation efforts. Such knowledge is urgently needed for Unio crassus, a highly imperilled species originally widespread throughout Europe and southwest Asia. Recent studies have resurrected several species from synonymy based on mitochondrial data, revealing U. crassus to be a complex of cryptic species. To address long-standing taxonomic uncertainties hindering effective conservation, we integrate morphometric, phylogenetic, and phylogeographic analyses to examine species diversity within the U. crassus complex across its entire range. Phylogenetic analyses were performed using cytochrome c oxidase subunit I (815 specimens from 182 populations) and, for selected specimens, whole mitogenome sequences and Anchored Hybrid Enrichment (AHE) data on ∼ 600 nuclear loci. Mito-nuclear discordance was detected, consistent with mitochondrial DNA gene flow between some species during the Pliocene and Pleistocene. Fossil-calibrated phylogenies based on AHE data support a Mediterranean origin for the U. crassus complex in the Early Miocene. The results of our integrative approach support 12 species in the group: the previously recognised Unio bruguierianus, Unio carneus, Unio crassus, Unio damascensis, Unio ionicus, Unio sesirmensis, and Unio tumidiformis, and the reinstatement of five nominal taxa: Unio desectusstat. rev., Unio gontieriistat. rev., Unio mardinensisstat. rev., Unio nanusstat. rev., and Unio vicariusstat. rev. Morphometric analyses of shell contours reveal important morphospace overlaps among these species, highlighting cryptic, but geographically structured, diversity. The distribution, taxonomy, phylogeography, and conservation of each species are succinctly described.

Genetic diversity of glutamate-rich protein (GLURP) in Plasmodium falciparum isolates from school-age children in Kinshasa, DRC

Malaria infections in school-age children further make it difficult to control the disease's spread. Moreover, the genetic diversity of glutamate-rich protein, potentially a candidate for vaccine development, has not yet been investigated in the Democratic Republic of Congo. Therefore, we aimed to assess the genetic diversity of the immunodominant C-terminal repetitive region (R2) of Plasmodium falciparum glutamate-rich protein gene (pfglurp) among school-age children living in Kinshasa, DRC. We conducted nested PCR targeting R2 of pfglurp and the amplicon were directly sequenced. We summarized the prevalence of mutations of bases and amino acids and indicated the amino acid repeat sequence in the R2 region by the unit code. We then statistically analyzed whether there was a relationship between the number of mutations in the pfglurp gene and attributes. In 221 samples, haplotype 1 was the most common (n = 137, 61.99%), with the same sequence as the 3D7 strain. Regarding the number of base mutations, it was higher in urban areas than rural areas (p = 0.0363). When genetic neutrality was tested using data from 171 samples of the single strain, Tajima's D was -1.857 (p = 0.0059). In addition, FST as the genetic distance between all attributes was very small and no significant difference was observed. This study clarified the genetic mutation status and relevant patient attributes among School-age children in the DRC. We found that urban areas are more likely to harbour pfglurp mutations. Future research needs to clarify the reason and mechanism involved.

Geomorphic impacts within Red River Fault and island shifting as witnessed by the phylogeography of the largest water strider

Palaeogeological events and climate oscillations profoundly impact the demographics and distributions of small-range species, increasing the extinction risk. The largest water strider worldwide, Gigantometra gigas (Hemiptera: Gerridae), exhibits restricted distributions in Vietnam and southern China. Herein, we generated three genomic datasets (mitogenomes, 146 nuclear protein-coding genes and single nucleotide polymorphisms) with ecological niche modelling (ENM) to explicitly test whether the present-day distribution of G. gigas actually resulted from geographical and climatic effects. We found that the origin of this largest water strider reached the divergence time of the genus within Gerridae, providing a greater opportunity to explore its response to geographic movements. The right-lateral motion of the Red River Fault facilitated the divergence of two phylogeographic lineages, resulting in the "north-south component" genetic pattern in G. gigas. The Hainan and southeast Vietnam populations of the southern linage were completely separated by the Beibu Gulf but exhibited similar genetic compositions, confirming that Hainan had a continental origin and that Hainan Island joined with the Indo-China Peninsula to promote gene exchange among populations. Additionally, we noticed the low genetic diversity but long demographic history of the northern lineage, which displayed population dynamics opposite to those of other organisms. Integrating the demographic changes and ENM findings revealed that suitable habitat contraction and rapid demographic decline during the Last Glacial Maximum (LGM) triggered the low genetic diversity of the northern lineage. Overall, the demographic history of the largest water strider was mainly shaped by geographical features, and first provided evidence from the phylogeographic perspective of aquatic insects to support the hypothesis of Hainan Island shifting.

Wide spectrum of novel and rare hemoglobin variants in the multi-ethnic Indian population: A review

The hemoglobin (Hb) variants are qualitative abnormalities due to production of structurally abnormal globin proteins. They are categorized based on the type of mutation present in the α1, α2, β, Gγ, Aγ and δ globin genes. So far, more than 1550 Hb variants are reported in the database. They could lead to Hb polymerization, Hb instability, altered oxygen affinity and decreased oxygen-carrying capacity of Hb or have no clinical manifestations. In India, ethnic diversity, consanguinity, regional variations and migration result in the presence of different Hb variants. We have compiled all the variants of α, β and δ globin chains in heterozygous, homozygous and in compound heterozygous forms reported from India in the last 52 years. Of the 63 rare and novel hemoglobin variants reported from India, 22 were α-globin chain variants, 37 were β-globin chain variants and 4 were δ-globin chain variants. Twelve novel Hb variants (Hb J Rajappan, Hb Koya Dora, Hb Rampa, Hb Godavari, Hb Chandigarh, Hb D Agri, Hb Lucknow, Hb Vellore, Hb Midnapore, Hb Bijnor, Hb A2Tianhe and Hb A2Saurashtra) were identified among persons of Indian origin. Majority of them were picked up on HPLC. Some of the variants like Hb Titusville, Hb Shimonoseki, Hb Chandigarh, Hb D Agri, Hb Yaizu and Hb Vellore eluted in the HbS window whereas variants like HbD Iran, Hb St. Louis, Hb G Coushata, HbM Saskatoon, Hb Lucknow, Hb Grange-Blanche and Hb Tianshui showed falsely elevated HbA2. Hence, careful and systematic investigations are required to identify them.

Distribution of human papillomavirus among Vietnamese women with cervical cancer and unusual genetic variability of HPV16

HPV16, with typical mutations that differ in geographical distribution and carcinogenic potency, has implications for cervical cancer screening, clinical diagnosis, and treatment. DNASTAR and MEGA were used to identify HPV16 variants and construct a phylogenetic tree. The most prevalent HPV genotypes were HPV16 (63.9%), HPV18 (26.7%), and other HPV (6.9%). HPV16 alterations were found in all E6, E7, and L1 genes, including 15 missense and 18 synonymous mutations. Missense mutations include R10G, Q14H, D25E, H78Y, L83V (E6); M29V, R35K, L78R, L95P (E7); H73Y, T176 N, N178T, T317P, T386S, L472F/I (L1). HPV16 sublineages include A1 (17.2%), A2 (0.9%), A3 (56.0%), A4 (19.0%), D1 (4.3%), and D3 (2.6%). Although several mutations in the oncoproteins E6, E7, and L1 have been detected, mutations known to be associated with cervical cancer risk, such as D25E and L83V, occur at a relatively low frequency. This suggests that HPV16 mutations are associated with cervical cancer through a complicated mechanism.

Identifying a potentially invasive population in the native range of a species: The enlightenment from the phylogeography of the yellow spotted stink bug, Erthesina fullo (Hemiptera: Pentatomidae)

The yellow spotted stink bug (YSSB), Erthesina fullo (Thunberg, 1783) is an important Asian pest that has recently successfully invaded Europe and an excellent material for research on the initial stage of biological invasion. Here, we reported the native evolutionary history, recent invasion history, and potential invasion threats of YSSB for the first time based on population genetic methods [using double digest restriction-site associated DNA (ddRAD) data and mitochondrial COI and CYTB] and ecological niche modelling. The results showed that four lineages (east, west, southwest, and Hainan Island) were established in the native range with a strong east-west differentiation phylogeographical structure, and the violent climate fluctuation might cause population divergence during the Middle and Upper Pleistocene. In addition, land bridges and monsoon promote dispersal and directional genetic exchanging between island populations and neighboring continental populations. The east lineage (EA) was identified as the source of invasion in Albania. EA had the widest geographical distribution among all other lineages, with a star-like haplotype network with the main haplotype as the core. It also had a rapid population expansion history, indicating that the source lineage might have stronger diffusion ability and adaptability. Our findings provided a significant biological basis for fine tracking of invasive source at the lineage or population level and promote early invasion warning of potential invasive species on a much subtler lineage level.

Shared phylogeographic patterns and environmental responses of co-distributed soybean pests: Insights from comparative phylogeographic studies of Riptortus pedestris and Riptortus linearis in the subtropics of East Asia

Comparative phylogeographic studies of closely related species sharing co-distribution areas can elucidate the role of shared historical factors and environmental changes in shaping their phylogeographic pattern. The bean bugs, Riptortus pedestris and Riptortus linearis, which both inhabit subtropical regions in East Asia, are recognized as highly destructive soybean pests. Many previous studies have investigated the biological characteristics, pheromones, chemicals and control mechanisms of these two pests, but few studies have explored their phylogeographic patterns and underlying factors. In this study, we generated a double-digest restriction site-associated DNA sequencing (ddRAD-seq) dataset to investigate phylogeographic patterns and construct ecological niche models (ENM) for both Riptortus species. Our findings revealed similar niche occupancies and population genetic structures between the two species, with each comprising two phylogeographic lineages (i.e., the mainland China and the Indochina Peninsula clades) that diverged approximately 0.1 and 0.3 million years ago, respectively. This divergence likely resulted from the combined effects of temperatures variation and geographical barriers in the mountainous regions of Southwest China. Further demographic history and ENM analyses suggested that both pests underwent rapid expansion prior to the Last Glacial Maximum (LGM). Furthermore, ENM predicts a northward shift of both pests into new soybean-producing regions due to global warming. Our study indicated that co-distribution soybean pests with overlapping ecological niches and similar life histories in subtropical regions of East Asia exhibit congruent phylogeographic and demographic patterns in response to shared historical biogeographic drivers.

The Poultry Red Mite, Dermanyssus gallinae, travels far but not frequently, and takes up permanent residence on farms

Management of Dermanyssus gallinae, a cosmopolitan hematophagous mite responsible for damage in layer poultry farming, is hampered by a lack of knowledge of its spatio-temporal population dynamics. Previous studies have shown that the circulation of this pest between farms is of strictly anthropogenic origin, that a mitochondrial haplogroup has been expanding on European farms since the beginning of the 21st century and that its local population growth may be particularly rapid. To refine our understanding of how D. gallinae spreads within and among farms, we characterized the genetic structure of mite populations at different spatial scales and sought to identify the main factors interrupting gene flow between poultry houses and between mitochondrial haplogroups. To this end, we selected and validated the first set of nuclear microsatellite markers for D. gallinae and sequenced a region of the CO1-encoding mitochondrial gene in a subsample of microsatellite-genotyped mites. We also tested certain conditions required for effective contamination of a poultry house through field experimentation, and conducted a survey of practices during poultry transfers. Our results confirm the role of poultry transport in the dissemination of mite populations, but the frequency of effective contamination after the introduction of contaminated material into poultry houses seems lower than expected. The high persistence of mites on farms, even during periods when poultry houses are empty and cleaned, and the very large number of nodes in the logistic network (large number of companies supplying pullets or transporting animals) undoubtedly explain the very high prevalence on farms. Substantial genetic diversity was measured in farm populations, probably as a result of the mite's known haplodiploid mode of sexual reproduction, coupled with the dense logistic network. The possibility of the occasional occurrence of asexual reproduction in this sexually reproducing mite was also revealed in our analyses, which could explain the extreme aggressiveness of its demographic dynamics under certain conditions.

Evolutionary analysis of LMP-1 genetic diversity in EBV-associated nasopharyngeal carcinoma: Bioinformatic insights into oncogenic potential

EBV latent membrane protein 1 (LMP-1) is an important oncogene involved in the induction and maintenance of EBV infection and the activation of several cell survival and proliferative pathways. The genetic diversity of LMP-1 has an important role in immunogenicity and tumorigenicity allowing escape from host cell immunity and more metastatic potential of LMP-1 variants. This study explored the evolutionary of LMP-1 in EBV-infected patients at an advanced stage of nasopharyngeal carcinoma (NPC). Detection of genetic variability in LMP-1 genes was carried out using Sanger sequencing. Bioinformatic analysis was conducted for translation and nucleotide alignment. Phylogenetic analysis was used to construct a Bayesian tree for a deeper understanding of the genetic relationships, evolutionary connections, and variations between sequences. Genetic characterization of LMP-1 in NPC patients revealed the detection of polymorphism in LMP-1 Sequences. Motifs were identified within three critical LMP-1 domains, such as PQQAT within CTAR1 and YYD within CTAR2. The presence of the JACK3 region at specific sites within CTAR3, as well as repeat regions at positions (122-132) and (133-143) within CTAR3, was also annotated. Additionally, several mutations were detected including 30 and 69 bp deletions, 33 bp repeats, and 15 bp insertion. Although LMP-1 strains appear to be genetically diverse, they are closely related to 3 reference strains: prototype B95.8, Med- 30 bp deletion, and Med + 30 bp deletion. In our study, one of the strains harboring the 30 bp deletion had both bone and bone marrow metastasis which could be attributed to the fact that LMP-1 is involved in tumor metastasis, evasion and migration of NPC cells. This study provided valuable insights into genetic variability in LMP-1 sequences of EBV in NPC patients. Further functional studies would provide a more comprehensive understanding of the molecular characteristics, epidemiology, and clinical implications of LMP-1 polymorphisms in EBV-related malignancies.

Genetic underpinnings of arthropod community distributions in Populus trichocarpa

Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. We surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.

Genetic Variation in Male Aggression Is Influenced by Genotype of Prior Social Partners in Drosophila melanogaster

AbstractSocial behaviors can be influenced by the genotypes of interacting individuals through indirect genetic effects (IGEs) and can also display developmental plasticity. We investigated how developmental IGEs, which describe the effects of a prior social partner's genotype on later behavior, can influence aggression in male Drosophila melanogaster. We predicted that developmental IGEs cannot be estimated by simply extending the effects of contextual IGEs over time and instead have their own unique effects on behavior. On day 1 of the experiment, we measured aggressive behavior in 15 genotypic pairings (n = 600 males). On day 2, each of the males was paired with a new opponent, and aggressive behavior was again measured. We found contextual IGEs on day 1 of the experiment and developmental IGEs on day 2 of the experiment: the influence of the day 1 partner's genotype on the focal individual's day 2 behavior depended on the genotypic identity of both the day 1 partner and the focal male. Importantly, the developmental IGEs in our system produced fundamentally different dynamics than the contextual IGEs, as the presence of IGEs was altered over time. These findings represent some of the first empirical evidence demonstrating developmental IGEs, a first step toward incorporating developmental IGEs into our understanding of behavioral evolution.

The role of historical biogeography in shaping colour morph diversity in the common wall lizard

The maintenance of polymorphisms often depends on multiple selective forces, but less is known on the role of stochastic or historical processes in maintaining variation. The common wall lizard (Podarcis muralis) is a colour polymorphic species in which local colour morph frequencies are thought to be modulated by natural and sexual selection. Here, we used genome-wide single-nucleotide polymorphism data to investigate the relationships between morph composition and population biogeography at a regional scale, by comparing morph composition with patterns of genetic variation of 54 populations sampled across the Pyrenees. We found that genetic divergence was explained by geographic distance but not by environmental features. Differences in morph composition were associated with genetic and environmental differentiation, as well as differences in sex ratio. Thus, variation in colour morph frequencies could have arisen via historical events and/or differences in the permeability to gene flow, possibly shaped by the complex topography and environment. In agreement with this hypothesis, colour morph diversity was positively correlated with genetic diversity and rates of gene flow and inversely correlated with the likelihood of the occurrence of bottlenecks. Concurrently, we did not find conclusive evidence for selection in the two colour loci. As an illustration of these effects, we observed that populations with higher proportions of the rarer yellow and yellow-orange morphs had higher genetic diversity. Our results suggest that processes involving a decay in overall genetic diversity, such as reduced gene flow and/or bottleneck events have an important role in shaping population-specific morph composition via non-selective processes.

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

New insight into genetic diversity of zoonotic-potential Ancylostoma ceylanicum in stray cats living in Bangkok, Thailand, based on deep amplicon sequencing

AIMS

This study aimed to characterize feline hookworms from stray cats living in Bangkok.

METHODS AND RESULTS

A total of 56 hookworm-positive faecal samples were identified for hookworm species by using PCR targeting the ITS1, 5.8S, and ITS2 fragment and qPCR targeting ITS2. Of 56 samples, 96.4% (54/56) were identified as Ancylostoma ceylanicum and 1.8% (1/56) as Ancylostoma caninum. With qPCR, 89.3% (50/56) were identified as single A. ceylanicum infection and 5.4% (3/56) as coinfection of A. ceylanicum and A. caninum. For genetic characterization of A. ceylanicum, 10 samples were pooled, and the partial COI gene was amplified, followed by deep amplicon sequencing. Five pooled samples were analysed, and 99.73% were identified with A. ceylanicum sequences, which were allocated into 19 haplotypes (AC01-AC19). Genetic diversity findings for A. ceylanicum in Asia revealed that three of eight haplotypes considered of zoonotic significance occurred in humans, dogs, and cats, including haplotypes H01, H20, and H21. The predominant haplotype in this study, AC01, was clustered with H01-a zoonotic haplotype.

CONCLUSIONS

The diversity obtained by deep amplicon sequencing supported that the A. ceylanicum community had high genetic variation. Deep amplicon sequencing was a useful method to determine source, zoonotic potential, and host-parasite relationship.

Subtypes of Blastocystis in Tibetan Antelope (Pantholops hodgsonii)

Blastocystis is a protist that is distributed in the gut tract of humans and animals. However, the reports about Blastocystis infection in Tibetan antelope are scarce. We collected 173 Tibetan antelope feces samples from Xinjiang, Qinghai and Xizang, and amplified the SSU rRNA gene of 600 bp region of Blastocystis in our research. Fifty-one samples in total were positive for Blastocystis, with all subtypes being ST31. The lowest prevalence of Blastocystis was observed in Xizang (2/20, 9.1%), followed by Qinghai (18/92, 16.4%), Xinjiang (31/61, 33.7%). The highest prevalence of Blastocystis in Tibetan antelope was detected during the summer was (19/30, 38.8%). This is the first research work regarding the Blastocystis subtypes ST31 in Tibetan antelope. Our research provides information for future researches on the distribution of this Blastocystis subtype and the control of Blastocystis infection.

Genetic characteristics of Blastocystis sp. in cattle from Hebei Province, China

Blastocystis sp. is a protozoan parasite that infects the intestines of humans and animals, causing chronic diseases such as skin rashes, abdominal pain, and irritable bowel syndrome. A survey was conducted to determine the prevalence and genetic diversity of Blastocystis sp. infection in cattle, in Hebei Province, China. 2746 cattle fecal samples were collected from 11 cities in Hebei Province and analyzed using polymerase chain reaction targeting the Blastocystis sp. barcoding gene. MEGA, PhyloSuite, and PopART were used to analyze the subtype, sequence signature, pairwise genetic distance, and genetic diversity indices. The results showed that the Blastocystis sp. detection rate was 12.60% (346/2746). The infection rate in different herds was affected by region, age, breeding mode, and variety; that is, the infection rates in areas of southern Hebei, cattle under one year old, intensive raising, and dairy cattle were higher than the infection rates in northern Hebei, cattle over one year old, scatter feeding, and beef cattle. Seven Blastocystis subtypes were identified, namely, ST1, ST2, ST5, ST10, ST14, ST21, and ST26; ST10 was the dominant subtype, and ST14 was the second most common subtype. A total of 374 polymorphic and conserved sites were obtained, including 273 invariable (monomorphic) sites and 101 variable (polymorphic) sites, accounting for 27.01% of all nucleotides. The nucleotide diversity index (Pi) was 0.07749, and the haplotype (gene) diversity index (Hd) was 0.946. This study provides the first comprehensive information on the epidemiological situation of Blastocystis sp. infection in cattle from Hebei Province, China, and revealed rich genetic diversity of Blastocystis sp.

Investigation of Chlamydia pecorum in livestock from Switzerland reveals a high degree of diversity in bovine strains

Chlamydia pecorum is a widespread veterinary chlamydial species causing endemic infections in livestock, such as ruminants and pigs, globally. However, there is limited contemporary knowledge on infecting strain diversity in various hosts. This study aimed to evaluate the genetic diversity of C. pecorum strains infecting Swiss livestock through C. pecorum genotyping and phylogenetic analyses in comparison to the global population, while also assessing chlamydial strains for plasmid carriage. A total of 263 C. pecorum positive samples from clinically healthy ruminant and pig herds (Bovines = 216, sheep = 25, pigs = 14) as well as placentae from eight C. pecorum positive ruminant abortion cases from other Swiss herds were investigated. The ompA and Multi-Locus sequence typing revealed novel C. pecorum genotypes, and bovine strains exhibited considerable genetic diversity, contrasting with lower diversity in sheep and pig strains. C. pecorum plasmid was detected in 100.0% of sheep (41/41) and pig (255/255) samples, and in 69.4% of bovine samples (150/216). In contrast, no plasmid was detected in the eight C. pecorum-positive ruminant abortion cases either representing plasmid-less strains or possibly escaping PCR detection due to autolysis of the placenta. This study supports the genetic diversity of C. pecorum strains, particularly in bovines, and identifies novel sequence types in Swiss livestock.

A comprehensive review of HVS-I mitochondrial DNA variation of 19 Iranian populations

Iran is located along the Central Asian corridor, a natural artery that has served as a cross-continental route since the first anatomically modern human populations migrated out of Africa. We compiled and reanalyzed the HVS-I (hypervariable segment-I) of 3840 mitochondrial DNA (mtDNA) sequences from 19 Iranian populations and from 26 groups from adjacent countries to give a comprehensive review of the maternal genetic variation and investigate the impact of historical events and cultural factors on the maternal genetic structure of modern Iranians. We conclude that Iranians have a high level of genetic diversity. Thirty-six haplogroups were observed in Iran's populations, and most of them belong to widespread West-Eurasian haplogroups, such as H, HV, J, N, T, and U. In contrast, the predominant haplogroups observed in most of the adjacent countries studied here are H, M, D, R, U, and C haplogroups. Using principal component analysis, clustering, and genetic distance-based calculations, we estimated moderate genetic relationships between Iranian and other Eurasian groups. Further, analyses of molecular variance and comparing geographic and genetic structures indicate that mtDNA HVS-I sequence diversity does not exhibit any sharp geographic structure in the country. Barring a few from some culturally distinct and naturally separated minorities, most Iranian populations have a homogenous maternal genetic structure.

Pedigree diversity and implications for genetic selection of Katahdin sheep

The Katahdin hair breed gained popularity in the United States as low input and prolific, with a propensity to exhibit parasite resistance. With the introduction of genomically enhanced estimated breeding values (GEBV) to the Katahdin genetic evaluation, defining the diversity present in the breed is pertinent. Utilizing pedigree records (n = 92,030) from 1984 to 2019 from the National Sheep Improvement Program, our objectives were to (i) estimate the completeness and quality of the pedigree, (ii) calculate diversity statistics for the whole pedigree and relevant reference subpopulations and (iii) assess the impact of current diversity on genomic selection. Reference 1 was Katahdins born from 2017 to 2019 (n = 23,494), while reference 2 was a subset with at least three generations of Katahdin ancestry (n = 9327). The completeness of the whole pedigree, and the pedigrees of reference 1 and reference 2, were above 50% through the fourth, fifth and seventh generation of ancestors, respectively. Effective population size (Ne) averaged 111 animals with a range from 42.2 to 451.0. The average generation interval was 2.9 years for the whole pedigree and reference 1, and 2.8 years for reference 2. The mean individual inbreeding and average relatedness coefficients were 1.62% and 0.91%, 1.74% and 0.90% and 2.94% and 1.46% for the whole pedigree, reference 1, and reference 2, respectively. There were over 300 effective founders in the whole pedigree and reference 1, with 169 in reference 2. Effective number of ancestors were over 150 for the whole pedigree and reference 1, while there were 67 for reference 2. Prediction accuracies increased as the reference population grew from 1k to 7.5k and plateaued at 15k animals. Given the large number of founders and ancestors contributing to the base genetic variation in the breed, the Ne is sufficient to maintain diversity while achieving progress with selection. Stable low rates of inbreeding and relatedness suggest that incorporating genetic conservation in breeding decisions is currently not of high priority. Current Ne suggests that with limited genotyping, high levels of accuracy for genomic prediction can be achieved. However, intense selection on GEBV may cause loss of genetic diversity long term.

Vertically inherited microbiota and environment modifying behaviours conceal genetic variation in dung beetle life history

Diverse organisms actively manipulate their (sym)biotic and physical environment in ways that feed back on their own development. However, the degree to which these processes affect microevolution remains poorly understood. The gazelle dung beetle both physically modifies its ontogenetic environment and structures its biotic interactions through vertical symbiont transmission. By experimentally eliminating (i) physical environmental modifications and (ii) the vertical inheritance of microbes, we assess how environment modifying behaviour and microbiome transmission shape heritable variation and evolutionary potential. We found that depriving larvae of symbionts and environment modifying behaviours increased additive genetic variance and heritability for development time but not body size. This suggests that larvae's ability to manipulate their environment has the potential to modify heritable variation and to facilitate the accumulation of cryptic genetic variation. This cryptic variation may become released and selectable when organisms encounter environments that are less amenable to organismal manipulation or restructuring. Our findings also suggest that intact microbiomes, which are commonly thought to increase genetic variation of their hosts, may instead reduce and conceal heritable variation. More broadly, our findings highlight that the ability of organisms to actively manipulate their environment may affect the potential of populations to evolve when encountering novel, stressful conditions.

Boom-bust population dynamics drive rapid genetic change

Increasing environmental threats and more extreme environmental perturbations place species at risk of population declines, with associated loss of genetic diversity and evolutionary potential. While theory shows that rapid population declines can cause loss of genetic diversity, populations in some environments, like Australia's arid zone, are repeatedly subject to major population fluctuations yet persist and appear able to maintain genetic diversity. Here, we use repeated population sampling over 13 y and genotype-by-sequencing of 1903 individuals to investigate the genetic consequences of repeated population fluctuations in two small mammals in the Australian arid zone. The sandy inland mouse (Pseudomys hermannsburgensis) experiences marked boom-bust population dynamics in response to the highly variable desert environment. We show that heterozygosity levels declined, and population differentiation (FST) increased, during bust periods when populations became small and isolated, but that heterozygosity was rapidly restored during episodic population booms. In contrast, the lesser hairy-footed dunnart (Sminthopsis youngsoni), a desert marsupial that maintains relatively stable population sizes, showed no linear declines in heterozygosity. These results reveal two contrasting ways in which genetic diversity is maintained in highly variable environments. In one species, diversity is conserved through the maintenance of stable population sizes across time. In the other species, diversity is conserved through rapid genetic mixing during population booms that restores heterozygosity lost during population busts.

Multiple genetic variants at the SLC30A8 locus affect local super-enhancer activity and influence pancreatic β-cell survival and function

Variants at the SLC30A8 locus are associated with type 2 diabetes (T2D) risk. The lead variant, rs13266634, encodes an amino acid change, Arg325Trp (R325W), at the C-terminus of the secretory granule-enriched zinc transporter, ZnT8. Although this protein-coding variant was previously thought to be the sole driver of T2D risk at this locus, recent studies have provided evidence for lowered expression of SLC30A8 mRNA in protective allele carriers. In the present study, we examined multiple variants that influence SLC30A8 allele-specific expression. Epigenomic mapping has previously identified an islet-selective enhancer cluster at the SLC30A8 locus, hosting multiple T2D risk and cASE associations, which is spatially associated with the SLC30A8 promoter and additional neighboring genes. Here, we show that deletion of variant-bearing enhancer regions using CRISPR-Cas9 in human-derived EndoC-βH3 cells lowers the expression of SLC30A8 and several neighboring genes and improves glucose-stimulated insulin secretion. While downregulation of SLC30A8 had no effect on beta cell survival, loss of UTP23, RAD21, or MED30 markedly reduced cell viability. Although eQTL or cASE analyses in human islets did not support the association between these additional genes and diabetes risk, the transcriptional regulator JQ1 lowered the expression of multiple genes at the SLC30A8 locus and enhanced stimulated insulin secretion.

Copy number and selection of MHC genes in ruminants are related to habitat, average life span and diet

The ruminants, as the main group of livestock, have been extensively studied in terms of their physiology, endocrinology, biochemistry, genetics, and nutrition. Despite the wide geographic distribution and habitat diversity of animals in this group, their ecology and evolution remain poorly understood. In this study, we analyzed the gene copy number, selection, and ecological and evolutionary processes that have affected the evolution of major histocompatibility complex (MHC) genes across ruminant lineages based on available genomic data. The 51 species analyzed represented all six families of ruminants. Our finding indicated that the architecture of the MHC region is conserved in ruminants, but with variable copy numbers of MHC-I, MHC-IIA, and MHC-IIB genes. No lineage-specific gene duplication was observed in the MHC genes. The phylogenetic generalized least squares regression (PGLS) model revealed association between ecological and biological factors (habitat and lifespan) and gene duplication in DQA and DQB, but not in DRB. The selection pressure of DQA and DQB were related with lifespan, diet, and the ratio of genetic repeat elements. These results suggest that the MHC evolution in ruminants, including copy number and selection, has been influenced by genetic repeat elements, pathogen exposure risk, and intrinsic cost of possessing multiple MHC genes.

[Characterization of sequences, expression profiling, and natural allelic variation analysis of the MYC gene family in sorghum (Sorghum bicolor)]

Sorghum aphid (Melanaphis sacchari) and head smut fungi (Sporisorium reilianum) infesting sorghum cause delayed growth and development, and reduce yield and quality. This study use bioinformatics and molecular biological approaches to profile the gene expression pattern during sorghum development and under pest infestation, and analyzed the natural allelic DNA variation of sorghum MYC gene family. The findings provide insights for potential application in breeding the stress resistant and high productivity sorghum varieties. The results indicated that there are 28 MYC genes identified in sorghum genome, distributed on 10 chromosomes. The bHLH_MYC_N and HLH domains are the conserved domains of the MYC gene in sorghum. Gene expression analysis showed that SbbHLH35.7g exhibited high expression levels in leaves, SbAbaIn showed strong expression in early grains, and SbMYC2.1g showed high expression levels in mature pollen. In anti-aphid strains at the 5-leaf stage, SbAbaIn, SbLHW.4g and SbLHW.2g were significantly induced in leaves, while SbbHLH35.7g displayed the highest expression level in panicle tissue, which was significantly induced by the infection of head smut. Promoter cis-element analysis identified methyl jasmonate (MJ), abscisic acid (ABA), salicylic acid (SA) and MYB-binding sites related to drought-stress inducibility. Furthermore, genomic resequencing data analysis revealed natural allelic DNA variations such as single nucleotide polymorphism (SNP) and insertion-deletion (INDEL) for the key SbMYCs. Protein interaction network analysis using STRING indicated that SbAbaIn interacts with TIFYdomain protein, and SbbHLH35.7g interacts with MDR and imporin. SbMYCs exhibited temporal and spatial expression patterns and played vital roles during the sorghum development. Infestation by sugarcane aphids and head smut fungi induced the expression of SbAbaIn and SbbHLH35.7g, respectively. SbAbaIn modulated the jasmonic acid (JA) pathway to regulate the expression of defensive genes, conferring resistance to insects. On the other hand, SbbHLH35.7g participated in detoxification reactions to defend against pathogens.

[Genetic background of lily germplasm resources based on SSR markers]

To study the genetic background of lily (Lilium spp.) germplasm resources, and accurately evaluate and select excellent germplasm for genetic improvement of lily, we analyzed the genetic background of 62 lily germplasm accessions from 11 provinces of China by using simple sequence repeat (SSR) molecular markers. The results showed that 15 out of 83 pairs of lily SSR primers were polymorphic. A total of 157 allelic loci were amplified, with the number of alleles per locus ranging from 5 to 19 and the average number of effective alleles per locus being 4.162 8. The average observed heterozygosity and expected heterozygosity were 0.228 2 and 0.694 1, respectively. The average polymorphic information content was 0.678 8. The average Nei's diversity index and Shannon's information index were 0.694 1 and 1.594 9, respectively, indicating that the tested lily germplasm had high genetic diversity. The 62 germplasm accessions were classified into 5 groups by the unweighted pair group method with arithmetic mean (UPGMA) and into 3 groups by the principal component analysis. The two analyses revealed a geographic correlation among different groups. The majority of lily germplasm accessions from the same source tended to cluster together. The population structure analysis classified the lily accessions into 4 populations and 1 mixed population. The above results provide a theoretical basis and genetic resources for the precise identification and breeding of lily germplasm resources.

Genetic diversity of KPC-2-producing Klebsiella pneumoniae complex from aquatic ecosystems

During the COVID-19 pandemic, the occurrence of carbapenem-resistant Klebsiella pneumoniae increased in human clinical settings worldwide. Impacted by this increase, international high-risk clones harboring carbapenemase-encoding genes have been circulating in different sources, including the environment. The blaKPC gene is the most commonly disseminated carbapenemase-encoding gene worldwide, whose transmission is carried out by different mobile genetic elements. In this study, blaKPC-2-positive Klebsiella pneumoniae complex strains were isolated from different anthropogenically affected aquatic ecosystems and characterized using phenotypic, molecular, and genomic methods. K. pneumoniae complex strains exhibited multidrug-resistant and extensively drug-resistant profiles, spotlighting the resistance to carbapenems, ceftazidime-avibactam, colistin, and tigecycline, which are recognized as last-line antimicrobial treatment options. Molecular analysis showed the presence of several antimicrobial resistance, virulence, and metal tolerance genes. In-depth analysis showed that the blaKPC-2 gene was associated with three different Tn4401 isoforms (i.e., Tn4401a, Tn4401b, and Tn4401i) and NTEKPC elements. Different plasmid replicons were detected and a conjugative IncN-pST15 plasmid harboring the blaKPC-2 gene associated with Tn4401i was highlighted. K. pneumoniae complex strains belonging to international high-risk (e.g., ST11 and ST340) and unusual clones (e.g., ST323, ST526, and ST4216) previously linked to clinical settings. In this context, some clones were reported for the first time in the environmental sector. Therefore, these findings evidence the occurrence of carbapenemase-producing K. pneumoniae complex strains in aquatic ecosystems and contribute to the monitoring of carbapenem resistance worldwide.

PPanG: a precision pangenome browser enabling nucleotide-level analysis of genomic variations in individual genomes and their graph-based pangenome

Graph-based pangenome is gaining more popularity than linear pangenome because it stores more comprehensive information of variations. However, traditional linear genome browser has its own advantages, especially the tremendous resources accumulated historically. With the fast-growing number of individual genomes and their annotations available, the demand for a genome browser to visualize genome annotation for many individuals together with a graph-based pangenome is getting higher and higher. Here we report a new pangenome browser PPanG, a precise pangenome browser enabling nucleotide-level comparison of individual genome annotations together with a graph-based pangenome. Nine rice genomes with annotations were provided by default as potential references, and any individual genome can be selected as the reference. Our pangenome browser provides unprecedented insights on genome variations at different levels from base to gene, and reveals how the structures of a gene could differ for individuals. PPanG can be applied to any species with multiple individual genomes available and it is available at https://cgm.sjtu.edu.cn/PPanG .

Treatment of ZC4H2 Variant-Associated Spastic Paraplegia with Selective Dorsal Rhizotomy and Intensive Postoperative Rehabilitation: A Case Report

Selective dorsal rhizotomy (SDR) has been used to treat children with spastic cerebral palsy (CP), and its beneficial effect on quality of life and ambulation has been confirmed in long-term follow-up studies. However, the role of SDR in the treatment of spasticity in patients with hereditary spastic paraplegia (HSP) and related disorders is not well-established. Here, we report the first patient with the ZC4H2 variant who underwent SDR to treat spastic paraplegia. Abnormal gait was discovered during a regular checkup at the age of 3 years and 9 months, and she was diagnosed with spastic paraplegia. She was heterozygous for the ZC4H2 variant and underwent SDR at the age of 5 years and 11 months, which alleviated the spasticity. The patient underwent inpatient postoperative rehabilitation for 4 months and continued outpatient physiotherapy after discharge. The Gross Motor Function Measure-88 score and maximum walking speed decreased transiently 1 month postoperatively, but gradually recovered, and continuously improved 6 months postoperatively. SDR and postoperative intensive rehabilitation were effective in improving motor and walking functions up to 6 months after surgery, although long-term follow-up is needed to draw conclusions.

Patrilineal segmentary systems provide a peaceful explanation for the post-Neolithic Y-chromosome bottleneck

Studies have found a pronounced decline in male effective population sizes worldwide around 3000-5000 years ago. This bottleneck was not observed for female effective population sizes, which continued to increase over time. Until now, this remarkable genetic pattern was interpreted as the result of an ancient structuring of human populations into patrilineal groups (gathering closely related males) violently competing with each other. In this scenario, violence is responsible for the repeated extinctions of patrilineal groups, leading to a significant reduction in male effective population size. Here, we propose an alternative hypothesis by modelling a segmentary patrilineal system based on anthropological literature. We show that variance in reproductive success between patrilineal groups, combined with lineal fission (i.e., the splitting of a group into two new groups of patrilineally related individuals), can lead to a substantial reduction in the male effective population size without resorting to the violence hypothesis. Thus, a peaceful explanation involving ancient changes in social structures, linked to global changes in subsistence systems, may be sufficient to explain the reported decline in Y-chromosome diversity.

Comprehensive genome assembly reveals genetic diversity and carcass consumption insights in critically endangered Asian king vultures

The Asian king vulture (AKV), a vital forest scavenger, is facing globally critical endangerment. This study aimed to construct a reference genome to unveil the mechanisms underlying its scavenger abilities and to assess the genetic relatedness of the captive population in Thailand. A reference genome of a female AKV was assembled from sequencing reads obtained from both PacBio long-read and MGI short-read sequencing platforms. Comparative genomics with New World vultures (NWVs) and other birds in the Family Accipitridae revealed unique gene families in AKV associated with retroviral genome integration and feather keratin, contrasting with NWVs' genes related to olfactory reception. Expanded gene families in AKV were linked to inflammatory response, iron regulation and spermatogenesis. Positively selected genes included those associated with anti-apoptosis, immune response and muscle cell development, shedding light on adaptations for carcass consumption and high-altitude soaring. Using restriction site-associated DNA sequencing (RADseq)-based genome-wide single nucleotide polymorphisms (SNPs), genetic relatedness and inbreeding status of five captive AKVs were determined, revealing high genomic inbreeding in two females. In conclusion, the AKV reference genome was established, providing insights into its unique characteristics. Additionally, the potential of RADseq-based genome-wide SNPs for selecting AKV breeders was demonstrated.

Identification and characterization of functionally relevant SSR markers in natural Dalbergia odorifera populations

BACKGROUND

Dalbergia odorifera is a rare and precious rosewood specie, which is valued for its amber tones, abstract figural patterns, and impermeability to water and insects. However, the information on genetic diversity and marker-assisted selection breeding of D. odorifera is still limited. Simple sequence repeat (SSR) markers are an ideal tool for genetic diversity analysis and marker-assisted molecular breeding for complex traits.

RESULTS

Here, we have developed SSR markers within candidate genes and used them to explore the genetic diversity among D. odorifera germplasm resources. A total of 635 SSR loci were identified. The proportions of mono-, di- and tri-nucleotide repeat motifs were 52.28%, 22.99% and 21.42%, respectively. From these, a total of 114 SSR primers were synthesized, of which 24 SSR markers displayed polymorphism (polymorphic information content (PIC) > 0.25). Subsequently, these polymorphic markers were used for the genetic diversity analysis of 106 D. odorifera individuals from 11 natural populations. According to the genetic diversity analysis of D. odorifera natural populations, the average observed heterozygosity (Ho) was 0.500, the average expected heterozygosity (He) was 0.524, and the average Shannon's information index (I) was 0.946. These indicated that the natural populations had moderate genetic diversity. AMOVA analysis showed that 5% of the total variation was within the individuals of a population, whereas 95% of the variation was among the individuals of the populations, indicating a high degree of genetic variation between populations. On the basis of their genetic structures, these populations could be divided into four groups.

CONCLUSIONS

Our study provides important experimental resources for genetic studies and assists in the program of molecular breeding of D. odorifera wood formation.

[Current status and prospects of genetic research on IgA nephropathy]

IgA nephropathy is the most common primary glomerulonephritis worldwide, and genetic factors may play an important role in its pathogenesis. Following candidate gene association analysis and genome-wide linkage study, genome-wide association studies (GWAS) have found multiple susceptibility genes related to the pathogenesis and clinical phenotype of IgA nephropathy. Meanwhile, structural variation and epigenetic changes are also closely related to IgA nephropathy. Genetic variants have been found to explain about 11% of its heritability. In the current era of genomic medicine, how to find more susceptible genes/loci, whole genome sequencing studies (WGS) provide clues to further understand the genetic variation of IgA nephropathy. How to find the cell type-specific susceptibility genes associated with IgA nephropathy, multi-omics studies will conduct comprehensive analysis via single-cell sequencing, expression quantitative trait locus (eQTL) and genomics to find the pathogenic genes and offer insights into the development of targeted drugs, which will be the trend and direction of future research.

The size diversity of the Pteridaceae family chloroplast genome is caused by overlong intergenic spacers

BACKGROUND

While the size of chloroplast genomes (cpDNAs) is often influenced by the expansion and contraction of inverted repeat regions and the enrichment of repeats, it is the intergenic spacers (IGSs) that appear to play a pivotal role in determining the size of Pteridaceae cpDNAs. This provides an opportunity to delve into the evolution of chloroplast genomic structures of the Pteridaceae family. This study added five Pteridaceae species, comparing them with 36 published counterparts.

RESULTS

Poor alignment in the non-coding regions of the Pteridaceae family was observed, and this was attributed to the widespread presence of overlong IGSs in Pteridaceae cpDNAs. These overlong IGSs were identified as a major factor influencing variations in cpDNA size. In comparison to non-expanded IGSs, overlong IGSs exhibited significantly higher GC content and were rich in repetitive sequences. Species divergence time estimations suggest that these overlong IGSs may have already existed during the early radiation of the Pteridaceae family.

CONCLUSIONS

This study reveals new insights into the genetic variation, evolutionary history, and dynamic changes in the cpDNA structure of the Pteridaceae family, providing a fundamental resource for further exploring its evolutionary research.

Analyzing genetic diversity in luffa and developing a Fusarium wilt-susceptible linked SNP marker through a single plant genome-wide association (sp-GWAS) study

BACKGROUND

Luffa (Luffa spp.) is an economically important crop of the Cucurbitaceae family, commonly known as sponge gourd or vegetable gourd. It is an annual cross-pollinated crop primarily found in the subtropical and tropical regions of Asia, Australia, Africa, and the Americas. Luffa serves not only as a vegetable but also exhibits medicinal properties, including anti-inflammatory, antidiabetic, and anticancer effects. Moreover, the fiber derived from luffa finds extensive applications in various fields such as biotechnology and construction. However, luffa Fusarium wilt poses a severe threat to its production, and existing control methods have proven ineffective in terms of cost-effectiveness and environmental considerations. Therefore, there is an urgent need to develop luffa varieties resistant to Fusarium wilt. Single-plant GWAS (sp-GWAS) has been demonstrated as a promising tool for the rapid and efficient identification of quantitative trait loci (QTLs) associated with target traits, as well as closely linked molecular markers.

RESULTS

In this study, a collection of 97 individuals from 73 luffa accessions including two major luffa species underwent single-plant GWAS to investigate luffa Fusarium wilt resistance. Utilizing the double digest restriction site associated DNA (ddRAD) method, a total of 8,919 high-quality single nucleotide polymorphisms (SNPs) were identified. The analysis revealed the potential for Fusarium wilt resistance in accessions from both luffa species. There are 6 QTLs identified from 3 traits, including the area under the disease progress curve (AUDPC), a putative disease-resistant QTL, was identified on the second chromosome of luffa. Within the region of linkage disequilibrium, a candidate gene homologous to LOC111009722, which encodes peroxidase 40 and is associated with disease resistance in Cucumis melo, was identified. Furthermore, to validate the applicability of the marker associated with resistance from sp-GWAS, an additional set of 21 individual luffa plants were tested, exhibiting 93.75% accuracy in detecting susceptible of luffa species L. aegyptiaca Mill.

CONCLUSION

In summary, these findings give a hint of genome position that may contribute to luffa wild resistance to Fusarium and can be utilized in the future luffa wilt resistant breeding programs aimed at developing wilt-resistant varieties by using the susceptible-linked SNP marker.

Fast and accurate variant identification tool for sequencing-based studies

BACKGROUND

Accurate identification of genetic variants, such as point mutations and insertions/deletions (indels), is crucial for various genetic studies into epidemic tracking, population genetics, and disease diagnosis. Genetic studies into microbiomes often require processing numerous sequencing datasets, necessitating variant identifiers with high speed, accuracy, and robustness.

RESULTS

We present QuickVariants, a bioinformatics tool that effectively summarizes variant information from read alignments and identifies variants. When tested on diverse bacterial sequencing data, QuickVariants demonstrates a ninefold higher median speed than bcftools, a widely used variant identifier, with higher accuracy in identifying both point mutations and indels. This accuracy extends to variant identification in virus samples, including SARS-CoV-2, particularly with significantly fewer false negative indels than bcftools. The high accuracy of QuickVariants is further demonstrated by its detection of a greater number of Omicron-specific indels (5 versus 0) and point mutations (61 versus 48-54) than bcftools in sewage metagenomes predominated by Omicron variants. Much of the reduced accuracy of bcftools was attributable to its misinterpretation of indels, often producing false negative indels and false positive point mutations at the same locations.

CONCLUSIONS

We introduce QuickVariants, a fast, accurate, and robust bioinformatics tool designed for identifying genetic variants for microbial studies. QuickVariants is available at https://github.com/caozhichongchong/QuickVariants .

Clinical implications of VUS reclassification in a single-centre series from application of ACMG/AMP classification rules specified for BRCA1/2

BACKGROUND

BRCA1/2 testing is crucial to guide clinical decisions in patients with hereditary breast/ovarian cancer, but detection of variants of uncertain significance (VUSs) prevents proper management of carriers. The ENIGMA (Evidence-based Network for the Interpretation of Germline Mutant Alleles) BRCA1/2 Variant Curation Expert Panel (VCEP) has recently developed BRCA1/2 variant classification guidelines consistent with ClinGen processes, specified against the ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular-Pathology) classification framework.

METHODS

The ClinGen-approved BRCA1/2-specified ACMG/AMP classification guidelines were applied to BRCA1/2 VUSs identified from 2011 to 2022 in a series of patients, retrieving information from the VCEP documentation, public databases, literature and ENIGMA unpublished data. Then, we critically re-evaluated carrier families based on new results and checked consistency of updated classification with main sources for clinical interpretation of BRCA1/2 variants.

RESULTS

Among 166 VUSs detected in 231 index cases, 135 (81.3%) found in 197 index cases were classified by applying BRCA1/2-specified ACMG/AMP criteria: 128 (94.8%) as Benign/Likely Benign and 7 (5.2%) as Pathogenic/Likely Pathogenic. The average time from the first report as 'VUS' to classification using this approach was 49.4 months. Considering that 15 of these variants found in 64 families had already been internally reclassified prior to this work, this study provided 121 new reclassifications among the 151 (80.1%) remaining VUSs, relevant to 133/167 (79.6%) families.

CONCLUSIONS

These results demonstrated the effectiveness of new BRCA1/2 ACMG/AMP classification guidelines for VUS classification within a clinical cohort, and their important clinical impact. Furthermore, they suggested a cadence of no more than 3 years for regular review of VUSs, which however requires time, expertise and resources.

Tissue-specific genetic variation suggests distinct molecular pathways between body shape phenotypes and colorectal cancer

It remains unknown whether adiposity subtypes are differentially associated with colorectal cancer (CRC). To move beyond single-trait anthropometric indicators, we derived four multi-trait body shape phenotypes reflecting adiposity subtypes from principal components analysis on body mass index, height, weight, waist-to-hip ratio, and waist and hip circumference. A generally obese (PC1) and a tall, centrally obese (PC3) body shape were both positively associated with CRC risk in observational analyses in 329,828 UK Biobank participants (3728 cases). In genome-wide association studies in 460,198 UK Biobank participants, we identified 3414 genetic variants across four body shapes and Mendelian randomization analyses confirmed positive associations of PC1 and PC3 with CRC risk (52,775 cases/45,940 controls from GECCO/CORECT/CCFR). Brain tissue-specific genetic instruments, mapped to PC1 through enrichment analysis, were responsible for the relationship between PC1 and CRC, while the relationship between PC3 and CRC was predominantly driven by adipose tissue-specific genetic instruments. This study suggests distinct putative causal pathways between adiposity subtypes and CRC.

Epidemiological and molecular surveillance of norovirus in the Brazilian Amazon: description of recombinant genotypes and improvement of evolutionary analysis

Noroviruses are highly infectious, genetically diverse viruses. Global outbreaks occur frequently, making molecular surveillance important for infection monitoring. This cross-sectional descriptive study aimed to monitor cases of norovirus gastroenteritis in the Brazilian Amazon. Fecal samples were tested by immunoenzymatic assay, RT-PCR and genetic sequencing for the ORF1/ORF2 and protease regions. Bayesian inference with a molecular clock was employed to construct the phylogeny. The norovirus prevalence was 25.8%, with a higher positivity rate among children aged 0-24 months. Genogroup GII accounted for 98.1% of the sequenced samples, while GI accounted for 1.9% of them. The GII.P16/GII.4 genotype was the most prevalent, with an evolution rate of 2.87x10-3 and TMRCA estimated in 2012. This study demonstrates that norovirus is a primary causative agent of gastroenteritis and provides data on viral genetic diversity that may facilitate infection surveillance and vaccine development.

From sub-Saharan Africa to China: Evolutionary history and adaptation of Drosophila melanogaster revealed by population genomics

Drosophila melanogaster is a widely used model organism for studying environmental adaptation. However, the genetic diversity of populations in Asia is poorly understood, leaving a notable gap in our knowledge of the global evolution and adaptation of this species. We sequenced genomes of 292 D. melanogaster strains from various ecological settings in China and analyzed them along with previously published genome sequences. We have identified six global genetic ancestry groups, despite the presence of widespread genetic admixture. The strains from China represent a unique ancestry group, although detectable differentiation exists among populations within China. We deciphered the global migration and demography of D. melanogaster, and identified widespread signals of adaptation, including genetic changes in response to insecticides. We validated the effects of insecticide resistance variants using population cage trials and deep sequencing. This work highlights the importance of population genomics in understanding the genetic underpinnings of adaptation, an effort that is particularly relevant given the deterioration of ecosystems.

Temporal analysis of genetic diversity and gene flow in the threatened catfish Pseudoplatystoma magdaleniatum from a dammed neotropical river

The striped catfish Pseudoplatystoma magdaleniatum is a large-sized migratory species from the north Andes region, endemic to Magdalena basin and one of the major fishery resources. Despite the estimated reduction of over 80% of the fisheries production of this species throughout the basin in recent decades, its population in the lower Magdalena-Cauca basin showed healthy genetics after molecular analyses. However, the current conservation status of this species and several habitat disturbances demand the re-evaluation of its population genetics to infer evolutionary risks and assess potential changes. This work analyzed a total of 164 samples from the Cauca River collected downstream the Ituango Dam between 2019-2021 using species-specific microsatellite markers to compare the genetic diversity and structure in samples collected between 2010-2014 from the lower Magdalena-Cauca basin, previously analyzed. Our results showed a relatively stable panmictic population over time (4 to 10 years), with high genetic diversity and evidence of recent bottleneck. Promoting habitat connectivity to conserve gene flow, characterizing diversity and genetic structure over the entire basin, and integrating the results with future monitoring are important aspects for the management planning for P. magdaleniatum in the Magdalena-Cauca basin.

Tracing genetic diversity captures the molecular basis of misfolding disease

Genetic variation in human populations can result in the misfolding and aggregation of proteins, giving rise to systemic and neurodegenerative diseases that require management by proteostasis. Here, we define the role of GRP94, the endoplasmic reticulum Hsp90 chaperone paralog, in managing alpha-1-antitrypsin deficiency on a residue-by-residue basis using Gaussian process regression-based machine learning to profile the spatial covariance relationships that dictate protein folding arising from sequence variants in the population. Covariance analysis suggests a role for the ATPase activity of GRP94 in controlling the N- to C-terminal cooperative folding of alpha-1-antitrypsin responsible for the correction of liver aggregation and lung-disease phenotypes of alpha-1-antitrypsin deficiency. Gaussian process-based spatial covariance profiling provides a standard model built on covariant principles to evaluate the role of proteostasis components in guiding information flow from genome to proteome in response to genetic variation, potentially allowing us to intervene in the onset and progression of complex multi-system human diseases.

Sequence characteristics, genetic diversity and phylogenetic analysis of the Cucurbita ficifolia (Cucurbitaceae) chloroplasts genome

BACKGROUND

Curcubita ficifolia Bouché (Cucurbitaceae) has high value as a food crop and medicinal plant, and also has horticultural value as rootstock for other melon species. China is home to many different cultivars, but the genetic diversity of these resources and the evolutionary relationships among them, as well as the differences between C. ficifolia and other Cucurbita species, remain unclear.

RESULTS

We investigated the chloroplast (cp) genomes of 160 C. ficifolia individuals from 31 populations in Yunnan, a major C. ficifolia production area in China. We found that the cp genome of C. ficifolia is ~151 kb and contains 128 genes, of which 86 are protein coding genes, 34 encode tRNA, and eight encode rRNAs. We also identified 64 SSRs, mainly AT repeats. The cp genome was found to contain a total of 204 SNP and 57 indels, and a total of 21 haplotypes were found in the 160 study individuals. The reverse repeat (IR) region of C. ficifolia contained a few differences compared with this region in the six other Cucurbita species. Sequence difference analysis demonstrated that most of the variable regions were concentrated in the single copy (SC) region. Moreover, the sequences of the coding regions were found to be more similar among species than those of the non-coding regions. The phylogenies reconstructed from the cp genomes of 61 representative species of Cucurbitaceae reflected the currently accepted classification, in which C. ficifolia is sister to the other Cucurbita species, however, different interspecific relationships were found between Cucurbita species.

CONCLUSIONS

These results will be valuable in the classification of C. ficifolia genetic resources and will contribute to our understanding of evolutionary relationships within the genus Cucurbita.

DArTseq-based SNP markers reveal high genetic diversity among early generation fall armyworm tolerant maize inbred lines

Diversity analysis using molecular markers serves as a powerful tool in unravelling the intricacies of inclusivity within various populations and is an initial step in the assessment of populations and the development of inbred lines for host plant resistance in maize. This study was conducted to assess the genetic diversity and population structure of 242 newly developed S3 inbred lines using 3,305 single nucleotide polymorphism (SNP) markers and to also assess the level of homozygosity achieved in each of the inbred lines. A total of 1,184 SNP markers were found highly informative, with a mean polymorphic information content (PIC) of 0.23. Gene diversity was high among the inbred lines, ranging from 0.04 to 0.50, with an average of 0.27. The residual heterozygosity of the 242 S3 inbred lines averaged 8.8%, indicating moderately low heterozygosity levels among the inbred lines. Eighty-four percent of the 58,322 pairwise kinship coefficients among the inbred lines were near zero (0.00-0.05), with only 0.3% of them above 0.50. These results revealed that many of the inbred lines were distantly related, but none were redundant, suggesting each inbred line had a unique genetic makeup with great potential to provide novel alleles for maize improvement. The admixture-based structure analysis, principal coordinate analysis, and neighbour-joining clustering were concordant in dividing the 242 inbred lines into three subgroups based on the pedigree and selection history of the inbred lines. These findings could guide the effective use of the newly developed inbred lines and their evaluation in quantitative genetics and molecular studies to identify candidate lines for breeding locally adapted fall armyworm tolerant varieties in Ghana and other countries in West and Central Africa.

Genetic barriers more than environmental associations explain Serratia marcescens population structure

Bacterial species often comprise well-separated lineages, likely emerged and maintained by genetic isolation and/or ecological divergence. How these two evolutionary actors interact in the shaping of bacterial population structure is currently not fully understood. In this study, we investigate the genetic and ecological drivers underlying the evolution of Serratia marcescens, an opportunistic pathogen with high genomic flexibility and able to colonise diverse environments. Comparative genomic analyses reveal a population structure composed of five deeply-demarcated genetic clusters with open pan-genome but limited inter-cluster gene flow, partially explained by Restriction-Modification (R-M) systems incompatibility. Furthermore, a large-scale research on hundred-thousands metagenomic datasets reveals only a partial habitat separation of the clusters. Globally, two clusters only show a separate gene composition coherent with ecological adaptations. These results suggest that genetic isolation has preceded ecological adaptations in the shaping of the species diversity, an evolutionary scenario coherent with the Evolutionary Extended Synthesis.

Haitian coffee agroforestry systems harbor complex arabica variety mixtures and under-recognized genetic diversity

Though facing significant challenges, coffee (Coffea arabica) grown in Haitian agroforestry systems are important contributors to rural livelihoods and provide several ecosystem services. However, little is known about their genetic diversity and the variety mixtures used. In light of this, there is a need to characterize Haitian coffee diversity to help inform revitalization of this sector. We sampled 28 diverse farms in historically important coffee growing regions of northern and southern Haiti. We performed KASP-genotyping of SNP markers and HiPlex multiplex amplicon sequencing for haplotype calling on our samples, as well as several Ethiopian and commercial accessions from international collections. This allowed us to assign Haitian samples to varietal groups. Our analyses revealed considerable genetic diversity in Haitian farms, higher in fact than many farmers realized. Notably, genetic structure analyses revealed the presence of clusters related to Typica, Bourbon, and Catimor groups, another group that was not represented in our reference accession panel, and several admixed individuals. Across the study areas, we found both mixed-variety farms and monovarietal farms with the historical and traditional Typica variety. This study is, to our knowledge, the first to genetically characterize Haitian C. arabica variety mixtures, and report the limited cultivation of C. canephora (Robusta coffee) in the study area. Our results show that some coffee farms are repositories of historical, widely-abandoned varieties while others are generators of new diversity through genetic mixing.

Genetic insights into Cyphocharax magdalenae (Characiformes: Curimatidae): Microsatellite loci development and population analysis in the Cauca River, Colombia

Cyphocharax magdalenae, a Colombian freshwater fish species, plays a vital role in nutrients distribution and serves as a significant food source for other fish species and local fishing communities. Considered a short-distance migratory species, C. magdalenae populations face substantial extinction risk due to human activities impacting their habitats. To address the lack of knowledge on genetic diversity and population structure, this study used next-generation sequencing technology to develop species-specific microsatellite loci and conducted a population genetics analysis of C. magdalenae in the middle and lower sections of the Cauca River, Colombia. Out of 30 pairs of microsatellite primers evaluated in 324 individuals, 14 loci were found to be polymorphic, at linkage equilibrium and, in at least one population, their genotypic frequencies were in Hardy-Weinberg equilibrium. Results showed high genetic diversity levels compared to other neotropical Characiformes, with inbreeding coefficients similar to those reported for phylogenetically related species. Moreover, C. magdalenae exhibits seasonal population structure (rainy-dry) consisting of two genetic stocks showing bottleneck signals and high effective population sizes. This information is essential for understanding the current species genetics and developing future management programs for this fishery resource.