Genome-wide Variation Research Articles

Identification of genetic fingerprint of type I interferon therapy in visceral metastases of melanoma

Malignant melanoma is a difficult-to-treat skin cancer with increasing incidence worldwide. Although type-I interferon (IFN) is no longer part of guidelines, several melanoma patients are treated with type-I interferon (IFN) at some point of the disease, potentially affecting its genetic progression. We run genome-wide copy number variation (CNV) analysis on previously type-I IFN-treated (n = 17) and control (n = 11) visceral metastases of melanoma patients. Results were completed with data from the TCGA and MM500 databases. We identified metastasis- and brain metastasis-specific gene signatures mostly affected by CN gains. Some cases were genetically resistant to IFN showing characteristic gene alterations (e.g. ABCA4 or ZEB2 gain and alterations of DNA repair genes). Analysis of a previously identified type-I IFN resistance gene set indicates that only a proportion of these genes was exclusive for the IFN-treated metastases reflecting a possible selective genomic pressure of endogenous IFNs during progression. Our data suggest that previous type-I IFN treatment and/or endogenous IFN production by immune response affect genomic progression of melanoma which may have clinical relevance, potentially influence immune checkpoint regulation in the tumor microenvironment.

Near complete assembly of Pyricularia penniseti infecting Cenchrus grass identified its eight core chromosomes

Fungi from the Pyricularia genus cause blast disease in many economically important crops and grasses, such as wheat, rice, and Cenchrus grass JUJUNCAO. Structure variation associated with the gain and loss of effectors contributes largely to the adaptive evolution of this fungus towards diverse host plants. A telomere-to-telomere genome assembly would facilitate the identification of genome-wide structural variations through comparative genomics. Here, we report a telomere-to-telomere, near-complete genome assembly of a Pyricularia penniseti isolate JC-1 infecting JUJUNCAO. The assembly consists of eight core chromosomes and two supernumerary chromosomes, named mini1 and mini2, spanning 42.1 Mb. We annotated 12,156 protein-coding genes and identified 4.54% of the genome as repetitive sequences. The two supernumerary chromosomes contained fewer genes and more repetitive sequences than the core chromosomes. Our genome and results provide valuable resources for the future study in genome evolution, structure variation and host adaptation of the Pyricularia fungus.

Genome-Wide Genetic Diversity and Population Structure of Charybdis feriata (Crustacea, Decapoda, and Portunidae) Along the Southeast Coast of China Inferred from Genotyping-by-Sequencing (GBS) Approach

Genome-Wide Genetic Diversity and Population Structure of Charybdis feriata (Crustacea, Decapoda, and Portunidae) Along the Southeast Coast of China Inferred from Genotyping-by-Sequencing (GBS) Approach

A novel functional IKBKE variant activating NFAT in a patient with polyarthritis and a remittent fever

BackgroundIKBKE is a negative regulator of T cell activation and one of the key activators of type I interferon (IFN) and NFκB signaling via non-classical pathways. The upstream single nucleotide polymorphism of IKBKE (rs2297550-G) is a genome-wide association study risk variant of systemic lupus erythematosus, and is associated with decreased IKBKE expression in T cells by expression quantitative trait locus analysis.Case presentationA 48-year-old female had a remittent fever, arthritis, and oral ulcers for 20 years. She had a poor response to corticosteroids or disease-modifying antirheumatic drugs, including the tumor necrosis factor-α antagonist, etanercept, and the anti-interleukin-6 receptor antibody, tocilizumab.MethodShe participated in the Initiative on Rare and Undiagnosed Disease (IRUD), and whole-exome sequencing (WES) was performed. Functional analyses were conducted by transfecting the identified variants into reporter cells to assess the activation of NFAT and NFκB signaling. Additionally, peripheral blood RNA- sequencing (RNA-seq) data were compared with those from healthy individuals to evaluate the gene expression profiles of immune cells.ResultWES identified a novel heterozygous c.1877G>A, p(Cys626Tyr) variant in IKBKE. Functional analysis indicated that this variant led to increased activity of NFAT (p = 0.015) and decreased activity of NFκB and type I IFN (p = 0.00068 and 0.00044, respectively). The patient had a remarkably low proportion of Naïve CD4 T cells. RNA-seq of peripheral blood immune cell subsets revealed significant differences in gene expression, especially in T cells.ConclusionA novel functional heterozygous variant in IKBKE is described in a patient with a remittent fever and arthritis. The data suggest that IKBKE is an important negative regulator of inflammation, particularly in T cells, and this IKBKE variant might be the underlying cause of a novel autoinflammatory pathology.

Contrasting Genomic Diversity and Inbreeding Levels Among Two Closely Related Falcon Species With Overlapping Geographic Distributions.

Genomic resources are valuable to examine historical demographic patterns and their effects to better inform management and conservation of threatened species. We evaluated population trends and genome-wide variation in the near-threatened Orange-breasted Falcon (Falco deiroleucus) and its more common sister species, the Bat Falcon (F. rufigularis), to explore how the two species differ in genomic diversity as influenced by their contrasting long-term demographic histories. We generated and aligned whole genome resequencing data for 12 Orange-breasted Falcons and 9 Bat Falcons to an annotated Gyrfalcon (F. rusticolus) reference genome that retained approximately 22.4 million biallelic autosomal SNPs (chromosomes 1-22). Our analyses indicated much lower genomic diversity in Orange-breasted Falcons compared to Bat Falcons. All sampled Orange-breasted Falcons were significantly more inbred than the sampled Bat Falcons, with values similar to those observed in island-mainland species comparisons. The distribution of runs of homozygosity showed variation suggesting long-term low population size and the possibility of bottlenecks in Orange-breasted Falcons contrasting with consistently larger populations in Bat Falcons. Analysis of genetic load suggests that Orange-breasted Falcons are less likely to experience inbreeding depression than Bat Falcons due to reduced inbreeding load but are at elevated risk from fixation of deleterious gene variants and perhaps a reduced adaptive potential. These genomic analyses highlight differences in the historical demography of two closely related species that have influenced their current genomic diversity and should result in differing strategies for their continued conservation.

Genome-wide diversity, linkage disequilibrium, and admixture in the main Colombian Creole pig breeds

Colombian Creole pigs have adapted to tropical conditions for over 500 years. They have been modified by natural and artificial selection in different regions. At present, the diversity and current introgression status are unknown. The objective was to estimate the genomic diversity, linkage disequilibrium, population structure, and admixture of four Colombian pig breeds and their relationship with other breeds worldwide. Three Colombian pig breeds (SPE-San Pedreño, 11 samples; ZUN-Zungo, 11 samples; CM-Casco de Mula, ten samples) from the conservation nucleus and one biotype not recognized as a breed (CCH-Criollo Chocoano, seven samples) were genotyped using the Illumina GGP-Porcine80K chip. Open-access data from seven international breeds were also included. Colombian Creole pigs showed moderate genetic differentiation (FST 0.14) globally, but several groups of animals separated, suggesting local clustering due to geographical isolation or different founding effects. Colombian Creole pigs showed breed imprinting and specific grouping in all analyses except for CCH, which, like the Ecuadorian Creole, was a cluster of admixtures. The Colombian Creole pigs revealed a significant relationship with the Iberian pig and some other breeds to varying degrees. However, good maintenance of the conservation nucleus was evidenced. Potential adaptive genes, mainly related to immunological functions, were found, according to FST and pcadapt analyses. This study provides a foundation and scientific data for policy decisions on zoogenetic resources.

Genome-wide association study and genotypic variation for the major tocopherol content in rice grain.

Rice tocopherols, vitamin E compounds with antioxidant activity, play essential roles in human health. Even though the key genes involved in vitamin E biosynthetic pathways have been identified in plants, the genetic architecture of vitamin E content in rice grain remains unclear. A genome-wide association study (GWAS) on 179 genotypically diverse rice accessions with 34,323 SNP markers was conducted to detect QTLs that define total and α- tocopherol contents in rice grains. Total and α-tocopherol contents had a strong positive correlation and varied greatly across the accessions, ranging from 0.230-31.76 and 0.011-30.83 (μg/g), respectively. A total of 13 QTLs were identified, which were spread across five of the rice chromosomes. Among the 13 QTLs, 11 were considered major with phenotypic variation explained (PVE) greater than 10%. Twelve transcription factor (TF) genes, one microprotein (miP), and a transposon were found to be associated with the QTLs with putative roles in controlling tocopherol contents. Moreover, intracellular transport proteins, ABC transporters, nonaspanins, and SNARE, were identified as associated genes on chromosomes 1 and 8. In the vicinity of seven QTLs, protein kinases were identified as key signaling factors. Haplotype analysis revealed the QTLs qAlph1.1, qTot1.1, qAlph2.1, qAlph6.1, qTot6.1, and qTot8.3 to have significant haplogroups. Quantitative RT-PCR validated the expression direction and magnitude of WRKY39 (Os02g0265200), PIP5Ks (Os08g0450800), and MADS59 (Os06g0347700) in defining the major tocopherol contents. This study provides insights for ongoing biofortification efforts to breed and/or engineer vitamin E and antioxidant levels in rice and other cereals.

Odontogenic Myxomas Harbor Recurrent Copy Number Alterations and a Distinct Methylation Signature.

Odontogenic myxoma is a rare, benign, and locally aggressive tumor that develops in the tooth-bearing areas of the jaw. The molecular mechanisms underlying odontogenic myxomas are unknown and no diagnostic markers are available to date. The aim of this study was to analyze DNA methylation and copy number variations in odontogenic myxomas to identify new molecular signatures for diagnostic decision-making. We collected a cohort of 16 odontogenic myxomas from 2006 to 2021 located in the mandible (n = 10) and maxilla (n = 6) with available formalin-fixed paraffin-embedded or fresh frozen tumor tissue from a biopsy or resection material. Genome-wide DNA methylation and copy number variation data were generated from 12 odontogenic myxomas using the Illumina Infinium Methylation EPIC array, interrogating >850,000 CpG sites. Unsupervised clustering and dimensionality reduction (Uniform Manifold Approximation and Projection) revealed that odontogenic myxomas formed a distinct DNA methylation class. Copy number profiling showed recurrent whole-chromosome gains (trisomies) of chromosomes 5, 8, and 20 in all cases, and of chromosomes 10, 12, and 17 in all except one case. In conclusion, odontogenic myxomas harbor recurrent copy number patterns and a distinct DNA methylation profile, which can be used as an additional diagnostic tool in the appropriate clinical and radiologic context. Further research is needed to explain the genetic mechanisms caused by these alterations that drive these locally aggressive neoplasms.

Genome-Wide Association Studies of Hair Whorl in Pigs.

In pigs, a hair whorl refers to hairs that form a ring of growth around the direction of the hair follicle at the dorsal hip. In China, a hair whorl is considered a negative trait that affects marketing, and no studies have been conducted to demonstrate whether hair whorl affects pig performance and provide an explanation for its genetic basis. Performance-measured traits and slaughter-measured traits of hair whorl and non-hair whorl pigs were differentially analyzed, followed by genome-wide association analysis (GWAS) and copy number variation (CNV) methods to investigate the genetic basis of hair whorl in pigs. Differential analysis of 2625 pigs (171 hair whorl and 2454 non-hair whorl) for performance measures showed that hair whorl and non-hair whorl pigs differed significantly (p < 0.05) in traits such as live births, total litter size, and healthy litter size (p < 0.05), while differential analysis of carcass and meat quality traits showed a significant difference only in the 45 min pH (p = 0.0265). GWAS identified 4 SNP loci significantly associated with the hair whorl trait, 2 of which reached genome-significant levels, and 23 candidate genes were obtained by annotation with the Ensembl database. KEGG and GO enrichment analyses showed that these genes were mainly enriched in the ErbB signaling, endothelial apoptosis regulation, and cell proliferation pathways. In addition, CNV analysis identified 652 differential genes between hair whorl and non-hair whorl pigs, which were mainly involved in the signal transduction, transcription factor activity, and nuclear and cytoplasmic-related pathways. The candidate genes and copy number variation differences identified in this study provide a new theoretical basis for pig breeding efforts.

An integrated strain-level analytic pipeline utilizing longitudinal metagenomic data.

The advancement in DNA-sequencing technology has enabled the high-resolution identification of microorganisms in microbial communities. Since different microbial strains within species may contain extreme phenotypic variability (e.g., nutrition metabolism, antibiotic resistance, and pathogen virulence), investigating within-species variations holds great scientific promise in understanding the underlying mechanism of microbial biological processes. To fully utilize the shared genomic variants across longitudinal metagenomics samples collected in microbiome studies, we develop an integrated analytic pipeline (LongStrain) for longitudinal metagenomics data. It concurrently leverages the information on proportions of mapped reads for individual strains and genome-wide SNVs to enhance the efficiency and accuracy of strain identification. Our method helps to understand strains' dynamic changes and their association with genome-wide variants. Given the fast-growing longitudinal studies of microbial communities, LongStrain which streamlines analyses of large-scale raw sequencing data should be of great value in microbiome research communities.

Evolution of unexpected diversity in a putative mating type locus and its correlation with genome variability reveals likely asexuality in the model mycorrhizal fungus Rhizophagus irregularis

BackgroundArbuscular mycorrhizal fungi (AMF) form mutualistic partnerships with approximately 80% of plant species. AMF, and their diversity, play a fundamental role in plant growth, driving plant diversity, and global carbon cycles. Knowing whether AMF are sexual or asexual has fundamental consequences for how they can be used in agricultural applications. Evidence for and against sexuality in the model AMF, Rhizophagus irregularis, has been proposed. The discovery of a putative mating-type locus (MAT locus) in R. irregularis, and the previously suggested recombination among nuclei of a dikaryon R. irregularis isolate, potentially suggested sexuality. Unless undergoing frequent sexual reproduction, evolution of MAT-locus diversity is expected to be very low. Additionally, in sexual species, MAT-locus evolution is decoupled from the evolution of arbitrary genome-wide loci.ResultsWe studied MAT-locus diversity of R. irregularis. This was then compared to diversification in a phosphate transporter gene (PTG), that is not involved in sex, and to genome-wide divergence, defined by 47,378 single nucleotide polymorphisms. Strikingly, we found unexpectedly high MAT-locus diversity indicating that either it is not involved in sex, or that AMF are highly active in sex. However, a strongly congruent evolutionary history of the MAT-locus, PTG and genome-wide arbitrary loci allows us to reject both the hypothesis that the MAT-locus is involved in mating and that the R. irregularis lineage is sexual.ConclusionOur finding shapes the approach to developing more effective AMF strains and is highly informative as it suggests that introduced strains applied in agriculture will not exchange DNA with native populations.

Ethnic and region-specific genetic risk variants of stroke and its comorbid conditions can define the variations in the burden of stroke and its phenotypic traits.

Burden of stroke differs by region, which could be attributed to differences in comorbid conditions and ethnicity. Genomewide variation acts as a proxy marker for ethnicity, and comorbid conditions. We present an integrated approach to understand this variation by considering prevalence and mortality rates of stroke and its comorbid risk for 204 countries from 2009 to 2019, and Genome-wide association studies (GWAS) risk variant for all these conditions. Global and regional trend analysis of rates using linear regression, correlation, and proportion analysis, signifies ethnogeographic differences. Interestingly, the comorbid conditions that act as risk drivers for stroke differed by regions, with more of metabolic risk in America and Europe, in contrast to high systolic blood pressure in Asian and African regions. GWAS risk loci of stroke and its comorbid conditions indicate distinct population stratification for each of these conditions, signifying for population-specific risk. Unique and shared genetic risk variants for stroke, and its comorbid and followed up with ethnic-specific variation can help in determining regional risk drivers for stroke. Unique ethnic-specific risk variants and their distinct patterns of linkage disequilibrium further uncover the drivers for phenotypic variation. Therefore, identifying population- and comorbidity-specific risk variants might help in defining the threshold for risk, and aid in developing population-specific prevention strategies for stroke.

Commonly used genomic arrays may lose information due to imperfect coverage of discovered variants for autism spectrum disorder

BackgroundCommon genetic variation has been shown to account for a large proportion of ASD heritability. Polygenic scores generated for autism spectrum disorder (ASD-PGS) using the most recent discovery data, however, explain less variance than expected, despite reporting significant associations with ASD and other ASD-related traits. Here, we investigate the extent to which information loss on the target study genome-wide microarray weakens the predictive power of the ASD-PGS.MethodsWe studied genotype data from three cohorts of individuals with high familial liability for ASD: The Early Autism Risk Longitudinal Investigation (EARLI), Markers of Autism Risk in Babies-Learning Early Signs (MARBLES), and the Infant Brain Imaging Study (IBIS), and one population-based sample, Study to Explore Early Development Phase I (SEED I). Individuals were genotyped on different microarrays ranging from 1 to 5 million sites. Coverage of the top 88 genome-wide suggestive variants implicated in the discovery was evaluated in all four studies before quality control (QC), after QC, and after imputation. We then created a novel method to assess coverage on the resulting ASD-PGS by correlating a PGS informed by a comprehensive list of variants to a PGS informed with only the available variants.ResultsPrior to imputations, None of the four cohorts directly or indirectly covered all 88 variants among the measured genotype data. After imputation, the two cohorts genotyped on 5-million arrays reached full coverage. Analysis of our novel metric showed generally high genome-wide coverage across all four studies, but a greater number of SNPs informing the ASD-PGS did not result in improved coverage according to our metric.Limitations.The studies we analyzed contained modest sample sizes. Our analyses included microarrays with more than 1-million sites, so smaller arrays such as Global Diversity and the PsychArray were not included. Our PGS metric for ASD is only generalizable to samples of European ancestries, though the coverage metric can be computed for traits that have sufficiently large-sized discovery findings in other ancestries.ConclusionsWe show that commonly used genotyping microarrays have incomplete coverage for common ASD variants, and imputation cannot always recover lost information. Our novel metric provides an intuitive approach to reporting information loss in PGS and an alternative to reporting the total number of SNPs included in the PGS. While applied only to ASD here, this metric can easily be used with other traits.

Influence of geographic isolation and the environment on gene flow among phenotypically diverse lizards

Lizards in the genus Anolis comprise hundreds of species that display a wide range of phenotypic variation closely related to their environment. One example is the Guadeloupean anole (Anolis marmoratus ssp.) that display extreme phenotypic variation, primarily in adult male color and pattern, with twelve described subspecies on the archipelago. Here we examine the relationship between phenotypic and genetic divergence among five subspecies on the two main islands and test the role of geographic isolation and the environment in reducing gene flow. We also examined two offshore island populations to assess the impact of complete geographic isolation on gene flow. We analyzed color phenotypes by measuring spectral reflectance and genomic diversity using SNPs. Genetic divergence was correlated with dorsolateral head and body color phenotypes, and slope and geographic distance were nearly equivalent at explaining this divergence. There was minimal genome-wide divergence at neutral loci among phenotypically disparate subspecies on the two main islands and their differentiation is consistent with a model of divergence with gene flow. Our spatial visualization of gene flow showed an impact of environmental features consistent with a hypothesis of ecologically driven divergence. Nonetheless, subspecies on the two main islands remain interconnected by substantial gene flow and their phenotypic variation is likely maintained at selection-gene flow equilibrium by divergent selection at loci associated with their color phenotypes. Greater isolation, such as inhabiting a remote island, may be required for reducing gene flow. Our findings highlight the role of the environment, adaptation, and geographic isolation on gene flow.

Clinical phenotype and genetic analysis of a rare case with 6p duplication and terminal deletion syndrome

To explore the genetic basis for a child with developmental delay and intellectual deficit (DD/ID). A child who was admitted to the Maternal and Child Health Care Hospital of Longhua District of Shenzhen City on June 3, 2023 due to DD/ID, craniofacial malformations, and recurrent infections of upper respiratory tract was selected as the study subject. G-banded chromosomal karyotyping was carried out for the child and her parents. Low-depth genome-wide copy number variation sequencing (CNV-seq) and chromosomal microarray analysis (CMA) were used to screen for genome-wide copy number variations (CNV), and fluorescence in situ hybridization (FISH) was used to verify the origin of candidate CNV. The child, an 8-year-old girl, had featured unexplained growth and intellectual development delay, multiple craniofacial malformations, and recurrent infections of the upper respiratory tract. She was found to have a karyotype of 46,XX,der(6)add(6)(q23), while both of her parents were normal. Both CNV-seq and CMA showed that the child has harbored a 21.38 Mb interstitial duplication at 6p25.3p22.3 and a 0.78 Mb terminal deletion at 6p25. FISH verified that both the duplication and deletion had occurred de novo. The abnormal phenotype of the child may be attributed to the 6p duplication and terminal deletion.

Meiosis I causes a high spontaneous mutation rate in a multicellular red alga (Pyropia yezoensis) with a complex life cycle

Mutations are the origin of genetic diversity and are fundamental parameters needed to understand the molecular evolution of species. Estimations of mutation rates have been conducted for many diverse taxa, although rates in several major eukaryotic lineages remain unexplored. Here, the first estimation is reported of the spontaneous mutation rate for the multicellular eukaryote red alga, Pyropia yezoensis, which exhibits a complex life cycle. An estimated mutation rate of 2.97 × 10−8 (95 % CI: 2.16 × 10−8–3.99 × 10−8) per site per generation was generated for the primary life cycle, the sexual cycle, which is the highest sexual mutation rate among published sexual plants. Combined with tetrad analysis, meiosis I was identified as the primary period responsible for the high mutation rate during the complex life cycle of P. yezoensis. This result provides direct evidence for the “meiosis is mutagenic” hypothesis for multicellular organisms. The accurate estimate of the mutation rate of P. yezoensis also informs several immediate applications. Based on the above estimate, the effective population size (Ne) of P. yezoensis was estimated at about 19,000, with extensive haploid phases and asexual reproduction through monospores possibly leading to linked selection that may reduce the genome-wide genetic diversity of P. yezoensis and consequently influence Ne estimation. Lastly, P. yezoensis was estimated to have diverged from P. haitanensis about 4.2 Ma, representing a more recent date than estimates from fossil-calibrated phylogenies. These findings provide valuable new information for understanding the evolution of red algae, in addition to the underlying mechanism of mutations.

Spatially exploring RNA biology in archival formalin-fixed paraffin-embedded tissues.

The capability to spatially explore RNA biology in formalin-fixed paraffin-embedded (FFPE) tissues holds transformative potential for histopathology research. Here, we present pathology-compatible deterministic barcoding in tissue (Patho-DBiT) by combining in situ polyadenylation and computational innovation for spatial whole transcriptome sequencing, tailored to probe the diverse RNA species in clinically archived FFPE samples. It permits spatial co-profiling of gene expression and RNA processing, unveiling region-specific splicing isoforms, and high-sensitivity transcriptomic mapping of clinical tumor FFPE tissues stored for 5 years. Furthermore, genome-wide single-nucleotide RNA variants can be captured to distinguish malignant subclones from non-malignant cells in human lymphomas. Patho-DBiT also maps microRNA regulatory networks and RNA splicing dynamics, decoding their roles in spatial tumorigenesis. Single-cell level Patho-DBiT dissects the spatiotemporal cellular dynamics driving tumor clonal architecture and progression. Patho-DBiT stands poised as a valuable platform to unravel rich RNA biology in FFPE tissues to aid in clinical pathology evaluation.

Mapping genome-wide diversity and population dynamics in Indian chicken breeds for targeted conservation and breeding

ABSTRACT 1. Genetic improvement and widespread use of artificial selection may have impacted the genetic make-up of Indian chicken breeds. The genetic architecture of contemporary chicken population of India needs to be assessed for future improvement and conservation programmes. This study utilised whole-genome sequences in 180 chicken samples from 16 indigenous breeds, along with the Red Jungle Fowl and the commercial White Leghorn. 2. A panel of 76 978 genome-wide single-nucleotide polymorphisms (SNP) was selected for comparative genome analysis after stringent screening. Breeds originating from the eastern regions of India exhibited higher genomic diversity, indicative of a rich repository of distinct germplasm. Conversely, the Uttara breed, from the northern hilly areas, display considerable genetic differentiation with diminished diversity compared to others, underscoring conservation concerns. The average coefficient (FIS) of 0.084 caution the need to mitigate risks associated with inbreeding. 3. The study revealed that the analysis of 76 978 genome-wide SNP will serve as a cornerstone in refining conservation strategies, to design interventions with greater precision. 4. The contribution of Red Jungle Fowl to the gene pool of all native breeds was supported by this study. Genetic structuring indicated a relationship among breeds based on geographical proximity, underscored by varying levels of admixture.

Comparative genomics and evolutionary analysis of dengue virus strains circulating in Pakistan.

Dengue fever virus (DENV) poses a significant public health risk in tropical and subtropical regions across the world. Although the dengue fever virus (DENV) exhibits significant genetic diversity and has the potential to evolve, there is a lack of comprehensive research on the comparative genomics and evolutionary dynamics of the virus in Pakistan. Phylogenetic analysis demonstrated the circulation of all four dengue virus serotypes (DENV-1, -2, -3, and -4) with prevalent genotypes III and V within DENV-1, cosmopolitan genotype within DENV-2, genotype III within DENV-3, and genotype I within DENV-4 during 2006-2014. Based on the complete envelope region, genome-wide residue signature and genetic diversity indicate that there is a high level of genetic diversity among DENV-1 strains, while DENV-3 strains exhibit the least genetic diversity. Comparative analysis of all four DENV serotypes revealed that certain codons in DENV-2 and -4 were subject to strong purifying selection, while a few codon sites in the envelope region showed evidence of positive selection. These findings provided valuable insights into the comparative genomics and evolutionary pattern of DENV strains reported from Pakistan. Whether those characteristics conferred a fitness advantage to DENV-1 genotypes within a specific geography and time interval warrants further investigations. The findings of the current study will contribute to tracking disease dynamics, understanding virus transmission and evolution, and formulating effective disease control strategies.

Genome-wide copy number variation regions in indigenous (Bos indicus) cattle breeds of Tamil Nadu, India.

Identification of large scale structural polymorphisms (Copy number variations) of more than 50 bp between the individuals of a species would help in knowing genetic diversity, phenotypic variability, adaptability to tropical environment and disease resistance. Read depth-based method implemented in CNVnator was used for calling copy number variant regions on sequenced data obtained from WGS from 15 pooled samples belonging to five draught cattle breeds of Tamil Nadu. A total of 11,605 CNV regions (CNVRs) were observed covering a genome size of 18.63 percent. Among these, 11,459 were restricted to autosomes, consisting of 11,013 deletions, 353 duplications and 93 complex events. These CNVRs were annotated to 4,989 candidate genes. A total of 8,291 numbers of CNVRs were shared among the five cattle breeds as also supported by PCA and STRUCTURE analyses and 1,172 CNVRs were breed-specific. Four out of five selected breed-specific CNVRs were validated using real-time PCR. Genes with CNVRs are related to milk production (BTN1A1, ABCA1 and LAP3), disease resistance (TLR4 and DNAH8), adaptability (SOD1, CAST and SMARCAL1), growth (EGFR, NKAIN3), reproduction (BRWD1 and PDE6D), meat and carcass traits (MAP3K5 and NCAM1) and exterior (ATRN and MITF) traits. Gene enrichment analysis based on the gene list retrieved from the CNVRs disclosed over-represented terms (p<0.01) associated with milk fat production. NETWORK analysis had identified 13 putative candidate genes involved in milk fat percentage, milk fat yield, lactation persistency, milk yield, heat tolerance, calving ease, growth and conformation traits. The genome-wide CNVRs identified in the present study produced genome-wide partial CNV map in indigenous cattle breeds of Tamil Nadu.