Browsing by Author "Ma, Li"
Now showing 1 - 20 of 20
Results Per Page
Sort Options
Item Analyses of inter-individual variations of sperm DNA methylation and their potential implications in cattle(Springer Nature, 2019-11-21) Liu, Shuli; Fang, Lingzhao; Zhou, Yang; Santos, Daniel J.A.; Xiang, Ruidong; Daetwyler, Hans D.; Chamberlain, Amanda J.; Cole, John B.; Li, Cong-jun; Yu, Ying; Ma, Li; Zhang, Shengli; Liu, George E.DNA methylation has been shown to be involved in many biological processes, including X chromosome inactivation in females, paternal genomic imprinting, and others.Item Array CGH-based detection of CNV regions and their potential association with reproduction and other economic traits in Holsteins(Springer Nature, 2019-03-07) Liu, Mei; Fang, Lingzhao; Liu, Shuli; Pan, Michael G.; Seroussi, Eyal; Cole, John B.; Ma, Li; Chen, Hong; Liu, George E.Copy number variations (CNVs) are structural variants consisting of large-scale insertions and deletions of genomic fragments. Exploring CNVs and estimating their effects on phenotypes are useful for genome selection but remain challenging in the livestock.Item Bos taurus–indicus hybridization correlates with intralocus sexual-conflict effects of PRDM9 on male and female fertility in Holstein cattle(Springer Nature, 2019-08-28) Seroussi, Eyal; Shirak, Andrey; Gershoni, Moran; Ezra, Ephraim; Santos, Daniel Jordan de Abreu; Ma, Li; Liu, George E.Crossover localization during meiotic recombination is mediated by the fast-evolving zinc-finger (ZnF) domain of gene PRDM9. To study its impact on dairy cattle performance, we compared its genetic variation between the relatively small Israeli (IL) Holsteins and the North American (US) Holsteins that count millions.Item Characterization of Copy Number Variation’s Potential Role in Marek’s Disease(MDPI, 2017-05-09) Xu, Lingyang; He, Yanghua; Ding, Yi; Sun, Guirong; Carrillo, Jose Adrian; Li, Yaokun; Ghaly, Mona M.; Ma, Li; Zhang, Huanmin; Liu, George E.; Song, JiuzhouMarek’s Disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. Chicken Lines 63 and 72, as well as their recombinant congenic strains (RCS) with varied susceptibility to MD, are ideal models to study the complex mechanisms of genetic resistance to MD. In this study, we investigated copy number variation (CNV) in these inbred chicken lines using the Affymetrix Axiom HD 600 K SNP genotyping array. We detected 393 CNV segments across all ten chicken lines, of which 12 CNVs were specifically identified in Line 72. We then assessed genetic structure based on CNV and observed markedly different patterns. Finally, we validated two deletion events in Line 72 and correlated them with genes expression using qPCR and RNA-seq, respectively. Our combined results indicated that these two CNV deletions were likely to contribute to MD susceptibility.Item Characterization of recombination features and the genetic basis in multiple cattle breeds(Springer Nature, 2018-04-27) Shen, Botong; Jiang, Jicai; Seroussi, Eyal; Liu, George E.; Ma, LiCrossover generated by meiotic recombination is a fundamental event that facilitates meiosis and sexual reproduction. Comparative studies have shown wide variation in recombination rate among species, but the characterization of recombination features between cattle breeds has not yet been performed. Cattle populations in North America count millions, and the dairy industry has genotyped millions of individuals with pedigree information that provide a unique opportunity to study breed-level variations in recombination. Based on large pedigrees of Jersey, Ayrshire and Brown Swiss cattle with genotype data, we identified over 3.4 million maternal and paternal crossover events from 161,309 three-generation families. We constructed six breed- and sex-specific genome-wide recombination maps using 58,982 autosomal SNPs for two sexes in the three dairy cattle breeds. A comparative analysis of the six recombination maps revealed similar global recombination patterns between cattle breeds but with significant differences between sexes. We confirmed that male recombination map is 10% longer than the female map in all three cattle breeds, consistent with previously reported results in Holstein cattle. When comparing recombination hotspot regions between cattle breeds, we found that 30% and 10% of the hotspots were shared between breeds in males and females, respectively, with each breed exhibiting some breed-specific hotspots. Finally, our multiple-breed GWAS found that SNPs in eight loci affected recombination rate and that the PRDM9 gene associated with hotspot usage in multiple cattle breeds, indicating a shared genetic basis for recombination across dairy cattle breeds. Collectively, our results generated breed- and sex-specific recombination maps for multiple cattle breeds, provided a comprehensive characterization and comparison of recombination patterns between breeds, and expanded our understanding of the breed-level variations in recombination features within an important livestock species.Item Comparative transcriptome in large-scale human and cattle populations(Springer Nature, 2022-08-22) Yao, Yuelin; Liu, Shuli; Xia, Charley; Gao, Yahui; Pan, Zhangyuan; Canela-Xandri, Oriol; Khamseh, Ava; Rawlik, Konrad; Wang, Sheng; Li, Bingjie; Zhang, Yi; Pairo-Castineira, Erola; D’Mellow, Kenton; Li, Xiujin; Yan, Ze; Li, Cong-jun; Yu, Ying; Zhang, Shengli; Ma, Li; Cole, John B.; Ross, Pablo J.; Zhou, Huaijun; Haley, Chris; Liu, George E.; Fang, Lingzhao; Tenesa, AlbertCross-species comparison of transcriptomes is important for elucidating evolutionary molecular mechanisms underpinning phenotypic variation between and within species, yet to date it has been essentially limited to model organisms with relatively small sample sizes. Here, we systematically analyze and compare 10,830 and 4866 publicly available RNA-seq samples in humans and cattle, respectively, representing 20 common tissues. Focusing on 17,315 orthologous genes, we demonstrate that mean/median gene expression, inter-individual variation of expression, expression quantitative trait loci, and gene co-expression networks are generally conserved between humans and cattle. By examining large-scale genome-wide association studies for 46 human traits (average n = 327,973) and 45 cattle traits (average n = 24,635), we reveal that the heritability of complex traits in both species is significantly more enriched in transcriptionally conserved than diverged genes across tissues. In summary, our study provides a comprehensive comparison of transcriptomes between humans and cattle, which might help decipher the genetic and evolutionary basis of complex traits in both species.Item Comparative whole genome DNA methylation profiling across cattle tissues reveals global and tissue-specific methylation patterns(Springer Nature, 2020-07-06) Zhou, Yang; Liu, Shuli; Hu, Yan; Fang, Lingzhao; Gao, Yahui; Xia, Han; Schroeder, Steven G.; Rosen, Benjamin D.; Connor, Erin E.; Li, Cong-jun; Baldwin, Ransom L.; Cole, John B.; Van Tassell, Curtis P.; Yang, Liguo; Ma, Li; Liu, George E.Efforts to improve animal health, and understand genetic bases for production, may benefit from a comprehensive analysis of animal genomes and epigenomes. Although DNA methylation has been well studied in humans and other model species, its distribution patterns and regulatory impacts in cattle are still largely unknown. Here, we present the largest collection of cattle DNA methylation epigenomic data to date. Using Holstein cattle, we generated 29 whole genome bisulfite sequencing (WGBS) datasets for 16 tissues, 47 corresponding RNA-seq datasets, and 2 whole genome sequencing datasets. We did read mapping and DNA methylation calling based on two different cattle assemblies, demonstrating the high quality of the long-read-based assembly markedly improved DNA methylation results. We observed large differences across cattle tissues in the methylation patterns of global CpG sites, partially methylated domains (PMDs), hypomethylated regions (HMRs), CG islands (CGIs), and common repeats. We detected that each tissue had a distinct set of PMDs, which showed tissue-specific patterns. Similar to human PMD, cattle PMDs were often linked to a general decrease of gene expression and a decrease in active histone marks and related to long-range chromatin organizations, like topologically associated domains (TADs). We tested a classification of the HMRs based on their distributions relative to transcription start sites (TSSs) and detected tissue-specific TSS-HMRs and genes that showed strong tissue effects. When performing cross-species comparisons of paired genes (two opposite strand genes with their TSS located in the same HMR), we found out they were more consistently co-expressed among human, mouse, sheep, goat, yak, pig, and chicken, but showed lower consistent ratios in more divergent species. We further used these WGBS data to detect 50,023 experimentally supported CGIs across bovine tissues and found that they might function as a guard against C-to-T mutations for TSS-HMRs. Although common repeats were often heavily methylated, some young Bov-A2 repeats were hypomethylated in sperm and could affect the promoter structures by exposing potential transcription factor binding sites. This study provides a comprehensive resource for bovine epigenomic research and enables new discoveries about DNA methylation and its role in complex traits.Item Dissection of additive, dominance, and imprinting effects for production and reproduction traits in Holstein cattle(Springer Nature, 2017-05-30) Jiang, Jicai; Shen, Botong; O’Connell, Jeffrey R.; VanRaden, Paul M.; Cole, John B.; Ma, LiAlthough genome-wide association and genomic selection studies have primarily focused on additive effects, dominance and imprinting effects play an important role in mammalian biology and development. The degree to which these non-additive genetic effects contribute to phenotypic variation and whether QTL acting in a non-additive manner can be detected in genetic association studies remain controversial. To empirically answer these questions, we analyzed a large cattle dataset that consisted of 42,701 genotyped Holstein cows with genotyped parents and phenotypic records for eight production and reproduction traits. SNP genotypes were phased in pedigree to determine the parent-of-origin of alleles, and a three-component GREML was applied to obtain variance decomposition for additive, dominance, and imprinting effects. The results showed a significant non-zero contribution from dominance to production traits but not to reproduction traits. Imprinting effects significantly contributed to both production and reproduction traits. Interestingly, imprinting effects contributed more to reproduction traits than to production traits. Using GWAS and imputation-based fine-mapping analyses, we identified and validated a dominance association signal with milk yield near RUNX2, a candidate gene that has been associated with milk production in mice. When adding non-additive effects into the prediction models, however, we observed little or no increase in prediction accuracy for the eight traits analyzed. Collectively, our results suggested that non-additive effects contributed a non-negligible amount (more for reproduction traits) to the total genetic variance of complex traits in cattle, and detection of QTLs with non-additive effect is possible in GWAS using a large dataset.Item Functional annotation of the cattle genome through systematic discovery and characterization of chromatin states and butyrate-induced variations(Springer Nature, 2019-08-16) Fang, Lingzhao; Liu, Shuli; Liu, Mei; Kang, Xiaolong; Lin, Shudai; Li, Bingjie; Connor, Erin E.; Baldwin, Ransom L. VI; Tenesa, Albert; Ma, Li; Liu, George E.; Li, Cong-junThe functional annotation of genomes, including chromatin accessibility and modifications, is important for understanding and effectively utilizing the increased amount of genome sequences reported. However, while such annotation has been well explored in a diverse set of tissues and cell types in human and model organisms, relatively little data are available for livestock genomes, hindering our understanding of complex trait variation, domestication, and adaptive evolution. Here, we present the first complete global landscape of regulatory elements in cattle and explore the dynamics of chromatin states in rumen epithelial cells induced by the rumen developmental regulator—butyrate.Item Genome changes due to artificial selection in U.S. Holstein cattle(Springer Nature, 2019-02-11) Ma, Li; Sonstegard, Tad S.; Cole, John B.; VanTassell, Curtis P.; Wiggans, George R.; Crooker, Brian A.; Tan, Cheng; Prakapenka, Dzianis; Liu, George E.; Da, YangThe availability of a unique unselected Holstein line since 1964 provided a direct comparison between selected and unselected Holstein genomes whereas large Holstein samples provided unprecedented statistical power for identifying high-confidence SNP effects. Utilizing these unique resources, we aimed to identify genome changes affected by selection since 1964.Item Genome-wide recombination map construction from single sperm sequencing in cattle(Springer Nature, 2022-03-05) Yang, Liu; Gao, Yahui; Li, Mingxun; Park, Ki-Eun; Liu, Shuli; Kang, Xiaolong; Liu, Mei; Oswalt, Adam; Fang, Lingzhao; Telugu, Bhanu P.; Sattler, Charles G.; Li, Cong-jun; Cole, John B.; Seroussi, Eyal; Xu, Lingyang; Yang, Lv; Zhou, Yang; Li, Li; Zhang, Hongping; Rosen, Benjamin D.; Van Tassell, Curtis P.; Ma, Li; Liu, George E.Meiotic recombination is one of the important phenomena contributing to gamete genome diversity. However, except for human and a few model organisms, it is not well studied in livestock, including cattle. To investigate their distributions in the cattle sperm genome, we sequenced 143 single sperms from two Holstein bulls. We mapped meiotic recombination events at high resolution based on phased heterozygous single nucleotide polymorphism (SNP). In the absence of evolutionary selection pressure in fertilization and survival, recombination events in sperm are enriched near distal chromosomal ends, revealing that such a pattern is intrinsic to the molecular mechanism of meiosis. Furthermore, we further validated these findings in single sperms with results derived from sequencing its family trio of diploid genomes and our previous studies of recombination in cattle. To our knowledge, this is the first large-scale single sperm whole-genome sequencing effort in livestock, which provided useful information for future studies of recombination, genome instability, and male infertility.Item Genomic predictions combining SNP markers and copy number variations in Nellore cattle(Springer Nature, 2018-06-05) Hay, El Hamidi A.; Utsunomiya, Yuri T.; Xu, Lingyang; Zhou, Yang; Neves, Haroldo H. R.; Carvalheiro, Roberto; Bickhart, Derek M.; Ma, Li; Garcia, Jose Fernando; Liu, George E.Due to the advancement in high throughput technology, single nucleotide polymorphism (SNP) is routinely being incorporated along with phenotypic information into genetic evaluation. However, this approach often cannot achieve high accuracy for some complex traits. It is possible that SNP markers are not sufficient to predict these traits due to the missing heritability caused by other genetic variations such as microsatellite and copy number variation (CNV), which have been shown to affect disease and complex traits in humans and other species. In this study, CNVs were included in a SNP based genomic selection framework. A Nellore cattle dataset consisting of 2230 animals genotyped on BovineHD SNP array was used, and 9 weight and carcass traits were analyzed. A total of six models were implemented and compared based on their prediction accuracy. For comparison, three models including only SNPs were implemented: 1) BayesA model, 2) Bayesian mixture model (BayesB), and 3) a GBLUP model without polygenic effects. The other three models incorporating both SNP and CNV included 4) a Bayesian model similar to BayesA (BayesA+CNV), 5) a Bayesian mixture model (BayesB+CNV), and 6) GBLUP with CNVs modeled as a covariable (GBLUP+CNV). Prediction accuracies were assessed based on Pearson’s correlation between de-regressed EBVs (dEBVs) and direct genomic values (DGVs) in the validation dataset. For BayesA, BayesB and GBLUP, accuracy ranged from 0.12 to 0.62 across the nine traits. A minimal increase in prediction accuracy for some traits was noticed when including CNVs in the model (BayesA+CNV, BayesB+CNV, GBLUP+CNV). This study presents the first genomic prediction study integrating CNVs and SNPs in livestock. Combining CNV and SNP marker information proved to be beneficial for genomic prediction of some traits in Nellore cattle.Item GWAS and fine-mapping of livability and six disease traits in Holstein cattle(Springer Nature, 2020-01-13) Freebern, Ellen; Santos, Daniel J. A.; Fang, Lingzhao; Jiang, Jicai; Parker Gaddis, Kristen L.; Liu, George E.; VanRaden, Paul M.; Maltecca, Christian; Cole, John B.; Ma, LiHealth traits are of significant economic importance to the dairy industry due to their effects on milk production and associated treatment costs. Genome-wide association studies (GWAS) provide a means to identify associated genomic variants and thus reveal insights into the genetic architecture of complex traits and diseases. The objective of this study is to investigate the genetic basis of seven health traits in dairy cattle and to identify potential candidate genes associated with cattle health using GWAS, fine mapping, and analyses of multi-tissue transcriptome data.Item Insights from Initial Variant Detection by Sequencing Single Sperm in Cattle(MDPI, 2021-11-15) Yang, Liu; Gao, Yahui; Boschiero, Clarissa; Li, Li; Zhang, Hongping; Ma, Li; Liu, George E.Meiotic de novo mutation (DNM) is one of the important phenomena contributing to gamete genome diversity. However, except for humans and a few model organisms, they are not well studied in livestock, including cattle. Moreover, bulk sperm samples have been routinely utilized in experiments, which include millions of single sperm cells and only report high-frequency variants. In this study, we isolated and sequenced 143 single sperms from two Holstein bulls and identified hundreds of candidate DNM events in ten sperms with deep sequencing coverage. We estimated DNM rates ranging from 1.08 × 10−8 to 3.78 × 10−8 per nucleotide per generation. We further validated 12 out of 14 selected DNM events using Sanger sequencing. To our knowledge, this is the first single sperm whole-genome sequencing effort in livestock, which provided useful information for future studies of point mutations and male fertility. Our preliminary results pointed out future research directions and highlighted the importance of uniform whole genome amplification, deep sequence coverage, and dedicated software pipelines for genetic variant detection using single-cell sequencing data.Item Integrating RNA-Seq with GWAS reveals novel insights into the molecular mechanism underpinning ketosis in cattle(Springer Nature, 2020-07-17) Yan, Ze; Huang, Hetian; Freebern, Ellen; Santos, Daniel J. A.; Dai, Dongmei; Si, Jingfang; Ma, Chong; Cao, Jie; Guo, Gang; Liu, George E.; Ma, Li; Fang, Lingzhao; Zhang, YiKetosis is a common metabolic disease during the transition period in dairy cattle, resulting in long-term economic loss to the dairy industry worldwide. While genetic selection of resistance to ketosis has been adopted by many countries, the genetic and biological basis underlying ketosis is poorly understood. We collected a total of 24 blood samples from 12 Holstein cows, including 4 healthy and 8 ketosis-diagnosed ones, before (2 weeks) and after (5 days) calving, respectively. We then generated RNA-Sequencing (RNA-Seq) data and seven blood biochemical indicators (bio-indicators) from leukocytes and plasma in each of these samples, respectively. By employing a weighted gene co-expression network analysis (WGCNA), we detected that 4 out of 16 gene-modules, which were significantly engaged in lipid metabolism and immune responses, were transcriptionally (FDR < 0.05) correlated with postpartum ketosis and several bio-indicators (e.g., high-density lipoprotein and low-density lipoprotein). By conducting genome-wide association signal (GWAS) enrichment analysis among six common health traits (ketosis, mastitis, displaced abomasum, metritis, hypocalcemia and livability), we found that 4 out of 16 modules were genetically (FDR < 0.05) associated with ketosis, among which three were correlated with postpartum ketosis based on WGCNA. We further identified five candidate genes for ketosis, including GRINA, MAF1, MAFA, C14H8orf82 and RECQL4. Our phenome-wide association analysis (Phe-WAS) demonstrated that human orthologues of these candidate genes were also significantly associated with many metabolic, endocrine, and immune traits in humans. For instance, MAFA, which is involved in insulin secretion, glucose response, and transcriptional regulation, showed a significantly higher association with metabolic and endocrine traits compared to other types of traits in humans. In summary, our study provides novel insights into the molecular mechanism underlying ketosis in cattle, and highlights that an integrative analysis of omics data and cross-species mapping are promising for illustrating the genetic architecture underpinning complex traits.Item Investigation of rumen long noncoding RNA before and after weaning in cattle(Springer Nature, 2022-07-20) Marceau, Alexis; Gao, Yahui; Baldwin VI, Ransom L.; Li, Cong-jun; Jiang, Jicai; Liu, George E.; Ma, LiThis study aimed to identify long non-coding RNA (lncRNA) from the rumen tissue in dairy cattle, explore their features including expression and conservation levels, and reveal potential links between lncRNA and complex traits that may indicate important functional impacts of rumen lncRNA during the transition to the weaning period. A total of six cattle rumen samples were taken with three replicates from before and after weaning periods, respectively. Total RNAs were extracted and sequenced with lncRNA discovered based on size, coding potential, sequence homology, and known protein domains. As a result, 404 and 234 rumen lncRNAs were identified before and after weaning, respectively. However, only nine of them were shared under two conditions, with 395 lncRNAs found only in pre-weaning tissues and 225 only in post-weaning samples. Interestingly, none of the nine common lncRNAs were differentially expressed between the two weaning conditions. LncRNA averaged shorter length, lower expression, and lower conservation scores than the genome overall, which is consistent with general lncRNA characteristics. By integrating rumen lncRNA before and after weaning with large-scale GWAS results in cattle, we reported significant enrichment of both pre- and after-weaning lncRNA with traits of economic importance including production, reproduction, health, and body conformation phenotypes. The majority of rumen lncRNAs are uniquely expressed in one of the two weaning conditions, indicating a functional role of lncRNA in rumen development and transition of weaning. Notably, both pre- and post-weaning lncRNA showed significant enrichment with a variety of complex traits in dairy cattle, suggesting the importance of rumen lncRNA for cattle performance in the adult stage. These relationships should be further investigated to better understand the specific roles lncRNAs are playing in rumen development and cow performance.Item The role of a lens survival pathway including sox2 and αA-crystallin in the evolution of cavefish eye degeneration(Springer Nature, 2014-08-28) Ma, Li; Parkhurst, Amy; Jeffery, William RThe teleost Astyanax mexicanus is a single species consisting of eyed surface-dwelling (surface fish) and blind cave-dwelling (cavefish) morphs. Cavefish eyes are lost through apoptosis of the lens, which in turn promotes the degeneration of other optic tissues. The αA-crystallin (αA-crys) gene is strongly downregulated in the cavefish lens and is located in a genomic region (QTL) responsible for eye loss. Therefore, αA-crys has been proposed as a candidate for regulating cavefish eye degeneration. The purpose of this study was to determine the mechanism of αA-crys downregulation and its role in cavefish eye degeneration. The involvement of αA-crys in eye degeneration was confirmed by knocking down its expression in surface fish, which led to apoptosis of the lens. The underlying reason for αA-crys downregulation in cavefish was investigated by comparing genomic αA-crys DNA sequences in surface fish and cavefish, however, no obvious cis-regulatory factors were discovered. Furthermore, the cavefish αA-crys allele is expressed in surface fish x cavefish F1 hybrids, indicating that evolutionary changes in upstream genes are most likely responsible for αA-crys downregulation. In other species, Sox2 is one of the transcription factors that regulate lens crystallin genes during eye development. Determination of sox2 expression patterns during surface fish and cavefish development showed that sox2 is specifically downregulated in the cavefish lens. The upstream regulatory function of Sox2 was demonstrated by knockdown in surface fish, which abolished αA-crys expression and induced lens apoptosis. The results suggest that αA-crys is required for normal eye development in cavefish via suppression of lens apoptosis. The regulatory changes involved in αA-crys downregulation in cavefish are in trans-acting factors rather than cis-acting mutations in the αA-crys gene. Therefore, αA-crys is unlikely to be the mutated gene(s) associated with an Astyanax eye QTL. The results reveal a genetic pathway leading from sox2 to αA-crys that is required for survival of the lens in Astyanax surface fish. Defects in this pathway may be involved in lens apoptosis and thus a cause of cavefish eye degeneration.Item Ruminant-specific multiple duplication events of PRDM9 before speciation(Springer Nature, 2017-03-14) Padhi, Abinash; Shen, Botong; Jiang, Jicai; Zhou, Yang; Liu, George E.; Ma, LiUnderstanding the genetic and evolutionary mechanisms of speciation genes in sexually reproducing organisms would provide important insights into mammalian reproduction and fitness. PRDM9, a widely known speciation gene, has recently gained attention for its important role in meiotic recombination and hybrid incompatibility. Despite the fact that PRDM9 is a key regulator of recombination and plays a dominant role in hybrid incompatibility, little is known about the underlying genetic and evolutionary mechanisms that generated multiple copies of PRDM9 in many metazoan lineages. The present study reports (1) evidence of ruminant-specific multiple gene duplication events, which likely have had occurred after the ancestral ruminant population diverged from its most recent common ancestor and before the ruminant speciation events, (2) presence of three copies of PRDM9, one copy (lineages I) in chromosome 1 (chr1) and two copies (lineages II & III) in chromosome X (chrX), thus indicating the possibility of ancient inter- and intra-chromosomal unequal crossing over and gene conversion events, (3) while lineages I and II are characterized by the presence of variable tandemly repeated C2H2 zinc finger (ZF) arrays, lineage III lost these arrays, and (4) C2H2 ZFs of lineages I and II, particularly the amino acid residues located at positions −1, 3, and 6 have evolved under strong positive selection. Our results demonstrated two gene duplication events of PRDM9 in ruminants: an inter-chromosomal duplication that occurred between chr1 and chrX, and an intra-chromosomal X-linked duplication, which resulted in two additional copies of PRDM9 in ruminants. The observation of such duplication between chrX and chr1 is rare and may possibly have happened due to unequal crossing-over millions of years ago when sex chromosomes were independently derived from a pair of ancestral autosomes. Two copies (lineages I & II) are characterized by the presence of variable sized tandem-repeated C2H2 ZFs and evolved under strong positive selection and concerted evolution, supporting the notion of well-established Red Queen hypothesis. Collectively, gene duplication, concerted evolution, and positive selection are the likely driving forces for the expansion of ruminant PRDM9 sub-family.Item Towards the detection of copy number variation from single sperm sequencing in cattle(Springer Nature, 2022-03-17) Yang, Liu; Gao, Yahui; Oswalt, Adam; Fang, Lingzhao; Boschiero, Clarissa; Neupane, Mahesh; Sattler, Charles G.; Li, Cong-jun; Seroussi, Eyal; Xu, Lingyang; Yang, Lv; Li, Li; Zhang, Hongping; Rosen, Benjamin D.; Van Tassell, Curtis P.; Zhou, Yang; Ma, Li; Liu, George E.Copy number variation (CNV) has been routinely studied using bulk-cell sequencing. However, CNV is not well studied on the single-cell level except for humans and a few model organisms. We sequenced 143 single sperms of two Holstein bulls, from which we predicted CNV events using 14 single sperms with deep sequencing. We then compared the CNV results derived from single sperms with the bulk-cell sequencing of one bull’s family trio of diploid genomes. As a known CNV hotspot, segmental duplications were also predicted using the bovine ARS-UCD1.2 genome. Although the trio CNVs validated only some single sperm CNVs, they still showed a distal chromosomal distribution pattern and significant associations with segmental duplications and satellite repeats. Our preliminary results pointed out future research directions and highlighted the importance of uniform whole genome amplification, deep sequence coverage, and dedicated software pipelines for CNV detection using single cell sequencing data.Item X-inactivation informs variance-based testing for X-linked association of a quantitative trait(Springer Nature, 2015) Ma, Li; Hoffman, Gabriel; Keinan, AlonThe X chromosome plays an important role in human diseases and traits. However, few X-linked associations have been reported in genome-wide association studies, partly due to analytical complications and low statistical power. In this study, we propose tests of X-linked association that capitalize on variance heterogeneity caused by various factors, predominantly the process of X-inactivation. In the presence of X-inactivation, the expression of one copy of the chromosome is randomly silenced. Due to the consequent elevated randomness of expressed variants, females that are heterozygotes for a quantitative trait locus might exhibit higher phenotypic variance for that trait. We propose three tests that build on this phenomenon: 1) A test for inflated variance in heterozygous females; 2) A weighted association test; and 3) A combined test. Test 1 captures the novel signal proposed herein by directly testing for higher phenotypic variance of heterozygous than homozygous females. As a test of variance it is generally less powerful than standard tests of association that consider means, which is supported by extensive simulations. Test 2 is similar to a standard association test in considering the phenotypic mean, but differs by accounting for (rather than testing) the variance heterogeneity. As expected in light of X-inactivation, this test is slightly more powerful than a standard association test. Finally, test 3 further improves power by combining the results of the first two tests. We applied the these tests to the ARIC cohort data and identified a novel X-linked association near gene AFF2 with blood pressure, which was not significant based on standard association testing of mean blood pressure. Variance-based tests examine overdispersion, thereby providing a complementary type of signal to a standard association test. Our results point to the potential to improve power of detecting X-linked associations in the presence of variance heterogeneity.