Single-molecule Long-read Sequencing Research Articles

HPV genotyping in clinical samples using long-read single-molecule real-time sequencing.

Human papillomavirus (HPV) genotyping is widely used, particularly in combination with high-risk (HR) HPV tests for cervical cancer screening. We developed a genotyping method using sequences of approximately 800 bp in the E6/E7 region obtained by PacBio single molecule real-time sequencing (SMRT) and evaluated its performance against MY09-11 L1 sequencing and after the APTIMA HPV genotyping assay. The levels of concordance of PacBio E6/E7 SMRT sequencing with MY09-11 L1 sequencing and APTIMA HPV genotyping were 100% and 90.8%, respectively. The sensitivity of PacBio E6/EA7 SMRT was slightly greater than that of L1 sequencing and, as expected, lower than that of HR-HPV tests. In the context of cervical cancer screening, PacBio E6/E7 SMRT is then best used after a positive HPV test. PacBio E6/E7 SMRT genotyping is an attractive alternative for HR and LR-HPV genotyping of clinical samples. PacBio SMRT sequencing provides unbiased genotyping and can detect multiple HPV infections and haplotypes within a genotype.

Pseudoexon activation by deep intronic variation in GNE myopathy with thrombocytopenia.

GNE myopathy is a rare autosomal recessive disorder caused by pathogenic variants in the GNE gene, which is essential for the sialic acid biosynthesis pathway. Although over 300 GNE variants have been reported, some patients remain undiagnosed with monoallelic pathogenic variants. This study aims to analyze the entire GNE genomic region to identify novel pathogenic variants. Patients with clinically compatible GNE myopathy and monoallelic pathogenic variants in the GNE gene were enrolled. The other GNE pathogenic variant was verified using comprehensive methods including exon 2 quantitative polymerase chain reaction and nanopore long-read single-molecule sequencing (LRS). A deep intronic GNE variant, c.862+870C>T, was identified in nine patients from eight unrelated families. This variant generates a cryptic splice site, resulting in the activation of a novel pseudoexon between exons 5 and 6. It results in the insertion of an extra 146 nucleotides into the messengerRNA (mRNA), which is predicted to result in a truncated humanGNE1(hGNE1) protein. Peanut agglutinin（PNA） lectin staining of muscle tissues showed reduced sialylation of mucin O-glycans on sarcolemmal glycoproteins. Notably, a third of patients with the c.862+870C>T variant exhibited thrombocytopenia. A common core haplotype harboring the deep intronic GNE variant was found in all these patients. The transcript with pseudoexon activation potentially affects sialic acid biosynthesis via nonsense-mediated mRNA decay, or resulting in a truncated hGNE1 protein, which interferes with normal enzyme function. LRS is expected to be more frequently incorporated in genetic analysis given its efficacy in detecting hard-to-find pathogenic variants.

Genome Sequence of Fusarium oxysporum f. sp. malus HS2 Causing Apple Replant Disease

Apple replant disease (ARD) has been observed to cause 50% mortality in fields where apples are replanted, which has seriously limited the sustainable development of China's apple industry. Our previous study revealed that Fusarium oxysporum is one of the causal agents of ARD; however, its pathogenic mechanism is still unclear. In this study, we obtained the draft genome of Fusarium oxysporum f. sp. malus HS2 using a combination of long-read single-molecule real-time sequencing technology (Pacific Biosciences [PacBio]) and MGISEQ sequencing. The total length of the assembled was 56,825,273 bp, and the GC content was 47.72%. The genome provides a valuable resource for facilitating our understanding the response of F. oxysporum to apple. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Novel missense mutation c.797T>C (p.Met266Thr) gives rise to the rare B(A) phenotype in a Chinese family.

B(A) phenotype is usually formed by nucleotide mutations in the ABO*B.01 allele, with their products exhibiting glycosyltransferases (GTs) A and B overlapping functionality. We herein report a B(A) allele found in a Chinese family. The entire ABO genes of the probands, including flanking regulatory regions, were sequenced through PacBio third-generation long-read single-molecule real-time sequencing. 3D molecular models of the wild-type and mutant GTB were generated using the DynaMut web server. The effect of the mutation on the enzyme function was predicted by PROVEAN and PolyPhen2. The predictions of stability changes were performed using DynaMut and SNPeffect. Based on serological and sequencing features, we concluded the two probands as possible cases of the B(A) phenotype. Crystallization analysis showed that Thr266 substitution does not disrupt the hydrogen bonds. However, some changes in interatomic contacts, such as loss of ionic interactions and hydrophobic contacts, and addition of weak hydrogen bonds, may have affected protein stability to some extent. This mutation was predicted to have a benign effect on enzyme function and slightly reduce protein stability. The probands had the same novel B(A) allele with a c.797T>C (p.Met266Thr) mutation on the ABO*B.01 backbone.

Structural variants involved in high-altitude adaptation detected using single-molecule long-read sequencing

Structural variants (SVs), accounting for a larger fraction of the genome than SNPs/InDels, are an important pool of genetic variation, enabling environmental adaptations. Here, we perform long-read sequencing data of 320 Tibetan and Han samples and show that SVs are highly involved in high-altitude adaptation. We expand the landscape of global SVs, apply robust models of selection and population differentiation combining SVs, SNPs and InDels, and use epigenomic analyses to predict enhancers, target genes and biological functions. We reveal diverse Tibetan-specific SVs affecting the regulatory circuitry of biological functions, including the hypoxia response, energy metabolism and pulmonary function. We find a Tibetan-specific deletion disrupts a super-enhancer and downregulates EPAS1 using enhancer reporter, cellular knock-out and DNA pull-down assays. Our study expands the global SV landscape, reveals the role of gene-regulatory circuitry rewiring in human adaptation, and illustrates the diverse functional roles of SVs in human biology.

Identification of Hb Lepore, Hb anti-Lepore, and α-globin gene triplications by long-read single-molecule real-time sequencing.

Hemoglobin (Hb) Lepore and Hb anti-Lepore are infrequent fusion gene variants that result from nonhomologous crossovers during meiosis. Conventional molecular testing methods may face challenges in identifying these variants. During Hb analysis using capillary electrophoresis, we encountered 6 cases with unusual Hb variants. Our aim was to identify the alterations in their globin genes. Gap-polymerase chain reaction (PCR), reverse dot-blot assay (RDB), Sanger sequencing, multiplex ligation-dependent probe amplification (MLPA), and long-read single-molecule real-time (SMRT) sequencing were used to confirm the presence of globin gene alterations. The routine thalassemia gene test kit using the gap-PCR and RDB techniques did not detect common gene variations. Direct sequencing failed to identify any known or unknown globin gene alterations. The MLPA analysis, however, revealed the possible presence of α-globin gene triplications as well as 2 types of fusion gene alterations. Further analysis using long-read SMRT sequencing accurately identified 3 rare gene variations: αααanti-3.7, Hb Lepore-Boston-Washington, and Hb anti-Lepore P-India. Conventional methods may overlook rare thalassemias or Hb variants. Long-read SMRT sequencing has the potential to identify breakpoints in fusion genes, demonstrating that it is a promising technique for detecting rare thalassemias.

Genome Sequencing of Strain Vm263, a Hypotoxic Strain of Valsa mali Causing Apple Valsa Canker

Apple valsa canker, caused by Valsa mali, has seriously endangered the sustainable development of the apple industry in China. Our previous study revealed that low-temperature induction increased the pathogenicity of strain Vm263, which is a hypotoxic strain of V. mali causing apple valsa canker. However, the regulatory mechanism of this phenomenon was still unclear. In this study, whole genome sequencing of strain Vm263 was performed using a combination of long-read single-molecule real-time sequencing technology (Pacific Biosciences) and MGISEQ sequencing. The genome size of strain Vm263 was 42.83 Mb and the assembly N50 was 5.09 Mb. The genome provides a valuable resource for facilitating our understanding of the response of V. mali to low temperature. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Viral Diversity in Benthic Abyssal Ecosystems: Ecological and Methodological Considerations.

Viruses are the most abundant 'biological entities' in the world's oceans. However, technical and methodological constraints limit our understanding of their diversity, particularly in benthic abyssal ecosystems (>4000 m depth). To verify advantages and limitations of analyzing virome DNA subjected either to random amplification or unamplified, we applied shotgun sequencing-by-synthesis to two sample pairs obtained from benthic abyssal sites located in the North-eastern Atlantic Ocean at ca. 4700 m depth. One amplified DNA sample was also subjected to single-molecule long-read sequencing for comparative purposes. Overall, we identified 24,828 viral Operational Taxonomic Units (vOTUs), belonging to 22 viral families. Viral reads were more abundant in the amplified DNA samples (38.5-49.9%) compared to the unamplified ones (4.4-5.8%), with the latter showing a greater viral diversity and 11-16% of dsDNA viruses almost undetectable in the amplified samples. From a procedural point of view, the viromes obtained by direct sequencing (without amplification step) provided a broader overview of both ss and dsDNA viral diversity. Nevertheless, our results suggest that the contextual use of random amplification of the same sample and long-read technology can improve the assessment of viral assemblages by reducing off-target reads.

Comprehensive Profiling of Alternative Splicing and Alternative Polyadenylation during Fruit Ripening in Watermelon (Citrullus lanatus).

Fruit ripening is a highly complicated process that is accompanied by the formation of fruit quality. In recent years, a series of studies have demonstrated post-transcriptional control play important roles in fruit ripening and fruit quality formation. Till now, the post-transcriptional mechanisms for watermelon fruit ripening have not been comprehensively studied. In this study, we conducted PacBio single-molecule long-read sequencing to identify genome-wide alternative splicing (AS), alternative polyadenylation (APA) and long non-coding RNAs (lncRNAs) in watermelon fruit. In total, 6,921,295 error-corrected and mapped full-length non-chimeric (FLNC) reads were obtained. Notably, more than 42,285 distinct splicing isoforms were derived from 5,891,183 intron-containing full-length FLNC reads, including a large number of AS events associated with fruit ripening. In addition, we characterized 21,506 polyadenylation sites from 11,611 genes, 8703 of which have APA sites. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that fructose and mannose metabolism, starch and sucrose metabolism and carotenoid biosynthesis were both enriched in genes undergoing AS and APA. These results suggest that post-transcriptional regulation might potentially have a key role in regulation of fruit ripening in watermelon. Taken together, our comprehensive PacBio long-read sequencing results offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of watermelon fruit ripening.

Establishment of tongue microbiota by 18 months of age and determinants of its microbial profile.

Understanding the development of oral microbiota early in life and the factors that influence it is important for preventing the establishment of dysbiotic oral microbiota later in life. This study demonstrates that the tongue microbiota undergoes early development from 4 to 18 months of age and converges into two types of microbiota showing indications of adult characteristics, with either S. salivarius or Neisseria-dominance. Interestingly, their divergence was strongly determined by their weaning status and the dietary frequencies of sweetened beverages, snacks, and fruits, suggesting that dietary habits during this period might influence the establishment of the oral microbiota. These findings may contribute to the development of novel preventive strategies against oral microbiota-related diseases.

Salivary microbiome profiles for different clinical phenotypes of pituitary adenomas by single-molecular long-read sequencing.

Pituitary adenomas (PAs) are common benign brain tumors. Although associations between the gastrointestinal microbiome and PA have been reported previously, studies on the distributions of salivary microbial species in PA patients and among their different clinical phenotypes are currently lacking. In this study, saliva samples from 42 patients and 20 healthy individuals were selected for third-generation sequencing. The PA group included four clinical phenotypes: adrenocorticotropic hormone-secreting PA (n = 6), growth hormone-secreting PA (n = 9), prolactin-secreting PA (n = 18), and nonfunctioning PA (n = 9). All samples were sequenced, and the data were clustered and de-chimerized to obtain information regarding the abundance of operational taxonomic units. We found that the species distributions in the saliva of PA patients were more abundant than those of healthy individuals. A total of 82 genera were identified across all samples, of which 14 and 17 genera were more abundant in the saliva samples of patients with PA and healthy individuals, respectively. In the phenotypic functional prediction, the phenotypes of anaerobic and Gram-positive organisms were more commonly seen in patients with PA than in healthy individuals. The bioinformatics prediction indicated that multiple metabolic pathways were involved in the pathogenesis of PA. In conclusion, this study highlighted the associations of salivary microbiome profiles with PA, which may improve the existing understanding of the pathogenesis of PA and provide diagnostic and therapeutic targets for PA.IMPORTANCEThe gut and salivary microbiomes have been widely reported to be significantly associated with a number of neurological disorders. The stability of the microbiome in the oral cavity makes it a potentially ideal sample that can be conveniently obtained for the investigation of microbiome-based pathogenesis in diseases. In the present study, we used a single-molecule long-read sequencing technique to study the distribution of the salivary microbiota in patients with pituitary adenoma (PA) and healthy individuals, as well as among four clinical phenotypes of PA. We found that the diversity of salivary microbes was more abundant in PA patients than in healthy individuals. We also observed some unique genera in different PA phenotypes. The bioinformatics-based functional predictions identified potential links between microbes and different clinical phenotypes of PA. This study improves the existing understanding of the pathogenesis of PA and may provide diagnostic and therapeutic targets for PA.

Long- and short-read RNA sequencing from five reproductive organs of boar

The production of semen in boars involves multiple reproductive glands, including the testis (Tes), epididymis (Epi), vesicular gland (VG), prostate gland (PG), and bulbourethral gland (BG). However, previous studies on boar reproduction primarily focused on the testis, with little attention paid to the other glands. Here, we integrated single-molecule long-read sequencing with short-read sequencing to characterize the RNA landscape from five glands of Banna mini-pig inbred line (BMI) and Diannan small-ear pigs (DSE). We identified 110,996 full-length isoforms from 22,298 genes, and classified the alternative splicing (AS) events in these five glands. Transcriptome-wide variation analysis indicated that the number of single nucleotide polymorphisms (SNPs) in five tissues of BMI was significantly lower than that in the non-inbred pig, DSE, revealing the effect of inbreeding on BMI. Additionally, we performed small-RNA sequencing and identified 299 novel miRNAs across all glands. Overall, our findings provide a comprehensive overview of the RNA landscape within these five glands, paving the path for future investigations on reproductive biology and the impact of inbreeding on pig transcriptome.

The sequencing of the key genes and end products in the TLR4 signaling pathway from the kidney of Rana dybowskii exposed to Aeromonas hydrophila.

Infectious diseases caused by Aeromonas hydrophila (AH) have reduced the populations of Rana dybowskii). However, little is known about the immune response of R. dybowskii against AH infections. The toll-like receptor (TLR) signaling pathway has been identified as a critical component in innate immunity, responsible for identifying pathogen-associated molecular patterns in pathogens. Our study used the next-generation sequencing technique and single-molecule long-read sequencing to determine the structures of transcript isoforms and functions of genes in the kidneys of R. dybowskii, as well as identify and validate the related genes in the TLR4 signaling pathway. In total, 628,774 reads of inserts were identified, including 300,053 full-length non-chimeric reads and 233,592 non-full-length reads. Among the transcriptome sequences, 124 genes were identified as homologs of known genes in the TLR4 pathway especially inflammatory cytokines and receptors. Our findings shed light on the structures and functions of R. dybowskii genes exposed to AH and confirm the presence of both MyD88-dependent and independent pathways in R. dybowskii. Our work reveals how various functional proteins in amphibians at the initial stage of immune response are activated and complete their corresponding functions in a short time.

Rapid detection of myeloid neoplasm fusions using single-molecule long-read sequencing.

Recurrent gene fusions are common drivers of disease pathophysiology in leukemias. Identifying these structural variants helps stratify disease by risk and assists with therapy choice. Precise molecular diagnosis in low-and-middle-income countries (LMIC) is challenging given the complexity of assays, trained technical support, and the availability of reliable electricity. Current fusion detection methods require a long turnaround time (7-10 days) or advance knowledge of the genes involved in the fusions. Recent technology developments have made sequencing possible without a sophisticated molecular laboratory, potentially making molecular diagnosis accessible to remote areas and low-income settings. We describe a long-read sequencing DNA assay designed with CRISPR guides to select and enrich for recurrent leukemia fusion genes, that does not need a priori knowledge of the abnormality present. By applying rapid sequencing technology based on nanopores, we sequenced long pieces of genomic DNA and successfully detected fusion genes in cell lines and primary specimens (e.g., BCR::ABL1, PML::RARA, CBFB::MYH11, KMT2A::AFF1) using cloud-based bioinformatics workflows with novel custom fusion finder software. We detected fusion genes in 100% of cell lines with the expected breakpoints and confirmed the presence or absence of a recurrent fusion gene in 12 of 14 patient cases. With our optimized assay and cloud-based bioinformatics workflow, these assays and analyses could be performed in under 8 hours. The platform's portability, potential for adaptation to lower-cost devices, and integrated cloud analysis make this assay a candidate to be placed in settings like LMIC to bridge the need of bedside rapid molecular diagnostics.

First Report of a Low-Frequency Mosaic Mutation in the Hydroxymethylbilane Synthase Gene Causing Acute Intermittent Porphyria.

Acute porphyrias are a group of monogenetic inborn errors of heme biosynthesis, characterized by acute and potentially life-threatening neurovisceral attacks upon exposure to certain triggering factors. Biochemical analyses can determine the type of acute porphyria, and subsequent genetic analysis allows for the identification of pathogenic variants in the specific gene, which provides information for family counselling. In 2017, a male Swiss patient was diagnosed with an acute porphyria while suffering from an acute attack. The pattern of porphyrin metabolite excretion in urine, faeces, and plasma was typical for an acute intermittent porphyria (AIP), which is caused by inherited autosomal dominant mutations in the gene for hydroxymethylbilane synthase (HMBS), the third enzyme in the heme biosynthetic pathway. However, the measurement of HMBS enzymatic activity in the erythrocytes was within the normal range and Sanger sequencing of the HMBS gene failed to detect any pathogenic variants. To explore the molecular basis of the apparent AIP in this patient, we performed third-generation long-read single-molecule sequencing (nanopore sequencing) on a PCR product spanning the entire HMBS gene, including the intronic sequences. We identified a known pathogenic variant, c.77G>A, p.(Arg26His), in exon 3 at an allelic frequency of ~22% in the patient's blood. The absence of the pathogenic variant in the DNA of the parents and the results of additional confirmatory studies supported the presence of a de novo mosaic mutation. To our knowledge, such a mutation has not been previously described in any acute porphyria. Therefore, de novo mosaic mutations should be considered as potential causes of acute porphyrias when no pathogenic genetic variant can be identified through routine molecular diagnostics.

Characterization of complex structural variation in the CYP2D6-CYP2D7-CYP2D8 gene loci using single-molecule long-read sequencing.

Complex regions in the human genome such as repeat motifs, pseudogenes and structural (SVs) and copy number variations (CNVs) present ongoing challenges to accurate genetic analysis, particularly for short-read Next-Generation-Sequencing (NGS) technologies. One such region is the highly polymorphic CYP2D loci, containing CYP2D6, a clinically relevant pharmacogene contributing to the metabolism of >20% of common drugs, and two highly similar pseudogenes, CYP2D7 and CYP2D8. Multiple complex SVs, including CYP2D6/CYP2D7-derived hybrid genes are known to occur in different configurations and frequencies across populations and are difficult to detect and characterize accurately. This can lead to incorrect enzyme activity assignment and impact drug dosing recommendations, often disproportionally affecting underrepresented populations. To improve CYP2D6 genotyping accuracy, we developed a PCR-free CRISPR-Cas9 based enrichment method for targeted long-read sequencing that fully characterizes the entire CYP2D6-CYP2D7-CYP2D8 loci. Clinically relevant sample types, including blood, saliva, and liver tissue were sequenced, generating high coverage sets of continuous single molecule reads spanning the entire targeted region of up to 52kb, regardless of SV present (n = 9). This allowed for fully phased dissection of the entire loci structure, including breakpoints, to accurately resolve complex CYP2D6 diplotypes with a single assay. Additionally, we identified three novel CYP2D6 suballeles, and fully characterized 17 CYP2D7 and 18 CYP2D8 unique haplotypes. This method for CYP2D6 genotyping has the potential to significantly improve accurate clinical phenotyping to inform drug therapy and can be adapted to overcome testing limitations of other clinically challenging genomic regions.

PacBio single-molecule long-read sequencing provides new insights into the complexity of full-length transcripts in oriental river prawn, macrobrachium nipponense

BackgroundOriental river prawn (Macrobrachium nipponense) is one of the most dominant species in shrimp farming in China, which is a rich source of protein and contributes to a significant impact on the quality of human life. Thus, more complete and accurate annotation of gene models are important for the breeding research of oriental river prawn.ResultsA full-length transcriptome of oriental river prawn muscle was obtained using the PacBio Sequel platform. Then, 37.99 Gb of subreads were sequenced, including 584,498 circular consensus sequences, among which 512,216 were full length non-chimeric sequences. After Illumina-based correction of long PacBio reads, 6,599 error-corrected isoforms were identified. Transcriptome structural analysis revealed 2,263 and 2,555 alternative splicing (AS) events and alternative polyadenylation (APA) sites, respectively. In total, 620 novel genes (NGs), 197 putative transcription factors (TFs), and 291 novel long non-coding RNAs (lncRNAs) were identified.ConclusionsIn summary, this study offers novel insights into the transcriptome complexity and diversity of this prawn species, and provides valuable information for understanding the genomic structure and improving the draft genome annotation of oriental river prawn.

SARS-CoV-2 Co-Infections and Recombinations Identified by Long-Read Single-Molecule Real-Time Sequencing.

Co-infection with at least 2 strains of virus is the prerequisite for recombination, one of the means of genetic diversification. Little is known about the prevalence of these events in SARS-CoV-2, partly because it is difficult to detect them. We used long-read PacBio single-molecule real-time (SMRT) sequencing technology to sequence whole genomes and targeted regions for haplotyping. We identified 17 co-infections with SARS-CoV-2 strains belonging to different clades in 6829 samples sequenced between January and October, 2022 (prevalence 0.25%). There were 3 Delta/Omicron co-infections and 14 Omicron/Omicron co-infections (4 cases of 21K/21L, 1 case of 21L/22A, 2 cases of 21L/22B, 4 cases of 22A/22B, 2 cases of 22B/22C and 1 case of 22B/22E). Four of these patients (24%) also harbored recombinant minor haplotypes, including one with a recombinant virus that was selected in the viral quasispecies over the course of his chronic infection. While co-infections remain rare among SARS-CoV-2-infected individuals, long-read SMRT sequencing is a useful tool for detecting them as well as recombinant events, providing the basis for assessing their clinical impact, and a precise indicator of epidemic evolution. IMPORTANCE SARS-CoV-2 variants have been responsible for the successive waves of infection over the 3 years of pandemic. While co-infection followed by recombination is one driver of virus evolution, there have been few reports of co-infections, mainly between Delta and Omicron variants or between the first 2 Omicron variants 21K_BA.1 and 21L_BA.2. The 17 co-infections we detected during 2022 included cases with the recent clades of Omicron 22A, 22B, 22C, and 22E; 24% harbored recombinant variants. This study shows that long-read SMRT sequencing is well suited to SARS-CoV-2 genomic surveillance.

High-Resolution Detection of Translocation of Oral Bacteria to the Gut

Ectopic enrichment of oral microbes in the gut is a notable alteration in gut microbial balance. These microbes are likely delivered from the oral cavity with saliva and food; however, evidence of oral–gut microbial transmission is insufficient and needs further investigation. In this observational study, we examined 144 pairs of saliva and stool samples collected from community-dwelling adults to verify the oral–gut microbial link and identify the relevant influencing factors on the increased abundance of oral microbes within the gut. The bacterial composition of each sample was determined using PacBio single-molecule long-read sequencing of the full-length 16S ribosomal RNA gene and amplicon sequence variant (ASV) analysis. Although the bacterial compositions of salivary and gut microbiota were distinctly different, at least 1 ASV was shared between salivary and gut microbiota in 72.9% of subjects. Shared ASVs accounted for 0.0% to 63.1% (median 0.14%) of the gut microbiota in each subject and frequently included abundant Streptococcus salivarius and Streptococcus parasanguinis. Their total relative abundance in the gut was significantly higher in older subjects or those with dental plaque accumulation. The gut microbiota with ≥5% of shared ASVs displayed a higher abundance of Streptococcus, Lactobacillus, and Klebsiella and a lower abundance of Faecalibacterium, Blautia, Megamonas, and Parabacteroides. Our study presents evidence for the translocation of oral bacteria to the gut in community-dwelling adults and suggests that aging and dental plaque accumulation contribute to an increased abundance of oral microbes in the gut, which might be relevant to the compositional shift in the gut commensals.

Comparative genomics of the Leukocyte Receptor Complex in carnivores.

The mammalian Leukocyte Receptor Complex (LRC) chromosomal region may contain gene families for the killer cell immunoglobulin-like receptor (KIR) and/or leukocyte immunoglobulin-like receptor (LILR) collections as well as various framing genes. This complex region is well described in humans, mice, and some domestic animals. Although single KIR genes are known in some Carnivora, their complements of LILR genes remain largely unknown due to obstacles in the assembly of regions of high homology in short-read based genomes. As part of the analysis of felid immunogenomes, this study focuses on the search for LRC genes in reference genomes and the annotation of LILR genes in Felidae. Chromosome-level genomes based on single-molecule long-read sequencing were preferentially sought and compared to representatives of the Carnivora. Seven putatively functional LILR genes were found across the Felidae and in the Californian sea lion, four to five genes in Canidae, and four to nine genes in Mustelidae. They form two lineages, as seen in the Bovidae. The ratio of functional genes for activating LILRs to inhibitory LILRs is slightly in favor of inhibitory genes in the Felidae and the Canidae; the reverse is seen in the Californian sea lion. This ratio is even in all of the Mustelidae except the Eurasian otter, which has a predominance of activating LILRs. Various numbers of LILR pseudogenes were identified. The structure of the LRC is rather conservative in felids and the other Carnivora studied. The LILR sub-region is conserved within the Felidae and has slight differences in the Canidae, but it has taken various evolutionary paths in the Mustelidae. Overall, the process of pseudogenization of LILR genes seems to be more frequent for activating receptors. Phylogenetic analysis found no direct orthologues across the Carnivora which corroborate the rapid evolution of LILRs seen in mammals.