Discovery of genomic variations by whole-genome resequencing of the North American Araucana chicken

Gallus gallus (chicken) is phenotypically diverse, with over 60 recognized breeds, among the myriad species within the Aves lineage. Domestic chickens have been under artificial selection by humans for thousands of years for agricultural purposes. The North American Araucana (NAA) breed arose as a cross between the Chilean “Collonocas” that laid blue eggs and was rumpless and the “Quetros” that had unusual tufts but with tail. NAAs were introduced from South America in the 1940s and have been kept as show birds by enthusiasts since then due to several distinctive traits: laying eggs with blue eggshells, characteristic ear-tufts, a pea comb, and rumplessness. The population has maintained variants for clean-faced and tufted, as well as tailed and rumplessness traits making it advantageous for genetic studies. Genome resequencing of six NAA chickens with a mixture of these traits was done to 71-fold coverage using Illumina HiSeq 2000 paired-end reads. Trimmed and concordant reads were mapped to the Gallus_gallus-5.0 reference genome (galGal5), generated from a female Red Junglefowl (UCD001). To identify candidate genes that are associated with traits of the NAA, their genome was compared with the Korean Araucana, Korean Domestic and White Leghorn breeds. Genomic regions with significantly reduced levels of heterogeneity were detected on five different chromosomes in NAA. The sequence data generated confirm the identity of variants responsible for the blue eggshells, pea comb, and rumplessness traits of NAA and propose one for ear-tufts

Analysis of pollen-specific alternative splicing in Arabidopsis thaliana via semi-quantitative PCR

Alternative splicing enables a single gene to produce multiple mRNA isoforms by varying splice site selection. In animals, alternative splicing of mRNA isoforms between cell types is widespread and supports cellular differentiation. In plants, at least 20% of multi-exon genes are alternatively spliced, but the extent and significance of tissue-specific splicing is less well understood, partly because it is difficult to isolate cells of a single type. Pollen is a useful model system to study tissue-specific splicing in higher plants because pollen grains contain only two cell types and can be collected in large amounts without damaging cells. Previously, we identified pollen-specific splicing patterns by comparing RNA-Seq data from Arabidopsis pollen and leaves. Here, we used semi-quantitative PCR to validate pollen-specific splicing patterns among genes where RNA-Seq data analysis indicated splicing was most different between pollen and leaves. PCR testing confirmed eight of nine alternative splicing patterns, and results from the ninth were inconclusive. In four genes, alternative transcriptional start sites coincided with alternative splicing. This study highlights the value of the low-cost PCR assay as a method of validating RNA-Seq results

Many rice genes are differentially spliced between roots and shoots but cytokinin has minimal effect on splicing

Abstract Alternatively spliced genes produce multiple spliced isoforms, called transcript variants. In differential alternative splicing, transcript variant abundance differs across sample types. Differential alternative splicing is common in animal systems and influences cellular development in many processes, but its extent and significance is not as well known in plants. To investigate differential alternative splicing in plants, we examined RNA‐Seq data from rice seedlings. The data included three biological replicates per sample type, approximately 30 million sequence alignments per replicate, and four sample types: roots and shoots treated with exogenous cytokinin delivered hydroponically or a mock treatment. Cytokinin treatment triggered expression changes in thousands of genes but had negligible effect on splicing patterns. However, many genes were differentially spliced between mock‐treated roots and shoots, indicating that our methods were sufficiently sensitive to detect differential splicing between data sets. Quantitative fragment analysis of reverse transcriptase‐PCR products made from newly prepared rice samples confirmed 9 of 10 differential splicing events between rice roots and shoots. Differential alternative splicing typically changed the relative abundance of splice variants that co‐occurred in a data set. Analysis of a similar (but less deeply sequenced) RNA‐Seq data set from Arabidopsis showed the same pattern. In both the Arabidopsis and rice RNA‐Seq data sets, most genes annotated as alternatively spliced had small minor variant frequencies. Of splicing choices with abundant support for minor forms, most alternative splicing events were located within the protein‐coding sequence and maintained the annotated reading frame. A tool for visualizing protein annotations in the context of genomic sequence (ProtAnnot) together with a genome browser (Integrated Genome Browser) were used to visualize and assess effects of differential splicing on gene function. In general, differentially spliced regions coincided with conserved protein domains, indicating that differential alternative splicing is likely to affect protein function between root and shoot tissue in rice

Discovery of genomic variations by whole-genome resequencing of the North American Araucana chicken.

Gallus gallus (chicken) is phenotypically diverse, with over 60 recognized breeds, among the myriad species within the Aves lineage. Domestic chickens have been under artificial selection by humans for thousands of years for agricultural purposes. The North American Araucana (NAA) breed arose as a cross between the Chilean "Collonocas" that laid blue eggs and was rumpless and the "Quetros" that had unusual tufts but with tail. NAAs were introduced from South America in the 1940s and have been kept as show birds by enthusiasts since then due to several distinctive traits: laying eggs with blue eggshells, characteristic ear-tufts, a pea comb, and rumplessness. The population has maintained variants for clean-faced and tufted, as well as tailed and rumplessness traits making it advantageous for genetic studies. Genome resequencing of six NAA chickens with a mixture of these traits was done to 71-fold coverage using Illumina HiSeq 2000 paired-end reads. Trimmed and concordant reads were mapped to the Gallus_gallus-5.0 reference genome (galGal5), generated from a female Red Junglefowl (UCD001). To identify candidate genes that are associated with traits of the NAA, their genome was compared with the Korean Araucana, Korean Domestic and White Leghorn breeds. Genomic regions with significantly reduced levels of heterogeneity were detected on five different chromosomes in NAA. The sequence data generated confirm the identity of variants responsible for the blue eggshells, pea comb, and rumplessness traits of NAA and propose one for ear-tufts

Genome-wide association mapping and identification of candidate genes for the rumpless and ear-tufted traits of the Araucana chicken.

Araucana chickens are known for their rounded, tailless rumps and tufted ears. Inheritance studies have shown that the rumpless (Rp) and ear-tufted (Et) loci each act in an autosomal dominant fashion, segregate independently, and are associated with an increased rate of embryonic mortality. To find genomic regions associated with Rp and Et, we generated genome-wide SNP profiles for a diverse population of 60 Araucana chickens using the 60 K chicken SNP BeadChip. Genome-wide association studies using 40 rumpless and 11 tailed birds showed a strong association with rumpless on Gga 2 (P(raw) = 2.45×10(-10), P(genome) = 0.00575), and analysis of genotypes revealed a 2.14 Mb haplotype shared by all rumpless birds. Within this haplotype, a 0.74 Mb critical interval containing two Iroquois homeobox genes, Irx1 and Irx2, was unique to rumpless Araucana chickens. Irx1 and Irx2 are central for developmental prepatterning, but neither gene is known to have a role in mechanisms leading to caudal development. A second genome-wide association analysis using 30 ear-tufted and 28 non-tufted birds revealed an association with tufted on Gga 15 (P(raw) = 6.61×10(-7), P(genome) = 0.0981). We identified a 0.58 Mb haplotype common to tufted birds and harboring 7 genes. Because homozygosity for Et is nearly 100% lethal, we employed a heterozygosity mapping approach to prioritize candidate gene selection. A 60 kb region heterozygous in all Araucana chickens contains the complete coding sequence for TBX1 and partial sequence for GNB1L. TBX1 is an important transcriptional regulator of embryonic development and a key genetic determinant of human DiGeorge syndrome. Herein, we describe localization of Rp and Et and identification of positional candidate genes

<i>IRX1</i> and <i>IRX2</i> are misexpressed in Araucana<i>^Rp</i> tailbud.

<p><i>IRX1</i> expression pattern (A–F,M–P). (A,C,E) Control embryos express <i>IRX1</i> in the neural tube. <i>IRX1</i> expression in controls matches Araucana<i>^Rp</i> at HH14 (A,B), but is misexpressed in Araucana<i>^Rp</i> tailbud at HH15 (D-arrowhead). Normal expression at level of somites is within neural tube (transverse section, M). Misexpression in Araucana<i>^Rp</i> can be seen at the level of the chordoneural hinge (transverse section-N and sagittal section-P) compared to expression in HH15 control (sagittal section-O). Misexpression is maintained through HH16 in Araucana<i>^Rp</i> (F-arrow) as compared to control (E). <i>IRX2</i> expression pattern (G–L,Q). (G,I,K) Control embryos do not express <i>IRX2</i> in the tailbud. No difference in expression between controls and Araucana<i>^Rp</i> seen at HH14 (G,H). <i>IRX2</i> is misexpressed in HH15 (J-arrowhead) and HH16 (L-arrowhead) Araucana<i>^Rp</i>. Transverse section of <i>IRX2</i> (Q) shows expression similar to <i>IRX1</i> at level of chordoneural hinge. A–L - anterior to top. M,N,Q - dorsal to top. O,P - dorsal to left, anterior to top. Abbreviations: nt-neural tube, s-somite, n-notochord.</p

A Novel Gain-Of-Function Mutation of the Proneural <i>IRX1</i> and <i>IRX2</i> Genes Disrupts Axis Elongation in the Araucana Rumpless Chicken

<div><p>Axis elongation of the vertebrate embryo involves the generation of cell lineages from posterior progenitor populations. We investigated the molecular mechanism governing axis elongation in vertebrates using the Araucana rumpless chicken. Araucana embryos exhibit a defect in axis elongation, failing to form the terminal somites and concomitant free caudal vertebrae, pygostyle, and associated tissues of the tail. Through whole genome sequencing of six Araucana we have identified a critical 130 kb region, containing two candidate causative SNPs. Both SNPs are proximal to the <i>IRX1</i> and <i>IRX2</i> genes, which are required for neural specification. We show that <i>IRX1</i> and <i>IRX2</i> are both misexpressed within the bipotential chordoneural hinge progenitor population of Araucana embryos. Expression analysis of <i>BRA</i> and <i>TBX6,</i> required for specification of mesoderm, shows that both are downregulated, whereas <i>SOX2</i>, required for neural patterning, is expressed in ectopic epithelial tissue. Finally, we show downregulation of genes required for the protection and maintenance of the tailbud progenitor population from the effects of retinoic acid. Our results support a model where the disruption in balance of mesoderm and neural fate results in early depletion of the progenitor population as excess neural tissue forms at the expense of mesoderm, leading to too few mesoderm cells to form the terminal somites. Together this cascade of events leads to axis truncation.</p></div