Genomic insights from the first chromosome-scale assemblies of oat (Avena spp.) diploid species.

BACKGROUND: Cultivated hexaploid oat (Common oat; Avena sativa) has held a significant place within the global crop community for centuries; although its cultivation has decreased over the past century, its nutritional benefits have garnered increased interest for human consumption. We report the development of fully annotated, chromosome-scale assemblies for the extant progenitor species of the As- and Cp-subgenomes, Avena atlantica and Avena eriantha respectively. The diploid Avena species serve as important genetic resources for improving common oat's adaptive and food quality characteristics. RESULTS: The A. atlantica and A. eriantha genome assemblies span 3.69 and 3.78 Gb with an N50 of 513 and 535 Mb, respectively. Annotation of the genomes, using sequenced transcriptomes, identified ~ 50,000 gene models in each species-including 2965 resistance gene analogs across both species. Analysis of these assemblies classified much of each genome as repetitive sequence (~ 83%), including species-specific, centromeric-specific, and telomeric-specific repeats. LTR retrotransposons make up most of the classified elements. Genome-wide syntenic comparisons with other members of the Pooideae revealed orthologous relationships, while comparisons with genetic maps from common oat clarified subgenome origins for each of the 21 hexaploid linkage groups. The utility of the diploid genomes was demonstrated by identifying putative candidate genes for flowering time (HD3A) and crown rust resistance (Pc91). We also investigate the phylogenetic relationships among other A- and C-genome Avena species. CONCLUSIONS: The genomes we report here are the first chromosome-scale assemblies for the tribe Poeae, subtribe Aveninae. Our analyses provide important insight into the evolution and complexity of common hexaploid oat, including subgenome origin, homoeologous relationships, and major intra- and intergenomic rearrangements. They also provide the annotation framework needed to accelerate gene discovery and plant breeding.

Paralog analyses reveal gene duplication events and genes under positive selection in Ixodes scapularis and other ixodid ticks.

BACKGROUND: Hard ticks (family Ixodidae) are obligatory hematophagous ectoparasites of worldwide medical and veterinary importance. The haploid genomes of multiple species of ixodid ticks exceed 1 Gbp, prompting questions regarding gene, segmental and whole genome duplication in this phyletic group. The availability of the genome assembly for the black legged tick, Ixodes scapularis, and transcriptome datasets for multiple species of ticks offers an opportunity to assess the contribution of gene duplication to the genome. Here we developed a bioinformatics pipeline to identify and analyze duplicated genes (paralogs) using gene models from the prostriate tick, I. scapularis IscaW1.1 annotation and expressed sequence tags (ESTs) from I. scapularis and the metastriate ticks, Rhipicephalus microplus (southern cattle tick), R. appendiculatus (brown ear tick) and Amblyomma variegatum (tropical bont tick). RESULTS: Approximately 1-2% of I. scapularis gene models and 2-14% of ESTs from the four species represent duplicated genes. The ratio of non-synonymous to synonymous nucleotide substitution rates suggests ~ 25% of duplicated genes are under positive selection pressure in each species. Analyses of synonymous substitution rates provide evidence for two duplication events in I. scapularis and R. microplus involving several hundred genes. Conservative molecular clock estimates based on synonymous substitution rates for species of Anopheles mosquitoes and the fruit fly, Drosophila melanogaster, suggest these events occurred within the last 50 MYA. Mapping of paralogs to the I. scapularis genome assembly supports tandem, or possibly segmental duplication events. CONCLUSIONS: The present study marks the first genome-level analyses of gene duplication for the Ixodidae and provides insights into mechanisms shaping genome evolution in this group. At least two duplication events involving hundreds of genes may have occurred independently in the lineages prostriata and metastriata, with tandem and segmental duplication the most likely mechanisms for paralog generation. Duplicated genes under positive selection pressure may be linked to emerging functions in the tick and represent important candidates for further study.

A Comparison of transgenic and wild type soybean seeds: analysis of transcriptome profiles using RNA-Seq.

BACKGROUND: Soybean (Glycine max) has been bred for thousands of years to produce seeds rich in protein for human and animal consumption, making them an appealing bioreactor for producing valuable recombinant proteins at high levels. However, the effects of expressing recombinant protein at high levels on bean physiology are not well understood. To address this, we investigated whether gene expression within transgenic soybean seed tissue is altered when large amounts of recombinant proteins are being produced and stored exclusively in the seeds. We used RNA-Seq to survey gene expression in three transgenic soybean lines expressing recombinant protein at levels representing up to 1.61 % of total protein in seed tissues. The three lines included: ST77, expressing human thyroglobulin protein (hTG), ST111, expressing human myelin basic protein (hMBP), and 764, expressing a mutant, nontoxic form of a staphylococcal subunit vaccine protein (mSEB). All lines selected for analysis were homozygous and contained a single copy of the transgene. METHODS: Each transgenic soybean seed was screened for transgene presence and recombinant protein expression via PCR and western blotting.  Whole seed mRNA was extracted and cDNA libraries constructed for Illumina sequencing.  Following alignment to the soybean reference genome, differential gene expression analysis was conducted using edgeR and cufflinks.  Functional analysis of differentially expressed genes was carried out using the gene ontology analysis tool AgriGO. RESULTS: The transcriptomes of nine seeds from each transgenic line were sequenced and compared with wild type seeds. Native soybean gene expression was significantly altered in line 764 (mSEB) with more than 3000 genes being upregulated or downregulated. ST77 (hTG) and ST111 (hMBP) had significantly less differences with 52 and 307 differentially expressed genes respectively. Gene ontology enrichment analysis found that the upregulated genes in the 764 line were annotated with functions related to endopeptidase inhibitors and protein synthesis, but suppressed expression of genes annotated to the nuclear pore and to protein transport. No significant gene ontology terms were detected in ST77, and only a few genes involved in photosynthesis and thylakoid functions were downregulated in ST111. Despite these differences, transgenic plants and seeds appeared phenotypically similar to non-transgenic controls. There was no correlation between recombinant protein expression level and the quantity of differentially expressed genes detected. CONCLUSIONS: Measurable unscripted gene expression changes were detected in the seed transcriptomes of all three transgenic soybean lines analyzed, with line 764 being substantially altered. Differences detected at the transcript level may be due to T-DNA insert locations, random mutations following transformation or direct effects of the recombinant protein itself, or a combination of these. The physiological consequences of such changes remain unknown.

Statistical issues associated with modeling of synonymous mutation data.

The explosion of data in evolutionary bioinformatics has led to sometimes ad hoc, incomplete and even inaccurate data analyses. Taking dS data, namely, data on synonymous substitutions per synonymous sites, we go through a statistical analysis for modeling the time since duplications of genes. We explore the shortcomings of previous analyses, especially with a view towards their effect on inference for the gene duplication process. We present a statistical analysis which respects the assumptions of the models and the integrity of the data, and emphasize that exploratory data analysis, formulation of a data model, its estimation and finally, assessment of the model are important steps in a complete data analysis. Furthermore, for dS data, we develop Bayesian discrete-continuous mixture models and present analyses using two genomes.

A gold standard dataset and evaluation of methods for lineage abundance estimation from wastewater.

During the SARS-CoV-2 pandemic, genome-based wastewater surveillance sequencing has been a powerful tool for public health to monitor circulating and emerging viral variants. As a medium, wastewater is very complex because of its mixed matrix nature, which makes the deconvolution of wastewater samples more difficult. Here we introduce a gold standard dataset constructed from synthetic viral control mixtures of known composition, spiked into a wastewater RNA matrix and sequenced on the Oxford Nanopore Technologies platform. We compare the performance of eight of the most commonly used deconvolution tools in identifying SARS-CoV-2 variants present in these mixtures. The software evaluated was primarily chosen for its relevance to the CDC wastewater surveillance reporting protocol, which until recently employed a pipeline that incorporates results from four deconvolution methods: Freyja, kallisto, Kraken 2/Bracken, and LCS. We also tested Lollipop, a deconvolution method used by the Swiss SARS-CoV-2 Sequencing Consortium, and three additional methods not used in the C-WAP pipeline: lineagespot, Alcov, and VaQuERo. We found that the commonly used software Freyja outperformed the other CDC pipeline tools in correct identification of lineages present in the control mixtures, and that the VaQuERo method was similarly accurate, with minor differences in the ability of the two methods to avoid false negatives and suppress false positives. Our results also provide insight into the effect of the tiling primer scheme and wastewater RNA extract matrix on viral sequencing and data deconvolution outcomes.

Co-expression of soybean Dicer-like genes in response to stress and development.

Regulation of gene transcription and post-transcriptional processes is critical for proper development, genome integrity, and stress responses in plants. Many genes involved in the key processes of transcriptional and post-transcriptional regulation have been well studied in model diploid organisms. However, gene and genome duplication may alter the function of the genes involved in these processes. To address this question, we assayed the stress-induced transcription patterns of duplicated gene pairs involved in RNAi and DNA methylation processes in the paleopolyploid soybean. Real-time quantitative PCR and Sequenom MassARRAY expression assays were used to profile the relative expression ratios of eight gene pairs across eight different biotic and abiotic stress conditions. The transcriptional responses to stress for genes involved in DNA methylation, RNAi processing, and miRNA processing were compared. The strongest evidence for pairwise co-expression in response to stresses was exhibited by non-paralogous Dicer-like (DCL) genes GmDCL2a-GmDCL3a and GmDCL1b-GmDCL2b, most profoundly in root tissues. Among homoeologous or paralogous DCL genes, the Dicer-like 2 (DCL2) gene pair exhibited the strongest response to stress and most conserved co-expression pattern. This was surprising because the DCL2 duplication event is more ancient than the other DCL duplications. Possible mechanisms that may be driving the DCL2 co-expression are discussed.

Large-scale development of cost-effective SNP marker assays for diversity assessment and genetic mapping in chickpea and comparative mapping in legumes.

A set of 2486 single nucleotide polymorphisms (SNPs) were compiled in chickpea using four approaches, namely (i) Solexa/Illumina sequencing (1409), (ii) amplicon sequencing of tentative orthologous genes (TOGs) (604), (iii) mining of expressed sequence tags (ESTs) (286) and (iv) sequencing of candidate genes (187). Conversion of these SNPs to the cost-effective and flexible throughput Competitive Allele Specific PCR (KASPar) assays generated successful assays for 2005 SNPs. These marker assays have been designated as Chickpea KASPar Assay Markers (CKAMs). Screening of 70 genotypes including 58 diverse chickpea accessions and 12 BC(3) F(2) lines showed 1341 CKAMs as being polymorphic. Genetic analysis of these data clustered chickpea accessions based on geographical origin. Genotyping data generated for 671 CKAMs on the reference mapping population (Cicer arietinum ICC 4958 × Cicer reticulatum PI 489777) were compiled with 317 unpublished TOG-SNPs and 396 published markers for developing the genetic map. As a result, a second-generation genetic map comprising 1328 marker loci including novel 625 CKAMs, 314 TOG-SNPs and 389 published marker loci with an average inter-marker distance of 0.59 cM was constructed. Detailed analyses of 1064 mapped loci of this second-generation chickpea genetic map showed a higher degree of synteny with genome of Medicago truncatula, followed by Glycine max, Lotus japonicus and least with Vigna unguiculata. Development of these cost-effective CKAMs for SNP genotyping will be useful not only for genetics research and breeding applications in chickpea, but also for utilizing genome information from other sequenced or model legumes.

Comparative evolution of photosynthetic genes in response to polyploid and nonpolyploid duplication.

The likelihood of duplicate gene retention following polyploidy varies by functional properties (e.g. gene ontologies or protein family domains), but little is known about the effects of whole-genome duplication on gene networks related by a common physiological process. Here, we examined the effects of both polyploid and nonpolyploid duplications on genes encoding the major functional groups of photosynthesis (photosystem I, photosystem II, the light-harvesting complex, and the Calvin cycle) in the cultivated soybean (Glycine max), which has experienced two rounds of whole-genome duplication. Photosystem gene families exhibit retention patterns consistent with dosage sensitivity (preferential retention of polyploid duplicates and elimination of nonpolyploid duplicates), whereas Calvin cycle and light-harvesting complex gene families do not. We observed similar patterns in barrel medic (Medicago truncatula), which shared the older genome duplication with soybean but has evolved independently for approximately 50 million years, and in Arabidopsis (Arabidopsis thaliana), which experienced two nested polyploidy events independent from the legume duplications. In both soybean and Arabidopsis, Calvin cycle gene duplicates exhibit a greater capacity for functional differentiation than do duplicates within the photosystems, which likely explains the greater retention of ancient, nonpolyploid duplicates and larger average gene family size for the Calvin cycle relative to the photosystems.

Integration of physical and genetic maps of common bean through BAC-derived microsatellite markers.

BACKGROUND: Common bean (Phaseolus vulgaris L.) is the most important legume for direct human consumption and the goal of this study was to integrate a recently constructed physical map for the species with a microsatellite based genetic map using a BAC library from the genotype G19833 and the recombinant inbred line population DOR364 x G19833. RESULTS: We searched for simple sequence repeats (SSRs) in the 89,017 BAC-end sequences (BES) from the physical map and genetically mapped any polymorphic BES-SSRs onto the genetic map. Among the BES it was possible to identify 623 contig-linked SSRs, most of which were highly AT-rich. A subgroup of 230 di-nucleotide and tri-nucleotide based SSR primer pairs from these BACs was tested on the mapping parents with 176 single copy loci and 114 found to be polymorphic markers. Of these, 99 were successfully integrated into the genetic map. The 99 linkages between the genetic and physical maps corresponded to an equal number of contigs containing a total of 5,055 BAC clones. CONCLUSIONS: Class II microsatellites were more common in the BES than longer class I microsatellites. Both types of markers proved to be valuable for linking BAC clones to the genetic map and were successfully placed across all 11 linkage groups. The integration of common bean physical and genetic maps is an important part of comparative genome analysis and a prelude to positional cloning of agronomically important genes for this crop.

Fractionation of synteny in a genomic region containing tandemly duplicated genes across glycine max, Medicago truncatula, and Arabidopsis thaliana.

Extended comparison of gene sequences found on homeologous soybean Bacterial Artificial Chromosomes to Medicago truncatula and Arabidopsis thaliana genomic sequences demonstrated a network of synteny within conserved regions interrupted by gene addition and/or deletions. Consolidation of gene order among all 3 species provides a picture of ancestral gene order. The observation supports a genome history of fractionation resulting from gene loss/addition and rearrangement. In all 3 species, clusters of N-hydroxycinnamoyl/benzoyltransferase genes were identified in tandemly duplicated clusters. Parsimony-based gene trees suggest that the genes within the arrays have independently undergone tandem duplication in each species.

Gene duplication and paleopolyploidy in soybean and the implications for whole genome sequencing.

BACKGROUND: Soybean, Glycine max (L.) Merr., is a well documented paleopolyploid. What remains relatively under characterized is the level of sequence identity in retained homeologous regions of the genome. Recently, the Department of Energy Joint Genome Institute and United States Department of Agriculture jointly announced the sequencing of the soybean genome. One of the initial concerns is to what extent sequence identity in homeologous regions would have on whole genome shotgun sequence assembly. RESULTS: Seventeen BACs representing approximately 2.03 Mb were sequenced as representative potential homeologous regions from the soybean genome. Genetic mapping of each BAC shows that 11 of the 20 chromosomes are represented. Sequence comparisons between homeologous BACs shows that the soybean genome is a mosaic of retained paleopolyploid regions. Some regions appear to be highly conserved while other regions have diverged significantly. Large-scale "batch" reassembly of all 17 BACs combined showed that even the most homeologous BACs with upwards of 95% sequence identity resolve into their respective homeologous sequences. Potential assembly errors were generated by tandemly duplicated pentatricopeptide repeat containing genes and long simple sequence repeats. Analysis of a whole-genome shotgun assembly of 80,000 randomly chosen JGI-DOE sequence traces reveals some new soybean-specific repeat sequences. CONCLUSION: This analysis investigated both the structure of the paleopolyploid soybean genome and the potential effects retained homeology will have on assembling the whole genome shotgun sequence. Based upon these results, homeologous regions similar to those characterized here will not cause major assembly issues.

Sequence conservation of homeologous bacterial artificial chromosomes and transcription of homeologous genes in soybean (Glycine max L. Merr.).

The paleopolyploid soybean genome was investigated by sequencing homeologous BAC clones anchored by duplicate N-hydroxycinnamoyl/benzoyltransferase (HCBT) genes. The homeologous BACs were genetically mapped to linkage groups C1 and C2. Annotation of the 173,747- and 98,760-bp BACs showed that gene conservation in both order and orientation is high between homeologous regions with only a single gene insertion/deletion and local tandem duplications differing between the regions. The nucleotide sequence conservation extends into intergenic regions as well, probably due to conserved regulatory sequences. Most of the homeologs appear to have a role in either transcription/DNA binding or cellular signaling, suggesting a potential preference for retention of duplicate genes with these functions. Reverse transcriptase-PCR analysis of homeologs showed that in the tissues sampled, most homeologs have not diverged greatly in their transcription profiles. However, four cases of changes in transcription were identified, primarily in the HCBT gene cluster. Because a mapped locus corresponds to a soybean cyst nematode (SCN) QTL, the potential role of HCBT genes in response to SCN is discussed. These results are the first sequenced-based analysis of homeologous BACs in soybean, a diploidized paleopolyploid.

Lung development in laminin gamma2 deficiency: abnormal tracheal hemidesmosomes with normal branching morphogenesis and epithelial differentiation.

BACKGROUND: Laminin gamma2 (Lamc2), one of the polypeptides in laminin-332 (laminin-5), is prominent in the basement membrane of alveolar walls and airways of developing and adult lung. Laminins are important for lung morphogenesis and based on its localization, a function for laminin gamma2 in lung development has been hypothesized. Targeted deletion of the laminin gamma2 gene in mice results in skin blistering and neonatal death at 3-5 days after birth due to failure to thrive. METHODS: Examination of lung development in Lamc2-/- mice through 1-2 days postnatal was accomplished by morphometric analysis, lung bud culture, electron microscopy, immunohistochemical and immunofluorescence staining. RESULTS: Compared to littermate controls, Lamc2-/- lungs were similar in morphology during embryonic life. At post-natal day 1-2, distal saccules were mildly dilated by chord length measurements. Epithelial differentiation as evaluated by immunohistochemical staining for markers of ciliated cells, Clara cells, alveolar type I cells and alveolar type II cells did not reveal a difference between Lamc2-/- and littermate control lungs. Likewise, vascular development, smooth muscle cell differentiation, and elastic fiber formation looked similar, as did airway basement membrane ultrastructure. Branching morphogenesis by lung bud culture was similar in Lamc2-/- and littermate control lungs. Since laminin-332 is important for hemidesmosome formation, we examined the structure of tracheal hemidesmosomes by transmission electron microscopy. Compared to littermate controls, Lamc2-/- tracheal hemidesmosomes were less organized and lacked the increased electron density associated with the basement membrane abutting the hemidesmosome. CONCLUSION: These findings indicate that laminin gamma2 and laminin-332, despite their prominence in the lung, have a minimal role in lung development through the saccular stage.

Transcript Polymorphism Rates in Soybean Seed Tissue Are Increased in a Single Transformant of Glycine max.

Transgenic crops have been utilized for decades to enhance agriculture and more recently have been applied as bioreactors for manufacturing pharmaceuticals. Recently, we investigated the gene expression profiles of several in-house transgenic soybean events, finding one transformant group to be consistently different from our controls. In the present study, we examined polymorphisms and sequence variations in the exomes of the same transgenic soybean events. We found that the previously dissimilar soybean line also exhibited markedly increased levels of polymorphisms within mRNA transcripts from seed tissue, many of which are classified as gene expression modifiers. The results from this work will direct future investigations to examine novel SNPs controlling traits of great interest for breeding and improving transgenic soybean crops. Further, this study marks the first work to investigate SNP rates in transgenic soybean seed tissues and demonstrates that while transgenesis may induce abundant unanticipated changes in gene expression and nucleotide variation, phenotypes and overall health of the plants examined remained unaltered.

A Consensus Map in Cultivated Hexaploid Oat Reveals Conserved Grass Synteny with Substantial Subgenome Rearrangement.

Hexaploid oat ( L., 2 = 6 = 42) is a member of the Poaceae family and has a large genome (∼12.5 Gb) containing 21 chromosome pairs from three ancestral genomes. Physical rearrangements among parental genomes have hindered the development of linkage maps in this species. The objective of this work was to develop a single high-density consensus linkage map that is representative of the majority of commonly grown oat varieties. Data from a cDNA-derived single-nucleotide polymorphism (SNP) array and genotyping-by-sequencing (GBS) were collected from the progeny of 12 biparental recombinant inbred line populations derived from 19 parents representing oat germplasm cultivated primarily in North America. Linkage groups from all mapping populations were compared to identify 21 clusters of conserved collinearity. Linkage groups within each cluster were then merged into 21 consensus chromosomes, generating a framework consensus map of 7202 markers spanning 2843 cM. An additional 9678 markers were placed on this map with a lower degree of certainty. Assignment to physical chromosomes with high confidence was made for nine chromosomes. Comparison of homeologous regions among oat chromosomes and matches to orthologous regions of rice ( L.) reveal that the hexaploid oat genome has been highly rearranged relative to its ancestral diploid genomes as a result of frequent translocations among chromosomes. Heterogeneous chromosome rearrangements among populations were also evident, probably accounting for the failure of some linkage groups to match the consensus. This work contributes to a further understanding of the organization and evolution of hexaploid grass genomes.

Population Genomics Related to Adaptation in Elite Oat Germplasm.

Six hundred thirty five oat ( L.) lines and 4561 single-nucleotide polymorphism (SNP) loci were used to evaluate population structure, linkage disequilibrium (LD), and genotype-phenotype association with heading date. The first five principal components (PCs) accounted for 25.3% of genetic variation. Neither the eigenvalues of the first 25 PCs nor the cross-validation errors from = 1 to 20 model-based analyses suggested a structured population. However, the PC and = 2 model-based analyses supported clustering of lines on spring oat vs. southern United States origin, accounting for 16% of genetic variation ( < 0.0001). Single-locus -statistic () in the highest 1% of the distribution suggested linkage groups that may be differentiated between the two population subgroups. Population structure and kinship-corrected LD of = 0.10 was observed at an average pairwise distance of 0.44 cM (0.71 and 2.64 cM within spring and southern oat, respectively). On most linkage groups LD decay was slower within southern lines than within the spring lines. A notable exception was found on linkage group Mrg28, where LD decay was substantially slower in the spring subpopulation. It is speculated that this may be caused by a heterogeneous translocation event on this chromosome. Association with heading date was most consistent across location-years on linkage groups Mrg02, Mrg12, Mrg13, and Mrg24.

CONTRAILS: A tool for rapid identification of transgene integration sites in complex, repetitive genomes using low-coverage paired-end sequencing.

Transgenic crops have become a staple in modern agriculture, and are typically characterized using a variety of molecular techniques involving proteomics and metabolomics. Characterization of the transgene insertion site is of great interest, as disruptions, deletions, and genomic location can affect product selection and fitness, and identification of these regions and their integrity is required for regulatory agencies. Here, we present CONTRAILS (Characterization of Transgene Insertion Locations with Sequencing), a straightforward, rapid and reproducible method for the identification of transgene insertion sites in highly complex and repetitive genomes using low coverage paired-end Illumina sequencing and traditional PCR. This pipeline requires little to no troubleshooting and is not restricted to any genome type, allowing use for many molecular applications. Using whole genome sequencing of in-house transgenic Glycine max, a legume with a highly repetitive and complex genome, we used CONTRAILS to successfully identify the location of a single T-DNA insertion to single base resolution.

Large-scale development of cost-effective single-nucleotide polymorphism marker assays for genetic mapping in pigeonpea and comparative mapping in legumes.

Single-nucleotide polymorphisms (SNPs, >2000) were discovered by using RNA-seq and allele-specific sequencing approaches in pigeonpea (Cajanus cajan). For making the SNP genotyping cost-effective, successful competitive allele-specific polymerase chain reaction (KASPar) assays were developed for 1616 SNPs and referred to as PKAMs (pigeonpea KASPar assay markers). Screening of PKAMs on 24 genotypes [23 from cultivated species and 1 wild species (Cajanus scarabaeoides)] defined a set of 1154 polymorphic markers (77.4%) with a polymorphism information content (PIC) value from 0.04 to 0.38. One thousand and ninety-four PKAMs showed polymorphisms between parental lines of the reference mapping population (C. cajan ICP 28 × C. scarabaeoides ICPW 94). By using high-quality marker genotyping data on 167 F(2) lines from the population, a comprehensive genetic map comprising 875 PKAMs with an average inter-marker distance of 1.11 cM was developed. Previously mapped 35 simple sequence repeat markers were integrated into the PKAM map and an integrated genetic map of 996.21 cM was constructed. Mapped PKAMs showed a higher degree of synteny with the genome of Glycine max followed by Medicago truncatula and Lotus japonicus and least with Vigna unguiculata. These PKAMs will be useful for genetics research and breeding applications in pigeonpea and for utilizing genome information from other legume species.

Draft genome sequence of pigeonpea (Cajanus cajan), an orphan legume crop of resource-poor farmers.

Pigeonpea is an important legume food crop grown primarily by smallholder farmers in many semi-arid tropical regions of the world. We used the Illumina next-generation sequencing platform to generate 237.2 Gb of sequence, which along with Sanger-based bacterial artificial chromosome end sequences and a genetic map, we assembled into scaffolds representing 72.7% (605.78 Mb) of the 833.07 Mb pigeonpea genome. Genome analysis predicted 48,680 genes for pigeonpea and also showed the potential role that certain gene families, for example, drought tolerance-related genes, have played throughout the domestication of pigeonpea and the evolution of its ancestors. Although we found a few segmental duplication events, we did not observe the recent genome-wide duplication events observed in soybean. This reference genome sequence will facilitate the identification of the genetic basis of agronomically important traits, and accelerate the development of improved pigeonpea varieties that could improve food security in many developing countries.

Epithelial laminin alpha5 is necessary for distal epithelial cell maturation, VEGF production, and alveolization in the developing murine lung.

Laminin alpha5 is prominent in the basement membrane of alveolar walls, airways, and pleura in developing and adult lung. Targeted deletion of laminin alpha5 in mice causes developmental defects in multiple organs, but embryonic lethality has precluded examination of the latter stages of lung development. To identify roles for laminin alpha5 in lung development, we have generated an inducible lung epithelial cell-specific Lama5 null (SP-CLama5(fl/-)) mouse through use of the Cre/loxP system, the human surfactant protein C promoter, and the reverse tetracycline transactivator. SP-CLama5(fl/-) embryos exposed to doxycycline from E6.5 died a few hours after birth. Compared to control littermates, SP-CLama5(fl/-) lungs had dilated, enlarged distal airspaces, but basement membrane ultrastructure was preserved. Distal epithelial cell differentiation was perturbed, with a marked reduction of alveolar type II cells and a virtual absence of type I cells. Cell proliferation was reduced and apoptosis was increased. Capillary density was diminished, and this was associated with a decrease in total lung VEGF production. Overall, these findings indicate that epithelial laminin alpha5, independent of its structural function, is necessary for murine lung development, and suggest a role for laminin alpha5 in signaling pathways that promote alveolar epithelial cell differentiation and VEGF expression.